Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,798 G→A 13.9% L359L (CTG→TTG)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7980GA13.9% 331.5 / 32.3 122L359L (CTG→TTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (39/66);  new base A (11/6);  total (50/72)
Fisher's exact test for biased strand distribution p-value = 3.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 124 total aligned reads displayed.

ACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGA  >  minE/311728‑311867
                                                                      |                                                                     
aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                                                       >  1:612465/1‑71 (MQ=255)
aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                                                       >  1:3481233/1‑71 (MQ=255)
  ggCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGa                                                                     <  1:3631523/71‑1 (MQ=255)
      aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                                                 >  1:2637157/1‑71 (MQ=11)
      aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                                                 >  1:5237249/1‑71 (MQ=11)
      aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                                                 >  1:5215938/1‑71 (MQ=11)
        ggTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGa                                                                >  1:1236659/1‑70 (MQ=255)
            ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                                                                  <  1:122985/64‑1 (MQ=255)
            ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACa                                                            <  1:3948426/70‑1 (MQ=255)
             cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAcatacata                                                          <  1:128063/71‑1 (MQ=255)
             cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                                                          <  1:2566942/71‑1 (MQ=255)
               aacAACACCACGGGTACCGCAATGGTGGGGATCAGAGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                        >  1:5222363/3‑71 (MQ=2)
               aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                                    >  1:2780819/1‑59 (MQ=255)
               aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                                                        <  1:2970212/71‑1 (MQ=255)
                gCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAcc                                                       <  1:1008852/71‑1 (MQ=255)
                acAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                        >  1:1661494/2‑70 (MQ=2)
                  aacaccacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAg                                                     >  1:4596896/4‑71 (MQ=2)
                  aaTACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                           <  1:511256/65‑1 (MQ=255)
                   acaccacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAgg                                                    <  1:4670260/69‑1 (MQ=9)
                     accacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAgg                                                    <  1:1106224/69‑1 (MQ=9)
                      ccaccGGTACGGCAATGGTCGGATTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                                  >  1:4568543/1‑70 (MQ=255)
                      ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                           >  1:4393492/1‑61 (MQ=255)
                      ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAac                                                 <  1:403968/71‑1 (MQ=255)
                      ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGAAGGAACAg                                                                 >  1:2907685/1‑55 (MQ=255)
                       caccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                                                   <  1:5398554/52‑1 (MQ=255)
                        acgggTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAg                                                     >  1:3124164/4‑65 (MQ=2)
                        acgggTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAacac                                               <  1:5493422/68‑1 (MQ=9)
                        accGGTACGGCAATGGTCGGAATCAACTTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                                               <  1:4566776/71‑1 (MQ=255)
                        accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                                    <  1:4882378/66‑1 (MQ=255)
                        accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                                               <  1:3584368/71‑1 (MQ=255)
                         ccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAg                                                                 >  1:266931/1‑52 (MQ=255)
                         ccGGTACGGCAATGGTAGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                                                                  >  1:4796944/1‑51 (MQ=255)
                           ggTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                            >  1:1448485/1‑71 (MQ=255)
                              aCGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgat                                         >  1:4248622/1‑71 (MQ=255)
                               ccgCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGcat                                         <  1:786931/68‑1 (MQ=9)
                               cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                        <  1:4147432/71‑1 (MQ=255)
                               cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                        <  1:109107/71‑1 (MQ=255)
                                 gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                                                   >  1:2497879/1‑42 (MQ=255)
                                 gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                                                   >  1:1191143/1‑42 (MQ=255)
                                 gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatgat                                      <  1:1848565/71‑1 (MQ=255)
                                 gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatgat                                      <  1:5089969/71‑1 (MQ=255)
                                  cAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAcc                                                       <  1:1768456/53‑1 (MQ=2)
                                   aaTGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGACCACCAGcataatg                                     >  1:5560352/1‑69 (MQ=2)
                                   aaTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATaa                                                         >  1:2928644/1‑50 (MQ=255)
                                    aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                           <  1:1387417/47‑1 (MQ=255)
                                    aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                           <  1:1821335/47‑1 (MQ=255)
                                    aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                                                          <  1:2646836/48‑1 (MQ=255)
                                    aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                                              <  1:2180989/60‑1 (MQ=255)
                                     tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                                    >  1:1701577/1‑37 (MQ=255)
                                     tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                            <  1:3283076/61‑1 (MQ=255)
                                     tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                            <  1:3730336/61‑1 (MQ=255)
                                     tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                            <  1:1524122/61‑1 (MQ=255)
                                     tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                            <  1:2015780/61‑1 (MQ=255)
                                      ggTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                            >  1:4168476/1‑60 (MQ=255)
                                       gTCGGATTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                            <  1:4376172/59‑1 (MQ=255)
                                         cGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                              >  1:1014996/1‑71 (MQ=255)
                                            aaTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                           <  1:5031656/71‑1 (MQ=255)
                                             aTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAg                          >  1:1071591/1‑71 (MQ=255)
                                              tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                                              <  1:1713518/50‑1 (MQ=255)
                                              tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                        <  1:1918509/56‑1 (MQ=255)
                                              tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                        <  1:2094437/56‑1 (MQ=255)
                                                aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                        <  1:3908644/70‑1 (MQ=255)
                                                  cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                                    <  1:4869205/40‑1 (MQ=255)
                                                  cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                     >  1:2748303/1‑71 (MQ=255)
                                                  cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                     >  1:2754883/1‑71 (MQ=255)
                                                  cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                     >  1:5371946/1‑71 (MQ=255)
                                                  cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                     >  1:1417996/1‑71 (MQ=255)
                                                   gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                    <  1:4426564/71‑1 (MQ=255)
                                                   gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                    <  1:950230/71‑1 (MQ=255)
                                                   gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                    <  1:4208913/71‑1 (MQ=255)
                                                   gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                    <  1:1484380/71‑1 (MQ=255)
                                                   gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                    <  1:3121986/71‑1 (MQ=255)
                                                   gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                    <  1:3234591/71‑1 (MQ=255)
                                                    tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                        <  1:1937636/66‑1 (MQ=255)
                                                    tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                   <  1:3918965/71‑1 (MQ=255)
                                                    tCGCGCGGAAGTTCTGAAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                        <  1:2248828/66‑1 (MQ=255)
                                                      gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                  >  1:1197818/1‑70 (MQ=255)
                                                      gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                  >  1:3370189/1‑70 (MQ=255)
                                                      gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                 <  1:2932654/71‑1 (MQ=255)
                                                         cgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCAc              <  1:4399059/71‑1 (MQ=255)
                                                          ggAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACt             >  1:3698759/1‑71 (MQ=255)
                                                           gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAgg                                           <  1:3898201/40‑1 (MQ=255)
                                                           gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAgg                                           <  1:1332773/40‑1 (MQ=255)
                                                             aGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                              >  1:3504734/1‑51 (MQ=255)
                                                             aGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                              >  1:2479405/1‑51 (MQ=255)
                                                              gTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt         >  1:3065788/1‑71 (MQ=255)
                                                               ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                           <  1:4895700/52‑1 (MQ=255)
                                                               ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                           <  1:5266851/52‑1 (MQ=255)
                                                               ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg        <  1:5126597/71‑1 (MQ=255)
                                                               ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg        <  1:4650009/71‑1 (MQ=255)
                                                               ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg        <  1:1448128/71‑1 (MQ=255)
                                                               ttCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCg                               >  1:4140710/1‑48 (MQ=12)
                                                               ttCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGt         >  1:2819597/1‑70 (MQ=18)
                                                               ttCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGt         >  1:4502896/1‑70 (MQ=18)
                                                                tCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg        >  1:4024720/1‑70 (MQ=255)
                                                                tCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTg        >  1:3099907/1‑70 (MQ=18)
                                                                 cTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGc                         >  1:2516665/1‑52 (MQ=255)
                                                                 cTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa      >  1:2686001/1‑71 (MQ=255)
                                                                    cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAAt     <  1:207952/69‑1 (MQ=255)
                                                                     aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga  <  1:836204/71‑1 (MQ=255)
                                                                      |                                                                     
ACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGA  >  minE/311728‑311867

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: