Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	884,298	G→A	G209D (GGT→GAT)	mdfA →	multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	884,298	0	G	A	100.0%	59.4 / NA	18	G209D (GGT→GAT)	mdfA	multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (4/14);  total (4/14)

CCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGA  >  NC_000913/884055‑884534                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             cctaaatgaaGGAATCCTTGGAAGAGGCGTTTTGTATTAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTCCTCCGCTACTTGTTCGGCTGGGGGGCGGGGCGTGGTTCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGAATTGGCAGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGAc                                                                                                                                                                                                                                       <  2:131986/242‑1 (MQ=255)                                       aaGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCggg                                                                                                                                                                                                 <  1:80994/251‑1 (MQ=255)                                        aGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCgggg                                                                                                                                                                                                >  2:307825/1‑251 (MQ=255)                                               cgcgCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCggggggggtgg                                                                                                                                                                                         >  2:33054/1‑245 (MQ=255)                                                gcgcTGATGGCGAGCGTGGCGCTTATTGCTCCGCTACTTGGTCCGCTGGTGGGGGCGGCGTGGTTCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTTTTTGCCGCATTGGCAGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGggcgctggc                                                                                                                                                                                        <  2:94892/251‑1 (MQ=255)                                                               gTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCACTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCAATGAAAGAACTCGATAGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCgggggggggggaggcgggagttgttcg                                                                                                                                                                         >  1:32613/1‑229 (MQ=255)                                                                                      tggttcgttggggggcggggcgTGGTTCCATGTGGTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCTGGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGAt                                                                                                                                                    <  2:125469/237‑1 (MQ=255)                                                                                         gccgtgggtggGCGCGGCGTGGATCCGTGTGCTGCCCGTGGAGGGGGTGTTTGTTTTGTGTGCCGCATTGGCAGCGATCTC‑TTTTGCGGGTCTGAAACGAG‑CCATGCCGGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCGGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGccc                                                                                                                                               <  2:456066/245‑1 (MQ=255)                                                                                                      gcggcgTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTat                                                                                                                                  <  1:307825/251‑1 (MQ=255)                                                                                                                     tttctgcCCTGGGTGGGGATGTTTTTGTTGTTTGCCGGATAGACAGTGATCTCCTTTTTCGA‑TCTCCAA‑GAGGCCATGCCTGAAACCGCTACGCGTATAGGCGAGAAACTGTCTCTGAAAGAACTCGATCCTGACTATAAGCTGGTGCTGAAGAACGGCCGCGTTGTGGCGGGGGCGCTGGTGCGTGGATTCGTTAGTTTGCCGGTGCTGGCGTGGCCCGCCCAGTCGCGGATTAGCATCATTACCGGct                                                                                                                    <  2:200523/247‑2 (MQ=255)                                                                                                                           cccTGGGAGGGTATGTTTGTTTTGTTTGCGGCATTGACATGGTTCCCCTTTTGCGG‑TCGGCAACGTG‑CTAGGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCCCTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGtt                                                                                                             <  2:384013/251‑1 (MQ=255)                                                                                                                                      gttttttgttttgtTTGCCGCATTGGCAGCGATCTCTTTTTTCGG‑TTTGCAAAGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATg                                                                                                  <  2:228135/247‑1 (MQ=255)                                                                                                                                           tggttttgtgtGCCGCATGGGCAGCGATCTCCTTTTTCGG‑TATGAAACGAG‑CCAGGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAata                                                                                              <  2:346852/248‑1 (MQ=255)                                                                                                                                                                cagggATCTCCTTTTGCGG‑TCTGCAACGGG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTAtt                                                                        <  2:435176/247‑1 (MQ=255)                                                                                                                                                                 aGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCTGGGGGTGGGGTTCGTTAGTCTGCCGGTGGTGGGGGGGGTCGCCCCCTCGCCGATTTTCATCATTACCCGCGGGGAGGTGAGGAGCTATGAATATGGGTttctgctcgtggcctttt                                                                       >  1:72340/1‑239 (MQ=255)                                                                                                                                                                                                                           gCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCAcc               <  2:100888/251‑1 (MQ=255)                                                                                                                                                                                                                           gCGAGAAACTGTCACTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCAcc               <  1:313982/251‑1 (MQ=255)                                                                                                                                                                                                                                        aCTGAAAGAACTCGATCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGa  <  1:67637/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             CCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGG‑TCTGCAACGAG‑CCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGA  >  NC_000913/884055‑884534

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