Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,854,927	G→A	G336S (GGC→AGC)	emrD →	multidrug efflux pump EmrD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,854,927	0	G	A	100.0%	31.8 / NA	14	G336S (GGC→AGC)	emrD	multidrug efflux pump EmrD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/7);  total (7/7)

GTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATC  >  NC_000913/3854704‑3855132                                                                                                                                                                                                                                |                                                                                                                                                                                                              gTTTGCCGGTCGTCCCAATAAACGCTTCTCCACGTTATGGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTGGCTGATGTGGATCCCCGACTGGTTGGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTAGCGCGCTGGTcggcgg                                                                                                                                                                                               <  1:41639/240‑1 (MQ=255)                                                                        actagCTGATGTGGATCCACGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTACCGCAGCGCTGATCTTTTTCGGTGCCGGGATGCTGTTTACGCTGGCGCCCAGCGGCACGATGGAGCCGTTCCCCTTCCTGACGGGCACGGCTAGCGCGCTGGTCGGCGGTTTGCCAACCATTGGTTCCGGCGTCCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCTGCCTGGGGttgttg                                                                                                             >  1:343284/5‑251 (MQ=255)                                                                            ncTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTAGCGCGCTGGTCGGCGGTCTGCCAACCATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGc                                                                                                                                             >  1:49483/2‑215 (MQ=255)                                                                                                       cGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGATGGCGAACAGCGACGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTAGCGCGCTGGTAGGCGGTCTGCAAACCATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGCCCTTAATGGGATTGTTGATAgtgctgt                                                                              >  1:377464/1‑251 (MQ=255)                                                                                                                          ctgctCGTGCCCGCCGCGCTGTTCTTTTTCGGGGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTAGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTggn                                                                                                                                                                 <  1:59119/149‑2 (MQ=255)                                                                                                                              tCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGAACGGCTAGCGCGCTGGTCGGCGGCCTGCAAAACATTGATTCCGGCGTGCTGGCGTCGCCCTCTGCGATGTTGCCGCAAACCGGTCAGGGCCGCCTGGGGTTGTTGATGACCTTAATGGGGTTGTTGATCGTGCTGTgctgg                                                                         >  1:206497/1‑233 (MQ=255)                                                                                                                                     cgccgcGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTAGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAgg                                                                                                                              <  1:189496/173‑1 (MQ=255)                                                                                                                                                   tttgcGGTGCCGGGATGCTGGTTCCTCGGTCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTAgcg                                                                                                                                                                                                            <  1:333172/77‑1 (MQ=255)                                                                                                                                                              gggATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTAGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCg                                                                                                                                           >  1:31417/1‑135 (MQ=255)                                                                                                                                                                                                           cccTTCCTGGCGGGCACGGCTAGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGNATTNNNNNTcnnnnnnnnnn                                                                           >  2:305976/1‑144 (MQ=255)                                                                                                                                                                                                           cccTTCCTGGCGGGCACGGCTAGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGtcatn  <  1:305976/227‑2 (MQ=255)                                                                                                                                                                                                                      gggCACGGCTAGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGc                                       >  1:85925/1‑179 (MQ=255)                                                                                                                                                                                                                        gCACGGCTAGCGCTCTGGTCGGCGGGCTGCAAAACATTGGTTCCTGCGTGCTGGCGTCGCTCTCTGCGATGTTGCGGCAAACCGGTCAGGGCAGCCTGGGGTTGTGGAGGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGAcgcg                                                        <  1:360071/160‑1 (MQ=255)                                                                                                                                                                                                                        gCACGGCTAGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCtn                                                                                                                                                                   <  1:94821/53‑2 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                              GTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGCGTGATGAATGTCTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGCACGTCACCGCAGCATCGTCATC  >  NC_000913/3854704‑3855132

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