Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,929,115 3 bp→CTT intergenic (‑201/‑218) 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,1150AC93.7% 922.4 / 34.9 366intergenic (‑201/‑220)SAUSA300_RS09560/SAUSA300_RS09565transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/5);  new base C (334/9);  total (352/14)
Fisher's exact test for biased strand distribution p-value = 8.62e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
*NC_0077931,929,1160GT94.0% 1018.3 / 37.3 366intergenic (‑202/‑219)SAUSA300_RS09560/SAUSA300_RS09565transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/5);  new base T (335/9);  total (352/14)
Fisher's exact test for biased strand distribution p-value = 6.89e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01
*NC_0077931,929,1170GT94.0% 1022.0 / 36.7 366intergenic (‑203/‑218)SAUSA300_RS09560/SAUSA300_RS09565transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/5);  new base T (335/9);  total (352/14)
Fisher's exact test for biased strand distribution p-value = 6.89e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 366 total aligned reads displayed.

TATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGT  >  NC_007793/1929042‑1929194
                                                                          |||                                                                             
tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCAcatcat                                                      <  1:4904910/92‑1 (MQ=255)
    cagagatgatATAATATACAAGTCGCCTATCTCTCAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAg                                                    >  1:5084312/5‑100 (MQ=255)
    aacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCa                                                  >  1:2143166/5‑100 (MQ=255)
          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAgat                                            >  1:804867/4‑98 (MQ=255)
               taatatacaagTCGCCTATCTCTCAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGtcatc                                         <  1:5932157/91‑1 (MQ=255)
                       ctgTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                               >  1:1283671/3‑100 (MQ=255)
                        ggTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                              >  1:188852/2‑100 (MQ=255)
                         gTCGCCTATCTCTCAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                               <  1:5431238/100‑1 (MQ=255)
                         gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                              >  1:3761504/1‑100 (MQ=255)
                           ggCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                           >  1:5152104/2‑97 (MQ=255)
                                aTCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                      >  1:5955027/1‑100 (MQ=255)
                                    ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                    >  1:1626455/1‑98 (MQ=255)
                                    ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                    >  1:6326574/1‑98 (MQ=255)
                                    ctcAGTCGTCAACTTAGGCGCAGAGAGGAGGCGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                    >  1:5258960/1‑98 (MQ=255)
                                    ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                  >  1:372821/1‑99 (MQ=255)
                                    ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                  >  1:5375106/1‑99 (MQ=255)
                                    ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                  >  1:2982364/1‑99 (MQ=255)
                                    ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                  >  1:4668648/1‑99 (MQ=255)
                                    ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                  >  1:4930798/1‑99 (MQ=255)
                                    ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                  >  1:891351/1‑99 (MQ=255)
                                    ctcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                  >  1:1348974/1‑99 (MQ=255)
                                    atcAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGCCATCAGTGGCTGTGCTGTTGCGtatta                  >  1:745744/2‑99 (MQ=255)
                                      cAGTCGTCAACTTAGGCGCAGAGAGGAGGGGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                  >  1:3097919/1‑99 (MQ=255)
                                      cAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                >  1:1235934/1‑97 (MQ=255)
                                      cAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                >  1:172567/1‑97 (MQ=255)
                                             cAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt          >  1:4977425/1‑100 (MQ=255)
                                             cAACTTAGGCGCAGAGAGGAGGTGCATTACTTTGCTAACCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt           >  1:1697949/1‑100 (MQ=255)
                                                tatAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:2226034/3‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:6139546/2‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:3142927/2‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:4443496/2‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:3429492/2‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:6158960/2‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:5480596/2‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:5708488/2‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:1936190/2‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:6428607/2‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:2400773/2‑100 (MQ=255)
                                                cttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:6398252/2‑100 (MQ=255)
                                                cctAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:3643356/3‑100 (MQ=255)
                                                agcaggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg        >  1:302699/7‑100 (MQ=255)
                                                 ttAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc       >  1:5115454/1‑100 (MQ=255)
                                                 tccggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc       >  1:1488371/6‑100 (MQ=255)
                                                 ctAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTGTTATACTTATTGGc       >  1:4788237/2‑100 (MQ=255)
                                                 ctAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc       >  1:2669283/2‑100 (MQ=255)
                                                 ctAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc       >  1:5605013/2‑100 (MQ=255)
                                                 ctAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc       >  1:412829/2‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCCGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:36584/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:5230322/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:1019973/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:4511624/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:4319301/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:4742075/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:4181759/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:111218/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:12028/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:1330037/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:1546472/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:6457756/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:2242411/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:2455327/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:2613816/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:2704795/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:2754098/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:6239627/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:3481056/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:3614825/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:6098223/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:3958533/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:4239735/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATAGGCt      >  1:3243895/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCACTCCTCAGTGGCTGCCTTGTTGCGTATTATACTTATTGGCt      >  1:951275/1‑100 (MQ=255)
                                                  tAGGCGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGGGCACATCATAGCACCAGTCATCAGTGACTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:5874497/1‑100 (MQ=255)
                                                  tAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGcca     >  1:5008233/1‑98 (MQ=255)
                                                  ccgACGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:2783522/3‑100 (MQ=255)
                                                  caggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:6033614/5‑100 (MQ=255)
                                                  caggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      >  1:3853334/5‑100 (MQ=255)
                                                   cggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:6247477/4‑99 (MQ=255)
                                                   cggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:3300255/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:3785737/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:1403364/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:1806917/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:6634788/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:6575168/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:2113465/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:2315742/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:2492484/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:5141524/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:4064534/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:4360347/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:5562578/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:4583123/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGCGCTGTTGCGTATTATACTTATTGGCtt     >  1:3922340/4‑99 (MQ=255)
                                                   aggcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGCGGCTGTGCTGTTGCGTATTATACTTATCGGCtc     >  1:4187909/4‑99 (MQ=255)
                                                    cccGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCCTCAGTGGCTGTGCTGTTGCGTATTGTACTTATTGGCtta    >  1:4811617/3‑98 (MQ=255)
                                                    accGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTGTACTTATTGGCtta    >  1:4868828/3‑98 (MQ=255)
                                                     tcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg   >  1:1078013/2‑100 (MQ=255)
                                                     gcGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg   >  1:1038479/2‑100 (MQ=255)
                                                     aCGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCACAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt  <  1:2896195/100‑1 (MQ=255)
                                                      cGCAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt  >  1:4378176/1‑100 (MQ=255)
                                                                          |||                                                                             
TATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGT  >  NC_007793/1929042‑1929194

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