Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,840 C→A 8.1% V345F (GTC→TTC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8400CA 8.1% 284.6 / 6.4 86V345F (GTC→TTC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (35/44);  new base A (2/5);  total (37/49)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTT  >  minE/311770‑311909
                                                                      |                                                                     
ggAATCAACGTCTCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGAc                                                                       >  1:3416269/1‑71 (MQ=255)
           cgcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                            >  1:3432227/1‑71 (MQ=255)
           cgcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                            >  1:706391/1‑71 (MQ=255)
            gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                            >  1:3644790/1‑70 (MQ=255)
            gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                           <  1:2472535/71‑1 (MQ=255)
            gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                           >  1:3917369/1‑71 (MQ=255)
            gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                           >  1:2901924/1‑71 (MQ=255)
               cgGAAGTTCTGCAGGAACAGATACATAACCCGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCAc                                                        <  1:3786173/71‑1 (MQ=255)
                ggAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACt                                                       >  1:2842000/1‑71 (MQ=255)
                  atattCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTc                                                     <  1:4132534/68‑1 (MQ=16)
                  atattCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTc                                                     <  1:1257570/68‑1 (MQ=16)
                    attCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGt                                                   >  1:3806065/2‑71 (MQ=18)
                     ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                  <  1:3051962/71‑1 (MQ=255)
                     ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                  <  1:3151805/71‑1 (MQ=255)
                       cTGCAGGAACAGAAACATAACCAGGAACACGAGGAAGATCGCTTCGACCAGCGTTTTAACCACt                                                       >  1:2082536/1‑64 (MQ=255)
                          caagAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGaa                                                <  1:2376885/65‑1 (MQ=21)
                          cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                                >  1:3926964/1‑68 (MQ=255)
                           aagAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGaa                                                <  1:3641581/65‑1 (MQ=21)
                            ggAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATagag                                           <  1:3727808/71‑1 (MQ=255)
                                cAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTAACCACTTCGTGAATAGAGAtttt                                      <  1:282269/71‑1 (MQ=255)
                                 aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                      <  1:1906737/71‑1 (MQ=255)
                                 aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAttt                                       >  1:3224490/1‑70 (MQ=255)
                                  gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                                               >  1:3809363/1‑45 (MQ=255)
                                  gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTc                                     <  1:1770351/71‑1 (MQ=255)
                                   aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa                                              >  1:3316468/1‑61 (MQ=255)
                                   aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                      <  1:360248/69‑1 (MQ=255)
                                    tACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc                                   <  1:3823302/71‑1 (MQ=255)
                                     aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                                                     <  1:2960730/52‑1 (MQ=255)
                                     aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                                                     <  1:2870579/52‑1 (MQ=255)
                                     aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                                                     <  1:1019695/52‑1 (MQ=255)
                                     aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt                                                   <  1:3108121/54‑1 (MQ=255)
                                     aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt                                                   <  1:683765/54‑1 (MQ=255)
                                     aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                  <  1:2335233/55‑1 (MQ=255)
                                     aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                      >  1:2972066/1‑67 (MQ=255)
                                     aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc                                   >  1:4121586/1‑70 (MQ=255)
                                       aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                  <  1:2657994/53‑1 (MQ=255)
                                       aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa                                >  1:906935/1‑71 (MQ=255)
                                       aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa                                >  1:442460/1‑71 (MQ=255)
                                       aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa                                >  1:66945/1‑71 (MQ=255)
                                        tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc                               >  1:2112705/1‑71 (MQ=255)
                                        tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc                               >  1:2195688/1‑71 (MQ=255)
                                         aaCCAGGAACACGTGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                              >  1:1061898/1‑71 (MQ=255)
                                          aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                           >  1:217126/1‑41 (MQ=255)
                                          aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                      >  1:3868383/1‑62 (MQ=255)
                                            cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                           <  1:618998/39‑1 (MQ=255)
                                            cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                           <  1:3437539/39‑1 (MQ=255)
                                            cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa                                              >  1:2920808/1‑52 (MQ=255)
                                            cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                                  <  1:1361033/64‑1 (MQ=255)
                                            cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCg                           <  1:2200701/71‑1 (MQ=255)
                                               gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATagag                                           <  1:86175/52‑1 (MQ=255)
                                               gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc                                   >  1:2062484/1‑60 (MQ=255)
                                               gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtgg                        >  1:2425559/1‑71 (MQ=255)
                                                aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACTACTTCGTGAATAGAGATTTTCACGAACGTCgtgg                        <  1:1653084/70‑1 (MQ=255)
                                                 acacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggtg                      <  1:2719486/71‑1 (MQ=255)
                                                  cacGAGGATGATCGCTTCGACCAGCGTTTTAACCAc                                                        >  1:2161213/1‑36 (MQ=255)
                                                  cacGAGGATGATCGCTTCGACCAGCGTTTTAACCAc                                                        >  1:2148393/1‑36 (MQ=255)
                                                  cacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                                >  1:3436375/1‑44 (MQ=255)
                                                  cacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggtg                      <  1:770770/70‑1 (MQ=255)
                                                    cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                                <  1:2813876/42‑1 (MQ=255)
                                                    cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                   <  1:3546563/71‑1 (MQ=255)
                                                    cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                   <  1:1569925/71‑1 (MQ=255)
                                                    cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                   <  1:3830764/71‑1 (MQ=255)
                                                     gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCg                           <  1:3222918/62‑1 (MQ=255)
                                                     gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                   <  1:2843283/70‑1 (MQ=255)
                                                      aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                 <  1:3447892/71‑1 (MQ=255)
                                                        gatgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATag                                             <  1:1701625/41‑1 (MQ=255)
                                                          tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa                                >  1:1592630/1‑52 (MQ=255)
                                                          tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                              <  1:3144583/54‑1 (MQ=255)
                                                          tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATggg             >  1:2882058/1‑71 (MQ=255)
                                                           gatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                         <  1:2270954/58‑1 (MQ=255)
                                                            atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGt            <  1:2409392/67‑1 (MQ=255)
                                                            atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                         <  1:2593067/57‑1 (MQ=255)
                                                            atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGt            >  1:72631/1‑70 (MQ=255)
                                                              ggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTa           >  1:1326359/2‑69 (MQ=255)
                                                              cGCTTCGACCAGCGTTTTAACCACTTCGTGAATAgaga                                          >  1:3200662/1‑38 (MQ=255)
                                                              cGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTGACGAACGGCgtggt                       >  1:3754235/1‑57 (MQ=255)
                                                              cGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt                       >  1:2418844/1‑57 (MQ=255)
                                                               gCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                          <  1:3901708/53‑1 (MQ=255)
                                                               gCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                          <  1:2217799/53‑1 (MQ=255)
                                                                cTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                                  <  1:474611/44‑1 (MQ=255)
                                                                cTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                         <  1:1477670/53‑1 (MQ=255)
                                                                 ttCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTaaa         >  1:437715/1‑68 (MQ=255)
                                                                   cGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAAGGCGTGGTGTCGTATGGGTAAACa       >  1:2040389/1‑68 (MQ=255)
                                                                    gACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtt    <  1:3340499/70‑1 (MQ=255)
                                                                     aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt  <  1:1053984/71‑1 (MQ=255)
                                                                     aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt  <  1:3171169/71‑1 (MQ=255)
                                                                     aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtt    <  1:1015302/69‑1 (MQ=255)
                                                                      |                                                                     
GGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTT  >  minE/311770‑311909

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: