Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,495,805 A→G 46.7% A173A (GCT→GCC emrD ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,495,8050AG46.7% 9.6 / 91.3 105A173A (GCT→GCCemrDmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (42/14);  new base G (29/20);  total (71/34)
Fisher's exact test for biased strand distribution p-value = 9.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 105 total aligned reads displayed.

GGCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATCCAGCAGACCGCCG  >  minE/2495738‑2495871
                                                                   |                                                                  
ggCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAg                                                                   <  1:5132/69‑1 (MQ=255)
ggCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAg                                                                   <  1:222689/69‑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:3226441/69‑1 (MQ=255)
                                    gCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGcc                               <  1:963020/69‑1 (MQ=255)
                                    gCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGcc                               <  1:3222478/69‑1 (MQ=255)
                                      aTCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAg                             <  1:2747729/69‑1 (MQ=255)
                                      aTCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAg                             <  1:156088/69‑1 (MQ=255)
                                      aTCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAg                             <  1:2052741/69‑1 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                                           >  1:2130457/1‑52 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                                           >  1:2193301/1‑52 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                                           >  1:2343101/1‑52 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                                           >  1:2632045/1‑52 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                                           >  1:2668477/1‑52 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                                           >  1:278711/1‑52 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                                           >  1:1777501/1‑52 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                                           >  1:1658166/1‑52 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                                           >  1:1431967/1‑52 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGAt                                           >  1:1224803/1‑52 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGtt                           >  1:2388689/1‑68 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:3121353/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:3058501/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:2299709/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:1334713/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:460342/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:1525424/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:2112938/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:3042678/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:1534756/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:499243/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:1573857/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:1899884/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:2374508/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:3289840/1‑69 (MQ=255)
                                         cAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTc                          >  1:957537/1‑69 (MQ=255)
                                           gCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCca                        >  1:3087124/1‑69 (MQ=255)
                                           gCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCca                        >  1:2878614/1‑69 (MQ=255)
                                           gCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCca                        >  1:809541/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:2912159/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:1085166/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:2846673/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:2995633/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:3010094/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:3065942/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:3178733/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:448314/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:3269060/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:435823/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:3346340/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:3347539/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:2657431/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:1068872/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:1178878/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:1490561/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:155532/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:1593942/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:1629266/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:1697568/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:1792592/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:1950237/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:2098681/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:2323394/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:2333141/1‑69 (MQ=255)
                                                  aTACTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAACAGGCACGCCAGTTCCACATCGTa                 >  1:2437442/1‑69 (MQ=255)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATGAc                                        <  1:2880041/44‑1 (MQ=37)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:2567216/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:2815050/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:1177856/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:1178601/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:2877012/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:473862/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:262165/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:449111/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:399925/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:3389413/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:2316075/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:1667151/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:1862474/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:3018771/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:3195899/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:1962366/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:1983385/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:3113728/44‑1 (MQ=38)
                                                    aCTGAAAGTCACACCGGCACAAAGAACCAACAAAAAGAGATAAc                                        <  1:2233629/44‑1 (MQ=38)
                                                     cTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATcc              >  1:1990015/1‑69 (MQ=255)
                                                             cacaCCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGcc                               <  1:1183929/44‑1 (MQ=255)
                                                             cacaCCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATCCAGCAGAcc      >  1:3154017/1‑69 (MQ=255)
                                                                 ccAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATCCAGCAGAccgccg  >  1:1920469/1‑69 (MQ=255)
                                                                 ccAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATCCAGCAGAccgccg  >  1:2599431/1‑69 (MQ=255)
                                                                 ccAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATCCAGCAGAccgccg  >  1:1198538/1‑69 (MQ=255)
                                                                 ccAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATCCAGCAGAccgccg  >  1:951705/1‑69 (MQ=255)
                                                                 ccAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATCCAGCAGAccgccg  >  1:1092748/1‑69 (MQ=255)
                                                                   |                                                                  
GGCGCGTGCGCGGTGCATCGACCGGACGCGTTTCCGGCATCCAGCGGGCCATACTGAAGGTCACACCAGCACAAAGAACCAACAAAAAGAGATAACAGGCGCGCCAGTTCCACATCGTATCCAGCAGACCGCCG  >  minE/2495738‑2495871

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