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

TTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGT  >  NC_000913/2507900‑2508368
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
ttATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGCCAGGCGATTCTGGATACCGTCCCGGCCCTGGGGTTCAATTTTCACt                                                                                                                                                                                                                                          >  2:709552/1‑239 (MQ=255)
    cccTTTCTAATTTTTCTTTTTTTCAGAGGCAATCCTCTCTTCGATTTTTTTGATGATGCCCGCAGCGTTTTTAATTGCTTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCAt                                                                                                                                                                                                                                      <  2:397102/239‑1 (MQ=255)
      ctttCGAATTTTATTTTTTTTCAGGGGCAATTTTTCTTTCGATTTTTTTGATGATTCCCGGAGGTTTTTTAATTGCGTTCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTTGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATaa                                                                                                                                                                                                                                    <  2:417004/236‑1 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAACCGTTCGTTATCTTCAGGGGTAACTGCCACGGAGTGGATGATGATGGTCGATCCGGCGATATCCTGCGCCGTCAGCCGCTTCTGGATACGGTCCGCGCCCTGGGCTTCAATTTTAGCTTCAGaaac                                                                                                                                                                                                                                  >  1:424720/1‑237 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGTTATCCTGCGCCGTCAGTCGATTCTGGTTACCGTCCGCGCCCTGGGTTTCAATTTTAACTTCCTac                                                                                                                                                                                                                                    >  1:957098/1‑236 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCTTAAcc                                                                                                                                                                                                                                  >  1:113865/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAAcc                                                                                                                                                                                                                                  >  1:828840/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAAcc                                                                                                                                                                                                                                  >  1:148244/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAAcc                                                                                                                                                                                                                                  >  1:180679/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCAAAAcc                                                                                                                                                                                                                                  >  2:841358/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGCCAGGCGATTCTGGATACAGTCCGCGCCCTGGGTTTCCATTTTCACTTCTTAAcc                                                                                                                                                                                                                                  >  1:44817/1‑239 (MQ=255)
        gTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCCCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGCCGCTTCGGCGATATCCTGAGCCGTCAGGCGATTTTGGATACCGTCCGCTCCCTGGGTTTCAATTTTCCCTTCTTAAcc                                                                                                                                                                                                                                  >  1:43607/1‑239 (MQ=255)
         tCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATAGCCTGCGCCGTCAGGCGATTCTGGTAACCGTCCGCGCCCTGGGTTTCAATTTTCCCTTCAAAACCg                                                                                                                                                                                                                                 >  2:392423/1‑239 (MQ=255)
         tCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGCCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCTTAACCg                                                                                                                                                                                                                                 >  1:397102/1‑239 (MQ=255)
         acagatatTACTGTTGGTCAGAGGCAATCATCGCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGGCGCGTTATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAAc                                                                                                                                                                                                                                   <  1:902466/233‑1 (MQ=255)
                   cTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCCATTTTCCCTTCATAACCGGCTTCTAcc                                                                                                                                                                                                                      >  1:222551/1‑238 (MQ=255)
                        gttCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAACCGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCCCTTCATAACCGGCTTCTACCCCCGct                                                                                                                                                                                                                >  1:253153/1‑239 (MQ=255)
                        gttCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGTTGGTCGCTTCGGCGATTTCCTGCGCCGTCAGGCGGTTCTGGTTACCGTCCGCGCCCTGGGTTTCAATTTTTCCTTCTTAAACGGCTTTTAccccgtc                                                                                                                                                                                                                 >  2:902466/1‑233 (MQ=255)
                        gttCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGct                                                                                                                                                                                                                >  1:70129/1‑239 (MQ=255)
                         ttCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGctt                                                                                                                                                                                                               >  1:117543/1‑239 (MQ=255)
                         ttCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGct                                                                                                                                                                                                                >  1:207262/1‑238 (MQ=255)
                           cAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACAGGCTTCTACCGCCGcttct                                                                                                                                                                                                             >  1:292423/1‑239 (MQ=255)
                              agGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTTCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCttcttc                                                                                                                                                                                                           >  1:555782/1‑238 (MQ=255)
                              agGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGGCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTAACTTCATAACCGGCTTCTACCGCCGCTCCTTcc                                                                                                                                                                                                          >  1:1088394/1‑239 (MQ=255)
                                   atcatcTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAgcg                                                                                                                                                                                                      >  1:392152/1‑238 (MQ=255)
                                   atcatcTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAgcc                                                                                                                                                                                                      >  1:423143/1‑237 (MQ=255)
                                       ttcttttgttttCTTTGATGGTGCCCGGAGCGTTTTTAATTGCGTCCTGAAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGa                                                                                                                                                                                                 <  1:353089/230‑1 (MQ=255)
                                        ctctTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCTTAACCGGCTTCTACCGCCGCTTCTTCCCGCGCCTGa                                                                                                                                                                                                 >  1:80504/1‑238 (MQ=255)
                                           ttCGATTTCTTTGATGATGCCCGCAGCTTTTTTAATTGCGTCCTGCAAAGTGATTTCATAACCGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCgg                                                                                                                                                                                             <  1:3326/239‑1 (MQ=255)
                                              gATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCCCCCCTTCTTCCACCCCCTGAGCGGccc                                                                                                                                                                                          >  2:545755/1‑238 (MQ=255)
                                                   cTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCGGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCCTCTTCCAGCGCCTGACCGCCCATaaaa                                                                                                                                                                                     >  2:284291/1‑238 (MQ=255)
                                                              cccGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCCTAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGCGTGACCCa                                                                                                                                                                           >  1:179071/1‑238 (MQ=255)
                                                                     cGTTTTTAATTGCGTCCTGCAAGGTGATTTCAAAAACGGCGCGTGATGCGAAACGTTCGTTATCTTCCGGGGTAACTGCCAGGGAGGGGATGATGATGGTCGCTGCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCAt                                                                                                                                                                    <  1:392423/238‑1 (MQ=255)
                                                                           tAATTGGGTCCTGAAAAGTGATTTCAAAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCCCTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGggcag                                                                                                                                                             <  1:709552/239‑1 (MQ=255)
                                                                               tGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGcac                                                                                                                                                         >  1:35346/1‑239 (MQ=255)
                                                                                                tttCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGtt                                                                                                                                          >  1:123825/1‑237 (MQ=255)
                                                                                                tttCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTCCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTTGGGGCCAGGCCCTTCGGGCAGGACCATAAGGCAATCAGgttt                                                                                                                                        >  2:242432/1‑238 (MQ=255)
                                                                                                      aaaGGTCGCGTGATTCAAAACTTTCTTTATCTTCGGGGGTAACTGCCGCGGGGTGGATGATGATGGTCGTTTCGGCGATTTCCGGCGCCGTCAGGCGATTCTGGATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTc                                                                                                                                  <  1:688477/239‑1 (MQ=255)
                                                                                                            cgcgTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGtt                                                                                                                            >  1:889313/1‑239 (MQ=255)
                                                                                                              cgTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATATAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTCCTCATCTTTaa                                                                                                                          >  1:704604/1‑239 (MQ=255)
                                                                                                                      gAAGCGTTCGTTAACTTCCGGGGTAACTGCCCGGGGGTGGATGATGATGGCCGTTTCGGGGATTCCCTGCGCCTTCAGGCGATTCTGGATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGt                                                                                                                  <  1:810811/239‑1 (MQ=255)
                                                                                                                             tcgttATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCCGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCCGTGa                                                                                                           >  1:897522/1‑239 (MQ=255)
                                                                                                                                      gcatgggAAACTGCCACGGACGGGATGATGATGGTCGCTTCGGCGATAGCCTGCGCCGTCAGGCGATTCTGGATAACGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCGATAAGGTGTGAACCAGACCCAACGGGCCGGCACCTAAGGCCACACGTTCCTTACTCATCGTTCATCACCCGTTCCAGGATAGGTGACg                                                                                                  >  2:579274/5‑239 (MQ=255)
                                                                                                                                                     aCGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCCTAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACCTAAGGCAATCAGTTTCTTACTCCTCGTTAATTCCTCGTCCAGGGATATTTTACGGGGGAAATTATGacc                                                                                   >  2:877624/1‑237 (MQ=255)
                                                                                                                                                           gtgatgaggatGTTCTCTTCGGCGATATCCTGCGCCGTCAGGCTATTCTGGATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCAACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCa                                                                              <  1:841358/236‑1 (MQ=255)
                                                                                                                                                                      ggTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCCGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCCGTGATAGTTGAAGAGGGAAATTATGAAAGGGCAAAACAAATCGcc                                                                  >  2:577374/1‑239 (MQ=255)
                                                                                                                                                                       gTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATcccct                                                                 >  1:722158/1‑238 (MQ=255)
                                                                                                                                                                        tCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGt                                                                >  1:1014440/1‑239 (MQ=255)
                                                                                                                                                                             tCGGGGATACCCTGCGCGGTCGGGCGTTTCGGGATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGTTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                                                           <  1:545755/239‑1 (MQ=255)
                                                                                                                                                                              ggggtttttccggctccggcAGGCGATCTGGTATCCCGCCCCCTCCCTGGTTTTAATGTTTCCCTAAAAAACCGCCTTCTCCCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAg                                                          <  2:70129/231‑1 (MQ=255)
                                                                                                                                                                                                  cggcttttCTGGTTTCCGTCCGCGCCCTGGGTTTAAATTTTCATTTCATACCCGGTTGCTACCGCCGCTTCTGCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTgg                                      <  2:180679/233‑1 (MQ=255)
                                                                                                                                                                                                       aTTCTGGTACCCGGCCCCGCCCGGGGTTACAATTTCCTCTTCAACGCGGCCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAAt                                 <  1:503285/239‑1 (MQ=255)
                                                                                                                                                                                                              gcctccGTCCGCGCCGGGGGTTCAAATTTTCCCTTATAAACCGGCTTCTCCCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGGGTGAGCCAGGCCCATCGTGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATcc                           <  2:1021030/235‑1 (MQ=255)
                                                                                                                                                                                                              gagcccGGCCGCGCCTGGGGTTACAATTTCCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCAGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATcc                           <  1:877624/235‑1 (MQ=255)
                                                                                                                                                                                                                tcccccccGCCCCGTGGGTTTAATTTTTCCCTATAAACCGGGTTTTTCCCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGGGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCg                         <  2:207262/233‑1 (MQ=255)
                                                                                                                                                                                                                ctccGTCCGCGCCCGGGGTTCCATTTTCCCCTTCATACCCGTCTTCTCCCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCg                         <  1:564067/237‑1 (MQ=255)
                                                                                                                                                                                                                ccagttCCCCGCCCTGGGTTCCAATTTGCGCTTCATCCCCGTCTTCTACCGTCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTTAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATcc                           <  1:519871/232‑1 (MQ=255)
                                                                                                                                                                                                                atcccccccgccccTGGGGTTTAAATTTCACCTCAAAACCCGCTTTTACCCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTTGGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATcc                           <  2:35346/227‑1 (MQ=255)
                                                                                                                                                                                                                 ccggTCCGCGCCCGGGGTTTCAATTTTCCTTTCATACCGGGTTACACCCGCCTCTGCTTCCAGCGCCTGAGCGGCCATAAAGGTGTTAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCg                         <  1:284291/235‑1 (MQ=255)
                                                                                                                                                                                                                 cccGGCCGCCCCTGGGGTTACATTTTTCCATTAAGAAGCGGTTTATCCCGCCCCTTCTTCCGGGGCCTGGGGGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTAGTTCAGTGATAGTTGACGAGGGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCg                         <  2:407056/237‑1 (MQ=255)
                                                                                                                                                                                                                   cGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACCTAAGGCAATCAGCTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGTACTTATGAACGTGCAAAACAATTCGCCGTTACCAGCCCAATAATGCATTTACTGCAGAATTAATTCCGCtt                      >  1:508680/1‑238 (MQ=255)
                                                                                                                                                                                                                   cGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGAATTTACTGGAGAATTAATccc                          >  1:692373/1‑235 (MQ=255)
                                                                                                                                                                                                                     tCCGCCCCCTGGGTTTCAATTTTCCCTTCAAAACCGGCTTCTCCCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTa                    <  1:880968/239‑1 (MQ=255)
                                                                                                                                                                                                                     ccccggcCCTGGTTTTCATTTTTCCCTTCATACCGGGTTTCACCCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTa                    <  2:123825/234‑1 (MQ=255)
                                                                                                                                                                                                                        gcccccTGGGTTTAAATTTTCCCTTCAGACCCGTCTTCCACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACgg                 <  1:242432/236‑1 (MQ=255)
                                                                                                                                                                                                                               gggTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGACGTTAACGGACAAAAAATGCATTTACTGGAAAATTAATCCCGCTGTACGGacggtc           >  2:519871/1‑233 (MQ=255)
                                                                                                                                                                                                                                ggTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCGGTACGGAAggggga         >  2:762121/1‑234 (MQ=255)
                                                                                                                                                                                                                                   ttCCATTTTTCCCTCAAAACCGGTATTCCACCCCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGGGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACg       <  2:44817/239‑1 (MQ=255)
                                                                                                                                                                                                                                       aTTTTCACTTCATACCCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGt  <  1:577374/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
TTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGT  >  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: