Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,631 T→G 6.5% N81T (AAC→ACC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,6310TG 6.5% 327.9 / 15.3 124N81T (AAC→ACC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (51/65);  new base G (6/2);  total (57/67)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 131 total aligned reads displayed.

CTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTG  >  minE/312561‑312701
                                                                      |                                                                      
cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                                                                        <  1:3169355/71‑1 (MQ=255)
cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                                                                        <  1:2635633/71‑1 (MQ=255)
     cgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTagag                                                                   <  1:2077407/71‑1 (MQ=255)
       atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa                                                                 >  1:10652/1‑71 (MQ=255)
       atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa                                                                 >  1:2873614/1‑71 (MQ=255)
         atCCGCATCAGTACGAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACa                                                               >  1:33773/1‑71 (MQ=255)
            cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg                                                             <  1:3094677/70‑1 (MQ=255)
             gCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa                                                           >  1:1490258/1‑71 (MQ=255)
              cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa                                                           >  1:1144665/1‑70 (MQ=255)
              cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa                                                           <  1:3566946/70‑1 (MQ=255)
              cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc                                                          <  1:2874417/71‑1 (MQ=255)
              cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc                                                          <  1:2717478/71‑1 (MQ=255)
               aTCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTaga                                                                    <  1:203686/60‑1 (MQ=255)
                  aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa                                                      <  1:2076713/71‑1 (MQ=255)
                   gTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGt                                                            <  1:1073532/64‑1 (MQ=255)
                    tACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAgg                                                    <  1:640487/71‑1 (MQ=255)
                      ccAGACTCAAAGGTCAGGGTGATCTGCACGGTACGCGTGGAGTCACTGTTAGAGGACATGTACATCAGGtt                                                  <  1:2155597/71‑1 (MQ=255)
                       cAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa                                                      <  1:2557084/66‑1 (MQ=255)
                       cAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa                                                      <  1:1891773/66‑1 (MQ=255)
                        aGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc                                                          >  1:1876565/1‑61 (MQ=255)
                        aGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt                                                <  1:2655081/71‑1 (MQ=255)
                        aGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt                                                >  1:51735/1‑71 (MQ=255)
                        aGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt                                                <  1:2420090/71‑1 (MQ=255)
                         gACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt                                                >  1:1523062/1‑70 (MQ=255)
                         gACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                               <  1:2758741/71‑1 (MQ=255)
                         gACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                               >  1:3319888/1‑71 (MQ=255)
                          aCTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                               >  1:3163857/1‑70 (MQ=255)
                          aCTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                               >  1:2039442/1‑70 (MQ=255)
                          aCTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                               >  1:93835/1‑70 (MQ=255)
                           cTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGa                                             <  1:2061349/71‑1 (MQ=255)
                           cTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGa                                             <  1:2119506/71‑1 (MQ=255)
                            tCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACa                                                         <  1:1175603/58‑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)
                              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)
                                    cAGTGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                                     <  1:830142/70‑1 (MQ=255)
                                    cAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCa                                    >  1:1235896/1‑71 (MQ=255)
                                     aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                               >  1:981753/1‑59 (MQ=255)
                                     aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                               >  1:3058917/1‑59 (MQ=255)
                                     aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                               >  1:3008188/1‑59 (MQ=255)
                                     aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                               >  1:2498347/1‑59 (MQ=255)
                                     aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                                     <  1:2170609/69‑1 (MQ=255)
                                      gggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAg                                                     >  1:3177005/1‑52 (MQ=255)
                                      gggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCata                                  >  1:1217826/1‑71 (MQ=255)
                                      gggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCata                                  <  1:941128/71‑1 (MQ=255)
                                       ggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAgg                                                                  >  1:2908011/1‑38 (MQ=255)
                                       ggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa                                                                 >  1:579410/1‑39 (MQ=255)
                                       ggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATa                                           <  1:386238/61‑1 (MQ=255)
                                           aTCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTg                             >  1:3101348/1‑71 (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:1137642/1‑59 (MQ=255)
                                               gCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTc                          >  1:3525639/1‑70 (MQ=255)
                                               gCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTc                          >  1:3421168/1‑70 (MQ=255)
                                                cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGa                                             >  1:2179881/1‑50 (MQ=255)
                                                cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATa                                           >  1:737938/1‑52 (MQ=255)
                                                cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCt                              >  1:2218531/1‑65 (MQ=255)
                                                cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                            <  1:2435987/67‑1 (MQ=255)
                                                cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                            <  1:2255692/67‑1 (MQ=255)
                                                cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                            <  1:1882583/67‑1 (MQ=255)
                                                cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTc                          <  1:3539866/69‑1 (MQ=255)
                                                cACGCTACCGGCGGAATCGCTGGTGGAGGACATATACATCAGGTTATCGATACCGtt                                      >  1:2158999/1‑57 (MQ=255)
                                                 aCGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGAt                       <  1:3470332/71‑1 (MQ=255)
                                                    gTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAAc                    >  1:540513/1‑71 (MQ=255)
                                                     tACCGGCGGAATCGCTGGTGGAGGACATATACATCAGGTTATCGATa                                           >  1:2835802/1‑47 (MQ=16)
                                                     tACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAAcc                   >  1:3726382/1‑71 (MQ=255)
                                                       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)
                                                          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)
                                                             gAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCAc           >  1:1773381/1‑71 (MQ=255)
                                                              aGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCAcc          <  1:802943/71‑1 (MQ=255)
                                                              aGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCAcc          <  1:3109743/71‑1 (MQ=255)
                                                                tCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                                     <  1:3598410/42‑1 (MQ=255)
                                                                 cACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGtt                           <  1:2643291/51‑1 (MQ=255)
                                                                  aCTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCgtg       <  1:1065346/70‑1 (MQ=255)
                                                                  aCTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCgtg       <  1:3068805/70‑1 (MQ=255)
                                                                    tGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                                     <  1:2060816/38‑1 (MQ=255)
                                                                    tGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                                     <  1:3591114/38‑1 (MQ=255)
                                                                    tGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTcc    >  1:1747006/1‑71 (MQ=255)
                                                                     gTTAGAGGACATGTACATCAGGTTATCGATACCGTTCata                                  <  1:729273/40‑1 (MQ=255)
                                                                     gTTAGAGGACATGTACATCAGGTTATCGATACCGTTCata                                  <  1:730413/40‑1 (MQ=255)
                                                                      ttAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTcc    <  1:3427591/69‑1 (MQ=255)
                                                                      ttAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTg  >  1:1333341/1‑71 (MQ=255)
                                                                      ttAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTg  <  1:2326700/71‑1 (MQ=255)
                                                                      |                                                                      
CTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTG  >  minE/312561‑312701

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: