Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,077,473 C→T Q90Q (CAG→CAA)  unknown ← DNA topoisomerase I (EC 5.99.1.2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,077,4730CT98.2% 144.5 / ‑5.0 56Q90Q (CAG→CAA) unknownDNA topoisomerase I (EC 5.99.1.2)
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (19/36);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGTGCGTTAACGCGATCCATATTCAAATGTTCGGGGTGTTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATCTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGCTGGATTTCACTACATAATTGTTTCCTAAATATTTATTAATGGTTTT  >  1/1077229‑1077700
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
gCTGTGCGTTAACGCGATCCATATTCAAATGTTCGGGGTGTTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGAGGAGAATTTGATAg                                                                                                                                                                                                                               >  1:662273/1‑251 (MQ=255)
 cTGTGCGTTAACGCGATCCATATTCAAATGTTCGGGGTGTTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAg                                                                                                                                                                                                                               <  1:824490/250‑1 (MQ=255)
     gCGTTAACGCGATCCATATTCAAATGTTCGGGGTGTTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAg                                                                                                                                                                                                                               <  1:250537/246‑1 (MQ=255)
                  ccATATTCAAATGTTCGGGGTGTTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGa                                                                                                                                                                                                             >  1:220589/1‑251 (MQ=255)
                    atatTCAAATGTTCGGGGTGTTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCttt                                                                                                                                                                                                                       >  1:676535/1‑239 (MQ=255)
                    atatTCAAATGTTCGGGGTGTTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATa                                                                                                                                                                                                           <  1:832304/251‑1 (MQ=255)
                     tatTCAAATGTTCGGGGTGTTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCAC‑‑‑GCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAggg                                                                                                                                                                                                       <  1:335938/251‑1 (MQ=255)
                               ttCGGGGTGTTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTCTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGc                                                                                                                                                                                                <  1:807633/251‑1 (MQ=255)
                                     gtgtTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAgggg                                                                                                                                                                                                       >  1:492993/1‑238 (MQ=255)
                                     gtgtTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATg                                                                                                                                                                                           <  1:657233/250‑1 (MQ=255)
                                      tgtTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGTGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGc                                                                                                                                                                                          <  1:596429/250‑1 (MQ=255)
                                      tgtTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCg                                                                                                                                                                                         <  1:488718/251‑1 (MQ=255)
                                       gtTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAAc                                                                                                                                                                                                                  <  1:185186/225‑1 (MQ=255)
                                              ggCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGt                                                                                                                                                                                                      <  1:230284/230‑1 (MQ=255)
                                               gCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACt                                                                                                                                                                                  <  1:110686/249‑1 (MQ=255)
                                                   gTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAg                                                                                                                                                                              >  1:728997/1‑249 (MQ=255)
                                                   gTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAg                                                                                                                                                                              >  1:692413/1‑249 (MQ=255)
                                                              tttttGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCttt                                                                                                                                                                 <  1:783441/251‑1 (MQ=255)
                                                                           cATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTcc                                                                                                                                                                         >  1:346437/1‑230 (MQ=255)
                                                                           cATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTAt                                                                                                                                                     >  1:449772/1‑250 (MQ=255)
                                                                            aTTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCg                                                                                                                                                   <  1:171046/251‑1 (MQ=255)
                                                                               aaaCACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGc                                                                                                                                                  >  1:115192/1‑249 (MQ=255)
                                                                                   cccacGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGttt                                                                                                                                            <  1:287046/250‑1 (MQ=255)
                                                                                     cacGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGttttt                                                                                                                                          <  1:461826/251‑1 (MQ=255)
                                                                                         cgGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTttcttc                                                                                                                                      <  1:116301/251‑1 (MQ=255)
                                                                                           gCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGc                                                                                                                                                                                                <  1:730889/191‑1 (MQ=255)
                                                                                              ataACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCAc                                                                                                                                 <  1:17123/251‑1 (MQ=255)
                                                                                                         catcGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTtg                                                                                                                      <  1:392274/251‑1 (MQ=255)
                                                                                                                gccgATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGaaa                                                                                                               >  1:232874/1‑251 (MQ=255)
                                                                                                                  cgATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTGTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAAtg                                                                                                             <  1:425916/251‑1 (MQ=255)
                                                                                                                       aCTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTtggtgg                                                                                                                  <  1:316974/241‑1 (MQ=255)
                                                                                                                        cTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAAtggctttg                                                                                                       >  1:319964/1‑251 (MQ=255)
                                                                                                                                         gCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATGAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTgg                                                                                       <  1:327785/250‑1 (MQ=255)
                                                                                                                                              cGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACtg                                                                                                                                                                                                                >  1:787568/1‑124 (MQ=255)
                                                                                                                                                   ctcCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTtggtgg                                                                                                                  <  1:59895/213‑1 (MQ=255)
                                                                                                                                                   ctcCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACg                                                                            >  1:899382/1‑251 (MQ=255)
                                                                                                                                                   ctcCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACg                                                                            >  1:182503/1‑251 (MQ=255)
                                                                                                                                                   ctcCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAA‑CACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACgg                                                                           <  1:1002/251‑1 (MQ=255)
                                                                                                                                                        tctcTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCg                                                                                                                                                                                         <  1:397361/137‑1 (MQ=255)
                                                                                                                                                        tctcTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCttttt                                                                                                                                                               <  1:495855/163‑1 (MQ=255)
                                                                                                                                                        tctcTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAAt                                                                                                              <  1:254852/212‑1 (MQ=255)
                                                                                                                                                             aTCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAAt                                                                  >  1:474805/1‑251 (MQ=255)
                                                                                                                                                             aTCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAAt                                                                  >  1:117410/1‑251 (MQ=255)
                                                                                                                                                                       gTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATCTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAatat                                                         <  1:613820/250‑1 (MQ=255)
                                                                                                                                                                              ggTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGc                                                       >  1:900640/1‑245 (MQ=255)
                                                                                                                                                                              ggTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACg                                                 >  1:47480/1‑251 (MQ=255)
                                                                                                                                                                                tAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGc                                                <  1:204251/250‑1 (MQ=255)
                                                                                                                                                                                tAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGc                                                <  1:584885/250‑1 (MQ=255)
                                                                                                                                                                                                        agagaTTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAA‑CCACGGTTGGCAAATCGCGAATATGCCCCAc                                                  <  1:237719/223‑1 (MQ=255)
                                                                                                                                                                                                          agaTTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGCTGGATTTCACTACATAATTg                            <  1:689833/244‑1 (MQ=255)
                                                                                                                                                                                                          agaTTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGCTGGATTTCACTACATAATTGTTTcc                       <  1:621427/249‑1 (MQ=255)
                                                                                                                                                                                                                   ttCGGCAACCCCTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGCTGGATTTCACTACATAATTGTTTCCTAAAtattta             <  1:391303/250‑1 (MQ=255)
                                                                                                                                                                                                                    tCGGCAACCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGCTGGATTTCACTACATAATTGTTTCCTAAAtatttatt           <  1:328043/251‑1 (MQ=255)
                                                                                                                                                                                                                         aaCCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGCTGGATTTCACTACATAATTGTTTCCTAAATATTTATTAATgg      <  1:27042/251‑1 (MQ=255)
                                                                                                                                                                                                                         aaCCACTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGCTGGATTTCACTACATAATTGTTTCCTAAATATTTATTAATgg      <  1:551258/251‑1 (MQ=255)
                                                                                                                                                                                                                              cTTTTTCTTTGCCGGGGAGAATTTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGCTGGATTTCACTACATAATTGTTTCCTAAATATTTATTAATGGtttt  >  1:337419/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
GCTGTGCGTTAACGCGATCCATATTCAAATGTTCGGGGTGTTCAAAGGCTTGTTTAATGGCGTTTTTGGTAATTTCATTAAACACCACGCGGCTATAACGTTCATCATCGCCGCCGATGACTTCACGCAAATGCCAGGCAATCGCCTCTCCTTCTCTATCCAAATCGGTGGCGAGGTAAATATGATCAACTTTTTTTGCGAGAGATTTCAGTTCGGCAACCACTTTTTCTTTGCCGGGGAGAATCTGATAGTTGGCTTTCCAACTGTGATAAGGGTCAATGCCCATGCGTTTGACTAAAGCTTCCTGTTCTTTTTTGTGTTTTATCGCCTGTTTTTCTTCCGCACTGAGCCCTTTGGTGGAAATGGCTTTGGCTTTTTCTCCAGTGGAGGAACCCACGGTTGGCAAATCGCGAATATGCCCCACGCTGGATTTCACTACATAATTGTTTCCTAAATATTTATTAATGGTTTT  >  1/1077229‑1077700

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑