Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	205,776	C→T	V217V (GTC→GTT)	dnaE →	DNA polymerase III subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	205,776	0	C	T	87.5%	14.9 / NA	8	V217V (GTC→GTT)	dnaE	DNA polymerase III subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (4/3);  minor base G (0/1);  total (4/4)
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.65e-01

CAGCGCGCTGGTAGATGAGTGTGTCGCGTTTTATGAAGAACACTTCCCGGATCGCTATTTTCTCGAGCTGATCCGCACCGGCAGGCCGGATGAAGAAAGCTATCTGCACGCGGCGGTGGAACTGGCGGAAGCGCGCGGTTTGCCCGTCGTGGCGACCAACGACGTGCGCTTTATCGACAGCAGCGACTTTGACGCACACGAAATCCGCGTCGCGATCCACGACGGCTTTACCCTCGACGATCCTAAACGCCCGCGTAACTATTCGCCGCAGCAATATATGCGTAGCGAAGAGGAGATGTGTGAGCTGTTTGCCGACATCCCCGAAGCCCTTGCCAACACCGTTGAGATCGCCAAACGCTGTAACGTAACCGTGCGTCTTGGTGAATACTTCCTGCCGC  >  NC_000913/205566‑205963                                                                                                                                                                                                                   |                                                                                                                                                                                            cAGCGAGCTGGTAGATGAGTGTGTCGCGTTTTATGAAGAACACTTCCCGGATCGCTATTTTCTCGAGCTGATCCGCACCGGCAGGCCGGATGAAGAAAGCTATCTGCACGCGGCGGTGGAACTGGCGGAAGCGCGCGGTTTGCCCGTCGTGGCGACCAACGACGTGCGCTTTATCGACAGCAGCGACTTTGACGCACACGAAATCCGCGTTGCGATCCACGACGGCTTTACCCTCGACGATCCTAAACGcc                                                                                                                                                     >  1:129477/1‑251 (MQ=255)                    ggtgtCGGGTTTTATGAAGAACACTTCCCGGATCGCTATTTTCTCGAGCTGATCCTCACCGACAGGCCGGATGAAGAAAGCTAGCTGCACGCGGCGGTGGAACTGGCGGAAGCGCGCGGTTTGCCCGTCGTGGCGTCCAACGACGTGCGCTTTATCGAAAGCAGCGACTTTTACTCACACGAAATCCGCGTGGCGATCCACGACGGCTTTGCCCTCGACTGTCCGAAACGCCCGCGTATCTATTCGCCgca                                                                                                                                  <  1:140518/250‑1 (MQ=255)                            ggtttATGAAGAACACTTCCCGGATCGCTATTTTCTCGAGCAGATCCGCAACGGAAGTCCGGATGACGAACGCAATATGCACGCGACGGAGGAACTGGCGGACGCGCGCGGTTTGCAAGTAGTGGCGACCAACGACGAGCGCTTTATCGAAAGCAGAGCCTTTGAAGCACACGACATCCGCGTTGCGATCCACGACGGCTTTACCCTCGCCGATCCTaaaagc                                                                                                                                                      >  1:162466/3‑220 (MQ=255)                                                                cGAGCTGAGCCGCACCGGCAGGCCGGATGAAGAAAGCTATCTGCACGCGGCGGTGGAACTGGCGGAAGCGCGCGGTTTGCCCGTCGTGGCGACCAACGACGTGCGCTTTATCGACAGCAGCGACTTTGACGCACACGAAATCCGCGTTGCGATCCACGACGGCTTTACCCTCGACGATCCTAAACGCCCGCGTAACTATTCGCCGCAGCAATATATGCGTAGCGAAGAGGAGATGTGTGAGCTGTTTGccn                                                                                      <  1:41501/251‑2 (MQ=255)                                                                                   aGGCCGGATGAAGAAAGCTATCTGCACGCGGCGGTGGAACTGGCGGAAGCGCGCGGTTTGCCCGTCGTGGCGACCAACGACGTGCGCTTTATCGACAGCAGCGACTTTGACGCACACGAAATCCGCGTTGCGATCCACGACGGCTTTACCCTCGACGATCCTAAACGCCCGCGTAACTATTCGCCGCAGCAATATATGCGTAGCGAAGAGGAGATGTGTGAGCTGTTTGCCGACATCCCCGAAGCCCTTgn                                                                   <  1:19884/251‑2 (MQ=255)                                                                                        ggATGAAGAAAGCTATCTGCACGCGGCGGTGGAACTGGCGGAAGCGCGCGGTTTGCCCGTCGTGGCGACCAACGACGTGCGCTTTATCGACAGCAGCGACTTTGACGCACACGAAATCCGCGTTGCGATCCACGACGGCTTTACCCTCGACGATCCTAAACGCCCGCGTAACTATTCGCCGCAGCAATATATGCGTAGCGAAGAGGAGATGTGTGAGc                                                                                               >  1:77351/1‑218 (MQ=255)                                                                                                            aCGCGGCGGTGGCACTGGCGGAAGCGCGCGGTTTTCCCGGCGTGGCGACCAACGCAGTGCGCTTTATCGACGGCAGCGACTTTGACGCACACGAAATCCGAGTTGCGGTCCACGACGGCTTTAACCACGACGGTCCTAAACGCCCGCGTAACTATTCGACGCCGAAATATATGCGCAGCGAAGAGGAGATGTGTGAGCTGTTTGCCGACATCCCCGACGCCCTTGCCAACCCCGTTGAGATCACCCAACGc                                          >  1:23175/1‑251 (MQ=255)                                                                                                                                                    cgtGGCGACCAACGTCGTGCGCTTTATCGACAGCAGCGACTTTGACGCACACGAAATCCGCGTTGCGATCCACGACGGCTTTACCCTCGACGATCCTAAACGCCCGCGTAACTATTCGCCGCAGCAATATATGCGTAGCGAAGAGGAGATGTGTGAGCTGTTTGCCGACATCCCCGAAGCCCTGGCCAACACCGTTGAGAGCCCCAAACGCTGTAACGTAACCGTGCGTCTTGGTGAATACTTCCTGCCGc  <  1:27845/251‑1 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                            CAGCGCGCTGGTAGATGAGTGTGTCGCGTTTTATGAAGAACACTTCCCGGATCGCTATTTTCTCGAGCTGATCCGCACCGGCAGGCCGGATGAAGAAAGCTATCTGCACGCGGCGGTGGAACTGGCGGAAGCGCGCGGTTTGCCCGTCGTGGCGACCAACGACGTGCGCTTTATCGACAGCAGCGACTTTGACGCACACGAAATCCGCGTCGCGATCCACGACGGCTTTACCCTCGACGATCCTAAACGCCCGCGTAACTATTCGCCGCAGCAATATATGCGTAGCGAAGAGGAGATGTGTGAGCTGTTTGCCGACATCCCCGAAGCCCTTGCCAACACCGTTGAGATCGCCAAACGCTGTAACGTAACCGTGCGTCTTGGTGAATACTTCCTGCCGC  >  NC_000913/205566‑205963

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