Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,624 G→A 6.5% D83D (GAC→GAT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,6240GA 6.5% 385.0 / 10.1 122D83D (GAC→GATacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (55/59);  new base A (6/2);  total (61/61)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 122 total aligned reads displayed.

CTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACC  >  minE/312555‑312693
                                                                     |                                                                     
cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                                                                      >  1:3256021/1‑71 (MQ=255)
 tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                                                                      >  1:1926290/1‑70 (MQ=255)
  gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc                                                                     >  1:2730804/1‑70 (MQ=255)
  gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa                                                                    >  1:1004204/1‑71 (MQ=255)
   aaCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc                                                                     <  1:1145281/69‑1 (MQ=255)
   aaCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCACGGTGATCTGCACGGTACCCGTGGAGTCa                                                                    <  1:2806951/70‑1 (MQ=255)
      cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                                                                <  1:3169355/71‑1 (MQ=255)
         cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                                                                >  1:3093525/1‑68 (MQ=255)
          gcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa                                                                    <  1:1074987/63‑1 (MQ=255)
           cgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTagag                                                           <  1:2077407/71‑1 (MQ=255)
             atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa                                                         >  1:2873614/1‑71 (MQ=255)
             atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa                                                         >  1:10652/1‑71 (MQ=255)
               atCCGCATCAGTACGAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACa                                                       >  1:33773/1‑71 (MQ=255)
                  cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg                                                     <  1:3094677/70‑1 (MQ=255)
                  cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg                                                     >  1:2575464/1‑70 (MQ=255)
                   gCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa                                                   >  1:1490258/1‑71 (MQ=255)
                    cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa                                                   <  1:3566946/70‑1 (MQ=255)
                    cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa                                                   >  1:1144665/1‑70 (MQ=255)
                    cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc                                                  <  1:2717478/71‑1 (MQ=255)
                    cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc                                                  <  1:2874417/71‑1 (MQ=255)
                     aTCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTaga                                                            <  1:203686/60‑1 (MQ=255)
                     aTCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTaga                                                            <  1:919007/60‑1 (MQ=255)
                        aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTagag                                                           >  1:1795913/1‑58 (MQ=255)
                        aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa                                              <  1:2076713/71‑1 (MQ=255)
                        aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa                                              <  1:2134927/71‑1 (MQ=255)
                         gTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGt                                                    <  1:1073532/64‑1 (MQ=255)
                          tACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAgg                                            <  1:1200755/71‑1 (MQ=255)
                            ccAGACTCAAAGGTCAGGGTGATCTGCACGGTACGCGTGGAGTCACTGTTAGAGGACATGTACATCAGGtt                                          <  1:2155597/71‑1 (MQ=255)
                             cAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa                                              <  1:1891773/66‑1 (MQ=255)
                             cAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa                                              <  1:2557084/66‑1 (MQ=255)
                              aGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt                                        <  1:2655081/71‑1 (MQ=255)
                               gACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                       >  1:3319888/1‑71 (MQ=255)
                                aCTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                       >  1:93835/1‑70 (MQ=255)
                                aCTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                       >  1:3163857/1‑70 (MQ=255)
                                 cTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACt                                                                  <  1:2529929/42‑1 (MQ=255)
                                 cTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGa                                     <  1:2061349/71‑1 (MQ=255)
                                  tCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACa                                                 <  1:1175603/58‑1 (MQ=255)
                                  tCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGtt                                          <  1:2767310/65‑1 (MQ=255)
                                    aaaGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt                                        >  1:1652470/1‑65 (MQ=255)
                                    aaaGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGa                                     <  1:1442003/68‑1 (MQ=255)
                                    aaaGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATa                                   <  1:1621565/70‑1 (MQ=255)
                                    aaaGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATa                                   <  1:1616987/70‑1 (MQ=255)
                                    aaaGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATa                                   <  1:1521077/70‑1 (MQ=255)
                                    aaaGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATAc                                  >  1:3368615/1‑71 (MQ=255)
                                    aaaGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATAc                                  <  1:3728057/71‑1 (MQ=255)
                                     taGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt                                        <  1:2371790/63‑1 (MQ=255)
                                     aaGGTAAGGGTAATTGTCACGCTACCGGCGGAATCGCTGGTGGAGGACATATACATCAGGTTATCGATAcc                                 <  1:2900353/71‑1 (MQ=255)
                                      aGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCg                                >  1:76873/1‑71 (MQ=255)
                                       ggTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGt                               >  1:2738827/1‑71 (MQ=255)
                                       ggTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGt                               <  1:1404292/71‑1 (MQ=255)
                                           aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                       >  1:3058917/1‑59 (MQ=255)
                                           aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                       >  1:3008188/1‑59 (MQ=255)
                                           aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                       >  1:981753/1‑59 (MQ=255)
                                           aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                       >  1:2498347/1‑59 (MQ=255)
                                           aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                             <  1:2170609/69‑1 (MQ=255)
                                           aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCa                            >  1:3488319/1‑70 (MQ=255)
                                            gggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAg                                             >  1:3177005/1‑52 (MQ=255)
                                            gggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCata                          <  1:941128/71‑1 (MQ=255)
                                            gggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCata                          >  1:1217826/1‑71 (MQ=255)
                                             ggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAgg                                                          >  1:2908011/1‑38 (MQ=255)
                                             ggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa                                                         >  1:579410/1‑39 (MQ=255)
                                             ggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa                                                         >  1:562294/1‑39 (MQ=255)
                                             ggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATa                                   <  1:386238/61‑1 (MQ=255)
                                              gTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGa                                     >  1:2172799/1‑58 (MQ=255)
                                                 aTCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTa                                         <  1:688140/51‑1 (MQ=255)
                                                 aTCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTa                                         <  1:3700416/51‑1 (MQ=255)
                                                  tCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCg                                >  1:104342/1‑59 (MQ=255)
                                                  tCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                    >  1:1834067/1‑71 (MQ=255)
                                                   cTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc                                                  >  1:431966/1‑40 (MQ=255)
                                                   cTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                    >  1:2190386/1‑70 (MQ=255)
                                                     gCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                             >  1:3578790/1‑59 (MQ=255)
                                                     gCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTc                  >  1:3421168/1‑70 (MQ=255)
                                                     gCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTc                  >  1:3525639/1‑70 (MQ=255)
                                                      cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATa                                   >  1:737938/1‑52 (MQ=255)
                                                      cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCt                      >  1:2218531/1‑65 (MQ=255)
                                                      cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTg                     <  1:738254/66‑1 (MQ=255)
                                                      cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                    <  1:1629868/67‑1 (MQ=255)
                                                      cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                    <  1:1882583/67‑1 (MQ=255)
                                                      cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                    <  1:2435987/67‑1 (MQ=255)
                                                      cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                    <  1:2255692/67‑1 (MQ=255)
                                                      cACGCTACCGGCGGAATCGCTGGTGGAGGACATATACATCAGGTTATCGATACCGtt                              >  1:3112617/1‑57 (MQ=255)
                                                      cACGCTACCGGCGGAATCGCTGGTGGAGGACATATACATCAGGTTATCGATACCGtt                              >  1:2158999/1‑57 (MQ=255)
                                                       aCGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGAt               <  1:3470332/71‑1 (MQ=255)
                                                       aCGCTACCGGCGGAATCGCTGGTGGAGGACATATACATCAGGTTATAGATACCGTTCATATTCTGTTCGAt               <  1:3909656/71‑1 (MQ=19)
                                                          gTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAAc            >  1:540513/1‑71 (MQ=255)
                                                          gTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAAc            >  1:2965456/1‑71 (MQ=255)
                                                           tACCGGCGGAATCGCTGGTGGAGGACATATACATCAGGTTATCGATa                                   >  1:2835802/1‑47 (MQ=16)
                                                             cccGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCtg         <  1:3601211/71‑1 (MQ=255)
                                                             cccGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCt          <  1:156924/70‑1 (MQ=255)
                                                              ccGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                    <  1:2803940/59‑1 (MQ=255)
                                                              ccGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCtgt        <  1:602086/71‑1 (MQ=255)
                                                              ccGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCtgt        <  1:2066195/71‑1 (MQ=255)
                                                                gcggAATCGCTGGTGGAGGACATATACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTgcgt      >  1:1941906/3‑71 (MQ=21)
                                                                gcggAATCGCTGGTGGAGGACATATACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTgcgt      >  1:408035/3‑71 (MQ=21)
                                                                gTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTgtgt      <  1:3464463/71‑1 (MQ=255)
                                                                gTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTgtgt      <  1:1813060/71‑1 (MQ=255)
                                                                 cggAATCGCTGGTGGAGGACATATACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTgcgt      >  1:1819625/2‑70 (MQ=255)
                                                                   gAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCa    <  1:2885434/70‑1 (MQ=255)
                                                                    aGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCAcc  <  1:3109743/71‑1 (MQ=255)
                                                                    aGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCAcc  <  1:802943/71‑1 (MQ=255)
                                                                     |                                                                     
CTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACC  >  minE/312555‑312693

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: