Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,900,530	T→C	E352E (GAA→GAG)	yoaE ←	putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,900,530	0	T	C	100.0%	52.4 / NA	16	E352E (GAA→GAG)	yoaE	putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/14);  total (2/14)

TCCCGCGAGCACGACGCAGCACGCCCAACAGGTTGATCGGATCGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCG  >  NC_000913/1900406‑1900629                                                                                                                             |                                                                                                    tCCCGCGAGCACGACGCAGCACGCCCAACAGGTTGATCGGATCGAGGGTTTCCGGGACGATAATCGCCGGAGCCGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGc                                                                                             <  2:464872/133‑1 (MQ=255) tCCCGCGAGCACGACGCAGCACGCCCAACAGGTTGATCGGATCGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGc                                                                                             <  1:401723/133‑1 (MQ=255) tCCCGCGAGCACGACGCAGCACGCCCAACAGGTTGATCGGATCGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGc                                                                                             <  1:449365/133‑1 (MQ=255) tCCCGCGAGCACGACGCAGCACGCCCAACAGGTTGATCGGATCGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGc                                                                                             <  1:92983/133‑1 (MQ=255)                    caCGCCCAACAGGTTGATCGGATCGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGAtt                                                                          >  2:76701/1‑133 (MQ=255)                                       ggatcgAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGAc                                                       <  1:275197/133‑1 (MQ=255)                                       ggatcgAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGAc                                                       <  2:275627/133‑1 (MQ=255)                                           cgAGGGTTTCCGGGAGGAGAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATAc                                                   <  1:223267/133‑1 (MQ=255)                                           cgAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATAc                                                   <  2:440456/133‑1 (MQ=255)                                           cgAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAACCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATAc                                                   <  1:96212/133‑1 (MQ=255)                                                   tCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACa                                           >  1:377479/1‑133 (MQ=255)                                                     cGGGACGATAATCGCAGGAGACGCCGAAGCAATCGCCGCCATATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGa                                         <  2:204625/133‑1 (MQ=255)                                                       ggACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACt                                       <  1:76701/133‑1 (MQ=255)                                                            aTAATCACCGTAAACGCCGAAGAAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCCACATACCGGGAACAGACTGTGCg                                  <  2:126755/133‑1 (MQ=255)                                                                   ccGGAGACGCCGAAGCAATCGCCGCCACGTCATCGCCCTCTCCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAaga                           <  2:183544/133‑1 (MQ=255)                                                                   ccGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAaga                           <  2:377479/133‑1 (MQ=255)                                                                                            cacaTCAACGCCCTCTTCCAGCGCCACCAGCAGCTCTTTAGTACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCg  <  2:341861/133‑1 (MQ=255)                                                                                                                             |                                                                                                    TCCCGCGAGCACGACGCAGCACGCCCAACAGGTTGATCGGATCGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCG  >  NC_000913/1900406‑1900629

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