Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 309,834 G→A 18.6% T1013T (ACC→ACT) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE309,8340GA18.6% 157.6 / 26.9 70T1013T (ACC→ACT) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base A (8/5);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

CTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGAT  >  minE/309765‑309902
                                                                     |                                                                    
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:3590208/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:3942667/1‑55 (MQ=255)
                      accaccACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCa                                                               >  1:1844696/1‑55 (MQ=255)
                      accaccACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCa                                                               >  1:2423986/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:3538927/71‑1 (MQ=255)
                                                               cGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTAc      >  1:1971779/1‑71 (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)
                                                                     |                                                                    
CTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGAT  >  minE/309765‑309902

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: