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_000913	2,414,829	0	A	G	96.3%	87.4 / NA	27	E28G (GAA→GGA)	pta	phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (15/11);  total (16/11)
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.55e-01

AAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCT  >  NC_000913/2414593‑2415063                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           aagGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTc                                                                                                                                                                                                                              <  1:3758/251‑1 (MQ=255) aagGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTAACGGAACCAGCGTCGGGCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGttt                                                                                                                                                                                                                              >  1:23350/1‑250 (MQ=255)       cccGTCCGGCGGGGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCTCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGa                                                                                                                                                                                                                        <  2:165972/251‑1 (MQ=255)                               ccaacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGGAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGGCGGTCTGACCCGCGTCCGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGCTCGTCTGAGCGTTTTCCAACCTATCGCTCAGc                                                                                                                                                                                                >  2:61078/1‑251 (MQ=255)                                           tGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCAGGGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTg                                                                                                                                                                                    >  1:123207/1‑251 (MQ=255)                                                                  cgccTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCCGCGGCAACGACAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACCCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGGTGCCCCCGCTCCGACTAcc                                                                                                                                                             >  2:75228/1‑250 (MQ=255)                                                                         cTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAAAGGATAAACCGTGTCCCGAATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCCAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCAATCAGACTACGACTATcg                                                                                                                                                      >  1:178210/1‑251 (MQ=255)                                                                                            cagctcagcgGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATGATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTccac                                                                                                                                   <  2:7770/251‑1 (MQ=255)                                                                                             agctcagctGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTccacc                                                                                                                                  <  1:10850/251‑1 (MQ=255)                                                                                                              tGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGccgctgaacc                                                                                                                 >  1:188447/1‑251 (MQ=255)                                                                                                                  ttGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGGTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTTGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGccgctgaaccgct                                                                                                              <  1:140168/250‑1 (MQ=255)                                                                                                                  ttGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAAACCTGACCCCTATTATTATGCTGATCCCTACCGGAACCCGCGACGGTCTGACCCCCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGccgctgaaccgctg                                                                                                             >  1:134630/1‑251 (MQ=255)                                                                                                                        cccGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGCTCAGACTACGACTATCGTGCGGGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATg                                                                                                       >  2:9636/1‑251 (MQ=255)                                                                                                                             cAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTa                                                                                                  <  2:4955/251‑1 (MQ=255)                                                                                                                               aaTCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACg                                                                                                >  1:60245/1‑251 (MQ=255)                                                                                                                                     cggTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGa                                                                                            >  2:104484/1‑249 (MQ=255)                                                                                                                                            gAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCt                                                                                   <  1:161244/251‑1 (MQ=255)                                                                                                                                              aaGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCttt                                                                                 >  1:221353/1‑251 (MQ=255)                                                                                                                                                    aTAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGCCCCAGCGTCGGTCTGAACAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAAACGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGc                                                                            >  1:22249/1‑250 (MQ=255)                                                                                                                                                                                  cccTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAag                                                   >  1:129966/1‑245 (MQ=255)                                                                                                                                                                                                    cGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGAGCGTCGCAAACTACCACGCt                            <  2:25486/250‑1 (MQ=255)                                                                                                                                                                                                              aGCGTCAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGa                 <  1:81656/251‑1 (MQ=255)                                                                                                                                                                                                                   cAGCCTTGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTg            <  2:170264/251‑1 (MQ=255)                                                                                                                                                                                                                         tGGCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTc       <  1:104484/250‑1 (MQ=255)                                                                                                                                                                                                                          ggCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGt     >  2:74527/1‑251 (MQ=255)                                                                                                                                                                                                                           gCGTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGtt    >  2:207570/1‑251 (MQ=255)                                                                                                                                                                                                                             gTGATCCGTGCAATGGGACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCt  >  1:165808/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           AAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGACTACGACTATCGTGCGTGCGAACTCTTCCACCACGACGGCCGCTGAACCGCTGAAAATGAGCTACGTTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCT  >  NC_000913/2414593‑2415063

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