Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,633 A→G 6.4% S80S (TCT→TCC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,6330AG 6.4% 360.7 / 10.4 125S80S (TCT→TCCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (52/65);  new base G (6/2);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

CGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCA  >  minE/312566‑312703
                                                                   |                                                                      
cgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTagag                                                                     <  1:2077407/71‑1 (MQ=255)
  atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa                                                                   >  1:10652/1‑71 (MQ=255)
    atCCGCATCAGTACGAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACa                                                                 >  1:33773/1‑71 (MQ=255)
        gCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa                                                             >  1:3178849/1‑71 (MQ=255)
         cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa                                                             >  1:1144665/1‑70 (MQ=255)
         cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa                                                             <  1:3566946/70‑1 (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:3408802/71‑1 (MQ=255)
             aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa                                                        <  1:2134927/71‑1 (MQ=255)
               tACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAgg                                                      <  1:1200755/71‑1 (MQ=255)
               tACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAgg                                                      <  1:640487/71‑1 (MQ=255)
                  cAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa                                                        <  1:1891773/66‑1 (MQ=255)
                  cAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa                                                        <  1:2557084/66‑1 (MQ=255)
                   aGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc                                                            >  1:1876565/1‑61 (MQ=255)
                   aGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt                                                  >  1:51735/1‑71 (MQ=255)
                   aGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt                                                  <  1:2655081/71‑1 (MQ=255)
                    gACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTAt                                                  >  1:1523062/1‑70 (MQ=255)
                    gACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                                 >  1:3319888/1‑71 (MQ=255)
                     aCTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                                 >  1:93835/1‑70 (MQ=255)
                      cTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGa                                               <  1:2119506/71‑1 (MQ=255)
                      cTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGa                                               <  1:2061349/71‑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:1521077/70‑1 (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)
                               cAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATa                                             >  1:2874364/1‑64 (MQ=255)
                               cAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCa                                      >  1:1235896/1‑71 (MQ=255)
                                aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                                 >  1:3058917/1‑59 (MQ=255)
                                aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                                 >  1:3008188/1‑59 (MQ=255)
                                aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                                 >  1:2498347/1‑59 (MQ=255)
                                aGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATc                                                 >  1:981753/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:1217826/1‑71 (MQ=255)
                                  ggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAgg                                                                    >  1:2908011/1‑38 (MQ=255)
                                  ggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa                                                                   >  1:562294/1‑39 (MQ=255)
                                  ggTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa                                                                   >  1:579410/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)
                                      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:3578790/1‑59 (MQ=255)
                                          gCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                                       >  1:1137642/1‑59 (MQ=255)
                                          gCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTc                            >  1:3421168/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:2435987/67‑1 (MQ=255)
                                           cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGt                              <  1:1629868/67‑1 (MQ=255)
                                           cACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTc                            <  1:3539866/69‑1 (MQ=255)
                                           cACGCTACCGGCGGAATCGCTGGTGGAGGACATATACATCAGGTTATCGATACCGtt                                        >  1:3112617/1‑57 (MQ=255)
                                           cACGCTACCGGCGGAATCGCTGGTGGAGGACATATACATCAGGTTATCGATACCGtt                                        >  1:2158999/1‑57 (MQ=255)
                                            aCGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGAt                         <  1:3470332/71‑1 (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: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)
                                                       ggAGTCACTGTTAGAGGACATGTACATCAGGTTATCGa                                               >  1:990760/1‑38 (MQ=255)
                                                        gAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCAc             >  1:1773381/1‑71 (MQ=255)
                                                         aGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCAcc            <  1:3109743/71‑1 (MQ=255)
                                                           tCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                                       <  1:3598410/42‑1 (MQ=255)
                                                           tCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATAtt                                  <  1:1856524/47‑1 (MQ=255)
                                                            cACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGtt                             <  1:2643291/51‑1 (MQ=255)
                                                             aCTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCgtg         <  1:1065346/70‑1 (MQ=255)
                                                             aCTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCgtg         <  1:3068805/70‑1 (MQ=255)
                                                             aCTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCgtg         <  1:1075688/70‑1 (MQ=255)
                                                               tGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                                       <  1:2060816/38‑1 (MQ=255)
                                                               tGTTAGAGGACATGTACATCAGGTTATCGATACCGTTc                                       <  1:1365873/38‑1 (MQ=255)
                                                                gTTAGAGGACATGTACATCAGGTTATCGATACCGTTCata                                    <  1:729273/40‑1 (MQ=255)
                                                                gTTAGAGGACATGTACATCAGGTTATCGATACCGTTCata                                    <  1:2252912/40‑1 (MQ=255)
                                                                 ttAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTg    <  1:2407123/71‑1 (MQ=255)
                                                                 ttAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTg    <  1:2326700/71‑1 (MQ=255)
                                                                  tAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTg    >  1:290465/1‑70 (MQ=255)
                                                                   agagGACATGTACATCAGGTTATCGATACCGTTCATATTCTGtt                             >  1:2662229/1‑44 (MQ=255)
                                                                   agagGACATGTACATCAGGTTATCGATACCGTTCATATTCTGtt                             >  1:381772/1‑44 (MQ=255)
                                                                   agagGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGc   <  1:2878758/70‑1 (MQ=255)
                                                                   agagGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCa  >  1:541248/1‑71 (MQ=255)
                                                                   |                                                                      
CGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCA  >  minE/312566‑312703

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: