Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,989,072	G→A	E191K (GAA→AAA)	yqeI →	putative transcriptional regulator YqeI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,989,072	0	G	A	100.0%	74.0 / NA	21	E191K (GAA→AAA)	yqeI	putative transcriptional regulator YqeI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/11);  total (10/11)

AGGCGCATATTCTGCATATTGGCTATGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGGAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATA  >  NC_000913/2988996‑2989141                                                                             |                                                                      aGGCGCATATTCTGCATATTGGCTATGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAAt                                                                   >  2:489420/1‑81 (MQ=255)  ggCGCATATTCTGCATATTGGCTATGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAAtt                                                                  <  1:489404/81‑1 (MQ=255)  ggCGCATATTCTGCATATTGGCTATGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAAt                                                                   >  2:489434/1‑80 (MQ=255)            cTGCATATTGGCTATGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGc                                                         <  1:489405/80‑1 (MQ=255)                    tGGCTATGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAAt                                                >  1:489423/1‑81 (MQ=255)                      gCTATGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGt                                              <  2:489405/81‑1 (MQ=255)                         aTGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTaa                                            >  1:489428/1‑80 (MQ=255)                         aTGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTaa                                            >  2:489425/1‑80 (MQ=255)                          tGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAAc                                           <  1:489407/80‑1 (MQ=255)                              ataataatCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACag                                         >  1:489427/1‑78 (MQ=255)                                 ataatCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATa                                   <  2:489407/81‑1 (MQ=255)                                     tCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAAt                               <  1:489409/81‑1 (MQ=255)                                              gCCATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGtt                      <  2:489409/81‑1 (MQ=255)                                                cATTTTTCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGTTg                     >  2:489424/1‑80 (MQ=255)                                                    tttCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGTTGAAag                 <  1:489411/80‑1 (MQ=255)                                                      tCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGTTGAAagag               <  2:489411/80‑1 (MQ=255)                                                      tCAAGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGTTGAAAGAGt              >  1:489431/1‑81 (MQ=255)                                                        aaGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGTTGAAAGAGtt             <  1:489413/80‑1 (MQ=255)                                                        aaGGATTACAAGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGTTGAAAGAGTTc            >  1:489429/1‑81 (MQ=255)                                                                 aaGACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAatat   <  1:489414/81‑1 (MQ=255)                                                                   gACTGTTGCGAAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAAtata  >  1:489434/1‑80 (MQ=255)                                                                             |                                                                      AGGCGCATATTCTGCATATTGGCTATGGAATAATAATCAACCGAAGCCATTTTTCAAGGATTACAAGACTGTTGCGGAAATTAATGGCTGCCATTTTAATGTAACAGAAGATACAATTGATGGGTTGAAAGAGTTCGATAAATATA  >  NC_000913/2988996‑2989141

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: ‑