Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,579 A→G 11.7% T98T (ACT→ACC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5790AG11.7% 272.5 / 23.0 102T98T (ACT→ACCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (50/40);  new base G (5/7);  total (55/47)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 109 total aligned reads displayed.

TCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCA  >  minE/312511‑312649
                                                                    |                                                                      
tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                                       >  1:1851380/1‑70 (MQ=255)
 cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                                       <  1:1706671/69‑1 (MQ=255)
 cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                                      >  1:2067582/1‑70 (MQ=255)
 cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                                      <  1:2299812/70‑1 (MQ=255)
 cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                                      <  1:1148549/70‑1 (MQ=255)
 cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                                      <  1:1343322/70‑1 (MQ=255)
 cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                                     >  1:10732/1‑71 (MQ=255)
  ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTTCACTTGCGCGATATCAGGATCGGt                                                                      <  1:94435/68‑1 (MQ=255)
  ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtcc                                                                    >  1:2054891/2‑69 (MQ=255)
  ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtc                                                                     <  1:1321715/69‑2 (MQ=16)
    gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                    >  1:1500614/1‑69 (MQ=255)
    gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                    >  1:2540908/1‑69 (MQ=255)
     cGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccgg                                                                 <  1:2303622/71‑3 (MQ=255)
     cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                                                  <  1:1935227/70‑1 (MQ=255)
     cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAg                                                                 <  1:3529200/71‑1 (MQ=255)
     cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAg                                                                 <  1:497714/71‑1 (MQ=255)
      ggCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                                                >  1:484397/1‑71 (MQ=16)
       gcagcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGAc                                                               <  1:352053/71‑1 (MQ=16)
       gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                               <  1:2401338/71‑1 (MQ=255)
       gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                               <  1:503449/71‑1 (MQ=255)
        cagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                               >  1:315622/1‑70 (MQ=255)
        cagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                                              <  1:2520448/71‑1 (MQ=255)
         agcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                                              >  1:3771861/1‑70 (MQ=255)
         agcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                                              >  1:2187233/1‑70 (MQ=255)
          gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                            <  1:1549774/71‑1 (MQ=255)
             aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                         >  1:684315/1‑71 (MQ=255)
             aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                         >  1:689156/1‑71 (MQ=255)
             aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                         >  1:962620/1‑71 (MQ=255)
              cGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAgg                                                        >  1:605973/1‑71 (MQ=255)
              cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                                        >  1:1005238/1‑71 (MQ=255)
              cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCACGATATCCGCATCAGTACCAGACTCAAAgg                                                        >  1:460897/1‑71 (MQ=255)
              cGGCATCGCCAGCTGCAGTTTGTTATGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                                        >  1:3402512/1‑71 (MQ=255)
               ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                       >  1:651667/1‑71 (MQ=255)
                gCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTc                                                      <  1:811430/71‑1 (MQ=255)
                 cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCa                                                     <  1:3712403/71‑1 (MQ=255)
                    cGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                                              >  1:598168/1‑59 (MQ=255)
                     gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaaa                                                          >  1:3284913/1‑62 (MQ=255)
                      cgaGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGta                                                <  1:2853999/69‑2 (MQ=12)
                      ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                                                <  1:1102180/71‑1 (MQ=255)
                       gaGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGtaa                                               <  1:2916913/70‑3 (MQ=255)
                       gaGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGtaa                                               <  1:639487/70‑3 (MQ=255)
                        aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                              <  1:2366295/71‑1 (MQ=255)
                          cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                                           <  1:3310785/56‑1 (MQ=255)
                            gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                         >  1:1511077/1‑56 (MQ=255)
                            gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAAACTCAAAGGTCAGGGTGATCTg                                           >  1:2604417/1‑70 (MQ=255)
                             gaGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAgg                                                        >  1:201974/2‑56 (MQ=12)
                             cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                                         >  1:1096764/1‑71 (MQ=255)
                              aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                                     >  1:2509207/1‑42 (MQ=255)
                              aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                                        >  1:1255720/1‑71 (MQ=255)
                              aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                                        >  1:1640108/1‑71 (MQ=255)
                                tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                       <  1:949551/54‑1 (MQ=255)
                                  tGTTCTGCACCTGAACCTGCGCGATATCCGCATCa                                                                        >  1:3421359/1‑35 (MQ=255)
                                  tGTTCTGCACCTGAACCTGCGCGATATCCGCATCa                                                                        >  1:1010098/1‑35 (MQ=255)
                                  tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                                    >  1:1017383/1‑71 (MQ=255)
                                   gTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                                    <  1:188869/70‑1 (MQ=255)
                                     tCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAGTCAAAGGTCAggg                                                  <  1:3450658/54‑1 (MQ=255)
                                      cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGc                                          >  1:1241334/1‑61 (MQ=255)
                                      cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc                                 >  1:3303864/1‑70 (MQ=255)
                                      cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc                                 >  1:927970/1‑70 (MQ=255)
                                      cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                                >  1:541457/1‑71 (MQ=255)
                                       tGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGt                               >  1:54975/1‑71 (MQ=255)
                                       tGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGt                               >  1:1424638/1‑71 (MQ=255)
                                        gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                                <  1:2697790/69‑1 (MQ=255)
                                          aCCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                               <  1:3729467/36‑1 (MQ=255)
                                          aCCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                              >  1:63853/1‑69 (MQ=255)
                                          aCCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                              >  1:1255845/1‑69 (MQ=255)
                                            cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                            <  1:1536096/37‑1 (MQ=255)
                                            cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                                    >  1:1609936/1‑61 (MQ=255)
                                            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)
                                                  cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                    <  1:2635633/71‑1 (MQ=255)
                                                    gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa                        >  1:2965798/1‑65 (MQ=255)
                                                     cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGa                            <  1:2305820/60‑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)
                                                          tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                              <  1:2708622/37‑1 (MQ=255)
                                                           atCCGCATCAGTACGAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACa           >  1:33773/1‑71 (MQ=255)
                                                           atCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                                        >  1:2586712/1‑42 (MQ=255)
                                                              cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg         <  1:3094677/70‑1 (MQ=255)
                                                              cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg         >  1:2575464/1‑70 (MQ=255)
                                                               gCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa       >  1:3178849/1‑71 (MQ=255)
                                                               gCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa       >  1:1490258/1‑71 (MQ=255)
                                                                cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa       >  1:1144665/1‑70 (MQ=255)
                                                                cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc      <  1:2874417/71‑1 (MQ=255)
                                                                cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTAc      <  1:2717478/71‑1 (MQ=255)
                                                                 aTCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTaga                <  1:919007/60‑1 (MQ=255)
                                                                 aTCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTaga                <  1:203686/60‑1 (MQ=255)
                                                                    aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                              >  1:1721736/1‑43 (MQ=255)
                                                                    aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTagag               >  1:1795913/1‑58 (MQ=255)
                                                                    aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa  <  1:2134927/71‑1 (MQ=255)
                                                                    aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa  <  1:2076713/71‑1 (MQ=255)
                                                                    aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa  <  1:3408802/71‑1 (MQ=255)
                                                                    |                                                                      
TCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCA  >  minE/312511‑312649

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: