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,3640TC96.8% 108.8 / ‑3.6 31V250A (GTC→GCC) emrBmultidrug efflux pump membrane subunit EmrB
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (18/12);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTTGTTT  >  NC_000913/2813317‑2813410
                                               |                                              
tATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCt                                              >  1:402794/1‑50 (MQ=255)
tATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCt                                              >  1:1675168/1‑50 (MQ=255)
     tcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAg                                         <  1:48706/50‑1 (MQ=255)
      ccTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGc                                        >  1:1860101/1‑50 (MQ=255)
      ccTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGc                                        >  1:2796390/1‑50 (MQ=255)
         tACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGa                                     >  1:1514612/1‑50 (MQ=255)
          aCCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGa                                     <  1:556984/49‑1 (MQ=255)
          aCCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGAc                                    >  1:1244599/1‑50 (MQ=255)
           ccGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGAcc                                   >  1:751964/1‑50 (MQ=255)
                 tggCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACcgacga                              >  1:313152/1‑49 (MQ=255)
                 tggCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACcgacga                              >  1:1936915/1‑49 (MQ=255)
                   gCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATaa                           >  1:2552956/1‑50 (MQ=255)
                     ggtggtGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAcc                         >  1:1137514/1‑50 (MQ=255)
                     ggtggtGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAcc                         >  1:2492844/1‑50 (MQ=255)
                        ggtggCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGa                      >  1:868109/1‑50 (MQ=255)
                         gtggCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGAt                     <  1:407010/50‑1 (MQ=255)
                          tggCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATa                    <  1:472651/50‑1 (MQ=255)
                            gCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGt                  <  1:2618773/50‑1 (MQ=255)
                            gCTATCTGCCTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGt                  <  1:46639/50‑1 (MQ=39)
                             cTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTc                 >  1:902517/1‑50 (MQ=255)
                                tCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCg                <  1:658653/48‑1 (MQ=255)
                                tCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGAt              >  1:1969063/1‑50 (MQ=255)
                                  tGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCt            >  1:452042/1‑50 (MQ=255)
                                   gCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTg           >  1:1244974/1‑50 (MQ=255)
                                    cTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATc             <  1:692063/47‑1 (MQ=255)
                                       ccTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTc         <  1:173783/48‑1 (MQ=255)
                                       ccTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTc         <  1:1227041/48‑1 (MQ=255)
                                        cTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgtt      >  1:2499873/1‑50 (MQ=255)
                                         tGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGTTCTGTCgttg     >  1:2643880/1‑50 (MQ=39)
                                           aTTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgttg     <  1:648701/48‑1 (MQ=255)
                                            ttGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTTGttt  <  1:445169/50‑1 (MQ=255)
                                               |                                              
TATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTTGTTT  >  NC_000913/2813317‑2813410

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