Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,993,754	T→A	D39E (GAT→GAA)	ygeI →	protein YgeI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,993,754	0	T	A	100.0%	45.4 / NA	14	D39E (GAT→GAA)	ygeI	protein YgeI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (7/7);  total (7/7)

TCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAACACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCT  >  NC_000913/2993677‑2993826                                                                              |                                                                         tCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAACACATCAACACCAGTGCTACTGATACAGAACa                                                                        >  2:546458/1‑80 (MQ=255)         aGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAACACATCAACACCAGTGCTACTGATACAGAACAAATACAAg                                                                >  1:546464/1‑80 (MQ=255)          gCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAACACATCAACACCAGTGCTACTGATACAGAACAAATACAAGc                                                               <  1:546437/80‑1 (MQ=255)                gCAACGGGTGATGAAGTCGTATCATTAGATAAACACATCAACACCAGTGCTACTGATACAGAACAAATACAAGCATTta                                                          >  2:546460/1‑79 (MQ=255)                       gTGATGAAGTCGTATCATTAGATAAACACATCAACACCAGTGCTACTGATACAGAACAAATACAAGCATTTATAGTTTCaa                                                 <  2:546438/81‑1 (MQ=255)                            gaAGTCGTATCATTAGATAAACACATCAACACCAGTGCTACTGATACAGAACAAATACAAGCATTTATAGTTTCAACatg                                             >  2:546459/1‑80 (MQ=255)                                  gTATCATTAGATAAACACATCAACACCAGTGCTACTGATACAGAACAAATACAAGCATTTATAGTTTCAACATGGATGgc                                       <  1:546439/80‑1 (MQ=255)                                    aTCATTAGATAAACACATCAACACCAGTGCTACTGATACAGAACAAATACAAGCATTTATAGTTTCAACATGGATGGCGcc                                    >  2:546457/1‑81 (MQ=255)                                     tCATTAGATAAACACATCAACACCAGTGCTACTGATACAGAACAAATACAAGCATTTATAGTTTCAACATGGATGGCGcc                                    <  1:546440/80‑1 (MQ=255)                                       aTTAGATAAACACATCAACACCAGTGCTACTGATACAGAACAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCtt                                  <  1:546441/80‑1 (MQ=255)                                                      cAACACCAGTGCTACTGATACAGAACAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTa                  >  1:546466/1‑81 (MQ=255)                                                                   aCTGATACAGAACAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAAt     >  1:546465/1‑81 (MQ=255)                                                                    cTGATACAGAACAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATc    <  1:546442/81‑1 (MQ=255)                                                                       aTACAGAACAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCt  <  2:546443/80‑1 (MQ=255)                                                                              |                                                                         TCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAACACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCT  >  NC_000913/2993677‑2993826

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