Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,813,364	T→C	V250A (GTC→GCC)	emrB →	multidrug efflux pump membrane subunit EmrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,813,364	0	T	C	100.0%	125.0 / NA	33	V250A (GTC→GCC)	emrB	multidrug efflux pump membrane subunit EmrB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/21);  total (12/21)

ATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTTGT  >  NC_000913/2813318‑2813408                                               |                                             atcatcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCt                                            >  1:1457319/1‑49 (MQ=255)  tcatcCTTACCGTGGTGGCGGTGGTGGCTCTCTTCTTCCTGATTGCCt                                            <  1:2190565/48‑1 (MQ=38)  tcatcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCt                                            <  1:1964363/48‑1 (MQ=255)  tcatcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCt                                            <  1:1294904/48‑1 (MQ=255)  tcatcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTg                                           >  1:1385606/1‑49 (MQ=255)   catcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                          >  1:1886059/1‑49 (MQ=255)    atcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                          >  1:1935871/1‑48 (MQ=255)       cTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGc                                      >  1:712514/1‑49 (MQ=255)       cTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGc                                      >  1:171111/1‑49 (MQ=255)       cTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGc                                      <  1:709100/49‑1 (MQ=255)          aCCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGa                                   >  1:697019/1‑49 (MQ=255)            cGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGAcc                                 >  1:1621460/1‑49 (MQ=255)            cGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGAcc                                 <  1:1812069/49‑1 (MQ=255)             gtggtgGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGAcc                                 >  1:705730/1‑48 (MQ=255)             gtggtgGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACcg                                >  1:1267185/1‑49 (MQ=255)                gtggCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACcgac                              <  1:238022/48‑1 (MQ=255)                    cGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATaa                         <  1:232800/49‑1 (MQ=255)                     ggtggtGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAc                        <  1:1387550/49‑1 (MQ=255)                     ggtggtGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAc                        <  1:402461/49‑1 (MQ=255)                      gtggtgGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAcc                       <  1:1808640/49‑1 (MQ=255)                          tggCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGAt                   <  1:1445689/49‑1 (MQ=255)                            gCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATa                  <  1:5619/48‑1 (MQ=255)                            gCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAg                 >  1:1256665/1‑49 (MQ=255)                             cTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGt                <  1:77007/49‑1 (MQ=255)                              tATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTc               >  1:2129239/1‑49 (MQ=255)                                tCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCg              <  1:203016/48‑1 (MQ=255)                                tCTGCTTCCTGATTGCCCGGGAGCTGACCGACGATAACCCGATAGTc               <  1:69676/47‑1 (MQ=38)                                    cTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGAt            <  1:1136499/46‑1 (MQ=255)                                    cTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCt          <  1:2350051/48‑1 (MQ=255)                                       ccTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGt        <  1:403015/47‑1 (MQ=255)                                        cTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCg      <  1:1352185/48‑1 (MQ=255)                                           aTTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgttgt  <  1:2134784/49‑1 (MQ=255)                                           aTTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgtt    <  1:1388124/47‑1 (MQ=255)                                               |                                             ATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTTGT  >  NC_000913/2813318‑2813408

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