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.7 / 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)

ACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGA  >  NC_000913/2589402‑2589658                                                                                                                                      |                                                                                                                            acaacaTCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCg                                                                                                                        <  2:571994/139‑1 (MQ=255)     catcatGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCa                                                                                                                    >  1:241476/1‑139 (MQ=255)      atcatGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCaa                                                                                                                   >  1:673576/1‑139 (MQ=255)          tGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGaa                                                                                                               >  1:901871/1‑139 (MQ=255)           gTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAg                                                                                                              >  2:400171/1‑139 (MQ=255)                tgtTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTg                                                                                                                      <  1:1017710/126‑1 (MQ=255)                tgtTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTg                                                                                                                      >  2:1017710/1‑126 (MQ=255)                    tGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                     <  1:400171/139‑1 (MQ=255)                     gCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                     <  2:987174/138‑1 (MQ=255)                     gCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                     >  1:987174/1‑138 (MQ=255)                      ccACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCa                                                                                                   <  2:901871/139‑1 (MQ=255)                                 tCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAaagaag                                                                                        >  2:288807/1‑139 (MQ=255)                                     gCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTcg                                                                                    <  2:731756/139‑1 (MQ=255)                                          gggggTAACGGGCAAAACGTGGCGCGAGTGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTGAAGAAGCTCGCGTTa                                                                               >  2:815632/1‑139 (MQ=255)                                               tAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGcc                                                                          >  2:1216410/1‑139 (MQ=255)                                                     gCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCagcagc                                                                    >  1:954928/1‑139 (MQ=255)                                                     gCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCagcagc                                                                    >  1:336472/1‑139 (MQ=255)                                                      cAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACTAAGGCGTTTAACCAAATTAAAGAAGCTc                                                                                     <  1:789562/121‑1 (MQ=255)                                                      cAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTc                                                                                     >  2:789562/1‑121 (MQ=255)                                                              ggCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCt                                                                                      <  1:24533/112‑1 (MQ=255)                                                              ggCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCt                                                                                      >  2:24533/1‑112 (MQ=255)                                                                  cgAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTa                                                       <  2:512152/139‑1 (MQ=255)                                                                   gAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAg                                                      <  2:954928/139‑1 (MQ=255)                                                                    agtGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGCAGCTCGCGTTATCGCCAGCAGACTGCCAGCAATTAGc                                                     <  1:815632/137‑1 (MQ=255)                                                                    aaTGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGc                                                     <  2:160944/139‑1 (MQ=255)                                                                     aTGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCg                                                    <  1:1035621/139‑1 (MQ=255)                                                                       gTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGt                                                  <  2:336472/139‑1 (MQ=255)                                                                                ccTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGa                                                                                                                           <  2:1099598/57‑1 (MQ=255)                                                                                 cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGtt                                        >  1:850541/1‑139 (MQ=255)                                                                                   gAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCt                                      >  1:37131/1‑139 (MQ=255)                                                                                              gCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGc                                                                >  2:398420/1‑102 (MQ=255)                                                                                              gCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGc                                                                <  1:398420/102‑1 (MQ=255)                                                                                                                  ggCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGc                                                     <  2:320799/93‑1 (MQ=255)                                                                                                                  ggCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGc                                                     >  1:320799/1‑93 (MQ=255)                                                                                                                   gCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGc                                                                                       >  2:1118945/1‑58 (MQ=255)                                                                                                                   gCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGc                                                                                       <  1:1118945/58‑1 (MQ=255)                                                                                                                   gCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGc      >  2:693578/1‑139 (MQ=255)                                                                                                                   gCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAACAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGc      <  1:792638/139‑1 (MQ=255)                                                                                                                       cTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGa  <  2:673576/139‑1 (MQ=255)                                                                                                                       cTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGa  <  1:288807/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                            ACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGA  >  NC_000913/2589402‑2589658

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