Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,589,535	C→T	I647I (ATC→ATT)	acrD →	aminoglycoside/multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,589,535	0	C	T	100.0%	127.1 / NA	40	I647I (ATC→ATT)	acrD	aminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/21);  total (19/21)

AACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCGGGTCACGAT  >  NC_000913/2589404‑2589670                                                                                                                                    |                                                                                                                                        aacaTCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTg                                                                                                                                  <  1:56354/139‑1 (MQ=255)   catcatGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCa                                                                                                                                >  2:217995/1‑139 (MQ=255)   catcatGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCa                                                                                                                                >  2:760044/1‑139 (MQ=255)        tGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGaa                                                                                                                           <  1:760044/139‑1 (MQ=255)         gTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAg                                                                                                                          >  2:1275631/1‑139 (MQ=255)         gTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAg                                                                                                                          >  2:69244/1‑139 (MQ=255)             gtgtTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGt                                                                                                                      >  2:519185/1‑139 (MQ=255)                 ttGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTaa                                                                                                                  <  2:1272483/139‑1 (MQ=255)                  tGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                                 <  1:1261873/139‑1 (MQ=255)                  tGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                                 <  1:69244/139‑1 (MQ=255)                                  ggCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTc                                                                                                 >  1:1034042/1‑139 (MQ=255)                                       tGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGtt                                                                                            >  2:291994/1‑139 (MQ=255)                                       tGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGtt                                                                                            <  1:393283/139‑1 (MQ=255)                                        gggggTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGt                                                                                                                      >  2:374588/1‑112 (MQ=255)                                        gggggTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGt                                                                                                                      <  1:374588/112‑1 (MQ=255)                                        gggggTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACTAAATTAAAGAAGCTCGCGTTa                                                                                           >  2:867980/1‑139 (MQ=255)                                                   gCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGCGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCagcagc                                                                                >  1:1409851/1‑139 (MQ=255)                                                       aaCGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGc                                                                            <  2:539914/139‑1 (MQ=255)                                                            ggCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                                 >  1:676663/1‑97 (MQ=255)                                                            ggCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                                 <  2:676663/97‑1 (MQ=255)                                                            ggCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCa                                                                       <  1:529015/139‑1 (MQ=255)                                                             gcgcGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCaa                                                                      >  2:1361405/1‑139 (MQ=255)                                                                 gAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAg                                                                  <  1:1275631/139‑1 (MQ=255)                                                                  aaTGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGc                                                                 <  1:1076513/139‑1 (MQ=255)                                                                     gTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGt                                                              <  1:217995/139‑1 (MQ=255)                                                                     gTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGt                                                              <  1:519185/139‑1 (MQ=255)                                                                               cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGtt                                                    >  1:1013378/1‑139 (MQ=255)                                                                               cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGtt                                                    >  2:450131/1‑139 (MQ=255)                                                                                tGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                                 >  1:1437906/1‑77 (MQ=255)                                                                                tGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                                 <  2:1437906/77‑1 (MQ=255)                                                                                  aaaGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTg                                                 <  2:1034042/139‑1 (MQ=255)                                                                                   aaGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGc                                                <  2:1013378/139‑1 (MQ=255)                                                                                        tGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGtt                                           >  2:312421/1‑139 (MQ=255)                                                                                        tGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGtt                                           <  1:450131/139‑1 (MQ=255)                                                                                         ggAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                                 >  2:650290/1‑68 (MQ=255)                                                                                         ggAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                                 <  1:650290/68‑1 (MQ=255)                                                                                           aGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTg                                        >  2:958278/1‑139 (MQ=255)                                                                                                                 gCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCcagcag                                                                                 <  2:1455707/76‑1 (MQ=255)                                                                                                                 gCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCcagcag                                                                                 >  1:1455707/1‑76 (MQ=255)                                                                                                                                 tATTGAGCGTGTAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCGGGTCACGAt  <  2:1409851/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        AACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCGGGTCACGAT  >  NC_000913/2589404‑2589670

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