Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,508,136 A→C 21.8% V36G (GTG→GGG)  fryB ← putative enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,508,1360AC21.8% 123.6 / 5.1 69V36G (GTG→GGG) fryBputative enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base A (36/14);  new base C (5/11);  total (42/27)
Fisher's exact test for biased strand distribution p-value = 6.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGT  >  NC_000913/2507900‑2508368
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
ttATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGCCAGGCGATTCTGG‑ATACCGTCCCGGCCCTGGGGTTCAATTTTCACt                                                                                                                                                                                                                                          >  2:709552/1‑239 (MQ=255)
    cccTTTCTAATTTTTCTTTTTTTCAGAGGCAATCCTCTCTTCGATTTTTTTGATGATGCCCGCAGCGTTTTTAATTGCTTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCAt                                                                                                                                                                                                                                      <  2:397102/239‑1 (MQ=255)
      ctttCGAATTTTATTTTTTTTCAGGGGCAATTTTTCTTTCGATTTTTTTGATGATTCCCGGAGGTTTTTTAATTGCGTTCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTTGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATaa                                                                                                                                                                                                                                    <  2:417004/236‑1 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAACCGTTCGTTATCTTCAGGGGTAACTGCCACGGAGTGGATGATGATGGTCGATCCGGCGATATCCTGCGCCGTCAGCCGCTTCTGG‑ATACGGTCCGCGCCCTGGGCTTCAATTTTAGCTTCAGaaac                                                                                                                                                                                                                                  >  1:424720/1‑237 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGTTATCCTGCGCCGTCAGTCGATTCTGG‑TTACCGTCCGCGCCCTGGGTTTCAATTTTAACTTCCTac                                                                                                                                                                                                                                    >  1:957098/1‑236 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCTTAAcc                                                                                                                                                                                                                                  >  1:113865/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAAcc                                                                                                                                                                                                                                  >  1:828840/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAAcc                                                                                                                                                                                                                                  >  1:148244/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAAcc                                                                                                                                                                                                                                  >  1:180679/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCAAAAcc                                                                                                                                                                                                                                  >  2:841358/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGCCAGGCGATTCTGG‑ATACAGTCCGCGCCCTGGGTTTCCATTTTCACTTCTTAAcc                                                                                                                                                                                                                                  >  1:44817/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCCCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGCCGCTTCGGCGATATCCTGAGCCGTCAGGCGATTTTGG‑ATACCGTCCGCTCCCTGGGTTTCAATTTTCCCTTCTTAAcc                                                                                                                                                                                                                                  >  1:43607/1‑239 (MQ=255)
         tCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATAGCCTGCGCCGTCAGGCGATTCTGG‑TAACCGTCCGCGCCCTGGGTTTCAATTTTCCCTTCAAAACCg                                                                                                                                                                                                                                 >  2:392423/1‑239 (MQ=255)
         tCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGCCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCTTAACCg                                                                                                                                                                                                                                 >  1:397102/1‑239 (MQ=255)
         acagatatTACTGTTGGTCAGAGGCAATCATCGCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGGCGCGTTATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAAc                                                                                                                                                                                                                                   <  1:902466/233‑1 (MQ=255)
                   cTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCCATTTTCCCTTCATAACCGG‑CTTCTA‑cc                                                                                                                                                                                                                      >  1:222551/1‑238 (MQ=255)
                        gttCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAACCGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCCCTTCATAACCGG‑CTTCTA‑CCCCCGct                                                                                                                                                                                                                >  1:253153/1‑239 (MQ=255)
                        gttCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGTTGGTCGCTTCGGCGATTTCCTGCGCCGTCAGGCGGTTCTGG‑TTACCGTCCGCGCCCTGGGTTTCAATTTTTCCTTCTTAAACGG‑CTTTTA‑ccccgtc                                                                                                                                                                                                                 >  2:902466/1‑233 (MQ=255)
                        gttCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGct                                                                                                                                                                                                                >  1:70129/1‑239 (MQ=255)
                         ttCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGctt                                                                                                                                                                                                               >  1:117543/1‑239 (MQ=255)
                         ttCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGct                                                                                                                                                                                                                >  1:207262/1‑238 (MQ=255)
                           cAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACAGG‑CTTCTA‑CCGCCGcttct                                                                                                                                                                                                             >  1:292423/1‑239 (MQ=255)
                              agGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTTCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCttcttc                                                                                                                                                                                                           >  1:555782/1‑238 (MQ=255)
                              agGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGGCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTAACTTCATAACCGG‑CTTCTA‑CCGCCGCTCCTTcc                                                                                                                                                                                                          >  1:1088394/1‑239 (MQ=255)
                                   atcatcTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAgcg                                                                                                                                                                                                      >  1:392152/1‑238 (MQ=255)
                                   atcatcTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAgcc                                                                                                                                                                                                      >  1:423143/1‑237 (MQ=255)
                                       ttcttttgttttCTTTGATGGTGCCCGGAGCGTTTTTAATTGCGTCCTGAAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGa                                                                                                                                                                                                 <  1:353089/230‑1 (MQ=255)
                                        ctctTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCTTAACCGG‑CTTCTA‑CCGCCGCTTCTTCCCGCGCCTGa                                                                                                                                                                                                 >  1:80504/1‑238 (MQ=255)
                                           ttCGATTTCTTTGATGATGCCCGCAGCTTTTTTAATTGCGTCCTGCAAAGTGATTTCATAACCGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCgg                                                                                                                                                                                             <  1:3326/239‑1 (MQ=255)
                                              gATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCCCCCCTTCTTCCACCCCCTGAGCGGccc                                                                                                                                                                                          >  2:545755/1‑238 (MQ=255)
                                                   cTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCGGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCCTCTTCCAGCGCCTGACCGCCCATaaaa                                                                                                                                                                                     >  2:284291/1‑238 (MQ=255)
                                                              cccGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCCTAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGCGTGACCCa                                                                                                                                                                           >  1:179071/1‑238 (MQ=255)
                                                                     cGTTTTTAATTGCGTCCTGCAAGGTGATTTCAAAAACGGCGCGTGATGCGAAACGTTCGTTATCTTCCGGGGTAACTGCCAGGGAGGGGATGATGATGGTCGCTGCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCAt                                                                                                                                                                    <  1:392423/238‑1 (MQ=255)
                                                                           tAATTGGGTCCTGAAAAGTGATTTCAAAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCCCTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGggcag                                                                                                                                                             <  1:709552/239‑1 (MQ=255)
                                                                               tGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGcac                                                                                                                                                         >  1:35346/1‑239 (MQ=255)
                                                                                                tttCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGtt                                                                                                                                          >  1:123825/1‑237 (MQ=255)
                                                                                                tttCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTCCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTTGGGGCCAGGCCCTTCGGGCAGGACCATAAGGCAATCAGgttt                                                                                                                                        >  2:242432/1‑238 (MQ=255)
                                                                                                      aaaGGTCGCGTGATTCAAAACTTTCTTTATCTTCGGGGGTAACTGCCGCGGGGTGGATGATGATGGTCGTTTCGGCGATTTCCGGCGCCGTCAGGCGATTCTGG‑ATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTc                                                                                                                                  <  1:688477/239‑1 (MQ=255)
                                                                                                            cgcgTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGtt                                                                                                                            >  1:889313/1‑239 (MQ=255)
                                                                                                              cgTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATATAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTCCTCATCTTTaa                                                                                                                          >  1:704604/1‑239 (MQ=255)
                                                                                                                      gAAGCGTTCGTTAACTTCCGGGGTAACTGCCCGGGGGTGGATGATGATGGCCGTTTCGGGGATTCCCTGCGCCTTCAGGCGATTCTGG‑ATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGt                                                                                                                  <  1:810811/239‑1 (MQ=255)
                                                                                                                             tcgttATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCCGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCCGTGa                                                                                                           >  1:897522/1‑239 (MQ=255)
                                                                                                                                      gcatgggAAACTGCCACGGACGGGATGATGATGGTCGCTTCGGCGATAGCCTGCGCCGTCAGGCGATTCTGG‑ATAACGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCGATAAGGTGTGAACCAGACCCAACGGGCCGGCACCTAAGGCCACACGTTCCTTACTCATCGTTCATCACCCGTTCCAGGATAGGTGACg                                                                                                  >  2:579274/5‑239 (MQ=255)
                                                                                                                                                     aCGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCCTAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACCTAAGGCAATCAGTTTCTTACTCCTCGTTAATTCCTCGTCCAGGGATATTTTACGGGGGAAATTATGacc                                                                                   >  2:877624/1‑237 (MQ=255)
                                                                                                                                                           gtgatgaggatGTTCTCTTCGGCGATATCCTGCGCCGTCAGGCTATTCTGG‑ATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCAA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCa                                                                              <  1:841358/236‑1 (MQ=255)
                                                                                                                                                                      ggTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCCGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCCGTGATAGTTGAAGAGGGAAATTATGAAAGGGCAAAACAAATCGcc                                                                  >  2:577374/1‑239 (MQ=255)
                                                                                                                                                                       gTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATcccct                                                                 >  1:722158/1‑238 (MQ=255)
                                                                                                                                                                        tCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGt                                                                >  1:1014440/1‑239 (MQ=255)
                                                                                                                                                                             tCGGGGATACCCTGCGCGGTCGGGCGTTTCGGG‑ATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑TTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                                                           <  1:545755/239‑1 (MQ=255)
                                                                                                                                                                              ggggtttttccggctccggcAGGCGA‑TCTGGTATCCCGCCCCCTCCCTGGTTTTAATGTTTCCCTAAAAAACCGC‑CTTCTC‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAg                                                          <  2:70129/231‑1 (MQ=255)
                                                                                                                                                                                                  cggcttttCTGG‑TTTCCGTCCGCGCCCTGGGTTTAAATTTTCATTTCATACCCGG‑TTGCTA‑CCGCCGCTTCTGCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTgg                                      <  2:180679/233‑1 (MQ=255)
                                                                                                                                                                                                       aTTCTGG‑TACCCGGCCCCGCCCGGGGTTACAATTTCCTCTTCAACGCGGC‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAAt                                 <  1:503285/239‑1 (MQ=255)
                                                                                                                                                                                                              gcctccGTCCGCGCCGGGGGTTCAAATTTTCCCTTATAAACCGG‑CTTCTC‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGGGTGAGCCAGGCCCATCGTGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATcc                           <  2:1021030/235‑1 (MQ=255)
                                                                                                                                                                                                              gagcccGGCCGCGCCTGGGGTTACAATTTCCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCAGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATcc                           <  1:877624/235‑1 (MQ=255)
                                                                                                                                                                                                                tcccccccGCCCCGTGGGTTTAATTTTTCCCTATAAACCGGG‑TTTTTC‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGGGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCg                         <  2:207262/233‑1 (MQ=255)
                                                                                                                                                                                                                ctccGTCCGCGCCCGGGGTTCCATTTTCCCCTTCATACCCGT‑CTTCTC‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCg                         <  1:564067/237‑1 (MQ=255)
                                                                                                                                                                                                                ccagttCCCCGCCCTGGGTTCCAATTTGCGCTTCATCCCCGT‑CTTCTA‑CCGTCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTTAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATcc                           <  1:519871/232‑1 (MQ=255)
                                                                                                                                                                                                                atcccccccgccccTGGGGTTTAAATTTCACCTCAAAACCCG‑CTTTTACCCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTTGGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATcc                           <  2:35346/227‑1 (MQ=255)
                                                                                                                                                                                                                 ccggTCCGCGCCCGGGGTTTCAATTTTCCTTTCATACCGGG‑TTACAC‑CCGCCTCTGCTTCCAGCGCCTGAGCGGCCATAAAGGTGTTAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCg                         <  1:284291/235‑1 (MQ=255)
                                                                                                                                                                                                                 cccGGCCGCCCCTGGGGTTACATTTTTCCATTAAGAAGCGG‑TTTATC‑CCGCCCCTTCTTCCGGGGCCTGGGGGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTAGTTCAGTGATAGTTGACGAGGGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCg                         <  2:407056/237‑1 (MQ=255)
                                                                                                                                                                                                                   cGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACCTAAGGCAATCAGCTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGTACTTATGAACGTGCAAAACAATTCGCCGTTACCAGCCCAATAATGCATTTACTGCAGAATTAATTCCGCtt                      >  1:508680/1‑238 (MQ=255)
                                                                                                                                                                                                                   cGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGAATTTACTGGAGAATTAATccc                          >  1:692373/1‑235 (MQ=255)
                                                                                                                                                                                                                     tCCGCCCCCTGGGTTTCAATTTTCCCTTCAAAACCGG‑CTTCTC‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTa                    <  1:880968/239‑1 (MQ=255)
                                                                                                                                                                                                                     ccccggcCCTGGTTTTCATTTTTCCCTTCATACCGGG‑TTTCAC‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTa                    <  2:123825/234‑1 (MQ=255)
                                                                                                                                                                                                                        gcccccTGGGTTTAAATTTTCCCTTCAGACCCGT‑CTTCCA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACgg                 <  1:242432/236‑1 (MQ=255)
                                                                                                                                                                                                                               gggTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGACGTTAACGGACAAAAAATGCATTTACTGGAAAATTAATCCCGCTGTACGGacggtc           >  2:519871/1‑233 (MQ=255)
                                                                                                                                                                                                                                ggTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCGGTACGGAAggggga         >  2:762121/1‑234 (MQ=255)
                                                                                                                                                                                                                                   ttCCATTTTTCCCTCAAAACCGGTATTCCA‑CCCCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGGGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACg       <  2:44817/239‑1 (MQ=255)
                                                                                                                                                                                                                                       aTTTTCACTTCATACCCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGt  <  1:577374/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
TTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGG‑ATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGG‑CTTCTA‑CCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGT  >  NC_000913/2507900‑2508368

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