Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 310,821 C→T 5.9% L684L (CTG→CTA) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE310,8210CT 5.9% 225.5 / 2.3 68L684L (CTG→CTA) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (35/29);  new base T (2/2);  total (37/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.98e-01

GCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCAT  >  minE/310753‑310889
                                                                    |                                                                    
gctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCa                                                                     <  1:1218891/70‑1 (MQ=255)
 ctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCa                                                                     >  1:1615559/1‑69 (MQ=255)
   gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                                                  >  1:2075878/1‑70 (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:1111326/2‑68 (MQ=255)
         ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaa                                                              >  1:476855/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)
                 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:1089315/1‑71 (MQ=255)
                                        tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGAtt                            >  1:3760848/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:884100/71‑1 (MQ=255)
                                                    ggCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTa                <  1:2447681/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)
                                                                 aaTCAGCTCAAAGTCAAAGCCGGTTGCAGTGCCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAAc     >  1:3924140/1‑69 (MQ=255)
                                                                 aaTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCa   >  1:3296466/1‑71 (MQ=255)
                                                                  aTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGaa        <  1:934726/65‑1 (MQ=255)
                                                                  aTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCAt  >  1:845469/1‑71 (MQ=255)
                                                                   tCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGa                              >  1:1229398/1‑42 (MQ=255)
                                                                    |                                                                    
GCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCAT  >  minE/310753‑310889

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: