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

AAGACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCGGGTCACG  >  NC_000913/2589399‑2589668                                                                                                                                         |                                                                                                                                      aaGACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGa                                                                                                                                     <  2:788146/139‑1 (MQ=255)  aGACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAg                                                                                                                                    <  2:1606011/139‑1 (MQ=255)   gACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGc                                                                                                                                   >  2:1572942/1‑139 (MQ=255)             tGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGaa                                                                                                                         <  2:1247353/139‑1 (MQ=255)                  gtgtTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGt                                                                                                                    >  1:518699/1‑139 (MQ=255)                   tgtTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGtt                                                                                                                   <  1:818592/139‑1 (MQ=255)                     tttGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTa                                                                                                                 >  2:1486578/1‑139 (MQ=255)                       tGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                               >  2:1095876/1‑139 (MQ=255)                       tGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                               <  2:518699/139‑1 (MQ=255)                       tGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                               <  1:1572942/139‑1 (MQ=255)                       tGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                               <  1:125568/139‑1 (MQ=255)                          cACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCaa                                                                                                            <  2:560348/139‑1 (MQ=255)                                   ttCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGaa                                                                                                                         <  2:370902/117‑1 (MQ=255)                                   ttCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGaa                                                                                                                         >  1:370902/1‑117 (MQ=255)                                      tGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCt                                                                                                >  1:1533365/1‑139 (MQ=255)                                       ggCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTc                                                                                               <  1:803215/139‑1 (MQ=255)                                             gggggTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTg                                                                                                                                >  1:983869/1‑100 (MQ=255)                                             gggggTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTg                                                                                                                                <  2:983869/100‑1 (MQ=255)                                                     cGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCcagc                                                                                 <  2:786873/139‑1 (MQ=255)                                                       ggCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCcagcag                                                                               <  2:1476685/139‑1 (MQ=255)                                                        gCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCagcagc                                                                              >  2:641934/1‑139 (MQ=255)                                                         cAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCcagc                                                                                 <  1:121530/135‑1 (MQ=255)                                                         cAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCcagc                                                                                 >  2:121530/1‑135 (MQ=255)                                                         cAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGcc                                                                             <  1:641934/139‑1 (MQ=255)                                                                     cgAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTa                                                                 <  2:1259668/139‑1 (MQ=255)                                                                        aTGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCg                                                              <  2:465143/139‑1 (MQ=255)                                                                        aTGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCg                                                              <  2:1533365/139‑1 (MQ=255)                                                                          gTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGt                                                            <  1:1486578/139‑1 (MQ=255)                                                                                   ccTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTaaa                                                                                                      <  1:1367410/88‑1 (MQ=255)                                                                                   ccTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTaaa                                                                                                      >  2:1367410/1‑88 (MQ=255)                                                                                    cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                               >  1:104968/1‑78 (MQ=255)                                                                                    cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAAc                                                                                                               <  2:104968/78‑1 (MQ=255)                                                                                    cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTcgcg                                                                                            >  2:1139670/1‑97 (MQ=255)                                                                                    cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTcgcg                                                                                            <  1:1139670/97‑1 (MQ=255)                                                                                     tGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGtt                                                  <  1:770615/138‑1 (MQ=255)                                                                                     tGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGtt                                                  >  2:770615/1‑138 (MQ=255)                                                                                              ggAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTGGTTCTGc                                              >  2:986569/1‑133 (MQ=255)                                                                                              ggAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTGGTTCTGc                                              <  1:986569/133‑1 (MQ=255)                                                                                                 gCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTa                                                                 >  1:636707/1‑111 (MQ=255)                                                                                                 gCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTa                                                                 <  2:636707/111‑1 (MQ=255)                                                                                                   gAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGata                                   <  1:334248/139‑1 (MQ=255)                                                                                                    aaCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGaa                                                                                                                         <  2:1592285/52‑1 (MQ=255)                                                                                                    aaCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGaa                                                                                                                         >  1:1592285/1‑52 (MQ=255)                                                                                                    aaCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTg                                               >  2:1049340/1‑126 (MQ=255)                                                                                                    aaCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTg                                               <  1:1049340/126‑1 (MQ=255)                                                                                                    aaCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGatat                                  >  1:431207/1‑139 (MQ=255)                                                                                                                aGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGc                                                               <  1:256534/98‑1 (MQ=255)                                                                                                                aGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGc                                                               >  2:256534/1‑98 (MQ=255)                                                                                                                aGACCGGCACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGa                      <  1:1095876/139‑1 (MQ=255)                                                                                                                       cACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGc                                                               >  2:575056/1‑91 (MQ=255)                                                                                                                       cACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGc                                                               <  1:575056/91‑1 (MQ=255)                                                                                                                       cACCTCGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCt               <  1:1316927/139‑1 (MQ=255)                                                                                                                            cGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGt                                                            >  1:1487008/1‑89 (MQ=255)                                                                                                                            cGTTTGCCATTATTGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGt                                                            <  2:1487008/89‑1 (MQ=255)                                                                                                                                    attattGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCGGGTCACg  <  2:194331/139‑1 (MQ=255)                                                                                                                                    attattGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCGGGTCACg  <  1:51806/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      AAGACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCGGGTCACG  >  NC_000913/2589399‑2589668

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