Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,110,975	C→T	G51E (GGA→GAA)	yghD ←	putative type II secretion system M‑type protein YghD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,110,975	0	C	T	100.0%	35.0 / NA	13	G51E (GGA→GAA)	yghD	putative type II secretion system M‑type protein YghD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/2);  total (11/2)

CCTGAATATTCTCCCCACGATCGGCTATCCGCCTGATGACCACCGAGTGCGCGGAAGCACTGTCACTGATGACCCGCGATACCTGCGCAGGCGTTGTTGTGGCAGATTTTCGCGCTGCAATAATATCCGGCGCGGCGTTCTTCAG‑TCGGGTGTTCATCGCCACCAGCTGCTGCAATATTGTCTCCTGTTGCTCAATCCGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTGCC  >  NC_000913/3110742‑3111212                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              ccTGAATATTCTCCCCACGATCGGCTATCCGCCTGATGACCACCGAGTGCGCGGACCCACTGTCACTGATGACCCCCGCTCCATGCGCAGGCGTTGATCTGGCAGATTTTCGCGCTGCAATAATATCCGGCGCGGCGTTCTTCAG‑TCGGGAGTTCATCGCCACCAGCTGCTGCAATATTGTCTCCTGTTGCTCAATCCGTTCGCTCAACGGCTGCCAGCTGAGGCCGTAATATTCGACGCTAAACAGGAac                                                                                                                                                                                                                              >  1:88072/1‑251 (MQ=255)                         cTATCCGCCTGATGACCACCGAGTGCGCGGAAGCACTGTCACTGATGACCCGCGATACCTGCGCAGGCGTTGTTGTGGCAGATTTTCGCGCTGCAATAATATCCGGCGCGGCGTTCTTCAG‑TCGGGTGTTCATCGCCACCAGCTGCTGCAATATTGTCTCCTGTTGCTCAATCCGTTCGCTCAACGGCTGCCAGATGAGGCCGTAATATTCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGccc                                                                                                                                                                                                      >  1:458606/1‑251 (MQ=255)                                                                                   ntGCGCAGGCGTTGTTGTGGCAGATTTTCGCGCTGCAATAATATCCGGCACGGCGTTTTTCAGGTGGGGGGTT‑ATCGACACCAGCTGCTGCAATATTGTCTACTGTTGCTCAATCCGTTCGTTCAACGGCTGCCAGGTGAGGACGTAATATTCGACGCTAAACAGGAACACCACCGCTGCCAGGAACATGCCCTTTTCACGCGGAGaa                                                                                                                                                                                       >  1:57944/2‑208 (MQ=255)                                                                                                gttgtGGCAGATTTTCGCGCTGCAATAATATCCGGCCCGGCGCTCTTCAG‑TCGGGTGTTCATCGCCACCAGCTGCTGCAATATTGTCTCCTGTTGCTCAATCCGTTCGCTCAACGGCTGACAGATGAGGACGTAATATTCGACGCTAAACAGGAATACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGGGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACg                                                                                                                               >  1:15896/1‑251 (MQ=255)                                                                                                      ncAGATTTTCGCGCTGCAATAATATCCGGCGCGGCGTTCTTCAG‑TCGGGTGTTCATCGCCACCAGCTGCTGCACTATTGTCTCCTGTTGCTCAATCCGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATTCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCg                                                                                                                                                      >  1:101356/2‑222 (MQ=255)                                                                                                                     gCAATAATATCCGGCGCGGCGTTCTTCAG‑TCGGGTGTTCATCGCCACCAGCTGCTGCAATATTGTCTCCTGTTGCTCAATCCGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATTCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAAt                                                                                                          >  1:35238/1‑251 (MQ=255)                                                                                                                                       gcgTTCTTCAG‑TCGGGTGTTCATCACCACCAGCTCCTGAACTATTGTCTCCTGTTGATAAATCCGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATTCGACGCTAAACAGGAAcaccac                                                                                                                                                                                                                         >  1:183333/1‑122 (MQ=255)                                                                                                                                        cgTTCTTCAG‑TCGGGTGTTCATCGCCACCAGCTGCTGCAATATTGTCTCCTGTTGCTCAATCCGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATTCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCg                                                                                                                                                      <  1:140819/188‑1 (MQ=255)                                                                                                                                                                                               gCTCAATCCGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATTCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTa                                  >  1:387291/1‑250 (MQ=255)                                                                                                                                                                                                             cTCAANGGCTGNCNNATGAGGACGTCATATTCGANGCTAAACAGGAACACNACCGNTGCCANTANc                                                                                                                                                                                                            >  2:207806/1‑66 (MQ=255)                                                                                                                                                                                                             cTCAACGGCTGCCAGATGAGGGCGTAATATTCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCAcg                                                                                                                                                                                              <  1:207806/80‑1 (MQ=255)                                                                                                                                                                                                                       gCCAGATGAGGACGTAATATTCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGNCATNNNNNaannnnnnnnnn                                                                                                          >  2:318738/1‑144 (MQ=255)                                                                                                                                                                                                                                gtaCGTAATATTCGACGCTAAACAGGAACACCACCGCTGCAAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCAAGCCAGTGTTCGCCAGGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTGcc  >  1:461460/3‑249 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              CCTGAATATTCTCCCCACGATCGGCTATCCGCCTGATGACCACCGAGTGCGCGGAAGCACTGTCACTGATGACCCGCGATACCTGCGCAGGCGTTGTTGTGGCAGATTTTCGCGCTGCAATAATATCCGGCGCGGCGTTCTTCAG‑TCGGGTGTTCATCGCCACCAGCTGCTGCAATATTGTCTCCTGTTGCTCAATCCGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTGCC  >  NC_000913/3110742‑3111212

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