Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,492,583 (G)6→7 intergenic (‑58/+4) idnT ← / ← idnO L‑idonate/5‑ketogluconate/gluconate transporter/5‑keto‑D‑gluconate 5‑reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,492,5831.G90.2% 112.3 / 3.0 41intergenic (‑58/+4)idnT/idnOL‑idonate/5‑ketogluconate/gluconate transporter/5‑keto‑D‑gluconate 5‑reductase
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base G (21/16);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGG‑‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGGCGCGCCAGCTCG  >  NC_000913/4492343‑4492830
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
cAGCCCGCCGAGCGTGCTGCCAGTACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGCCGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCGCTTATGAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAAc                                                                                                                                                                                                                                                 <  2:204694/251‑1 (MQ=255)
    ccgccgAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGcc                                                                                                                                                                                                                                             >  1:149714/1‑251 (MQ=255)
          ggCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCACCCCTTCGGCAAATCCGACGACGGCAGCTACAAGAACGAGGGAAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAAc                                                                                                                                                                                                                                       <  2:437508/250‑1 (MQ=255)
                     gATACTATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGCCGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGAAAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATc                                                                                                                                                                                                                            <  1:2053/251‑1 (MQ=255)
                      aTACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTGCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATca                                                                                                                                                                                                                           <  1:296241/251‑1 (MQ=255)
                        cccATTTTGTAAGGAGGGCAGGACGGCCTGGGCACCCATCCCTTCGGCAAATCCGACGACGGCAGCTCCCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATcaac                                                                                                                                                                                                                         <  2:140841/250‑1 (MQ=255)
                                  ataGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGAACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAgg                                                                                                                                                                                                               >  1:325310/1‑251 (MQ=255)
                                   taGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGt                                                                                                                                                                                                              >  1:59093/1‑251 (MQ=255)
                                      agggCAGGACGGCTGGTGCATCCACCCCCTCGGAAAATCCGACGACGGCAGCTACCAAAACGAGGGCAATAAAGCCTTTAATTTTAAAGACGTTCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGAGTAATGGCAGAGTCTTTATCTCTTATTAGGGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGc                                                                                                                                                                                                           <  2:534747/248‑1 (MQ=255)
                                             gACGGCCTGTGCACCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTAc                                                                                                                                                                                                    <  1:240076/251‑1 (MQ=255)
                                                ggCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAAAAGGTGGCCGTTTACGa                                                                                                                                                                                                  >  1:591587/1‑250 (MQ=255)
                                                                        aaaTCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCCGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGATTTCGCGCTGCACGGAGGGGGGGCTGTAAAACAGCCACTAACATGCCTCCATCAACAAACAACAGTTGGCCGTTTACGAAATCAGAGGCTTTTGAAAAAAGGa                                                                                                                                                                          >  1:1320/1‑250 (MQ=255)
                                                                                    ggCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAcc                                                                                                                                                             >  2:162806/1‑251 (MQ=255)
                                                                                     gCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAcc                                                                                                                                                             >  1:16498/1‑250 (MQ=255)
                                                                                      cAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCAGCCGCAATAATGATTAATGGCATAGTCCTCATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCAACAACAAACAACCGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcagcaaca                                                                                                                                                            >  1:380126/1‑247 (MQ=255)
                                                                                                         aaTAAAGCCGTTAACTTTTAAGACGATCATCAGGATAAGAAGCAACGCGACGCCTGCCGCAATAATTATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGCGGCCGTTTACGAAATCAGAGTCTTTTGCAGAAAGGAACACCGCATCCCCAATCAGTTCCTCCGGATCTcc                                                                                                                                        >  2:73470/1‑251 (MQ=255)
                                                                                                            aaaGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTcccc                                                                                                                                      >  1:66814/1‑250 (MQ=255)
                                                                                                             aaGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGAGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTcccc                                                                                                                                      >  2:358413/1‑250 (MQ=255)
                                                                                                                                        aGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAAc                                                                                                         >  1:23590/1‑251 (MQ=255)
                                                                                                                                                  gcagcGCGACGCCTGCCGCAATACTGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTAGCGCTGCACGTAGGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGCTCTCCCCAGCGGGCGGCGGGGGTCCGTTTGCACCACCAGGCGGTg                                                                                                >  2:142589/1‑251 (MQ=255)
                                                                                                                                                         gACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGCGCCAGAAAAGCAGCGATTTCGCGCTGACCGCAGGCGGGG‑‑ATGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGaaagactc                                                                                        >  1:612210/1‑245 (MQ=255)
                                                                                                                                                          aCGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAAAAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACACCCAGGCGTGGAAGGcctct                                                                                       >  1:629797/1‑247 (MQ=255)
                                                                                                                                                                   gCAATAATGATTAAGGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAAcc                                                                              <  2:603644/251‑1 (MQ=255)
                                                                                                                                                                      ataatGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTTAAGGCTTCGTCCTCAACCAGt                                                                           <  1:504023/251‑1 (MQ=255)
                                                                                                                                                                               aaTGGCATAGTCTTTATCTCTTATTAGTTCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAAGCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTc                                                                  <  1:383908/251‑1 (MQ=255)
                                                                                                                                                                               aaTGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGGCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCCACCCGTGCTTTAGTc                                                                  >  1:21206/1‑251 (MQ=255)
                                                                                                                                                                                  ggtaTAGTCTTTTTCGCTTATTAGGTCCAGAAAAGCACGGCTGTCGCGCTCCACGTAGGTGGGG‑‑CTGTTAAACAGCCACGAACATGCCGCCATCAACAAACAACAGGTGGCCGGTTACGAAATCAAAGGCTTTTGAAGAAAGGAACATCTCAGCTCCAAGCAGTTCCTGCGGATCTCCCCAGCGTGCGTCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGGCCTCAACCAGTGCTTTAGTCagt                                                               <  2:562317/248‑3 (MQ=255)
                                                                                                                                                                                  ggCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCAt                                                                <  2:105125/250‑1 (MQ=255)
                                                                                                                                                                                        gTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCGGCCCGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAAGCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGtt                                                         <  2:226052/251‑1 (MQ=255)
                                                                                                                                                                                                  tttattaGTGCCAGAAAAGCAGCTCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGaaa                                                   <  1:603569/246‑1 (MQ=255)
                                                                                                                                                                                                   ttattaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCCGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTTAAATAGcc                                             >  1:91158/1‑251 (MQ=255)
                                                                                                                                                                                                    tattaGTGCCAGAACAGCAGCGCTTTCGCGCTGCAAGTAGGGGGGGCTGTTAAAAAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCCGGCGGTGACGGCTTCTTACTCAAACAGTGCTTTTATCCTTTTTGTTTTGAAATAGccc                                            >  2:257529/1‑251 (MQ=255)
                                                                                                                                                                                                             cAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGgcgccat                                    >  1:420161/1‑247 (MQ=255)
                                                                                                                                                                                                                 aaaGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACCTGCCGCCATCAACAAACAACAGGTGGCCGATTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAAACCGTGCTTTTGGCCTTTCTGTTTTTAAATAGCCCGGCGCAATAccct                               >  1:230073/1‑249 (MQ=255)
                                                                                                                                                                                                                  aaGCAGCGCTTTCGCGCTTGACGTAGGGGGGGCTGTTAAACAGCCACTAACATGACGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTGGAAGAAAGGAACACCGCAGCACCAAGCAGGTTCGGCTGATCTCCCTATCGTGCGGCGGGGTTCCGTTTGCACGACCAGGCGGTTAAGGCTTCGTCCTCAACCAGTGCTGTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGtt                              <  1:359733/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaGCAGCGCTTTCGCGCTGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGtt                              <  2:361306/251‑1 (MQ=255)
                                                                                                                                                                                                                     caccgctttggcccgtcACGTAGGGGGG‑‑GGTGTTAAACAGCCATTAACATCGCGCCATCAACAAACAACAGGTGGCCGTTTACGTAATCAGAGGCGTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATGTCCCCAGCGTGGGGCGGGGGTGGGTTTGCACAACCAGGGGGTGAAGGCTTGGTCCTCAGTCAGTGCTGTTGTCATTTCTGTTGTGAGATAGCCCGGCGCAATGCCTTTGAc                           <  2:380126/235‑1 (MQ=255)
                                                                                                                                                                                                                                cgcTGCACGTAGGGGGGGCTGTTAAACAGCAACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGAGGGGGTCCGGTTGCACACCCAGGCGGGGAAGGCTTCTGCCTCCACCAGTGGCTTAGTCCTTTCCGTTTTTAAATGAGCCAGACCAAAACCGTTGACCTGAattttt                 >  2:98673/1‑246 (MQ=255)
                                                                                                                                                                                                                                agcTGCAGGTATGGGGGG‑‑CTGTTAAGCAGACACTACCATGCCGCCATCAACCAACAACAGGTGGCCATTTACGAAATCAGAGGCATTTGAAGAAAGGAACACCGAAGCAACAATCAGTTCCTGCGGATCTCACCAGCGTGCGGCGGGGGTCCGTTTGCACACCCAGGAAGTGAAGGCTTTGTCCTCAACCATTGCTTTAGTAATTTTTGTTTTGGAATAGGCCGGGACAATACCGTTTAcccgtattattt                >  2:30768/2‑241 (MQ=255)
                                                                                                                                                                                                                                 ccgatcggtaggggggg‑‑CGGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATtgt                <  2:162011/239‑1 (MQ=255)
                                                                                                                                                                                                                                   tGCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCCACCAGTGCTTTAGTCATTTCTGTTTTGGAATAGCCCGGCGCAATACCGTTGACCTTGATATTTTGgc             >  1:62642/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gCACGTAGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGgcg            <  1:187874/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gtaggggggg‑‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCCGTGCTTTAGTCATTTCTGTTTTTAAATAGCCCGGCGCCATACCGTTGACCTGAATATtgtagatcggg        >  1:162011/5‑243 (MQ=255)
                                                                                                                                                                                                                                             ggtggg‑‑CTGTTACACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCTCCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCATTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGTAATACCGTTGACCTGAATATTGTGGCGCGCCAGCt    <  2:62683/248‑1 (MQ=255)
                                                                                                                                                                                                                                               gggg‑‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGGGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTTTGTTTTGAAATAGCCCGGCGCCATACCGGTGAACTGAATATTGTGGCGCGCCAGCTCg  >  1:187102/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
CAGCCCGCCGAGCGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGG‑‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATTGTGGCGCGCCAGCTCG  >  NC_000913/4492343‑4492830

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