Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 309,825 T→C 19.7% V1016V (GTA→GTG) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE309,8250TC19.7% 147.1 / 30.9 71V1016V (GTA→GTG) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (32/25);  new base C (9/5);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

ATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCAC  >  minE/309756‑309892
                                                                     |                                                                   
atCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTa                                                                    >  1:204623/1‑71 (MQ=255)
    tCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGt                                                                >  1:2942161/1‑71 (MQ=255)
    tCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGt                                                                >  1:2778805/1‑71 (MQ=255)
    tCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGt                                                                >  1:725538/1‑71 (MQ=255)
     cATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCg                                                                 <  1:3114024/69‑1 (MQ=255)
     cATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGtt                                                               >  1:859975/1‑71 (MQ=255)
      aTTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGtt                                                               <  1:3180454/70‑1 (MQ=255)
      aTTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTg                                                              >  1:3272525/1‑71 (MQ=255)
      aTTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTg                                                              >  1:2099800/1‑71 (MQ=255)
         cTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCgg                                                           >  1:3006742/1‑71 (MQ=255)
            gCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCa                                                        >  1:1647239/1‑71 (MQ=255)
                 tAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCa                                                        >  1:2418527/1‑66 (MQ=255)
                 tAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt                                                    <  1:1489815/70‑1 (MQ=255)
                   aaGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATcc                                                  <  1:749825/70‑1 (MQ=255)
                    aGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCg                                                >  1:1225334/1‑71 (MQ=255)
                    aGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCg                                                >  1:3948018/1‑71 (MQ=255)
                      cggcggcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGcc                                              >  1:1194655/1‑71 (MQ=255)
                      caacggcggATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCATtcctcc                                              >  1:2110854/4‑71 (MQ=255)
                         cggcggATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCATtcct                                                >  1:1561204/1‑65 (MQ=255)
                          ggcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCAtt                                          <  1:1048865/71‑1 (MQ=255)
                          ggcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCAt                                           >  1:1330923/1‑70 (MQ=255)
                           gcgAACCACCACAAAGACTACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                         <  1:2209098/71‑1 (MQ=255)
                           gcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                         <  1:3528139/71‑1 (MQ=255)
                           gcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                         <  1:1784113/71‑1 (MQ=255)
                           gcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                         <  1:3590208/71‑1 (MQ=255)
                           gcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                         <  1:2805596/71‑1 (MQ=255)
                            cgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCg                                                                 <  1:457684/46‑1 (MQ=255)
                            cgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                         >  1:2941047/1‑70 (MQ=255)
                             gAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTaca                                       <  1:324198/71‑1 (MQ=255)
                               accaccACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCa                                                     >  1:1844696/1‑55 (MQ=255)
                               accaccACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCa                                                     >  1:2423986/1‑55 (MQ=255)
                               accaccACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCa                                                     >  1:3942667/1‑55 (MQ=255)
                                ccaccaCAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATccc                                                 >  1:2098850/1‑58 (MQ=255)
                                ccaccaCAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCg                                    >  1:834336/1‑71 (MQ=255)
                                   ccacAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                         <  1:2752933/63‑1 (MQ=255)
                                   ccacAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                         <  1:460194/63‑1 (MQ=255)
                                    cacaAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTAc                                >  1:968756/1‑71 (MQ=255)
                                     acaAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAttc                                                  >  1:366737/1‑50 (MQ=16)
                                     acaAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCATTCCTCCCATTACCCCGATAcc                               >  1:3566960/1‑71 (MQ=255)
                                      caAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGccc                                             <  1:3565250/56‑1 (MQ=255)
                                       aaaGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTa                             >  1:1611431/1‑71 (MQ=255)
                                       aaaGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGcag                                                           >  1:3773625/1‑39 (MQ=19)
                                       aaaGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAtt                                                   >  1:1296639/1‑48 (MQ=16)
                                         aGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCATTCCTCCCATTACCCCGATAccca                             >  1:185324/1‑67 (MQ=255)
                                          gAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTg                          >  1:1474655/1‑71 (MQ=255)
                                            acaCCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCATTCCTCCCATTACCCCGATACCCACAGCg                        <  1:3082222/69‑1 (MQ=255)
                                            aTACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCg                        <  1:350623/71‑1 (MQ=255)
                                              aCCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCATTCCTCCCATTACCCCGATACCCACAGCGtt                      <  1:1227198/71‑1 (MQ=255)
                                               ccGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCATTCCTCCCATTACCCCGATACCCACAGCGTTc                     <  1:5099/71‑1 (MQ=255)
                                               ccGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGtt                      <  1:912952/70‑1 (MQ=255)
                                                cGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCATTCCTCCCATTACCCCGATACCCACAGCGTTCt                    <  1:756989/71‑1 (MQ=255)
                                                  gAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgc                  >  1:3107272/1‑71 (MQ=255)
                                                   taCGAAGAAGATTGCCAGCAACGTTGCAGAGACCATTCCTCCCATTACCCCGATACCCACAGCGTTCTgcg                 >  1:350509/2‑71 (MQ=255)
                                                    aCGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgcgc                <  1:3699715/71‑1 (MQ=255)
                                                    aCGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgcgc                <  1:3761026/71‑1 (MQ=255)
                                                    aCGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgcgc                >  1:1502455/1‑71 (MQ=255)
                                                    aCGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgcg                 >  1:2198179/1‑70 (MQ=255)
                                                    aCGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGGTCTgcgc                <  1:1468783/71‑1 (MQ=255)
                                                     cGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgcgcg               <  1:2391592/71‑1 (MQ=255)
                                                      gacgaagaTTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgcgcgc              <  1:3693684/68‑1 (MQ=255)
                                                      gaagaaGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgcgcgc              <  1:844117/71‑1 (MQ=255)
                                                      gaagaaGATTGCCAGCAACGTTGCAGAGACCATTCCTCCCATTACCCCGATACCCACAGCGTTCTgcgcgc              >  1:2518292/1‑71 (MQ=255)
                                                         gaagaTTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCgg           <  1:1516309/71‑1 (MQ=255)
                                                         gaagaTTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCg            >  1:432992/1‑70 (MQ=255)
                                                          aagaTTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCgg           >  1:343039/1‑70 (MQ=255)
                                                          aagaTTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCgg           >  1:1533154/1‑70 (MQ=255)
                                                            gaTTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAAc        <  1:1058779/71‑1 (MQ=255)
                                                              ttGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCa      <  1:1277295/71‑1 (MQ=255)
                                                               tGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAg     >  1:2240384/1‑71 (MQ=255)
                                                                    gTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCAc  >  1:3017790/1‑69 (MQ=255)
                                                                    gTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCAc  >  1:2411339/1‑69 (MQ=255)
                                                                     tACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGcc             >  1:3569415/1‑57 (MQ=255)
                                                                     tACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGcc             >  1:2374236/1‑57 (MQ=255)
                                                                     |                                                                   
ATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCAC  >  minE/309756‑309892

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: