Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,745,039	G→T	25.0%	D191Y (GAC→TAC)	gltA →	citrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,745,039	0	G	T	25.0%	36.6 / 4.7	16	D191Y (GAC→TAC)	gltA	citrate synthase
Reads supporting (aligned to +/- strand):  ref base G (8/4);  new base T (2/2);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTGGATATCAATAACCCGCAACACCGCGAAATCTCCGCGGTGCGCCTGGTCGCCAAGATGCCGACCCTGGCTGCGATGGTTTACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGACCGGATCTTCATCCTCCACGCCGACCACGAGCAGAACGCCT  >  NC_002947/4744914‑4745171                                                                                                                              |                                                                                                                                     cccTGGATATCAATAACCCGCAACACCGCGAAATCTCCGCGGTGCGCCTGGTCGCCAAGATGCCGACCCTGGCTGCGATGGTTTACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTa                                                                                                                            >  1:18088/1‑136 (MQ=255)     ggATATCAATAACCCGCAACACCGCGAAATCTCCGCGGTGCGCCTGGTCGCCAAGATGCCGACCCTGGCTGCGATGGTTTACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACTACCTGTCGTAcgcg                                                                                                                        >  2:189916/1‑136 (MQ=255)     ggATATCAATAACCCGCAACACCGCGAAATCTCCGCGGTGCGCCTGGTCGCCAAGATGCCGACCCTGGCTGCGATGGTTTACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACTACCTGTCGTAcgcg                                                                                                                        >  2:289698/1‑136 (MQ=255)                        cacCGCGAAATCTCCGCGGTGCGCCTGGTCGCCAAGATGCCGACCCTGGCTGCGATGGTTTACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACatga                                                                                                     >  1:112793/1‑136 (MQ=255)                                     ccGCGGTGCGCCTGGTCGCCAAGATGCCGACCCTGGCTGCGATGGTTTACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACAcccc                                                                                        <  1:176497/136‑1 (MQ=255)                                              gcCTGGTCGCCAAGATGCCGACCCTGGCTGCGATGGTTTACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGAt                                                                               <  2:150117/136‑1 (MQ=255)                                                                          tGCGATGGTTTACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGTTCAGCCCGGTGCTGGCCAAg                                                   <  1:347840/136‑1 (MQ=255)                                                                                   ttACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACTACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGAc                                          <  1:189916/136‑1 (MQ=255)                                                                                   ttACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACTACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGAc                                          <  1:289698/136‑1 (MQ=255)                                                                                     aCAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGACCg                                        >  1:72375/1‑136 (MQ=255)                                                                                               cATGGGCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGACCGGATCTTCAtc                              >  2:316427/1‑136 (MQ=255)                                                                                                   ggCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGACCGGATCTTCAtcctcc                          <  2:18088/136‑1 (MQ=255)                                                                                                    gCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGACCGGATCTTCATCCTCCa                         >  1:307545/1‑136 (MQ=255)                                                                                                    gCCAGCAGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGACCGGATCTTCATCCTCCa                         >  1:18060/1‑136 (MQ=255)                                                                                                            atgatgTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGACCGGATCTTCATCCTCCACGCCGAcc                 >  1:53025/1‑136 (MQ=255)                                                                                                                           aaCGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGACCGGATCTTCATCCTCCACGCCGACCACGAGCAGAACGcct  >  1:299854/1‑136 (MQ=255)                                                                                                                              |                                                                                                                                     CCCTGGATATCAATAACCCGCAACACCGCGAAATCTCCGCGGTGCGCCTGGTCGCCAAGATGCCGACCCTGGCTGCGATGGTTTACAAGTACTCCATGGGCCAGCCGATGATGTACCCGCGCAACGACCTGTCGTACGCGGAAAACTTCCTGCACATGATGTTCAACACCCCGTGCGAGATCAAACCGATCAGCCCGGTGCTGGCCAAGGCAATGGACCGGATCTTCATCCTCCACGCCGACCACGAGCAGAACGCCT  >  NC_002947/4744914‑4745171

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