Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 310,785 A→G 7.3% T696T (ACT→ACC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE310,7850AG 7.3% 193.2 / 4.6 68T696T (ACT→ACCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (26/37);  new base G (3/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTT  >  minE/310715‑310855
                                                                      |                                                                      
cGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGa                                                                        <  1:4411969/71‑1 (MQ=255)
    tGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGt                                                                      >  1:1446769/1‑69 (MQ=255)
        cGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                                <  1:4757332/71‑1 (MQ=255)
        cGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                                <  1:3711164/71‑1 (MQ=255)
         gTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                                               <  1:2909431/71‑1 (MQ=255)
         gTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                                               <  1:5607459/71‑1 (MQ=255)
            cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                                               <  1:4161674/68‑1 (MQ=255)
            cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                            >  1:3162654/1‑71 (MQ=255)
            cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                            >  1:3586773/1‑71 (MQ=255)
             gCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                           <  1:1106219/71‑1 (MQ=255)
             gCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                           <  1:4482940/71‑1 (MQ=255)
             gCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                           <  1:5023481/71‑1 (MQ=255)
                 gTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                       <  1:295943/71‑1 (MQ=255)
                   cAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCGTGAGTCAGTTTTTCGTGAcc                                                       >  1:280878/1‑69 (MQ=255)
                   cAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                       <  1:3969444/69‑1 (MQ=255)
                   cAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                       >  1:4977237/1‑69 (MQ=255)
                      cATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGc                                                  <  1:5482814/71‑1 (MQ=255)
                      cATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGc                                                  <  1:83119/71‑1 (MQ=255)
                        tatCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                                <  1:194227/71‑1 (MQ=255)
                            aGGGTGCTTCGCTGCTTCTGCAAGCAACTTGTTACGCGCCTGa                                                                        <  1:3392592/43‑1 (MQ=255)
                               gTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                       <  1:5345261/57‑1 (MQ=255)
                                tGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTc                                                                     >  1:160094/1‑42 (MQ=255)
                                tGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTc                                                                     >  1:3843480/1‑42 (MQ=255)
                                  cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                                      >  1:5678351/1‑71 (MQ=255)
                                  cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                                      <  1:4248127/71‑1 (MQ=255)
                                  cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                                       <  1:4569373/70‑1 (MQ=255)
                                  cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                                       <  1:2166150/70‑1 (MQ=255)
                                    tCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                                <  1:1571864/43‑1 (MQ=255)
                                    tCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                                <  1:3740370/43‑1 (MQ=255)
                                    tCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                                <  1:3167609/43‑1 (MQ=255)
                                     cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                     >  1:2705732/1‑69 (MQ=255)
                                     cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                     >  1:1054075/1‑69 (MQ=255)
                                     cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                     >  1:1785531/1‑69 (MQ=255)
                                       ctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                             <  1:4555066/43‑1 (MQ=255)
                                        tgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                >  1:3442614/1‑71 (MQ=255)
                                         gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                           >  1:1755266/1‑43 (MQ=255)
                                         gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                >  1:1410835/1‑70 (MQ=255)
                                           ttCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                              >  1:2339927/1‑70 (MQ=255)
                                           ttCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                              >  1:4169904/1‑70 (MQ=255)
                                            tCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGcc                                             >  1:1967700/1‑54 (MQ=255)
                                             cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAc                                                        <  1:1587239/42‑1 (MQ=255)
                                             cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa                            <  1:1220595/70‑1 (MQ=255)
                                             cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                             <  1:817166/69‑1 (MQ=255)
                                              tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                            >  1:3649713/1‑37 (MQ=255)
                                              tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                            >  1:1733112/1‑37 (MQ=255)
                                              tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                       <  1:1351816/42‑1 (MQ=255)
                                              tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                       <  1:70925/42‑1 (MQ=255)
                                              tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                       <  1:5243633/42‑1 (MQ=255)
                                               gcaagcaaCTTGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                                      >  1:1807289/1‑42 (MQ=255)
                                               gcaagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                                      >  1:2489477/1‑42 (MQ=255)
                                               ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc                         <  1:1626494/70‑1 (MQ=255)
                                                 aagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga                       <  1:999838/71‑3 (MQ=255)
                                                 aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa                       <  1:1283897/71‑1 (MQ=255)
                                                 aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa                       <  1:2247475/71‑1 (MQ=255)
                                                  agcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga                       >  1:1523076/1‑68 (MQ=255)
                                                   gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcgaaa                     >  1:1711755/1‑70 (MQ=255)
                                                   gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcgaaa                     >  1:4088067/1‑70 (MQ=255)
                                                   gcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag                     >  1:425769/1‑71 (MQ=255)
                                                    caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc                    >  1:4575506/1‑71 (MQ=255)
                                                    caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc                    >  1:2708401/1‑71 (MQ=255)
                                                     aaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                   <  1:5341281/71‑1 (MQ=255)
                                                       cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                           <  1:5583754/45‑1 (MQ=255)
                                                       cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                     >  1:436058/1‑51 (MQ=255)
                                                          gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                             <  1:2018286/56‑1 (MQ=255)
                                                          gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                             <  1:835662/56‑1 (MQ=255)
                                                          gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg              <  1:4700272/71‑1 (MQ=255)
                                                          gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg              <  1:4797019/71‑1 (MQ=255)
                                                            tACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg              >  1:4778008/1‑69 (MQ=255)
                                                                  ccTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTAccc      <  1:3974914/71‑1 (MQ=255)
                                                                  ccTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTAccc      <  1:577185/71‑1 (MQ=255)
                                                                      aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGtt  <  1:2198278/71‑1 (MQ=255)
                                                                      |                                                                      
CGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTT  >  minE/310715‑310855

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: