Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,744,337	(A)7→6	coding (10/474 nt)	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,331	0	A	.	100.0%	94.5 / NA	23	coding (4/474 nt)	yiaM	2,3‑diketo‑L‑gulonate:Na(+) symporter ‑ membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (11/12);  total (11/12)

CTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAG  >  NC_000913/3744253‑3744407                                                                               |                                                                             ctctAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGaaa                                                                            <  1:665872/81‑1 (MQ=255) ctctAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGaa                                                                             <  1:665873/80‑1 (MQ=255)   ctAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGaaaa                                                                           >  2:665911/1‑80 (MQ=255)      aCTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGaaaaaat                                                                        <  2:665873/80‑2 (MQ=255)         atatTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTc                                                                   >  2:665904/1‑81 (MQ=255)                 aatTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCaa                                                            >  2:665912/1‑80 (MQ=255)                      aGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTgg                                                      >  1:665920/1‑81 (MQ=255)                        gTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGc                                                     >  2:665914/1‑80 (MQ=255)                                acaTGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAAtct                                            <  2:665875/81‑1 (MQ=255)                                acaTGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAAtct                                            >  1:665914/1‑81 (MQ=255)                                acaTGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAAtct                                            <  1:665875/81‑1 (MQ=255)                                 caTGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAAtctc                                           >  2:665899/1‑81 (MQ=255)                                     cGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGcc                                        <  1:665877/80‑1 (MQ=255)                                      gACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGt                                      <  2:665877/81‑1 (MQ=255)                                      gACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGt                                      <  1:665879/81‑1 (MQ=255)                                         ttGAATTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACt                                   >  2:665917/1‑81 (MQ=255)                                            aaTTATGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACttt                                 <  1:665880/80‑1 (MQ=255)                                                aTGACAGACATAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCtg                             <  2:665880/80‑1 (MQ=255)                                                   acagacaTAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCtgtatt                         <  2:665881/81‑1 (MQ=255)                                                     agacaTAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCtgtattgt                       <  2:665882/81‑1 (MQ=255)                                                      gacaTAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCtgtattgta                      >  1:665919/1‑81 (MQ=255)                                                       acaTAGCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCtgtattgtat                     <  1:665884/81‑1 (MQ=255)                                                            gCCTCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTtata                <  1:665885/81‑1 (MQ=255)                                                               tCAAGGAATAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTTATaaa              >  1:665917/1‑80 (MQ=255)                                                                      aTAGCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTTATAAACATTAttt      >  2:665922/1‑81 (MQ=255)                                                                         gCTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTTATAAACATTAtttt     >  2:665923/1‑79 (MQ=255)                                                                          cTATG‑AAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAg  <  2:665885/81‑1 (MQ=255)                                                                               |                                                                             CTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAG  >  NC_000913/3744253‑3744407

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