Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,496 A→G 8.8% F459F (TTT→TTC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4960AG 8.8% 257.5 / 14.0 91F459F (TTT→TTCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base G (3/5);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 101 total aligned reads displayed.

GTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCC  >  minE/311426‑311566
                                                                      |                                                                      
gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                                                                        <  1:798959/71‑1 (MQ=255)
   ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                                                                        >  1:4080458/1‑68 (MQ=255)
   ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag                                                                     <  1:2055908/71‑1 (MQ=255)
   ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag                                                                     <  1:4008394/71‑1 (MQ=255)
      aCAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaag                                                                  >  1:1475278/1‑71 (MQ=255)
           gcCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGaag                                                                     <  1:3311858/63‑1 (MQ=12)
              aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                                                                      <  1:217482/59‑1 (MQ=255)
               ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                                          >  1:1240749/1‑70 (MQ=255)
               ttGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGaag                                                                     <  1:2683586/59‑1 (MQ=9)
                tGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                                                                      <  1:4143381/57‑1 (MQ=255)
                 gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                                       >  1:884535/1‑71 (MQ=255)
                 gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                                       >  1:3313949/1‑71 (MQ=255)
                 gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                                       >  1:3593586/1‑71 (MQ=255)
                  cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                                       >  1:1075745/1‑70 (MQ=255)
                  cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                                       >  1:1052156/1‑70 (MQ=255)
                  cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:2634234/71‑1 (MQ=255)
                  cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:3604915/71‑1 (MQ=255)
                  cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:445949/71‑1 (MQ=255)
                  cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:2696760/71‑1 (MQ=255)
                    gAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGaaga                                                                    <  1:612709/55‑1 (MQ=12)
                     aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                                                    >  1:1214635/1‑70 (MQ=255)
                     aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                                                    <  1:3294276/70‑1 (MQ=255)
                     aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                                                    >  1:2967883/1‑70 (MQ=255)
                     aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                                                   <  1:456851/71‑1 (MQ=255)
                     aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                                                   <  1:726494/71‑1 (MQ=255)
                      aaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                                                                        <  1:3808073/49‑1 (MQ=255)
                       aCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:928449/66‑1 (MQ=255)
                        cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAcc                                                <  1:603014/71‑1 (MQ=255)
                        cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAcc                                                <  1:2620980/71‑1 (MQ=255)
                           tGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAccgcc                                             <  1:1802132/71‑1 (MQ=255)
                             gTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGa                                           <  1:73839/71‑1 (MQ=255)
                                aTAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                                                                      <  1:1344794/41‑1 (MQ=255)
                                  agagAACTGACGATAGATAGCACCAGTAGGACCGCCAAAGAAGGCCAt                                                             >  1:2305752/1‑48 (MQ=255)
                                   gagaACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCAc                                     >  1:2879698/1‑71 (MQ=9)
                                   gagaACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                                                                   >  1:2737850/1‑41 (MQ=255)
                                   gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                                        >  1:2858637/1‑52 (MQ=255)
                                   gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                                                   <  1:3472467/57‑1 (MQ=255)
                                   gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                                                   <  1:54172/57‑1 (MQ=255)
                                   gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                                                   <  1:69656/57‑1 (MQ=255)
                                   gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAAt                                                   <  1:2858723/57‑1 (MQ=255)
                                    agaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:497705/53‑1 (MQ=255)
                                    agaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAcc                                    <  1:2908535/71‑1 (MQ=255)
                                     gaACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGa                                           >  1:3653763/1‑63 (MQ=9)
                                      aaCTGACGATAGATAGCACCAGTAGAACCGCCaaa                                                                      >  1:1190044/1‑35 (MQ=255)
                                      aaCTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAccgcc                                             <  1:1774125/60‑1 (MQ=255)
                                        cTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATcg                                >  1:3074943/1‑71 (MQ=255)
                                        cTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                                                    >  1:2611265/1‑51 (MQ=255)
                                          gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                               >  1:2217854/1‑70 (MQ=255)
                                          gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                               >  1:2164353/1‑70 (MQ=255)
                                           gcGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGa                             <  1:57355/70‑1 (MQ=255)
                                           aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                                                    >  1:1042407/1‑48 (MQ=255)
                                           aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                               <  1:135907/69‑1 (MQ=255)
                                           aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                               <  1:3675328/69‑1 (MQ=255)
                                           aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                               <  1:2122659/69‑1 (MQ=255)
                                           aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                               <  1:2649993/69‑1 (MQ=255)
                                            cGATAAATTGCCCCAGTAGAACCGCCGGAGAATGCCATCGGAATAAATAcc                                                <  1:2682735/51‑1 (MQ=255)
                                              atagatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                                        <  1:2630602/57‑1 (MQ=255)
                                              atagatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                                        <  1:566775/57‑1 (MQ=255)
                                               tagatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:510594/42‑1 (MQ=255)
                                                agatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATAcc                         >  1:203348/1‑70 (MQ=255)
                                                 gatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:2847626/40‑1 (MQ=255)
                                                  atagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                               <  1:440967/62‑1 (MQ=255)
                                                  atagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                               <  1:136807/62‑1 (MQ=255)
                                                    agCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGAc                      >  1:1004187/1‑69 (MQ=255)
                                                     gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa                             >  1:2320443/1‑61 (MQ=255)
                                                     gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCag                   <  1:304715/71‑1 (MQ=255)
                                                     gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCag                   <  1:1857996/71‑1 (MQ=255)
                                                     gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCag                   <  1:3208573/71‑1 (MQ=255)
                                                      cACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAccgcc                                             <  1:1104317/44‑1 (MQ=255)
                                                       aCCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAcc                                    <  1:1521534/52‑1 (MQ=255)
                                                          aGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAcc                                    >  1:3144201/1‑49 (MQ=255)
                                                          aGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                                     >  1:81596/1‑48 (MQ=255)
                                                          aGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAgcgc              >  1:3139594/1‑71 (MQ=255)
                                                           gTAGAACCGCCAAAGAAGGCCATCGGTACGAATAc                                                 >  1:821858/1‑35 (MQ=255)
                                                            tAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                               <  1:3356690/52‑1 (MQ=255)
                                                             aGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                                        >  1:2280982/1‑42 (MQ=255)
                                                             aGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                                        >  1:3104664/1‑42 (MQ=255)
                                                             aGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGccc            >  1:3841881/1‑70 (MQ=255)
                                                             aGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt           >  1:1238850/1‑71 (MQ=255)
                                                              gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                               <  1:3635217/50‑1 (MQ=255)
                                                              gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa                             >  1:2131698/1‑52 (MQ=255)
                                                              gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCa                    >  1:551963/1‑61 (MQ=255)
                                                              gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGccc            >  1:3752865/1‑69 (MQ=255)
                                                              gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt           >  1:3295307/1‑70 (MQ=255)
                                                              gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt           >  1:3172011/1‑70 (MQ=255)
                                                              gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt           >  1:1065092/1‑70 (MQ=255)
                                                               aaCCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa         <  1:2121909/71‑1 (MQ=255)
                                                                aCCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa         <  1:907978/70‑1 (MQ=255)
                                                                aCCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa         <  1:2506293/70‑1 (MQ=255)
                                                                aCCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa         <  1:307462/70‑1 (MQ=255)
                                                                 ccGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAgcgc              >  1:3389224/1‑64 (MQ=255)
                                                                 ccGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAgcgc              >  1:2925843/1‑64 (MQ=255)
                                                                 ccGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAgcgc              >  1:3394083/1‑64 (MQ=255)
                                                                  cGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATc      <  1:2007500/71‑1 (MQ=255)
                                                                   gCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                                     >  1:2684275/1‑39 (MQ=255)
                                                                   gCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGaa        >  1:1976328/1‑68 (MQ=255)
                                                                     cAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCAt                                  >  1:859337/1‑40 (MQ=255)
                                                                      aaaGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt           >  1:917748/1‑62 (MQ=255)
                                                                      aaaGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa         <  1:1683033/64‑1 (MQ=255)
                                                                      aaaGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTccc  <  1:1776318/71‑1 (MQ=255)
                                                                      |                                                                      
GTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCC  >  minE/311426‑311566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: