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_0009132,813,3640TC100.0% 107.6 / NA 29V250A (GTC→GCC) emrBmultidrug efflux pump membrane subunit EmrB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/17);  total (12/17)

ATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTT  >  NC_000913/2813318‑2813406
                                              |                                          
atcatcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTg                                         >  1:2098823/1‑50 (MQ=255)
 tcatcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                        >  1:1134427/1‑50 (MQ=255)
 tcatcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                        >  1:577277/1‑50 (MQ=255)
 tcatcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                        <  1:2076115/50‑1 (MQ=255)
   atcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                        >  1:841012/1‑48 (MQ=255)
         aCCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGAc                                >  1:1509129/1‑50 (MQ=255)
            gtggtgGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACcga                             >  1:1372142/1‑50 (MQ=255)
               gtggCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACcgacga                          <  1:1714017/50‑1 (MQ=255)
               gtggCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACcgacg                           <  1:1520735/49‑1 (MQ=255)
                  gCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATaa                       <  1:2166873/50‑1 (MQ=255)
                   cGGTGGTGGCTATCTGCTTCTTGATTGCCTGGGAGCTGACCGACGATAAc                      >  1:17301/1‑50 (MQ=255)
                   cGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAc                      >  1:1582385/1‑50 (MQ=255)
                     gtggtgGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAccc                    <  1:1617304/50‑1 (MQ=255)
                     gtggtgGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAccc                    <  1:1801596/50‑1 (MQ=255)
                     gtggtgGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAccc                    <  1:1487202/50‑1 (MQ=255)
                         tggCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATa                <  1:829891/50‑1 (MQ=255)
                         tggCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATa                <  1:714506/50‑1 (MQ=255)
                           gCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGt              <  1:521622/50‑1 (MQ=255)
                           gCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGt              <  1:690369/50‑1 (MQ=255)
                            cTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTc             <  1:101349/50‑1 (MQ=255)
                            cTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTc             <  1:1389274/50‑1 (MQ=255)
                            cTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACACGATAGTc             >  1:1431587/1‑50 (MQ=255)
                              aTCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGa           >  1:2188105/1‑50 (MQ=255)
                              aTCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGa           >  1:261575/1‑50 (MQ=255)
                                   cTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATc         <  1:1888187/47‑1 (MQ=255)
                                   cTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGt      <  1:8324/50‑1 (MQ=255)
                                     tCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCg    >  1:828051/1‑50 (MQ=255)
                                      ccTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCg    <  1:1085155/49‑1 (MQ=255)
                                       cTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgtt  <  1:2030877/50‑1 (MQ=255)
                                              |                                          
ATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTT  >  NC_000913/2813318‑2813406

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