Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,729 C→T 9.9% V382I (GTC→ATC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7290CT 9.9% 264.6 / 13.1 91V382I (GTC→ATC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (44/38);  new base T (5/4);  total (49/42)
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.90e-01

AACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAG  >  minE/311659‑311798
                                                                      |                                                                     
aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                                       >  1:5483622/1‑71 (MQ=255)
aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                                       >  1:469348/1‑71 (MQ=255)
aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                                       >  1:4240022/1‑71 (MQ=255)
aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                                       >  1:3244469/1‑71 (MQ=255)
   aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                                                                    <  1:841289/71‑1 (MQ=255)
    ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                                   <  1:572928/71‑1 (MQ=255)
    ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                                   <  1:4191323/71‑1 (MQ=255)
     cccGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCaa                                                                  <  1:1437865/70‑1 (MQ=16)
       cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                                <  1:1094316/71‑1 (MQ=255)
         aTGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGt                                                              <  1:1344227/71‑1 (MQ=12)
          tGGCGAGCACCATCCCTAACATTGTTATCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                                   <  1:4981848/65‑1 (MQ=255)
           ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                                 >  1:5263517/1‑66 (MQ=255)
            gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                           <  1:3276581/71‑1 (MQ=255)
            gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                           <  1:4340837/71‑1 (MQ=255)
            gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                           <  1:1676517/71‑1 (MQ=255)
            gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                           <  1:1865527/71‑1 (MQ=255)
              aaGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTccct                                                          >  1:85032/2‑69 (MQ=12)
               aGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAg                                                        >  1:2721745/1‑71 (MQ=255)
               aGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTccctaa                                                        >  1:1227868/1‑68 (MQ=12)
                gCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAg                                                        >  1:3962273/1‑70 (MQ=255)
                 cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                                                            >  1:4627673/1‑65 (MQ=255)
                 cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                                                       >  1:1881803/1‑70 (MQ=255)
                  aCCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaa                                                     <  1:1302844/71‑1 (MQ=12)
                   ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCa                                                      >  1:2964362/1‑69 (MQ=255)
                   ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCa                                                      >  1:2161714/1‑69 (MQ=255)
                   ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAAt                                                    <  1:4201267/71‑1 (MQ=255)
                   ccATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaac                                                    <  1:1962639/71‑2 (MQ=12)
                    cATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                                                       <  1:1733292/67‑1 (MQ=255)
                    cATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTc                                                             >  1:128890/1‑61 (MQ=12)
                    cATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaaca                                                   >  1:4166017/1‑69 (MQ=12)
                      tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                                                  <  1:4952290/54‑1 (MQ=255)
                      tCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaaca                                                   >  1:91188/1‑67 (MQ=12)
                        ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                                                       <  1:4266106/63‑1 (MQ=255)
                          gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                                             >  1:4877978/1‑71 (MQ=255)
                           aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                                                       <  1:2224202/60‑1 (MQ=255)
                           aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                            <  1:1693640/71‑1 (MQ=255)
                            aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                            >  1:4653480/1‑70 (MQ=255)
                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                                   <  1:1352518/46‑1 (MQ=255)
                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                                             >  1:1447142/1‑68 (MQ=255)
                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGAGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                                             >  1:1376820/1‑68 (MQ=255)
                              aTTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTAc                                         <  1:4654934/71‑1 (MQ=255)
                                 gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                           >  1:4456264/1‑50 (MQ=255)
                                 gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                           >  1:1202656/1‑50 (MQ=255)
                                 gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                           >  1:1032302/1‑50 (MQ=255)
                                 gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                           >  1:727980/1‑50 (MQ=255)
                                 gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACgg                                       >  1:1724965/1‑70 (MQ=255)
                                   tAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCaa                                    >  1:196454/1‑71 (MQ=255)
                                      cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                                                  <  1:800495/38‑1 (MQ=255)
                                      cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                                                  <  1:4855722/38‑1 (MQ=255)
                                       gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                                >  1:4259160/1‑39 (MQ=255)
                                       gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                                >  1:2100296/1‑39 (MQ=255)
                                       gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                                >  1:5462723/1‑39 (MQ=255)
                                       gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                                >  1:753075/1‑39 (MQ=255)
                                        tgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                                                       <  1:1681563/47‑1 (MQ=255)
                                         gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgt                                            <  1:4398025/57‑2 (MQ=255)
                                         gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                            <  1:3540056/57‑1 (MQ=255)
                                         gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                            <  1:2718617/57‑1 (MQ=255)
                                         gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCg                              >  1:3594483/1‑71 (MQ=255)
                                          tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                                 >  1:2446163/1‑35 (MQ=255)
                                          tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                                 >  1:2547344/1‑35 (MQ=255)
                                          tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                                 >  1:1771258/1‑35 (MQ=255)
                                          tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                                 >  1:5407435/1‑35 (MQ=255)
                                          tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                                 >  1:5392864/1‑35 (MQ=255)
                                          tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                                 >  1:1336855/1‑35 (MQ=255)
                                          tttATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                            <  1:141612/56‑1 (MQ=255)
                                          tttATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCgg                             <  1:494426/71‑1 (MQ=255)
                                            tATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                           <  1:4050109/39‑1 (MQ=255)
                                             aTCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAAt                          >  1:2996679/1‑71 (MQ=255)
                                              tCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATc                         <  1:2303522/71‑1 (MQ=255)
                                                 agaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCa                        <  1:208985/69‑1 (MQ=255)
                                                  gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCATCGGTACGGCaa                                    >  1:893220/1‑56 (MQ=255)
                                                  gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                     <  1:113161/71‑1 (MQ=255)
                                                  gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg                     <  1:2329970/71‑1 (MQ=255)
                                                   aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAAc                      >  1:2228641/1‑69 (MQ=255)
                                                   aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGt                    >  1:2290489/1‑71 (MQ=255)
                                                   aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAAAACCACCGGTACGGCAATGGTCGGAATCaa                       >  1:2715058/1‑68 (MQ=255)
                                                       cAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcg                >  1:2071402/1‑71 (MQ=255)
                                                          aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCaa                                    <  1:3393648/48‑1 (MQ=255)
                                                          aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCaa                                    <  1:3895571/48‑1 (MQ=255)
                                                          aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcgcg              >  1:4846186/1‑70 (MQ=255)
                                                          aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCgg             <  1:3663923/71‑1 (MQ=255)
                                                           ggcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGa            <  1:2655864/71‑1 (MQ=255)
                                                            gcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATg                                  >  1:2210611/1‑48 (MQ=255)
                                                            gcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATg                                  >  1:778627/1‑48 (MQ=255)
                                                            gcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCACCGTCGCGCgg             >  1:3497152/1‑69 (MQ=255)
                                                            gcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCgg             >  1:991608/1‑69 (MQ=255)
                                                                cAAGGACGGCAAAGGGCCCGAGCAATACCACCGGTACGGCAATGGTCGGa                            <  1:4347077/50‑1 (MQ=255)
                                                                 aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt      <  1:1388741/71‑1 (MQ=255)
                                                                  aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg     <  1:1954891/71‑1 (MQ=255)
                                                                  aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg     >  1:2082518/1‑71 (MQ=255)
                                                                  aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg     >  1:2930044/1‑71 (MQ=255)
                                                                  aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGGCGGAATCAACGTCGCGCGGAAGTTCTg     <  1:1104062/71‑1 (MQ=255)
                                                                  aGGACGGCAAAGGTCCCGAGAAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg     >  1:3748501/1‑71 (MQ=255)
                                                                     aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg  >  1:612465/1‑71 (MQ=255)
                                                                     aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg  >  1:3481233/1‑71 (MQ=255)
                                                                      |                                                                     
AACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAG  >  minE/311659‑311798

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: