Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,790 G→A 12.6% N361N (AAC→AAT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7900GA12.6% 314.4 / 25.1 111N361N (AAC→AATacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (40/57);  new base A (7/7);  total (47/64)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 118 total aligned reads displayed.

AAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGT  >  minE/311724‑311860
                                                                  |                                                                      
aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt                                                                    <  1:1388741/71‑1 (MQ=255)
 aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                                   <  1:1954891/71‑1 (MQ=255)
 aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                                   >  1:2930044/1‑71 (MQ=255)
 aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                                   >  1:2082518/1‑71 (MQ=255)
 aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGGCGGAATCAACGTCGCGCGGAAGTTCTg                                                                   <  1:1104062/71‑1 (MQ=255)
    aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                                                >  1:612465/1‑71 (MQ=255)
    aCGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAg                                                                >  1:3481233/1‑71 (MQ=255)
      ggCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGa                                                              <  1:3631523/71‑1 (MQ=255)
          aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                                          >  1:5215938/1‑71 (MQ=11)
          aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                                          >  1:2637157/1‑71 (MQ=11)
          aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                                          >  1:5237249/1‑71 (MQ=11)
           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                                        <  1:4473877/56‑1 (MQ=255)
           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                                        <  1:4773430/56‑1 (MQ=255)
           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                                        <  1:2961535/56‑1 (MQ=255)
           aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                                        <  1:2566969/56‑1 (MQ=255)
            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)
                   aacAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                 >  1:4922309/3‑71 (MQ=2)
                   aacAACACCACGGGTACCGCAATGGTGGGGATCAGAGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                 >  1:5222363/3‑71 (MQ=2)
                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                             >  1:2780819/1‑59 (MQ=255)
                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                             >  1:2609659/1‑59 (MQ=255)
                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                             >  1:2727410/1‑59 (MQ=255)
                   aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                                                 <  1:2970212/71‑1 (MQ=255)
                    gCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAcc                                                <  1:1008852/71‑1 (MQ=255)
                    acAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                                 >  1:1975921/2‑70 (MQ=2)
                    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)
                          ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                           >  1:4377344/1‑70 (MQ=255)
                          ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAac                                          <  1:403968/71‑1 (MQ=255)
                           caccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                                            <  1:5398554/52‑1 (MQ=255)
                            acgggTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAg                                              >  1:3124164/4‑65 (MQ=2)
                            accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                             <  1:4882378/66‑1 (MQ=255)
                            accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                                        <  1:3584368/71‑1 (MQ=255)
                            accGGCACCGCGATCGTCGGGATCAGCGTCGCCCGCCAATTCTGCAGGAACAGGAACATCACCAGGAagac                                        <  1:4785712/71‑1 (MQ=2)
                             ccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAg                                                          >  1:266931/1‑52 (MQ=255)
                             ccGGTACGGCAATGGTAGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                                                           >  1:4796944/1‑51 (MQ=255)
                               ggTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                     >  1:1448485/1‑71 (MQ=255)
                                 tACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGga                                   >  1:2878751/1‑71 (MQ=255)
                                  aCGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgat                                  >  1:4248622/1‑71 (MQ=255)
                                   ccgCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGcat                                  <  1:786931/68‑1 (MQ=9)
                                   ccgCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGcat                                  <  1:3775086/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)
                                        aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGa                                                         >  1:1807280/1‑42 (MQ=255)
                                        aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                                    <  1:1387417/47‑1 (MQ=255)
                                        aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                                       <  1:2180989/60‑1 (MQ=255)
                                         tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                             >  1:1701577/1‑37 (MQ=255)
                                         tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                             <  1:1423131/53‑1 (MQ=255)
                                         tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                     <  1:1524122/61‑1 (MQ=255)
                                         tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                     <  1:3730336/61‑1 (MQ=255)
                                         tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                     <  1:3283076/61‑1 (MQ=255)
                                         tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                     <  1:3154050/61‑1 (MQ=255)
                                         tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                     <  1:5267957/61‑1 (MQ=255)
                                          ggTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                                     >  1:4168476/1‑60 (MQ=255)
                                          ggTCGGAATCAACGTCGCGAGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCtt                          >  1:5645169/1‑71 (MQ=255)
                                             cGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                       >  1:1014996/1‑71 (MQ=255)
                                                aaTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                    <  1:5031656/71‑1 (MQ=255)
                                                 aTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAg                   >  1:2390574/1‑71 (MQ=255)
                                                 aTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAg                   >  1:1071591/1‑71 (MQ=255)
                                                  tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                                       <  1:1713518/50‑1 (MQ=255)
                                                  tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                 <  1:2094437/56‑1 (MQ=255)
                                                  tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                                 <  1:1918509/56‑1 (MQ=255)
                                                    aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                 <  1:3908644/70‑1 (MQ=255)
                                                      cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                             <  1:4149973/40‑1 (MQ=255)
                                                      cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt              >  1:5371946/1‑71 (MQ=255)
                                                      cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt              >  1:2754883/1‑71 (MQ=255)
                                                      cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt              >  1:2748303/1‑71 (MQ=255)
                                                      cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt              >  1:1417996/1‑71 (MQ=255)
                                                       gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt             <  1:950230/71‑1 (MQ=255)
                                                       gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt             <  1:4426564/71‑1 (MQ=255)
                                                       gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt             <  1:4208913/71‑1 (MQ=255)
                                                       gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt             <  1:3234591/71‑1 (MQ=255)
                                                       gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt             <  1:3121986/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)
                                                          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)
                                                                  |                                                                      
AAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGT  >  minE/311724‑311860

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: