Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,508,134 A→C 21.9% 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,1340AC21.9% 126.4 / 5.9 73V36G (GTG→GGG) fryBputative enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base A (36/17);  new base C (5/12);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 9.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGTTATCTTCCGGGGTAACTGCCACGGAGTGGATGATGATGGTCGCTTCGGCGATATCCTGCGCCGTCAGGCGATTCTGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGT  >  NC_000913/2507898‑2508366
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
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)
                               ggtattcttcttttggtttttttTGATGGTGCCCGAAGGTTTTTTATTGGGTTCCTGAAAGGTGATTCCATAAACGGCGCGTGATTCGAAGCGTTCTTTATCTTCCGGGGTACCTGCCACGGAGGGGATGAGGATGGTCGCTTCGGCTATACCCTGCGCCGTCAGGCGATTCTGGATCCCGTCCGCGCCCTGGGTTTCAATTTTCACTTCATAACCGGCTTCTACCGCCGCTTCTTcc                                                                                                                                                                                                          <  2:148244/221‑1 (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)
                                                      ggttgttcgcgctgggttttttAATTGCGTCCGGAAAAGTGGTTTCATAATGTGCGGGTGATTAGAACGGTCGGTTACCTCCCGGGGTAACTGCCACGGAGGGGATGATGATGGGCGCTTGGGCGTTATCTGGCGCCGTCAGGCGATTCGGGATACCGTCCGCGCCCTGGGTTTCAATTTTCACTTAATAACCGGTTTCTACCGCCGCTTCTTCCAGCGCCTGAGCGGCCATAAAGgtg                                                                                                                                                                                  <  2:828840/222‑1 (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)
                                                                                                                                                                                                               gcccgcccccctgtcttGGGTTTCTGTTTTCACCTCAAACCCGTTTTCAACCGCCGTTTCTTCCCGCGCCTGACCGGCCATAAAGGTGTGAGCCAGGCTCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATcc                           <  2:424720/223‑1 (MQ=255)
                                                                                                                                                                                                               gagcccGGCCGCGCCTGGGGTTACAATTTCCACTTCATAACCGGCTTCTACCGCCGCTTCTTCCAGCGCCAGAGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATcc                           <  1:877624/235‑1 (MQ=255)
                                                                                                                                                                                                                 tccccccccccccttGGTTTTCATTTTTCCCTTCAAACCCGGTTTCACCCGCCGCTTCTTCCAGCGCCTGGGCGGCCATAAAGGTGTGAGCCAGGCCCATCGGGCAGGCACATAAGGCAATCAGTTTCTTACTCATCGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCg                         <  2:1014440/225‑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/2507898‑2508366

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