Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,291,428 G→A D242N (GAT→AAT)  SAUSA300_RS06340 → zinc protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,291,4280GA98.9% 540.3 / ‑4.6 179D242N (GAT→AAT) SAUSA300_RS06340zinc protease
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (177/0);  total (179/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

Only 100 of 180 total aligned reads displayed.

CTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACCGGAGGATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAGCGCCTC  >  NC_007793/1291333‑1291521
                                                                                               |                                                                                             
cTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGt                                                                                           >  1:3097273/1‑100 (MQ=255)
 tGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGtt                                                                                          >  1:2414838/1‑100 (MQ=255)
 tGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGtt                                                                                          >  1:1630207/1‑100 (MQ=255)
 tGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGtt                                                                                          >  1:1467448/1‑100 (MQ=255)
              cGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTg                                                                             >  1:1881840/1‑100 (MQ=255)
                       aaaCAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAAt                                                                    >  1:1885522/1‑100 (MQ=255)
                                 aGGATGCTCGTAATAAAGTTAACCAACCCAAACTCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAAt                                                          >  1:1271271/1‑100 (MQ=255)
                                               aaaGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGAc                                            >  1:2515152/1‑100 (MQ=255)
                                                        ttaccCAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTg                                   >  1:5463232/3‑100 (MQ=255)
                                                         aaccCAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTgg                                  >  1:4275477/1‑100 (MQ=255)
                                                         aaccCAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTgg                                  >  1:5566874/1‑100 (MQ=255)
                                                            ccAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGCCTAATGCTTGGttt                               >  1:736988/1‑100 (MQ=255)
                                                             cAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGtttt                              >  1:4727207/1‑100 (MQ=255)
                                                             cAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGtttt                              >  1:2994342/1‑100 (MQ=255)
                                                             cAAAATCGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGtttt                              >  1:921703/1‑100 (MQ=255)
                                                                 gccGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTaaaa                          >  1:2780541/3‑100 (MQ=255)
                                                                 cccGAACGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTaaaa                          >  1:4718083/3‑100 (MQ=255)
                                                                  gtgcaCGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTaaaaa                         >  1:3317786/5‑100 (MQ=255)
                                                                     aaCGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATaa                      >  1:5415655/1‑100 (MQ=255)
                                                                       ctaGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAAc                    >  1:1648747/3‑100 (MQ=255)
                                                                       cGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAAc                    >  1:2695001/1‑100 (MQ=255)
                                                                       cGAGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAAc                    >  1:1132232/1‑100 (MQ=255)
                                                                         aGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCa                  >  1:3769961/1‑100 (MQ=255)
                                                                         aGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCa                  >  1:28818/1‑100 (MQ=255)
                                                                         aGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCa                  >  1:5226107/1‑100 (MQ=255)
                                                                         aGGACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCa                  >  1:2924998/1‑100 (MQ=255)
                                                                          ggACTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCAt                 >  1:4746705/1‑100 (MQ=255)
                                                                             cTTGTTGATGAACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTAc              >  1:170358/1‑100 (MQ=255)
                                                                                agtgatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaag           >  1:2931613/3‑100 (MQ=255)
                                                                                 ctgatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaaga          >  1:3924988/2‑100 (MQ=255)
                                                                                 atgatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaaga          >  1:1115696/2‑100 (MQ=255)
                                                                                  tgatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaa         >  1:1717081/1‑100 (MQ=255)
                                                                                  tgatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaa         >  1:1327936/1‑100 (MQ=255)
                                                                                  tgatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaa         >  1:3821439/1‑100 (MQ=255)
                                                                                   gatgaATCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaag        >  1:2455471/1‑100 (MQ=255)
                                                                                   gatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaag        >  1:2929436/1‑100 (MQ=255)
                                                                                   gatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaag        >  1:4302719/1‑100 (MQ=255)
                                                                                   gatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaag        >  1:3629877/1‑100 (MQ=255)
                                                                                   gatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaag        >  1:4598493/1‑100 (MQ=255)
                                                                                   gatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaag        >  1:2512662/1‑100 (MQ=255)
                                                                                   gatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaag        >  1:1216808/1‑100 (MQ=255)
                                                                                   gatgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACaagaag        >  1:5718271/1‑100 (MQ=255)
                                                                                   gatgaACCGGAGAATGTTAAAGAAGCATTTGGTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCagatcggaaga        >  1:4465621/1‑90 (MQ=255)
                                                                                   cctgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTAAAAAATAAACCATTACaagaag        >  1:5018327/3‑100 (MQ=255)
                                                                                    ttgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:1921323/2‑100 (MQ=255)
                                                                                    gtgaACCGGAGAATGTTAAAGAGGCATTTGTTACTGAACCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:17735/2‑100 (MQ=255)
                                                                                    gtgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:2641073/2‑100 (MQ=255)
                                                                                    gtgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:4626504/2‑100 (MQ=255)
                                                                                    ctgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:4416281/2‑100 (MQ=255)
                                                                                    ctgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:4139775/2‑100 (MQ=255)
                                                                                    ctgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:5776959/2‑100 (MQ=255)
                                                                                    ctgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:1027123/2‑100 (MQ=255)
                                                                                    ctgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:1025119/2‑100 (MQ=255)
                                                                                    ctgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCCTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:4499280/2‑100 (MQ=255)
                                                                                    atgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:2220612/1‑100 (MQ=255)
                                                                                    atgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:1408937/1‑100 (MQ=255)
                                                                                    atgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgc       >  1:5450213/1‑100 (MQ=255)
                                                                                     tgaACCGTAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:804745/1‑100 (MQ=255)
                                                                                     tgaACCGGAGGATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:521775/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:128330/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:5052183/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:2280881/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:2479834/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:2554877/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:5264517/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:2125881/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:2124151/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:3040123/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:2070365/1‑100 (MQ=255)
                                                                                     tgaACCGGAGAATGCTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcg      >  1:1192449/1‑100 (MQ=255)
                                                                                      tagccGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:1693857/4‑100 (MQ=255)
                                                                                      gaACCGGGGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:3977687/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:3589014/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:2138315/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:5114016/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:2002528/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:615428/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:5075797/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:1060977/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:853049/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:4999855/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:3358660/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:364094/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:3538230/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:3690535/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:3900580/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:3945353/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:4263062/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:2970814/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:2583035/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:4485453/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:4550711/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:4789102/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATAAACAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgagc     >  1:2057676/1‑100 (MQ=255)
                                                                                      gaACCGGAGAATGTTAAAGAAGCATTTGGTACTGAATCTATGAAAATTCAATCACCGAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:2260059/1‑100 (MQ=255)
                                                                                      aaaCCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAgcgc     >  1:4240917/2‑100 (MQ=255)
                                                                                       aaCCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAGCGcc    >  1:3293688/1‑100 (MQ=255)
                                                                                        gccGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAGCGCCt   >  1:1886520/2‑100 (MQ=255)
                                                                                        aCCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAGCGCCt   >  1:725316/1‑100 (MQ=255)
                                                                                        aCCGGAGAATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAGCGCCt   >  1:4823163/1‑100 (MQ=255)
                                                                                               |                                                                                             
CTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACCGGAGGATGTTAAAGAAGCATTTGTTACTGAATCTATGAAAATTCAATCACCAAGACTAATGCTTGGTTTTAAAAATAAACCATTACAAGAAGCGCCTC  >  NC_007793/1291333‑1291521

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