Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,664 G→C 17.6% G403G (GGC→GGG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,6640GC17.6% 218.0 / 32.4 85G403G (GGC→GGGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (36/34);  new base C (9/6);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAA  >  minE/311595‑311734
                                                                     |                                                                      
cggcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAgg                                                                        >  1:96511/1‑70 (MQ=255)
cggcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGc                                                                       >  1:2992948/1‑71 (MQ=255)
    aaaCCTTCTTCCGCCATAACACGCTAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                                  >  1:3383171/1‑71 (MQ=255)
     aaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                                  <  1:2664213/71‑1 (MQ=255)
     aaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                                  <  1:290721/71‑1 (MQ=255)
     aaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                                  >  1:1955805/1‑71 (MQ=255)
      aCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATgg                                                                 >  1:1359380/1‑71 (MQ=255)
        cttcttCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCg                                                               <  1:123282/71‑1 (MQ=255)
         ttcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGc                                                                >  1:2761354/1‑69 (MQ=255)
         ttcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGa                                                              >  1:3515405/1‑71 (MQ=255)
          tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGc                                                                       >  1:1427650/1‑61 (MQ=255)
          tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                                             >  1:488728/1‑71 (MQ=255)
          tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                                             >  1:3675685/1‑71 (MQ=255)
              ccGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCa                                                           <  1:1651762/69‑1 (MQ=255)
               cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCAcc                                                         <  1:3713951/70‑1 (MQ=255)
               cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCAcc                                                         <  1:274408/70‑1 (MQ=255)
               cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCAcc                                                         <  1:903897/70‑1 (MQ=255)
                  cATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATcc                                                     >  1:333152/1‑71 (MQ=255)
                      acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                                 >  1:2133125/1‑71 (MQ=255)
                      acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGa                                                  >  1:3355033/1‑70 (MQ=255)
                       cacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATgg                                                                 <  1:1873226/54‑1 (MQ=255)
                       cacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                                 >  1:3239778/1‑70 (MQ=255)
                        acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATc                                                      <  1:3904518/64‑1 (MQ=255)
                        acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATc                                                      <  1:3175583/64‑1 (MQ=255)
                        acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                               <  1:1348260/71‑1 (MQ=255)
                         cGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAt                                              <  1:1507946/71‑1 (MQ=255)
                         cGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAt                                              >  1:3003479/1‑71 (MQ=255)
                           cTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAtt                                             <  1:1341490/70‑1 (MQ=255)
                            tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGc                                                                       >  1:2432467/1‑43 (MQ=255)
                              aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGa                                                              >  1:1735466/1‑50 (MQ=255)
                              aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGtt                                          <  1:2023695/70‑1 (MQ=255)
                              aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTa                                         <  1:692119/71‑1 (MQ=255)
                                cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                       <  1:848997/71‑1 (MQ=255)
                                cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                       >  1:1053744/1‑71 (MQ=255)
                                 gTTTTCTACCACAACGATGGCGTCATCCACCAACAGGc                                                                       >  1:3020226/1‑38 (MQ=255)
                                 gTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCg                                      >  1:992196/1‑71 (MQ=255)
                                 gTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCg                                      <  1:2724515/71‑1 (MQ=255)
                                  ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAg                                        <  1:1369796/68‑1 (MQ=255)
                                    ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                                   >  1:2238548/1‑71 (MQ=255)
                                    ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                                   >  1:282127/1‑71 (MQ=255)
                                     tCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                  >  1:3860039/1‑71 (MQ=255)
                                       tACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                >  1:3897047/1‑71 (MQ=255)
                                       tACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                >  1:1699142/1‑71 (MQ=255)
                                       caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                >  1:93297/2‑69 (MQ=9)
                                       caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                >  1:482721/2‑69 (MQ=9)
                                       caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                >  1:3350701/2‑69 (MQ=9)
                                       caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttg                                 <  1:1540903/69‑2 (MQ=9)
                                       caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttg                                 <  1:193640/69‑2 (MQ=9)
                                        accaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAt                               >  1:2544236/1‑71 (MQ=9)
                                        accaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAt                               <  1:1366210/71‑1 (MQ=9)
                                        aCCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                               >  1:3452488/1‑71 (MQ=255)
                                         ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                               <  1:1304492/70‑1 (MQ=255)
                                         ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATc                              <  1:1727597/71‑1 (MQ=255)
                                          cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                               <  1:3109005/53‑1 (MQ=255)
                                          cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                <  1:1515133/68‑1 (MQ=255)
                                          cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                             <  1:564412/71‑1 (MQ=255)
                                          cacaACGATCGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                             <  1:1458078/71‑1 (MQ=255)
                                            caACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag                           <  1:2228932/71‑1 (MQ=255)
                                             aaCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATccc                                                    >  1:1344730/1‑45 (MQ=255)
                                             aaCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga                          <  1:2189356/71‑1 (MQ=255)
                                                  tGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                      <  1:2971964/70‑1 (MQ=255)
                                                     cGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAg                  <  1:3797414/71‑1 (MQ=255)
                                                      atCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaaaag                 >  1:3276913/2‑69 (MQ=9)
                                                      atCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaaaag                 >  1:2778838/2‑69 (MQ=9)
                                                       tCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGtt                                  <  1:1540298/53‑1 (MQ=16)
                                                       tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg                             <  1:1899985/58‑1 (MQ=255)
                                                        cATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAg                        >  1:227123/1‑62 (MQ=255)
                                                         aTCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                   >  1:563275/1‑66 (MQ=255)
                                                          tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc             <  1:1824706/71‑1 (MQ=255)
                                                          tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc             <  1:2575103/71‑1 (MQ=255)
                                                              ccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                  >  1:3401248/1‑46 (MQ=255)
                                                               caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa        <  1:2470315/71‑1 (MQ=255)
                                                                agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagcggc             <  1:1648596/63‑1 (MQ=12)
                                                                agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat       >  1:2319320/3‑70 (MQ=12)
                                                                agcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat       <  1:2299073/69‑2 (MQ=12)
                                                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                       >  1:1262751/1‑39 (MQ=255)
                                                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                       >  1:426267/1‑39 (MQ=255)
                                                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt                               <  1:3787862/47‑1 (MQ=255)
                                                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGaa                         >  1:2471585/1‑53 (MQ=255)
                                                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc       <  1:2742190/71‑1 (MQ=255)
                                                                 gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGatg      >  1:3815727/2‑69 (MQ=12)
                                                                 gcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGatg      >  1:576654/2‑69 (MQ=12)
                                                                 aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                  >  1:2683269/1‑43 (MQ=255)
                                                                 aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                  >  1:1569027/1‑43 (MQ=255)
                                                                 aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg      >  1:3497243/1‑71 (MQ=255)
                                                                   aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACgg     >  1:3534739/1‑70 (MQ=255)
                                                                    ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGc                       <  1:3191325/51‑1 (MQ=255)
                                                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa  <  1:2124814/71‑1 (MQ=255)
                                                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa  <  1:1796073/71‑1 (MQ=255)
                                                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa  <  1:3864917/71‑1 (MQ=255)
                                                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa  <  1:2363271/71‑1 (MQ=255)
                                                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa   <  1:3006533/70‑1 (MQ=255)
                                                                     |                                                                      
CGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAA  >  minE/311595‑311734

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: