Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,427 T→A 5.1% V482V (GTA→GTT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4270TA 5.1% 231.7 / 3.5 79V482V (GTA→GTTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (40/35);  new base A (2/2);  total (42/37)
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.95e-01

TGATCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCA  >  minE/311359‑311496
                                                                    |                                                                     
tGATCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAc                                                                     >  1:2016078/1‑71 (MQ=255)
tGATCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAc                                                                     >  1:442773/1‑71 (MQ=255)
  aTCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAcc                                                                    >  1:3283865/1‑70 (MQ=255)
  aTCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAcc                                                                    >  1:1667070/1‑70 (MQ=255)
  aTCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAcc                                                                    <  1:2097881/70‑1 (MQ=255)
  aTCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAcc                                                                    >  1:1611210/1‑70 (MQ=255)
  aTCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAcc                                                                    >  1:2741941/1‑70 (MQ=255)
  aTCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAcc                                                                    >  1:3489411/1‑70 (MQ=255)
  aTCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAcc                                                                    >  1:382858/1‑70 (MQ=255)
  aTCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAcc                                                                    >  1:441084/1‑70 (MQ=255)
    cGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGa                                                                  >  1:2675018/1‑70 (MQ=255)
     gCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTTGAGTCAGGATCAACGCCACCAGTACCGACa                                                                <  1:837090/71‑1 (MQ=255)
     gCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTa                                                                      <  1:2726375/65‑1 (MQ=255)
     gCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGa                                                                  >  1:359258/1‑69 (MQ=255)
        tttGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAgcg                                                             >  1:1413060/1‑71 (MQ=255)
           ggCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGcc                                                           <  1:3797385/70‑1 (MQ=255)
                tCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCATGATCAACGCCACCAGTACCGACAGCGCCa                                                          >  1:2586122/1‑66 (MQ=255)
                tCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCCCCAGTAc                                                                     >  1:4108672/1‑55 (MQ=255)
                tCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTAcc                                                                    <  1:3076650/56‑1 (MQ=255)
                tCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACATCGcc                                                           >  1:1030651/1‑65 (MQ=255)
                 cGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGa                                                                  >  1:100972/1‑57 (MQ=255)
                  ggTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTg                                                    >  1:876342/1‑70 (MQ=255)
                  ggTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTg                                                    >  1:1291885/1‑70 (MQ=255)
                      tCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                <  1:1169174/70‑1 (MQ=255)
                      tCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                <  1:1858629/70‑1 (MQ=255)
                      tCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                >  1:2111685/1‑70 (MQ=255)
                      tCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                >  1:3133498/1‑70 (MQ=255)
                      tCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                <  1:4109940/70‑1 (MQ=255)
                        aGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAAt                                             <  1:1121023/71‑1 (MQ=255)
                        aGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAAt                                             >  1:1570549/1‑71 (MQ=255)
                        aGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAAt                                             >  1:268024/1‑71 (MQ=255)
                           aTGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATgg                                           >  1:792499/1‑70 (MQ=255)
                           aTGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGt                                          <  1:3270665/71‑1 (MQ=255)
                              gtggCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTaa                                        <  1:3944650/70‑1 (MQ=255)
                                 gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAAt                                             >  1:1696424/1‑62 (MQ=255)
                                 gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                      >  1:2363516/1‑69 (MQ=255)
                                 gCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                      >  1:1001928/1‑69 (MQ=255)
                                  cacaAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTaa                                        <  1:2640482/66‑1 (MQ=255)
                                   acaAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAgcgc                                                            <  1:2569232/45‑1 (MQ=255)
                                   acaAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGa                                                   <  1:1627369/54‑1 (MQ=255)
                                   acaAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGa                                                   <  1:1246194/54‑1 (MQ=255)
                                    caAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTg                                                       <  1:1773358/49‑1 (MQ=255)
                                     aaaGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGa                                                   <  1:3397926/52‑1 (MQ=255)
                                     aaaGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGa                                                   <  1:3044367/52‑1 (MQ=255)
                                     aaaGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAAc                                                <  1:958860/55‑1 (MQ=255)
                                        gcgCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGa                             >  1:181615/3‑71 (MQ=255)
                                         agCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                            <  1:1056014/71‑1 (MQ=255)
                                              gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCa                                                          >  1:1183668/1‑36 (MQ=255)
                                              gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGa                             >  1:342043/1‑65 (MQ=255)
                                              gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgatag                       <  1:115343/71‑1 (MQ=255)
                                              gAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgat                         <  1:3740230/69‑1 (MQ=255)
                                               aGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAAt                                             <  1:10949/48‑1 (MQ=255)
                                                   aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag                   >  1:3684563/1‑70 (MQ=255)
                                                   aGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGatagatag                   >  1:362022/1‑70 (MQ=255)
                                                     gATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATa                                      <  1:2757028/49‑1 (MQ=255)
                                                       tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGt                                          <  1:543723/43‑1 (MQ=255)
                                                       tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATaga                                    >  1:425825/1‑49 (MQ=255)
                                                       tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                >  1:1124834/1‑69 (MQ=255)
                                                       tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCAc                >  1:3731864/1‑69 (MQ=255)
                                                       tCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa              <  1:3324/71‑1 (MQ=255)
                                                        cAATGCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAgaga                                  <  1:324753/50‑1 (MQ=255)
                                                        cAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                          >  1:941713/1‑58 (MQ=255)
                                                         aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAg             >  1:1952394/4‑70 (MQ=9)
                                                         aatgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGAGATAAATTGCCCCAg             >  1:239808/4‑70 (MQ=9)
                                                         aaCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATagag                                   <  1:1619759/48‑1 (MQ=255)
                                                         aaCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAg             <  1:3192083/70‑1 (MQ=255)
                                                          aCGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACg                           <  1:2226893/55‑1 (MQ=255)
                                                           tgCCACCAGAACAGAAAGCGCCATTGCCGAAACGATGGTGATAgaga                                  <  1:3478613/46‑1 (MQ=19)
                                                           cGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg          <  1:2575960/71‑1 (MQ=255)
                                                            gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCa              >  1:398186/1‑66 (MQ=255)
                                                            gCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTa           >  1:3591795/1‑69 (MQ=255)
                                                                ccaGTACCGACAGCGCCATTGCTGAAACAATGGTAATagag                                   <  1:1045963/41‑1 (MQ=255)
                                                                 cagaacagaaAGCGCCATTGCCGAAACGATGGTGATAGAGAAc                                >  1:76246/1‑43 (MQ=16)
                                                                 caGTACCGACAGCGCCATTGCTGAAACAATGGTAATATAGAACTGACgat                         <  1:2672565/50‑1 (MQ=255)
                                                                 caGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACgat                         <  1:2252326/50‑1 (MQ=255)
                                                                 caGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACga                          >  1:976084/1‑49 (MQ=255)
                                                                  aGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                            >  1:3109767/1‑46 (MQ=255)
                                                                  aGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                            >  1:1847989/1‑46 (MQ=255)
                                                                   gTACCGACAGCGCCATTGCTGAAACAATGGTAATag                                     <  1:3117052/36‑1 (MQ=255)
                                                                   gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa  <  1:798959/71‑1 (MQ=255)
                                                                    tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGAc                            <  1:742265/44‑1 (MQ=255)
                                                                    tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg          <  1:1599318/62‑1 (MQ=255)
                                                                    tACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAg          <  1:3384108/62‑1 (MQ=255)
                                                                    |                                                                     
TGATCGCCTTTGGCAATCGGTTTCAGCATGGTGGCACAAAGAGCTGGAGTCAGGATCAACGCCACCAGTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCA  >  minE/311359‑311496

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: