Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,744,369	C→T	A14A (GCC→GCT)	yiaM →	2,3‑diketo‑L‑gulonate:Na(+) symporter ‑ membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,744,369	0	C	T	100.0%	65.5 / NA	20	A14A (GCC→GCT)	yiaM	2,3‑diketo‑L‑gulonate:Na(+) symporter ‑ membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/6);  total (14/6)

GGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTCCTTAATCTTATTTATTTGTGGCG  >  NC_000913/3744123‑3744607                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               ggATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTAc                                                                                                                                                                                                                                            >  2:167323/1‑251 (MQ=255)                                       ctGAGATACGAAAAATTGTGGATAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTATAAGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGatat                                                                                                                                                                                                      >  2:96975/1‑251 (MQ=255)                                         gagaTACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGata                                                                                                                                                                                                       <  1:82793/248‑1 (MQ=255)                                              aCGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTc                                                                                                                                                                                               <  2:77914/251‑1 (MQ=255)                                                    aaTTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACa                                                                                                                                                                                          >  2:223734/1‑250 (MQ=255)                                                                                  gAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGtt                                                                                                                                                           >  2:152041/1‑251 (MQ=255)                                                                                                                         aaTGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTa                                                                                                                     <  1:182674/250‑1 (MQ=255)                                                                                                                             ttATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCttt                                                                                                                >  1:253161/1‑251 (MQ=255)                                                                                                                                 tGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATg                                                                                                            >  1:143787/1‑251 (MQ=255)                                                                                                                                   ctctAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATgg                                                                                                           >  1:241772/1‑250 (MQ=255)                                                                                                                                           atatTAAGAATTCAGGTTGCTGCACATTCGACTTGAATTAGGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTAGTTTAAGATATGGTTTTCAGACAAGCATTTTAGCTGTTGATGAATTGTCACGTTATTTATGTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGc                                                                                                    <  2:145350/249‑1 (MQ=255)                                                                                                                                                                                     acagacaTAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAAATCTCCCCCGCAt                                                        >  2:289166/1‑251 (MQ=255)                                                                                                                                                                                      cagacaTAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATg                                                       <  1:238199/251‑1 (MQ=255)                                                                                                                                                                                           aTAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTCATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGc                                                  >  2:193772/1‑251 (MQ=255)                                                                                                                                                                                               ccTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGcgacg                                              >  1:110373/1‑251 (MQ=255)                                                                                                                                                                                                   aaGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTc                                          >  2:239581/1‑251 (MQ=255)                                                                                                                                                                                                        aTAGCTATGAAAAAAATACTCGCAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGctct                                     >  2:278740/1‑251 (MQ=255)                                                                                                                                                                                                                                 cAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTACTTAATCttattt            >  2:164844/1‑251 (MQ=255)                                                                                                                                                                                                                                  aaTACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTCCTTAATCttattt            >  2:269856/1‑250 (MQ=255)                                                                                                                                                                                                                                           gATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTCCTTAATCTTATTTATTTGTGgcg  <  1:152041/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               GGATGTTCACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTCCTTAATCTTATTTATTTGTGGCG  >  NC_000913/3744123‑3744607

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