Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,076,209	C→T	L22L (CTA→TTA)	yihY →	PF03631 family membrane protein YihY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,076,209	0	C	T	100.0%	18.2 / NA	9	L22L (CTA→TTA)	yihY	PF03631 family membrane protein YihY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (3/6);  total (3/6)

ATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAGGCACCGTACCCGTCCACTATGGGCCTGGCTAAAACTACTCTGGCAACGCATTGATGAGGACAACATGACAACCCTGGCAGGTAACCTTGCCTATGTGTCGTTGCTCTCATTAGTGCCGCTGGTTGCCGTTGTTTTTGCGCTTTTCGCCGCTTTTCCCATGTTTTCCGACGTCAGCATTCAGTTACGTCACTTTATTTTTGCCAACTTTCTGCCTGCTACTGGCGATGTTATTCAGCGGTATATCGAACAATTTGTTGCC  >  NC_000913/4075967‑4076433                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 aTGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATGACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGTACTATTTCGCGAAGGGGTTATGCTAAAATCCACTCAGGACAAAGCCAGGCACCGTACCCGTCCACTATGGGCCTGGCTAAAATTACTCTgg                                                                                                                                                                                                                          <  1:118635/251‑1 (MQ=255)                gCACAAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAGGCACCGTACCCGTCCACTATGGGCCTGGCTAAAATTACTCTGGCAACGCATtgatgan                                                                                                                                                                                                           <  1:31553/251‑2 (MQ=255)                                                         aCCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGATGGGCATTACCAGTATTCTGGTGAAAGATAAAAACACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAGGCACCGTACCCGTCCACTATGGGCCTGGCTAAAATTACTCTGGCAACGCATTGATGAGGACAACATGACAACCCTGGCAGGTAACCTTGCCTATGTGTc                                                                                                                                                                  >  1:83638/1‑251 (MQ=255)                                                                          ttttCGACTAGAACGTCGATAAGATAGAATGAGCCAATCATCTGGGCATTACCAGTATTCTGTTTAAAGATAAAACCACCATCCCGGGCTATTTCGCGAAGGTGTTATGCTAAAAACCATTCATGACAAAGCCATGCACCGTACCCGTCCACAATGGGCCTGGCTAAATTTACTCTCGCAACACATTGGTTAGGACAACATGACAACTCTATCAGGTAACCTTGCCGATGTTTCGTTGCTCTCATTAg                                                                                                                                                    <  1:109209/248‑1 (MQ=255)                                                                                                                           tACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAGGCACCGTACCCGTCCACTATGGGCCTGGCTAAAATTACTCTGGCAACGCATTGATGAGGACAACATGACAACCCTGGCAGGTAACCATGCCTATGTGTCGTTGCTCTCATTAGTGCCGCTGGTTGCCGTTGTTTTTGCGCTTTTcg                                                                                                                   <  1:110018/232‑1 (MQ=255)                                                                                                                                                                                                    aGGACAAAGCCAGGCAACGTACCCGTCCACTATGGGCCTGGCTAAAATTACTCTGGCAACGCATTGATGAGGACAACATGACAACCCTGGCAGGTAACCTTGCCTATGTGTCGTTGCTCTCATTAGTGCCGCTGGTTGCCGTTGTTTTTGCGCTTTTCGCCGCTTTTCCCATGTTTTCCGACGTCAGCATTCAGTTACGTCACTTTATTTTTGCCAACTTTCTGCCTGCTACTGGCGATGTTATTCAGCgn                       <  1:19014/251‑2 (MQ=255)                                                                                                                                                                                                          aaGCCAGGCACCGTACCCGTCCACTATGGGCCTGGCTAAAATTACTCTGGCAACGCATTGATGAGGACAACATGACAACCCTGGCAGGTAACCTTGCCTATGAGTCGTTGCTCTCATTAGTGCCGCTGGTTGCCGTTGTTTTTGCGCTTTTCGCCGCTTTTCCCATGTTTTCCGACGTCAGCAGTCAGTGACGTCACGTTATTTTTGCCAACTTTCTGCCTGCTACTGGCGATGTTATTCAGCGGTATATc                 <  1:169338/251‑1 (MQ=255)                                                                                                                                                                                                            gCCAGGCACCGTACCCGTCCACTATGGGCCTGGCTAAAATTACTCTGGCAACGCATTGATGAGGACAACATGACAACCCTGGCAGGTAACCTTGCCTATGTGTCGTTGCTCTCATTAGTGCCGCTGGTTGCCGTTGTTTTTGCGCTTTTCGCCGCTTTTCCCATGTTTTCCGACGTAAGCATTCAGTTACGTCACTTTATTTTTGCCAACTTTCTGCCTGCTACTGGCGATGTTATTCAg                          >  1:61098/1‑240 (MQ=255)                                                                                                                                                                                                                         cccGTCCACTATGGGCCTGGCTAAAATTACTCTGGCAACGCATTGATGAGGACCACATGCCAACCCTGGCAGGTAACCTTGCCTATGTGTCGTTGATCTCATTAGTGCCGCTGGTTGCCGTTGTTTTTGCGCTTTTCGCCGCTTTTCCCATGTTTTCCGACGTCAGCATTCAGTTACGTCACTTTATTTTTGCCAACTTTCTGCCTGCTACTGGCGATGTTATTCAGCGGTCTATCGACCAATTTGTTGcc  >  1:73037/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 ATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAGGCACCGTACCCGTCCACTATGGGCCTGGCTAAAACTACTCTGGCAACGCATTGATGAGGACAACATGACAACCCTGGCAGGTAACCTTGCCTATGTGTCGTTGCTCTCATTAGTGCCGCTGGTTGCCGTTGTTTTTGCGCTTTTCGCCGCTTTTCCCATGTTTTCCGACGTCAGCATTCAGTTACGTCACTTTATTTTTGCCAACTTTCTGCCTGCTACTGGCGATGTTATTCAGCGGTATATCGAACAATTTGTTGCC  >  NC_000913/4075967‑4076433

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