Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,513 T→C 7.5% Q120Q (CAA→CAG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5130TC 7.5% 191.8 / 6.4 67Q120Q (CAA→CAGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base C (3/2);  total (35/32)
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.42e-01

GATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTA  >  minE/312444‑312582
                                                                     |                                                                     
gATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtctt                                                                      >  1:3118476/1‑71 (MQ=255)
   aaCGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCg                                                                   >  1:3789940/1‑71 (MQ=255)
     cGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCgg                                                                  >  1:631382/1‑70 (MQ=255)
     cGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgc                                                                 <  1:3168788/71‑1 (MQ=255)
      gCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCgg                                                                  >  1:1275268/1‑69 (MQ=255)
        cGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgcagc                                                              <  1:2261467/71‑1 (MQ=255)
          aCAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCaa                                                            >  1:2210745/1‑71 (MQ=255)
          aCAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCaa                                                            >  1:3850119/1‑71 (MQ=255)
          aCAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCaa                                                            >  1:2944373/1‑71 (MQ=255)
          aCAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCaa                                                            <  1:3189533/71‑1 (MQ=255)
             aCCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACgg                                                         >  1:2683621/1‑71 (MQ=255)
             aCCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACgg                                                         >  1:3066913/1‑71 (MQ=255)
              ccATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGc                                                        <  1:1405669/71‑1 (MQ=255)
               cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa                                                       <  1:129090/71‑1 (MQ=255)
               cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa                                                       >  1:3013521/1‑71 (MQ=255)
               cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa                                                       <  1:3217671/71‑1 (MQ=255)
                 tCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc                                                     <  1:4006811/71‑1 (MQ=255)
                 tCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc                                                     <  1:339736/71‑1 (MQ=255)
                     gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                   >  1:3920429/1‑69 (MQ=255)
                     gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                   >  1:2771831/1‑69 (MQ=255)
                     gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                   >  1:2661605/1‑69 (MQ=255)
                     gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                   >  1:2656583/1‑69 (MQ=255)
                     gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                   >  1:3842665/1‑69 (MQ=255)
                     gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                   >  1:1051993/1‑69 (MQ=255)
                       aGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGccct                                                >  1:2937306/1‑68 (MQ=255)
                       aGCTGCAGGATGATTTCTCAACGATCACCCCTTGCTGCTGAActtctt                                                                      >  1:1505627/1‑48 (MQ=255)
                         ctgctgTATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                             <  1:2824893/71‑1 (MQ=255)
                         ctgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                             <  1:1574427/71‑1 (MQ=255)
                          tgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgcagc                                                              <  1:1488798/53‑1 (MQ=255)
                          tgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgcagc                                                              <  1:3339193/53‑1 (MQ=255)
                                 tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                      <  1:2469988/70‑1 (MQ=255)
                                 tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                      <  1:46484/70‑1 (MQ=255)
                                 tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                     <  1:2443771/71‑1 (MQ=255)
                                 tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                     <  1:263968/71‑1 (MQ=255)
                                   atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt                                    <  1:3458726/70‑1 (MQ=255)
                                      tctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCAt                                                      >  1:1570324/1‑49 (MQ=255)
                                      tctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTg                                 <  1:3633194/70‑1 (MQ=255)
                                       ctcAACGCTCACCCCTTGCTGCTGGACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCa                               >  1:3768052/1‑71 (MQ=255)
                                       ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                                <  1:1285477/70‑1 (MQ=255)
                                       ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                                <  1:2309462/70‑1 (MQ=255)
                                       ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCa                               >  1:3244722/1‑71 (MQ=255)
                                          aaCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCt                            >  1:737035/1‑71 (MQ=255)
                                           aCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTg                           <  1:2484024/71‑1 (MQ=255)
                                            cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                             >  1:1001038/1‑52 (MQ=255)
                                            cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                          <  1:3679459/71‑1 (MQ=255)
                                                gatccccTGCTGCTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCAcc                             >  1:2370368/4‑64 (MQ=255)
                                                cACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt                                    >  1:2061855/1‑57 (MQ=255)
                                                  ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                    <  1:3909276/71‑1 (MQ=255)
                                                  ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                    >  1:2161056/1‑71 (MQ=255)
                                                  ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                    >  1:3584377/1‑71 (MQ=255)
                                                     cTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                  >  1:144054/1‑70 (MQ=255)
                                                      ttGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGttt                                       <  1:2385447/48‑1 (MQ=255)
                                                       tgctgcTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                      <  1:983296/48‑1 (MQ=255)
                                                       tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTg                                      <  1:2691196/48‑1 (MQ=19)
                                                        gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCt                                  >  1:2949376/1‑51 (MQ=255)
                                                        gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat               <  1:1439362/70‑1 (MQ=255)
                                                        gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACATGCGCGata              >  1:4121072/1‑71 (MQ=255)
                                                        gctgctGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGt                                         >  1:4055741/1‑44 (MQ=16)
                                                        gctgctGAACATCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGt                                         >  1:1610271/1‑44 (MQ=16)
                                                         ctgctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGat               <  1:3652408/69‑1 (MQ=255)
                                                           gctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATcc           >  1:3016428/1‑71 (MQ=255)
                                                            ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                        >  1:3629191/1‑57 (MQ=255)
                                                            ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat               <  1:882196/66‑1 (MQ=255)
                                                            ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg          >  1:2171676/1‑71 (MQ=255)
                                                              gAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                     <  1:510482/58‑1 (MQ=255)
                                                               aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc         <  1:1552607/69‑1 (MQ=255)
                                                               aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt       <  1:3180559/71‑1 (MQ=255)
                                                                aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCAc                              >  1:3997405/1‑47 (MQ=255)
                                                                aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt       >  1:4026884/1‑70 (MQ=255)
                                                                 cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat               >  1:1272280/1‑61 (MQ=255)
                                                                 cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat               >  1:1338591/1‑61 (MQ=255)
                                                                 cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat               >  1:2182012/1‑61 (MQ=255)
                                                                  ttcttGCGGCAGCAACGGCATCGCCAGCTGCAGtt                                        <  1:3660103/35‑1 (MQ=255)
                                                                   tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg    >  1:35390/1‑70 (MQ=255)
                                                                    cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                        >  1:1235457/1‑49 (MQ=255)
                                                                    cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                        >  1:1240221/1‑49 (MQ=255)
                                                                    cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa  >  1:92574/1‑71 (MQ=255)
                                                                     |                                                                     
GATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTA  >  minE/312444‑312582

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: