Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,414,829 A→G E28G (GAA→GGA)  pta → phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,414,8290AG93.6% 151.3 / 0.2 48E28G (GAA→GGA) ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (22/22);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGG  >  NC_000913/2414586‑2415072
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
aacaaAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGAAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAAACAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCaa                                                                                                                                                                                                                                               >  1:497596/1‑250 (MQ=255)
 aaaaaGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAg                                                                                                                                                                                                                                             <  1:369339/249‑1 (MQ=255)
              ccGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAAAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTTATCCGGGCAATGGGACGGAAAGGCGTTCCTCTg                                                                                                                                                                                                                                  >  1:169942/1‑249 (MQ=255)
                              aTCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTAt                                                                                                                                                                                                                 <  1:346331/250‑1 (MQ=255)
                                                             aaGACGCGGGCCGCCTGACTGCCGGATTTCCCACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCAGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGAGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCg                                                                                                                                                                                 <  1:186623/251‑1 (MQ=255)
                                                                 cgcgAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAAACGTGTCCCGGATTATTATGCTGATCCCTAACGGAAACAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTTGGCTGAGCGTTTTTAAACCTATCGCTCAGCCGCGGACCGGTGGCGATACGCCCgctc                                                                                                                                                                              >  1:173470/1‑247 (MQ=255)
                                                                             gACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGGAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTTGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTTAGCGTTTTCAAACCTAACGCTCAGCCGCGTACCCGTGGCGGTGCGCCCGATCAAGCTACGAatat                                                                                                                                                                 >  1:523977/1‑251 (MQ=255)
                                                                               ctgcctgcTTTCACACCGCAAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTCACCAACGAGGGTAAAACCTGTCCCGTATTATTCTGCTGATACCTACCGGGAACAGCCTCGGGCTGAACACAGTCCACCTTGGCGGTATCCGTGGAATGGGAAGCACAAGCGTTCGTCTGCGCGATTTCAAACCTCTCGCTCCGCCGCGCACCGGGGGGGATGGGCCCGaaacggctccacccat                                                                                                                                                                 >  2:69803/1‑234 (MQ=255)
                                                                                             cccGCCAGCGCAGGTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGGATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGATctct                                                                                                                                                 <  2:524362/250‑1 (MQ=255)
                                                                                             acCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGTCCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCTTGCGTGCGAActct                                                                                                                                                 <  1:408847/251‑1 (MQ=255)
                                                                                                 ccAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAACCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTcc                                                                                                                                              >  1:229892/1‑250 (MQ=255)
                                                                                                      ccagcgGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCacg                                                                                                                                        <  2:92996/250‑1 (MQ=255)
                                                                                                        agctGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTAGTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCacgac                                                                                                                                      <  1:379499/251‑1 (MQ=255)
                                                                                                                         tttAACCCTCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCTGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGGCTGAGCGTTTTCAAACCTCTCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTGCCACCACGACGGccgctgaaccgctg                                                                                                                      <  2:323815/248‑1 (MQ=255)
                                                                                                                                 gCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGc                                                                                                             <  2:32711/251‑1 (MQ=255)
                                                                                                                                  ccAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGAGTGCGAACTCTTCCACCACGAAGGCCGCTGAACCGCTGAAAATGAGcc                                                                                                            >  2:19311/1‑250 (MQ=255)
                                                                                                                                    aaaTCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTAAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTAc                                                                                                          <  1:56712/251‑1 (MQ=255)
                                                                                                                                         ggcggTAACGAAAGAGGGTAAACCGTGTCCCGTTTTATAATGCTGATCGCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCGTGGCGTGATCCGTGCAATGGGACGCAAAGGCGGTCGTCGGAGCGTTTGCAAACCTATCGCTCAGCCGCGGACCGGTGGCGATGCGCCCGAGCAGACGACGACTAGCGTGCGTGCGAACTCGTCCACCACGGCGGCCGCTGAACCGCTGAAAATGAGCGACGTTGa                                                                                                     <  2:243098/251‑1 (MQ=255)
                                                                                                                                          gcggTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGa                                                                                                     >  1:204129/1‑250 (MQ=255)
                                                                                                                                              tAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGGCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTAc                                                                                                          <  1:397736/241‑1 (MQ=255)
                                                                                                                                              tAACGAAAGAGGATAAACAGTGTCCCGTATTATTATGCTGATCCCTACCGGACCCAGCGACGGTCTGACCAGCTTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGACGAGGCGCCCGATCAGACTCCGACTATCGTGCGTGCGCACTCTTCCACCCCGACGGCCGCTGAACCGCTGAAACTGAGCtcccgatcggaa                                                                                                 >  2:397736/1‑239 (MQ=255)
                                                                                                                                               aaCGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATc                                                                                                                                                                <  1:76691/186‑1 (MQ=255)
                                                                                                                                               aaCGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATc                                                                                                                                                                >  2:76691/1‑186 (MQ=255)
                                                                                                                                                aCGAAAGAGTATAAACCTGTGCCCGTCTTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACTAGCTTCAGCCTTGGCGTGATCCGTTCACTGGTACGTAAATGCGGTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGTCCGAGCAGGCTACGGCTATCGTGCGTGCGAGTTCTTGCGCCACGACGGCCGCTGAACCGCTGAAAATTAGCTACGTTGAAGGTCTg                                                                                              <  1:231067/251‑1 (MQ=255)
                                                                                                                                                 cGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCACTGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGc                                                                                             >  1:411391/1‑251 (MQ=255)
                                                                                                                                                  gAAAGGGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACGCTTCCACCGCGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCt                                                                                            <  2:278516/251‑1 (MQ=255)
                                                                                                                                                    aaGAGGTTAAACCGTGTCCCGGATTATTATGCTGGTCCCTACCGGAACCAGCGTGGGTCTGACCAGCTGCAGCCTTGGTGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGCACCGCTGAAAATGAGCTACGTTGAAGGTCTGCttg                                                                                          <  2:277187/251‑2 (MQ=255)
                                                                                                                                                       agGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGAGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCa                                                                                       <  1:103844/251‑1 (MQ=255)
                                                                                                                                                             aaCCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATc                                                                                 <  1:457042/251‑1 (MQ=255)
                                                                                                                                                              aCCGTGTCCCGTATTATTATGCTGACCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCa                                                                                <  2:173470/251‑1 (MQ=255)
                                                                                                                                                                  tgtCCCGTATTATTATGCTGATCCCTACAGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGaa                                                                             >  2:168444/1‑250 (MQ=255)
                                                                                                                                                                  tgtCCAGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGaa                                                                             >  1:203733/1‑250 (MQ=255)
                                                                                                                                                                   gtCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGCTGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGaa                                                                             >  1:336678/1‑249 (MQ=255)
                                                                                                                                                                       ccGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGAtg                                                                        <  1:271818/251‑1 (MQ=255)
                                                                                                                                                                           attattatGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCt                                                                    <  1:30847/251‑1 (MQ=255)
                                                                                                                                                                                tatGCTGAACCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATg                                                                <  1:489240/250‑1 (MQ=255)
                                                                                                                                                                                        tCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCCGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCATTCAGAAAGATGTGCGGATGGAag                                                            >  2:495894/1‑246 (MQ=255)
                                                                                                                                                                                             aCCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGAtcgtcg                                                   <  1:441062/250‑1 (MQ=255)
                                                                                                                                                                                                tgAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCaa                                                >  1:490642/2‑250 (MQ=255)
                                                                                                                                                                                                      aGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTAcc                                          >  2:247089/1‑250 (MQ=255)
                                                                                                                                                                                                                tGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACAcc                               >  1:73223/1‑251 (MQ=255)
                                                                                                                                                                                                                    catggTCAGCCTGGGCTTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCGGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAg                           <  2:447927/247‑1 (MQ=255)
                                                                                                                                                                                                                     aGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCaa                             >  2:25709/1‑248 (MQ=255)
                                                                                                                                                                                                                            gCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGa                    >  1:10115/1‑250 (MQ=255)
                                                                                                                                                                                                                                tGGCGTGATCCGTGCAATGGGCCGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCg               >  1:333866/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cgATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTg          >  2:123374/2‑251 (MQ=255)
                                                                                                                                                                                                                                         ccGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTg      <  1:12525/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCAAGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAgg  >  1:441257/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
AACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGG  >  NC_000913/2414586‑2415072

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