Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,718 G→A 11.5% A385A (GCC→GCT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7180GA11.5% 266.6 / 17.9 95A385A (GCC→GCTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (46/38);  new base A (6/5);  total (52/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

GTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGA  >  minE/311649‑311788
                                                                     |                                                                      
gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                                                       <  1:2435204/71‑1 (MQ=255)
gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                                                       >  1:5049028/1‑71 (MQ=255)
 tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggc                                                                      <  1:4854868/71‑1 (MQ=255)
 tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggc                                                                      <  1:1420939/71‑1 (MQ=255)
   aTCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcgg                                                                    >  1:539611/1‑71 (MQ=255)
    tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                                   <  1:5510265/71‑1 (MQ=255)
    tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcgg                                                                    >  1:4525407/1‑70 (MQ=255)
      caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTGTCGAGAAGCCAAAg                                                                        >  1:4383075/1‑64 (MQ=255)
      caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                                   >  1:65230/1‑69 (MQ=255)
      caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAg                                                                        >  1:2174257/1‑64 (MQ=255)
       accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                                                <  1:3755918/71‑1 (MQ=255)
       accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                                                <  1:2964428/71‑1 (MQ=255)
         gagcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagc                                                                      <  1:600496/60‑1 (MQ=9)
         gagcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGa                                                              <  1:3877372/68‑1 (MQ=12)
         caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                              >  1:2970675/1‑71 (MQ=255)
          aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                             >  1:5483622/1‑71 (MQ=255)
          aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                             >  1:4240022/1‑71 (MQ=255)
          aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                             >  1:3244469/1‑71 (MQ=255)
          aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                             >  1:469348/1‑71 (MQ=255)
           aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                                                <  1:5028315/67‑1 (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)
                ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                                   <  1:3557547/59‑1 (MQ=255)
                 cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                      <  1:1094316/71‑1 (MQ=255)
                   aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                                   <  1:3707842/56‑1 (MQ=255)
                   aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                                   <  1:219989/56‑1 (MQ=255)
                   aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                                   <  1:4838550/56‑1 (MQ=255)
                    tGGCGAGCACCATCCCTAACATTGTTATCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                         <  1:4981848/65‑1 (MQ=255)
                     ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                       >  1:5263517/1‑66 (MQ=255)
                      gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                 <  1:1676517/71‑1 (MQ=255)
                      gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                 <  1:3276581/71‑1 (MQ=255)
                      gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                 <  1:4340837/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)
                         aGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                                             >  1:4923939/1‑55 (MQ=12)
                         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)
                            aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                              >  1:68116/1‑52 (MQ=255)
                            aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                              >  1:785974/1‑52 (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)
                                tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCaa                                                                 >  1:2579569/1‑45 (MQ=255)
                                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)
                                      aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                  >  1:4653480/1‑70 (MQ=255)
                                       cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                              >  1:3814848/1‑41 (MQ=255)
                                       cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                              >  1:2130522/1‑41 (MQ=255)
                                       cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                         <  1:1352518/46‑1 (MQ=255)
                                       cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                                   >  1:1447142/1‑68 (MQ=255)
                                       cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGAGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                                   >  1:1376820/1‑68 (MQ=255)
                                           gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                 >  1:1032302/1‑50 (MQ=255)
                                           gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                 >  1:1202656/1‑50 (MQ=255)
                                           gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                 >  1:4456264/1‑50 (MQ=255)
                                           gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACgg                             >  1:1724965/1‑70 (MQ=255)
                                            ttAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                              <  1:4969280/36‑1 (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:2100296/1‑39 (MQ=255)
                                                 gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                      >  1:753075/1‑39 (MQ=255)
                                                 gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                                      >  1:5462723/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)
                                                   gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCg                    >  1:3594483/1‑71 (MQ=255)
                                                    tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                       >  1:2446163/1‑35 (MQ=255)
                                                    tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                       >  1:1336855/1‑35 (MQ=255)
                                                    tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                       >  1:2547344/1‑35 (MQ=255)
                                                    tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                       >  1:5407435/1‑35 (MQ=255)
                                                    tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                       >  1:1771258/1‑35 (MQ=255)
                                                    tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                                       >  1:5392864/1‑35 (MQ=255)
                                                    tttATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                                  <  1:141612/56‑1 (MQ=255)
                                                      tATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                                 <  1:4050109/39‑1 (MQ=255)
                                                       aTCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAAt                >  1:2996679/1‑71 (MQ=255)
                                                        tCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATc               <  1:2303522/71‑1 (MQ=255)
                                                            gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCATCGGTACGGCaa                          >  1:893220/1‑56 (MQ=255)
                                                            gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg           <  1:2329970/71‑1 (MQ=255)
                                                            gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg           <  1:113161/71‑1 (MQ=255)
                                                             aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAAc            >  1:2228641/1‑69 (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)
                                                                     |                                                                      
GTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGA  >  minE/311649‑311788

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: