Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,018,884	G→A	G942D (GGC→GAC)	ygfK →	putative oxidoreductase, Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,018,884	0	G	A	100.0%	21.4 / NA	12	G942D (GGC→GAC)	ygfK	putative oxidoreductase, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/4);  total (8/4)

GTCAATCCAGCGGAAATCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCTGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCAGAACCGTTTCCAGACGCTGCACCTCGACGCTTACTGTAACGAATGCGGCAACTGCGCTCAGTTCTGTCCGTGGAACGGTAAACCGTACAAAGACAAAATCACCGTCTTCAGCCTGGCGCAAGACTTTGATAACAGCAGCAACCCAGGCTTCCTTGTGGAAGATTGCCGGGTACGAGTACGTCTGAATAACCAAAGCTGGGTGTTAAACATCGACAGCAAAGGTCAGTTTAACAACGTACCACCGGAGCTGAACGATATGTGCCGCATCAT  >  NC_000913/3018652‑3019106                                                                                                                                                                                                                                         |                                                                                                                                                                                                                               ntCAATCCAGCGGAAATCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCTGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCAGAACCGTTTCCAGACGCTGCACCTCGACGCTTACTGTAACGAATGCGACAACTGCGCTCAGTTCTg                                                                                                                                                                                                              >  1:155776/2‑251 (MQ=255)                                                                              tGTCGCCCAGGAAGCCGCTCGCTGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCAGAACCGTTTCCAGACGCTGCACCTCGACGCTTACTGTAACGAATGCGACAACTGCGCTCAGTTCTGTCCGTGGAACGGTAAACCGTACAAAGACAAAATCACCGTCTGCAGCCTGgc                                                                                                                                                           <  1:57828/225‑1 (MQ=255)                                                                              tGTCGCCCAGGAAGCCGCTCGCTGCCNCGAATGNNACNACGTTTGCAGCAAGTGTGTNGANNNCNNNCNGAACCNCGNNAACGTCTCCNTTGCGGTCCCAGGCTTCCAGAACCGTTTCCAGACGCTGCACCTCGACGCTTACNGTAACGAATNNGACAACTGCGCTCAGTTCTGTNCGTNNANNGNTAANccn                                                                                                                                                                                           >  2:57828/1‑192 (MQ=255)                                                                                                                              caAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCAGAACCGTTTCCAGACGCTGCACCTCGACGCTTACTGTAACGAATGCGACAACTGCGCTCAGTTCTGTCCGTGGAACGGTAAACCGGACAAAGACAAAATCACCGTCTTCAGCCTGGCGCAAGACTTTGATAACAGCATTAACCCAGGCTTCCTTGTGGAAGATTGCCGGGTACGAGTACGTCTGAATAAcc                                                                                 <  1:147974/251‑1 (MQ=255)                                                                                                                                                cccGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCAGAACCGTTTACAGACGCTGCACCTCGACGCTTACTGTAAAGAATGCGACAAATGCGCTCAGCTCTGTACGTGGACCGGTAAACCGTCCAAAGACAAAATCACCGTCTTCAGCCTGGCGCAAGACTTTGATAACAGCAGCAACCCAGGCTTCCTTGTGGac                                                                                                                >  1:145126/1‑201 (MQ=255)                                                                                                                                                        gcgcCNNCGNNNCANNTGCGGTCCCAGGCTTCCAGCACCNTTTCCAGACGCTGCACCTAGACGCTTACTNTNANGCATNNGACNACTGCGCTCANTTCTGTCCGTg                                                                                                                                                                                                        >  2:190472/1‑106 (MQ=255)                                                                                                                                                        gcgcCAACGTCTCCATTGCGGTCCCAGGCTGCCAGAACCGTTTCCAGACGCTGCACCTCGACGCTTACTGTAACGAATGCGACAACTGCGCTCAGTTCTGTCCGTg                                                                                                                                                                                                        <  1:190472/106‑1 (MQ=255)                                                                                                                                                               cGTCTCCATTGCGGTCCCAGGCTTCCAGAACCGTTTCCAGACGCTGCACCTCGACGCTTACTGTAACGAATGCGACAACTGCGCTCAGTTCTGTCCGTGGAACGGTAAACCGTACAAAGACAAAATCACCGTCTTCAGCCTGGCGCAAGACttt                                                                                                                                                 >  1:143551/1‑154 (MQ=255)                                                                                                                                                                                                             cACCTCGACGCTTACTGTAACGAATGCGACAACTGCGCTCAGTTCTGTCCGTGGAACGGTAAACCGTACAAAGACAAAATCACCGTCTTCAGCCTGGCGCAAGACTTTGATAACAGCAGCAACCCAGGCTTCCTTGTGGAAGATTGCCGGGTACGAGTACGTCTGAATAACCAAAGCTGGGTGTTAAACATCGACAGCAAAGGTCAGTTTAACAACGTACCACCGGAGCTGAACGATATGTGCCGcatcat  >  1:208621/1‑251 (MQ=255)                                                                                                                                                                                                                      gCTTACTGTAACGAATGCGACAACTGCGCTCAGTTCTGTCCGTGGAACGGTAAACCGTACAAAGACAAAATCACCGTCTTCAGCCTGGCGCAAGACTTTGATAACAGCAGCAACCCAGGCTTCCTTGTGGAAGATTGCCggg                                                                                                      >  1:103551/1‑142 (MQ=255)                                                                                                                                                                                                                               naCGAATGCGACAACTGCGCTCAGTTCTGTCCGTGGAACGGTAAACCGTACAAAGACAAAATCACCGTCTTCAGCCTGGCGCAAGACTTTGATAACAGCAGCAACCCAGGCTTCCTTGTGGAAGATTGCCGGGTAc                                                                                                   >  1:150650/2‑136 (MQ=255)                                                                                                                                                                                                                               aaNGAATGCGACAACTGCGCTCAGTTCTGTCCGTGGAACGGTAAACCGTACAAAGACAAAATCACCGTCTTCAGCCTGGCGCAAGACTTTGATAACAGCAGCAACCCAGGCTTCCTTGTGGAAGATTGCCGGGTAc                                                                                                   <  2:150650/136‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                               GTCAATCCAGCGGAAATCTATCAACGTAAAGGCGATATCTCTATCACGCTGGTGAACAGTGACGATCGTGACGCGTTTGTCGCCCAGGAAGCCGCTCGCTGCCTCGAATGTAACTACGTTTGCAGCAAGTGTGTGGATGTCTGCCCGAACCGCGCCAACGTCTCCATTGCGGTCCCAGGCTTCCAGAACCGTTTCCAGACGCTGCACCTCGACGCTTACTGTAACGAATGCGGCAACTGCGCTCAGTTCTGTCCGTGGAACGGTAAACCGTACAAAGACAAAATCACCGTCTTCAGCCTGGCGCAAGACTTTGATAACAGCAGCAACCCAGGCTTCCTTGTGGAAGATTGCCGGGTACGAGTACGTCTGAATAACCAAAGCTGGGTGTTAAACATCGACAGCAAAGGTCAGTTTAACAACGTACCACCGGAGCTGAACGATATGTGCCGCATCAT  >  NC_000913/3018652‑3019106

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