Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,929,167 2 bp→TG intergenic (‑253/‑167) SAUSA300_RS09560 ← / → SAUSA300_RS09565 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,929,1670CT96.8% 1685.8 / 11.8 600intergenic (‑253/‑168)SAUSA300_RS09560/SAUSA300_RS09565transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (9/10);  new base T (574/7);  total (583/17)
Fisher's exact test for biased strand distribution p-value = 1.05e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
*NC_0077931,929,1680AG97.7% 1685.9 / 8.5 598intergenic (‑254/‑167)SAUSA300_RS09560/SAUSA300_RS09565transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/9);  new base G (576/7);  total (581/17)
Fisher's exact test for biased strand distribution p-value = 8.68e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 601 total aligned reads displayed.

CGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAT  >  NC_007793/1929069‑1929257
                                                                                                   ||                                                                                         
cGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCa                                                                                           >  1:2972408/1‑100 (MQ=255)
       tgtgcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                   >  1:5444003/5‑98 (MQ=255)
        tctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttt                                                                                   <  1:1692458/100‑1 (MQ=255)
         ttcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCTGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:2920090/2‑99 (MQ=255)
         ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                   >  1:3380702/1‑98 (MQ=255)
         ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                   >  1:4639282/1‑98 (MQ=255)
         ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                   >  1:3081317/1‑98 (MQ=255)
         ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                   >  1:3697334/1‑98 (MQ=255)
         ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                   >  1:5181569/1‑98 (MQ=255)
         ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                   >  1:5343827/1‑98 (MQ=255)
         ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTCCACATCATAGCACCAGTCACCAGTGGCTGTGCTGTTGCGtat                                                                                   >  1:2692583/1‑98 (MQ=255)
         ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTAGTACACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                   >  1:5483723/1‑98 (MQ=255)
         ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:3737120/1‑99 (MQ=255)
         ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:2115822/1‑99 (MQ=255)
         ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:217007/1‑99 (MQ=255)
         ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:2305763/1‑99 (MQ=255)
         ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:4204269/1‑99 (MQ=255)
         ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:2429145/1‑99 (MQ=255)
         ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:2632339/1‑99 (MQ=255)
         ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:3595092/1‑99 (MQ=255)
         ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:4272303/1‑99 (MQ=255)
         ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                 >  1:624282/1‑99 (MQ=255)
          tcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                                >  1:1496129/1‑98 (MQ=255)
           cAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAGTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                               >  1:3474488/1‑97 (MQ=255)
               cGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                             >  1:3834516/1‑95 (MQ=255)
               cGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                             >  1:68008/1‑95 (MQ=255)
                 tCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                           >  1:1864623/1‑100 (MQ=255)
                 ccAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                           >  1:4731026/2‑100 (MQ=255)
                     tatAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                       >  1:3023834/3‑100 (MQ=255)
                     cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTGTTgg                                                                       >  1:3569196/2‑100 (MQ=255)
                     cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                       >  1:4946161/2‑100 (MQ=255)
                     cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                       >  1:4422787/2‑100 (MQ=255)
                     cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                       >  1:1640852/2‑100 (MQ=255)
                     cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                       >  1:3884007/2‑100 (MQ=255)
                     cgtAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                       >  1:2597264/3‑100 (MQ=255)
                     aaaaggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                       >  1:425544/7‑100 (MQ=255)
                      ttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:2206382/1‑100 (MQ=255)
                      ttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:5233072/1‑100 (MQ=255)
                      ttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:4343258/1‑100 (MQ=255)
                      ttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:1147149/1‑100 (MQ=255)
                      tctggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:143035/6‑100 (MQ=255)
                      tcaggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:1994124/6‑100 (MQ=255)
                      gcaggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:1312625/6‑100 (MQ=255)
                      ctAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:4108412/2‑100 (MQ=255)
                      ctAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:1006017/2‑100 (MQ=255)
                      ctAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:2374186/2‑100 (MQ=255)
                      ctAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:5040911/2‑100 (MQ=255)
                      ctAGGCGCAGAGAGGAGGTGCATTACGTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:4984931/2‑100 (MQ=255)
                      atAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                      >  1:2830830/2‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTTCTGTTGCGTATTATACTTATTGTCt                                                                     >  1:4593085/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTTGCt                                                                     >  1:4554442/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:4317173/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:4888987/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:4820018/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:550321/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:5596535/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:800081/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:906856/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:449597/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:4170249/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:941956/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:3998396/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:1064497/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:1247442/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:1338929/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:155958/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:1756912/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:1811120/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:1917035/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:202941/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     <  1:2512845/100‑1 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:2718848/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:2754068/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:2766369/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:289393/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:3206887/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:3248096/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:3949441/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:3783957/1‑100 (MQ=255)
                       tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:3654950/1‑100 (MQ=255)
                       caggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                     >  1:3127682/5‑100 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTGTTATACTTATTGGCtt                                                                    >  1:5157367/4‑99 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    >  1:4464279/4‑99 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    >  1:5762244/4‑99 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    >  1:5720334/4‑99 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    >  1:5667533/4‑99 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    >  1:1282800/4‑99 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    >  1:3340493/4‑99 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    >  1:5282465/4‑99 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    >  1:746910/4‑99 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    >  1:3000597/4‑99 (MQ=255)
                        aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCACAGCACCAGTCATCAGCGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                    >  1:3435575/4‑99 (MQ=255)
                         ggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                   >  1:1214177/3‑98 (MQ=255)
                         agcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                   >  1:4058510/3‑98 (MQ=255)
                          gcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                  >  1:1612362/2‑100 (MQ=255)
                          aCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                  >  1:3303132/1‑100 (MQ=255)
                           cGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                 >  1:3937619/1‑100 (MQ=255)
                           cGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                 >  1:2668076/1‑100 (MQ=255)
                           cGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                 >  1:4858823/1‑100 (MQ=255)
                                                                                        aGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaa    <  1:5070764/100‑1 (MQ=255)
                                                                                                   ||                                                                                         
CGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAT  >  NC_007793/1929069‑1929257

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