Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 309,859 G→A 13.0% T1005I (ACC→ATC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE309,8590GA13.0% 200.9 / 17.8 77T1005I (ACC→ATC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (37/30);  new base A (6/4);  total (43/34)
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.66e-01

CACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAA  >  minE/309792‑309928
                                                                   |                                                                     
cacaAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTAc                                                                    >  1:968756/1‑71 (MQ=255)
 acaAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCATTCCTCCCATTACCCCGATAcc                                                                   >  1:3566960/1‑71 (MQ=255)
   aaaGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTa                                                                 >  1:1611431/1‑71 (MQ=255)
     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:1502455/1‑71 (MQ=255)
                aCGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgcgc                                                    <  1:3761026/71‑1 (MQ=255)
                aCGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTgcgc                                                    <  1:3699715/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:1533154/1‑70 (MQ=255)
                      aagaTTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCgg                                               >  1:343039/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:291112/1‑71 (MQ=255)
                                    cGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTAc                                >  1:1971779/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:3047510/39‑1 (MQ=255)
                                            tCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCt                                                        <  1:1327094/39‑1 (MQ=255)
                                            tCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAAc                         >  1:328478/1‑70 (MQ=255)
                                            tCACCATCCCCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCAcc                                     <  1:155341/57‑1 (MQ=255)
                                             cACCATCCCGCCCATTACACCGGTACCTACTGCGTTc                                                         >  1:1549396/1‑37 (MQ=255)
                                             cACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAg                                   >  1:2778931/1‑59 (MQ=255)
                                             cACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAAc                         >  1:2954398/1‑69 (MQ=255)
                                             cACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAAc                         >  1:3761830/1‑69 (MQ=255)
                                              aCCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACt                               <  1:3247268/62‑1 (MQ=255)
                                              aCCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCa                       <  1:3657933/70‑1 (MQ=255)
                                              aCCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCa                       >  1:927027/1‑70 (MQ=255)
                                              aCCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAg                      >  1:1082288/1‑71 (MQ=255)
                                              aCCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAg                      >  1:18761/1‑71 (MQ=255)
                                                cATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCTGAACCAGCACCAGTACTGATAACCAGCg                    <  1:3407678/71‑1 (MQ=255)
                                                cATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCg                    >  1:1425086/1‑71 (MQ=255)
                                                cATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCg                    >  1:3805356/1‑71 (MQ=255)
                                                  tCCCGCCCATTACACCGGTACCTACTGCGTTCTgcg                                                     >  1:2616681/1‑36 (MQ=255)
                                                  tCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCgg                   >  1:2289047/1‑70 (MQ=255)
                                                   cccGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCa                 >  1:1801455/1‑71 (MQ=255)
                                                    ccGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCg                    <  1:409043/67‑1 (MQ=255)
                                                       cccATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATaa              >  1:1930723/1‑70 (MQ=255)
                                                       cccATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAAc             >  1:2299343/1‑71 (MQ=255)
                                                       cccATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAAc             >  1:3348179/1‑71 (MQ=255)
                                                        ccATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGAt                            <  1:1395343/55‑1 (MQ=255)
                                                          aTTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGcc          >  1:1702184/1‑71 (MQ=255)
                                                             acacCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCAc                                      >  1:3937843/1‑40 (MQ=255)
                                                             acacCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTa                                 >  1:1573387/1‑45 (MQ=255)
                                                             acacCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGa        <  1:607219/70‑1 (MQ=255)
                                                               acCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGa     <  1:3418941/71‑1 (MQ=255)
                                                               acCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGa     <  1:771483/71‑1 (MQ=255)
                                                                ccGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGAt    <  1:2760635/71‑1 (MQ=255)
                                                                ccGATACCCACAGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAg       >  1:2664856/1‑68 (MQ=255)
                                                                ccGATACCCACAGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAg       >  1:1643547/1‑68 (MQ=255)
                                                                   gTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATaa  <  1:83574/70‑1 (MQ=255)
                                                                   |                                                                     
CACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAA  >  minE/309792‑309928

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: