Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 310,815 G→A 5.8% D686D (GAC→GAT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE310,8150GA 5.8% 226.7 / 2.9 69D686D (GAC→GATacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base A (2/2);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAA  >  minE/310746‑310883
                                                                     |                                                                    
gTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGt                                                                     >  1:453074/1‑71 (MQ=255)
   cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                                                                  >  1:1298517/1‑71 (MQ=255)
     tCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                                                 >  1:1500400/1‑70 (MQ=255)
     tCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTCACCAAGGCCAGCCTGGTCAAt                                                                 >  1:1631682/1‑70 (MQ=255)
      cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                                                 >  1:1303049/1‑69 (MQ=255)
       gctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCa                                                               <  1:1218891/70‑1 (MQ=255)
        ctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCa                                                               >  1:1615559/1‑69 (MQ=255)
          gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                                            >  1:2075878/1‑70 (MQ=255)
           ctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTc                                                                    <  1:377604/61‑1 (MQ=255)
           ctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                                                          <  1:498104/71‑1 (MQ=255)
           ctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                                                          <  1:3699435/71‑1 (MQ=255)
           ctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtc                                                           >  1:4019568/1‑70 (MQ=255)
             tCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                                                          >  1:2032543/1‑69 (MQ=255)
                gcaagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagt                                                      >  1:3668354/1‑70 (MQ=255)
                gcaagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCAAGCCTGGTCAATCAGCtcaaagt                                                      >  1:2797037/1‑70 (MQ=255)
                ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaa                                                        >  1:476855/2‑68 (MQ=255)
                ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaa                                                        >  1:1111326/2‑68 (MQ=255)
                 caagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagt                                                      <  1:3593032/69‑1 (MQ=255)
                  aagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga                                                   <  1:3619062/71‑3 (MQ=255)
                  aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca                                                    >  1:2451158/1‑70 (MQ=255)
                    gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCTATTAACtcaaa                                                        <  1:2812089/64‑1 (MQ=255)
                    gcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag                                                 >  1:3793891/1‑71 (MQ=255)
                     caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc                                                >  1:28976/1‑71 (MQ=255)
                      aaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                               >  1:976190/1‑71 (MQ=255)
                        cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                                                                   <  1:1189699/49‑1 (MQ=255)
                        cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                               <  1:3286078/69‑1 (MQ=255)
                         tGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGt                                            >  1:4017314/1‑71 (MQ=255)
                          ggTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATc                                                                >  1:823528/1‑50 (MQ=255)
                           gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGt                                            <  1:2653576/69‑1 (MQ=255)
                           gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                          <  1:1415522/71‑1 (MQ=255)
                            ttACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCATCCTGGTCAATCAGCTCAAAGTCAAAGCGGGTTg                                          <  1:1118833/70‑1 (MQ=255)
                            ttACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                          >  1:817641/1‑70 (MQ=255)
                               cgcgCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                          <  1:263900/67‑1 (MQ=255)
                               cgcgCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAg                                       <  1:2273834/70‑1 (MQ=255)
                               cgcgCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAg                                       <  1:2561521/70‑1 (MQ=255)
                                       aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGtt                              >  1:985914/1‑71 (MQ=255)
                                        gTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTc                             >  1:1171458/1‑71 (MQ=255)
                                         tgaGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTc                             >  1:2442801/3‑70 (MQ=255)
                                         tCAGTTTTTCGTGACCAAGGCCGGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                               <  1:2581810/52‑1 (MQ=255)
                                         tCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTcc                            <  1:2049929/71‑1 (MQ=255)
                                           aGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCa                           <  1:3308177/70‑1 (MQ=255)
                                           aGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCAc                          >  1:3860641/1‑71 (MQ=255)
                                             tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                                                         <  1:930832/38‑1 (MQ=255)
                                             tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCg                                              <  1:2477937/49‑1 (MQ=255)
                                             tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCa                           >  1:24111/1‑68 (MQ=255)
                                             tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCAc                          >  1:2252873/1‑69 (MQ=255)
                                              ttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGAt                       >  1:2747150/1‑71 (MQ=255)
                                               tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg                                             >  1:2260991/1‑48 (MQ=255)
                                               tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAg                                       >  1:424045/1‑54 (MQ=255)
                                               tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGAtt                      >  1:3760848/1‑71 (MQ=255)
                                               tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGAtt                      >  1:1089315/1‑71 (MQ=255)
                                               tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGAt                       <  1:1618945/70‑1 (MQ=255)
                                                    tGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc                                                <  1:3858891/40‑1 (MQ=255)
                                                    tGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTc                             <  1:3734477/59‑1 (MQ=255)
                                                      aCCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCggg                 >  1:657521/1‑69 (MQ=255)
                                                      aCCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGc                <  1:2443875/70‑1 (MQ=255)
                                                        cAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCg                                              >  1:1704989/1‑38 (MQ=255)
                                                        cAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTc                             <  1:4000403/55‑1 (MQ=255)
                                                         aaGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGc                                         <  1:3362337/42‑1 (MQ=255)
                                                         aaGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGc                                         <  1:129647/42‑1 (MQ=255)
                                                          aGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGt            >  1:2075322/1‑70 (MQ=255)
                                                           ggCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGtt           >  1:587665/1‑70 (MQ=255)
                                                           ggCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTa          <  1:2447681/71‑1 (MQ=255)
                                                           ggCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTa          <  1:884100/71‑1 (MQ=255)
                                                            gCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg                                             >  1:2974904/1‑35 (MQ=255)
                                                            gCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTcc                            <  1:1814448/52‑1 (MQ=255)
                                                                 ccTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGc     >  1:703462/1‑70 (MQ=255)
                                                                 ccTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCg    <  1:2063356/71‑1 (MQ=255)
                                                                   tGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGa   >  1:3218036/1‑70 (MQ=255)
                                                                    ggTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGaa  >  1:2977729/1‑70 (MQ=255)
                                                                     gTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCa               <  1:3636614/56‑1 (MQ=255)
                                                                     |                                                                    
GTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAA  >  minE/310746‑310883

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