Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,582 A→C 7.3% G97G (GGT→GGG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5820AC 7.3% 250.5 / 13.2 96G97G (GGT→GGGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (47/42);  new base C (3/4);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

CTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGG  >  minE/312512‑312651
                                                                      |                                                                     
cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                                       >  1:10732/1‑71 (MQ=255)
   gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                      >  1:2540908/1‑69 (MQ=255)
   gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                                      >  1:1500614/1‑69 (MQ=255)
    cGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccgg                                                                   <  1:2303622/71‑3 (MQ=255)
    cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                                                    <  1:1935227/70‑1 (MQ=255)
    cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAg                                                                   <  1:497714/71‑1 (MQ=255)
    cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAg                                                                   <  1:3529200/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:2187233/1‑70 (MQ=255)
        agcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                                                >  1:3771861/1‑70 (MQ=255)
         gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTc                                                               >  1:1205334/1‑70 (MQ=255)
         gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                              <  1:1549774/71‑1 (MQ=255)
         gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                              <  1:3253030/71‑1 (MQ=255)
            aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                           >  1:962620/1‑71 (MQ=255)
            aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                           >  1:689156/1‑71 (MQ=255)
            aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                           >  1:684315/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:639487/70‑3 (MQ=255)
                      gaGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGtaa                                                 <  1:2916913/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)
                           gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                             >  1:633201/1‑70 (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:1640108/1‑71 (MQ=255)
                             aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                                          >  1:1255720/1‑71 (MQ=255)
                               tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                         <  1:998822/54‑1 (MQ=255)
                               tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                         <  1:949551/54‑1 (MQ=255)
                               tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                         <  1:2747186/54‑1 (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:1424638/1‑71 (MQ=255)
                                      tGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGt                                 >  1:54975/1‑71 (MQ=255)
                                       gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                                  <  1:2500555/69‑1 (MQ=255)
                                       gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                                  <  1:2697790/69‑1 (MQ=255)
                                         aCCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                                 <  1:3729467/36‑1 (MQ=255)
                                         aCCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                                >  1:1255845/1‑69 (MQ=255)
                                         aCCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                                >  1:63853/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:2635633/71‑1 (MQ=255)
                                                 cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt                      <  1:3169355/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:10652/1‑71 (MQ=255)
                                                        atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa               >  1:2873614/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:2575464/1‑70 (MQ=255)
                                                             cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg           <  1:3094677/70‑1 (MQ=255)
                                                              gCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa         >  1:1490258/1‑71 (MQ=255)
                                                              gCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa         >  1:3178849/1‑71 (MQ=255)
                                                               cATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTa         <  1:3566946/70‑1 (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:203686/60‑1 (MQ=255)
                                                                aTCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTaga                  <  1:919007/60‑1 (MQ=255)
                                                                   aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                                >  1:1721736/1‑43 (MQ=255)
                                                                   aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTagag                 >  1:1795913/1‑58 (MQ=255)
                                                                   aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa    <  1:2076713/71‑1 (MQ=255)
                                                                   aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa    <  1:3408802/71‑1 (MQ=255)
                                                                   aGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCa    <  1:2134927/71‑1 (MQ=255)
                                                                    gTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGt          <  1:1073532/64‑1 (MQ=255)
                                                                     tACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAgg  <  1:640487/71‑1 (MQ=255)
                                                                     tACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAgg  <  1:1200755/71‑1 (MQ=255)
                                                                      |                                                                     
CTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATGTACATCAGG  >  minE/312512‑312651

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: