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%	188.1 / NA	58	I647I (ATC→ATT)	acrD	aminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/31);  total (27/31)

AAGACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAG  >  NC_000913/2589399‑2589659                                                                                                                                         |                                                                                                                             aaGACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGa                                                                                                                            <  2:1749216/139‑1 (MQ=255)             tGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGaa                                                                                                                <  1:95959/139‑1 (MQ=255)                  ttgtTTGCCACCGTTGGTCCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGt                                                                                                           >  2:277518/2‑139 (MQ=255)                   tgtTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGtt                                                                                                          >  1:1364060/1‑139 (MQ=255)                   tgtTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGtt                                                                                                          >  1:1416288/1‑139 (MQ=255)                       tGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      <  2:713874/139‑1 (MQ=255)                       tGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      <  2:1416288/139‑1 (MQ=255)                       tGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      <  1:1398905/139‑1 (MQ=255)                               ttGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      <  2:1487744/131‑1 (MQ=255)                               ttGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      >  1:1487744/1‑131 (MQ=255)                               ttGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTaa                                                                                              >  1:24185/1‑139 (MQ=255)                                 ggTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAaag                                                                                            >  2:911571/1‑139 (MQ=255)                                            tGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGtt                                                                                 <  2:601814/139‑1 (MQ=255)                                             gggggTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAg                                                                                                                           <  2:1707268/96‑1 (MQ=255)                                             gggggTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAg                                                                                                                           >  1:1707268/1‑96 (MQ=255)                                             gggggTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      >  1:1492746/1‑117 (MQ=255)                                             gggggTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      <  2:1492746/117‑1 (MQ=255)                                             gggggTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTa                                                                                >  2:76436/1‑139 (MQ=255)                                                                       aaTGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGc                                                      >  2:1079566/1‑139 (MQ=255)                                                                        aTGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCg                                                     <  1:76436/139‑1 (MQ=255)                                                                               tCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTgg                                              <  1:1079566/139‑1 (MQ=255)                                                                                 cGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTcgcg                                                                                   >  2:902206/1‑100 (MQ=255)                                                                                 cGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTcgcg                                                                                   <  1:902206/100‑1 (MQ=255)                                                                                  gCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      <  1:1710232/80‑1 (MQ=255)                                                                                  gCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      >  2:1710232/1‑80 (MQ=255)                                                                                   ccTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTgg                                              >  2:474534/1‑135 (MQ=255)                                                                                   ccTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTgg                                              <  1:474534/135‑1 (MQ=255)                                                                                   ccTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGt                                          >  2:1193284/1‑139 (MQ=255)                                                                                    cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGTAACGAAGGCGTTTAAc                                                                                                      <  1:581603/78‑1 (MQ=255)                                                                                    cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGTAACGAAGGCGTTTAAc                                                                                                      >  2:581603/1‑78 (MQ=255)                                                                                    cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGt                                                   <  2:1732844/129‑1 (MQ=255)                                                                                    cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGt                                                   >  1:1732844/1‑129 (MQ=255)                                                                                    cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGtt                                         >  1:1150817/1‑139 (MQ=255)                                                                                             tGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAAt                                                                                                 <  2:784854/74‑1 (MQ=255)                                                                                             tGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAAt                                                                                                 >  1:784854/1‑74 (MQ=255)                                                                                              ggAGCGAACTCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGttt                               <  2:1150817/139‑1 (MQ=255)                                                                                               gAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCa                                                               >  2:1637984/1‑106 (MQ=255)                                                                                               gAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCa                                                               <  1:1637984/106‑1 (MQ=255)                                                                                                    aaCGTGACAGTAAGAACGGCACCTCGTTTGCCATTATTGAGCGTGCAACGaa                                                                                                                <  2:1143680/52‑1 (MQ=255)                                                                                                    aaCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGaa                                                                                                                >  1:1143680/1‑52 (MQ=255)                                                                                                     aCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTaa                                                                                                       <  1:436437/60‑1 (MQ=255)                                                                                                     aCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTaa                                                                                                       >  2:436437/1‑60 (MQ=255)                                                                                                      cgcgACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGcc                                                                           <  1:1405368/87‑1 (MQ=255)                                                                                                      cgcgACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGcc                                                                           >  2:1405368/1‑87 (MQ=255)                                                                                                      cgcgACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATgg                       <  1:1193284/139‑1 (MQ=255)                                                                                                          aCAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGtt                   >  1:1284135/1‑139 (MQ=255)                                                                                                           cAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTg                  <  2:1364060/139‑1 (MQ=255)                                                                                                                aGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGa             <  2:24185/139‑1 (MQ=255)                                                                                                                     ggCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      >  1:647976/1‑45 (MQ=255)                                                                                                                     ggCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                      <  2:647976/45‑1 (MQ=255)                                                                                                                       caaacgaggtGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAg                                           >  1:491848/10‑102 (MQ=255)                                                                                                                       caaacgaggtGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAg                                           <  2:491848/93‑1 (MQ=255)                                                                                                                          cTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGa   <  2:1394698/139‑1 (MQ=255)                                                                                                                           tCGTTTGCCATTATTGAGCGTGCAACGAATGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAg  <  1:1158783/139‑1 (MQ=255)                                                                                                                            cGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTa                                                        <  2:1667338/84‑1 (MQ=255)                                                                                                                            cGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTa                                                        >  1:1667338/1‑84 (MQ=255)                                                                                                                                tGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAg               >  2:1110073/1‑121 (MQ=255)                                                                                                                                tGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAg               <  1:1110073/121‑1 (MQ=255)                                                                                                                                    attattGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGa   <  2:300987/129‑1 (MQ=255)                                                                                                                                    attattGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGa   >  1:300987/1‑129 (MQ=255)                                                                                                                                         |                                                                                                                             AAGACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAG  >  NC_000913/2589399‑2589659

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