Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,664,866	C→T	Q121* (CAA→TAA)	dmsD →	redox enzyme maturation protein DmsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,664,866	0	C	T	100.0%	27.5 / NA	12	Q121* (CAA→TAA)	dmsD	redox enzyme maturation protein DmsD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/7);  total (5/7)

CGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTG‑ATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGC  >  NC_000913/1664625‑1665067                                                                                                                                                                                                                                                  |                                                                                                                                                                                                           nggCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCAATAGGGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTAACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGGCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCCGTCTGGCTGGCTCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGCGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGTAAAACg                                                                                                                                                                                                      >  1:64149/2‑248 (MQ=255)            ntCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGCTCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGGGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGTAAAACGAACCGGAAGATCAt                                                                                                                                                                                        >  1:37224/2‑251 (MQ=255)                                 cGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCCCGCCAAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGTAAAACGAACCGGAAGATCATTTTGGTTCActg                                                                                                                                                                            >  1:197507/1‑242 (MQ=255)                                                             tCAGACACAGTGTGAAGAGACGCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGTAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTG‑ATGGCTGCGTGGTTGGCGGAGAATGGTCGc                                                                                                                                       <  1:127796/250‑1 (MQ=255)                                                               aGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGTAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTG‑ATGGCTGCGTGGTTGGCGGAGAATGGTCGCcag                                                                                                                                    <  1:114960/251‑1 (MQ=255)                                                                                           ggCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGGCGTCTCCGCCGCGGGGGGCGGTCTGGCTGGATCGCGACTCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGACAGGCATTCAGTTTGTAATGTAGTAAAACGACCCGGAAGATCATTTTGGTTCACTGCTGTTG‑ATGGCTGCGTGGTTGGCGGAGAATGGTCGCCCGACAGAATGCGCAGAATTGCtggcatggc                                                                                                        >  1:347944/1‑251 (MQ=255)                                                                                             ccTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGTAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTG‑ATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTTCTGGCATGGCAc                                                                                                      <  1:271393/251‑1 (MQ=255)                                                                                                                 ggCCCGTGGGCCCTGTCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCTAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGTAAAACTAACCGGAAGATCATTTTGGTTCACTGCTGTTG‑ATGGCTGCGTGGTTGGCGGATAATGGTCTCCAGACAGAATGCGAATAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACg                                                                                  <  1:265418/251‑1 (MQ=255)                                                                                                                                      ccgccgTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGTAAAACGAACCGGAAGATCAGTTTGGGTCACTGCTGTTG‑ATGTCTGGGGGGTTGGCGGAGAATGGTCGCCAGACATAATGCGAATAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTAt                                                                <  1:196771/248‑1 (MQ=255)                                                                                                                                                                                                 ggCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGTAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTG‑ATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGGGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGgcgc  <  1:255808/251‑1 (MQ=255)                                                                                                                                                                                                                       gagaAAGGCATTCATTTTTAAATGAAGTAAAACGAATCGGAAGATCATGTGGGTTCACTG‑TGTTGGATTGCTGCTTGTGGGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCCCCtttt                                                                                                 <  1:232487/134‑1 (MQ=255)                                                                                                                                                                                                                                              gAAGTAAAACGAACCGGAAGATCATTTTGGTTCGCTGCTGTTG‑ATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTNATGNNNNNAtnnnnnnnnnn                                                      >  2:216903/1‑144 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                           CGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTG‑ATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGC  >  NC_000913/1664625‑1665067

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