Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,495,796 G→A 48.0% T176T (ACC→ACT emrD ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,495,7960GA48.0% 7.0 / 89.4 102T176T (ACC→ACTemrDmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (37/16);  new base A (29/20);  total (66/36)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

GGTAAGCAGGCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATCC  >  minE/2495730‑2495859
                                                                  |                                                               
ggTAAGCAGGCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGt                                                               <  1:215951/69‑1 (MQ=255)
  tAAGCAGGCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTca                                                             >  1:303725/1‑69 (MQ=255)
      cAGGCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCacac                                                          <  1:507309/68‑1 (MQ=255)
      cAGGCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCacac                                                          <  1:389577/68‑1 (MQ=255)
        ggCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAg                                                       <  1:222689/69‑1 (MQ=255)
        ggCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAg                                                       <  1:5132/69‑1 (MQ=255)
         gcgcGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAg                                                       <  1:2849153/68‑1 (MQ=255)
                     tGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACca                                           >  1:306041/1‑68 (MQ=255)
                              cGGACGCGTTTCCGGCATCCAGCGGTCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAagag                                 <  1:666969/69‑1 (MQ=255)
                              cGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAagag                                 <  1:3015757/69‑1 (MQ=255)
                               ggACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAagag                                 <  1:1929263/68‑1 (MQ=255)
                                            gCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGcc                   <  1:963020/69‑1 (MQ=255)
                                            gCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGcc                   <  1:2220907/69‑1 (MQ=255)
                                            gCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGcc                   <  1:3222478/69‑1 (MQ=255)
                                            gCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGcc                   <  1:3226441/69‑1 (MQ=255)
                                              aTCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAg                 <  1:2052741/69‑1 (MQ=255)
                                              aTCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAg                 <  1:2747729/69‑1 (MQ=255)
                                              aTCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAg                 <  1:156088/69‑1 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                               >  1:2668477/1‑52 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                               >  1:2130457/1‑52 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                               >  1:1431967/1‑52 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                               >  1:2193301/1‑52 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                               >  1:1224803/1‑52 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                               >  1:2343101/1‑52 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                               >  1:1777501/1‑52 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                               >  1:278711/1‑52 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                               >  1:1658166/1‑52 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                               >  1:2632045/1‑52 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGtt               >  1:2388689/1‑68 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:2374508/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:957537/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:2299709/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:1334713/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:3042678/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:499243/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:460342/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:1573857/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:3289840/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:1899884/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:1534756/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:3121353/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:1525424/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:3058501/1‑69 (MQ=255)
                                                 cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc              >  1:2112938/1‑69 (MQ=255)
                                                   gCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCca            >  1:809541/1‑69 (MQ=255)
                                                   gCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCca            >  1:3087124/1‑69 (MQ=255)
                                                   gCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCca            >  1:2878614/1‑69 (MQ=255)
                                                   gCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCca            >  1:1722535/1‑69 (MQ=255)
                                                   gCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCca            >  1:3352675/1‑69 (MQ=255)
                                                        ccATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCg       >  1:2079763/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCg       >  1:1555726/1‑67 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGt      >  1:9663/1‑68 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:1593942/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:1657540/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:435823/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:3347539/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:3346340/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:1697568/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:3269060/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:448314/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:1792592/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:2657431/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:3178733/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:1950237/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:1629266/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:3065942/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:155532/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:2912159/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:1068872/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:2437442/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:1085166/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:2846673/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:2333141/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:2323394/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:2098681/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:2995633/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:3010094/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:1178878/1‑69 (MQ=255)
                                                          aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa     >  1:1490561/1‑69 (MQ=255)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATGAc                            <  1:2880041/44‑1 (MQ=37)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:449111/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:1177856/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:1178601/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:473862/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:1862474/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:399925/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:3389413/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:1667151/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:3195899/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:3113728/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:1962366/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:1983385/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:2233629/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:3018771/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:2316075/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:2877012/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:2815050/44‑1 (MQ=38)
                                                            aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                            <  1:2567216/44‑1 (MQ=38)
                                                             cTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATcc  >  1:3233897/1‑69 (MQ=255)
                                                             cTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATcc  >  1:1990015/1‑69 (MQ=255)
                                                                  |                                                               
GGTAAGCAGGCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATCC  >  minE/2495730‑2495859

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