Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	481,421	C→A	A995S (GCT→TCT)	acrB ←	multidrug efflux pump RND permease AcrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	481,421	0	C	A	100.0%	101.2 / NA	28	A995S (GCT→TCT)	acrB	multidrug efflux pump RND permease AcrB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/12);  total (16/12)

CGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAGGTCATCAGGATC  >  NC_000913/481359‑481482                                                               |                                                              cGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAAc                                                             >  1:2189363/1‑65 (MQ=255)    tGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGaaccagaaccag                                                          <  1:51066/65‑1 (MQ=255)     gCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGt                                                         >  1:2526311/1‑65 (MQ=255)       ggTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTAc                                                       >  1:1674935/1‑65 (MQ=255)       ggTCACCATCCCGCCCATTACACCGGCACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTAc                                                       >  1:3860984/1‑65 (MQ=37)        gTCACCATCCCGCCCATTTCACCGTTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACt                                                      <  1:1665520/65‑1 (MQ=37)        gTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACt                                                      <  1:2736072/65‑1 (MQ=255)        gTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACt                                                      <  1:1654079/65‑1 (MQ=255)        gTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACt                                                      <  1:1707036/65‑1 (MQ=255)        gTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACt                                                      <  1:1798202/65‑1 (MQ=255)          cACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGa                                                    <  1:1528388/65‑1 (MQ=255)               tCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAAcc                                               <  1:1559488/65‑1 (MQ=255)                cccGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCa                                              >  1:4023588/1‑65 (MQ=255)                     ccATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGc                                         <  1:2348026/65‑1 (MQ=255)                     ccATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGc                                         >  1:1377614/1‑65 (MQ=255)                          acacCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAAc                                    >  1:904310/1‑65 (MQ=255)                            acCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGc                                  >  1:2249537/1‑65 (MQ=255)                             ccGGTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGCCATAACGcc                                 >  1:764759/1‑65 (MQ=255)                               ggTACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGa                               >  1:3396119/1‑65 (MQ=255)                                 tACCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAgg                             <  1:2278752/65‑1 (MQ=255)                                 tACCTACTGCGGTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAgg                             <  1:1496124/65‑1 (MQ=255)                                  aCCTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAGGa                            >  1:2532687/1‑65 (MQ=255)                                    cTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAGGATa                          >  1:3704167/1‑65 (MQ=255)                                    cTACTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAGGATa                          >  1:2986562/1‑65 (MQ=255)                                      aCTGCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAGGATaaa                        >  1:90312/1‑65 (MQ=255)                                         gCGTTCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGc                     >  1:675631/1‑65 (MQ=255)                                             tCTGCGCGCCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGc                 >  1:1524347/1‑65 (MQ=255)                                                  gcgcCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAGGATCAACGCCAGCGAggc            >  1:1812812/1‑64 (MQ=25)                                                  gcgcCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAGGt            <  1:208494/65‑1 (MQ=255)                                                  gcgcCGGAACCAGAACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAGGt            <  1:983439/65‑1 (MQ=255)                                                            cAGAACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAGGTCATCAGGATc  >  1:3476205/1‑65 (MQ=255)                                                               |                                                              CGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAGGTCATCAGGATC  >  NC_000913/481359‑481482

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