Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 309,836 T→A 18.0% T1013S (ACC→TCC) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE309,8360TA18.0% 163.7 / 28.6 72T1013S (ACC→TCC) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base A (8/5);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

GCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAAC  >  minE/309768‑309905
                                                                    |                                                                     
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:3590208/71‑1 (MQ=255)
               gcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                                      <  1:3528139/71‑1 (MQ=255)
               gcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                                      <  1:2805596/71‑1 (MQ=255)
               gcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                                      <  1:1784113/71‑1 (MQ=255)
                cgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa                                                      >  1:2941047/1‑70 (MQ=255)
                 gAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTaca                                                    <  1:324198/71‑1 (MQ=255)
                   accaccACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCa                                                                  >  1:2423986/1‑55 (MQ=255)
                   accaccACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCa                                                                  >  1:3942667/1‑55 (MQ=255)
                   accaccACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCa                                                                  >  1:1844696/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)
                           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:1502455/1‑71 (MQ=255)
                                        aCGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgcgc                             <  1:3761026/71‑1 (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)
                                                          aCCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCa             >  1:467581/1‑69 (MQ=255)
                                                          aCCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAg            <  1:447574/70‑1 (MQ=255)
                                                            cGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTc                                  <  1:2388672/46‑1 (MQ=255)
                                                            cGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTAc         >  1:1971779/1‑71 (MQ=255)
                                                            cGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTAc         <  1:3538927/71‑1 (MQ=255)
                                                            cGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTAc         >  1:291112/1‑71 (MQ=255)
                                                              ttGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCAc               <  1:1483380/63‑1 (MQ=255)
                                                                gCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGa      >  1:2202386/1‑70 (MQ=255)
                                                                 cGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGAt     <  1:1876620/70‑1 (MQ=255)
                                                                  ggTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTg                                <  1:1274505/42‑1 (MQ=255)
                                                                  ggTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTg                                <  1:2328175/42‑1 (MQ=255)
                                                                    tCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCt                                 <  1:1327094/39‑1 (MQ=255)
                                                                    tCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCt                                 <  1:3047510/39‑1 (MQ=255)
                                                                    tCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAAc  >  1:328478/1‑70 (MQ=255)
                                                                    tCACCATCCCCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCAcc              <  1:155341/57‑1 (MQ=255)
                                                                    |                                                                     
GCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAAC  >  minE/309768‑309905

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: