Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,177,346	G→A	intergenic (‑142/+166)	gatY ← / ← fbaB	tagatose‑1,6‑bisphosphate aldolase 2/fructose‑bisphosphate aldolase class I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,177,346	0	G	A	95.8%	75.6 / NA	24	intergenic (‑142/+166)	gatY/fbaB	tagatose‑1,6‑bisphosphate aldolase 2/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/8);  minor base T (0/1);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTC  >  NC_000913/2177102‑2177590                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     cgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAg                                                                                                                                                                                                                                                <  1:58666‑M2/110‑1 (MQ=255)     gatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCa                                                                                                                                                                                                                                              <  1:63208‑M2/112‑1 (MQ=255)         tccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTgc                                                                                                                                                                                                                                        >  1:130047‑M2/134‑251 (MQ=255)                            ttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                                                                                                                                                                                     >  1:226005‑M2/115‑251 (MQ=255)                             tacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                                                                                                                                                                                     <  2:146367‑M2/137‑1 (MQ=255)                                                                catctgatagaaacagaagcccctggagcatctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAATCGACAATATGAAAATTATTCGAGGAGTGAAAGTCAAAAAAACGGCCTCCCGATAGGGAATCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCAt                                                                                                                                                                                 <  1:262033‑M2/173‑1 (MQ=255)                                                                                                             cagacactaaatcagtaagttggcagcatcaccATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCa                                                                                                                                    <  1:130191‑M2/218‑1 (MQ=255)                                                                                                                                             cATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAATTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGACACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCctacggggaggg                                                                                                    >  2:154638‑M2/2‑247 (MQ=255)                                                                                                                                              aTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGAGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGgtggtg                                                                                                   >  1:262332/1‑251 (MQ=255)                                                                                                                                                      attTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGGGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCGTATCCGTCCTACGGGGTGGTGTTGTTCGt                                                                                           >  2:193822/1‑251 (MQ=255)                                                                                                                                                           tttatAAATATCGAGATAACGCTCACAACAACGACAATATGATTATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGAGCGTGTTTTTATCCCGGATGCGGGGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAGCGCCTTATCCGTCCTCCGGGGGGGTGCTGTACGtata                                                                                        >  2:56652/5‑249 (MQ=255)                                                                                                                                                                   tatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATc                                                                               >  2:89880/1‑250 (MQ=255)                                                                                                                                                                          aTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGAt                                                                        >  2:167578/1‑250 (MQ=255)                                                                                                                                                                                 tCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAAttt                                                                >  1:6753/1‑251 (MQ=255)                                                                                                                                                                                    caAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGc                                                             <  2:80361/251‑1 (MQ=255)                                                                                                                                                                                                 atGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCGCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATaaa                                                 <  2:50480/250‑1 (MQ=255)                                                                                                                                                                                                                   gTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAACTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGAAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATCAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAAt                              >  1:117127/1‑251 (MQ=255)                                                                                                                                                                                                                      aaaGGCAAAAAAACGGCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAAta                             >  1:37701/1‑248 (MQ=255)                                                                                                                                                                                                                           cAAAAAAACGTCCTCACGATAGGGAAAACGTAGCAAAGTGCGTGTGTTTTTATGCCGGAGGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCTGCGTCGCAGATGGCATAGGTGGCCGGATGCAGCGTAAACGCCTTAACCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCa                      <  1:223384/251‑1 (MQ=255)                                                                                                                                                                                                                                   cGGCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCAGGCTGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAAGCTTCAGCGTCGCATCAGGCATCGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGCGATAAAAGTCCTGCACGGCGTTAATCAGTTTCACGCCGTc               >  1:235366/1‑250 (MQ=255)                                                                                                                                                                                                                                    ggCCTCCCGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAg             <  1:108631/251‑1 (MQ=255)                                                                                                                                                                                                                                           cGATAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCg       >  2:210897/1‑250 (MQ=255)                                                                                                                                                                                                                                             aTAGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGAtt    >  1:227410/1‑251 (MQ=255)                                                                                                                                                                                                                                               aGGGAAACCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTc  >  1:8612/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     CCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTC  >  NC_000913/2177102‑2177590

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