Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,697 C→T 7.1% T392T (ACG→ACA acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,6970CT 7.1% 279.5 / 7.8 84T392T (ACG→ACAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (42/36);  new base T (2/4);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGG  >  minE/311627‑311767
                                                                      |                                                                      
cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                                                        <  1:1421849/71‑1 (MQ=255)
  ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCTAGCACCATCCCGAACATTGTTAGc                                                                        <  1:866525/69‑1 (MQ=255)
  ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                                                                        <  1:2887343/69‑1 (MQ=255)
     tCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                                                                   >  1:4134395/1‑71 (MQ=255)
       caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                                                                 <  1:2056628/70‑3 (MQ=9)
         ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGttt                                                                  <  1:1844710/68‑1 (MQ=255)
         ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                                                 <  1:1766885/69‑1 (MQ=255)
           acaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCga                                                             <  1:3972220/71‑1 (MQ=255)
            caACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATc                                                               >  1:2141354/1‑68 (MQ=255)
            caACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCACGAACATTGTTAGCGTGttt                                                                  >  1:1351879/1‑65 (MQ=255)
             aaCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga                                                           <  1:2470879/71‑1 (MQ=255)
              aCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag                                                            >  1:2723166/1‑69 (MQ=255)
              aCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag                                                            >  1:988806/1‑69 (MQ=255)
                gATGGCGTCATCCACCATCAGGCCGATGGCGAGCTCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAg                                                         <  1:4105831/70‑1 (MQ=255)
                    gcatCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                                     >  1:3728942/4‑70 (MQ=9)
                    gcatCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                                     <  1:2358204/67‑1 (MQ=9)
                        cATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                                                 <  1:2515136/54‑1 (MQ=255)
                            caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtg                                                                     <  1:78444/46‑1 (MQ=255)
                            caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag                                                            >  1:1643518/1‑55 (MQ=255)
                            caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                              >  1:1610436/1‑69 (MQ=255)
                            caccaACAGGCCGATGGCAAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                              >  1:2142968/1‑69 (MQ=255)
                             accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                           >  1:3044013/1‑71 (MQ=255)
                                aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                        >  1:1331522/1‑71 (MQ=255)
                                 aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag                                                            <  1:2436159/50‑1 (MQ=255)
                                 aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                                       >  1:176189/1‑71 (MQ=255)
                                  cAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtg                                                                     >  1:1848640/1‑40 (MQ=255)
                                  cAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGc                                                        >  1:420902/1‑53 (MQ=255)
                                   aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag                                                            <  1:779333/48‑1 (MQ=255)
                                   aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                         <  1:1791344/67‑1 (MQ=255)
                                    ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                                     >  1:3299922/1‑70 (MQ=255)
                                    ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                                     >  1:190568/1‑70 (MQ=255)
                                    gcccGATGGCAAGCACCATCCCGAACATCGTTAGTGTGtt                                                                   <  1:648351/38‑1 (MQ=21)
                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                                  >  1:93436/1‑56 (MQ=255)
                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                                  >  1:3879823/1‑56 (MQ=255)
                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                   <  1:1182417/71‑1 (MQ=255)
                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                   <  1:2469561/71‑1 (MQ=255)
                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                   <  1:1963441/71‑1 (MQ=255)
                                     gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                   <  1:2540081/71‑1 (MQ=255)
                                     cccGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaaa                                                    <  1:3136104/53‑1 (MQ=12)
                                      ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                    >  1:2581942/1‑69 (MQ=255)
                                      ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                    <  1:2602027/69‑1 (MQ=255)
                                       cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                    >  1:2716595/1‑68 (MQ=255)
                                       cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                 >  1:774543/1‑71 (MQ=255)
                                       cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                 <  1:934981/71‑1 (MQ=255)
                                       cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                 <  1:2415528/71‑1 (MQ=255)
                                        gATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                           <  1:2622223/60‑1 (MQ=255)
                                         aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGt                               <  1:2688022/71‑1 (MQ=255)
                                           ggggAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                             >  1:1622186/4‑71 (MQ=255)
                                           ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                 >  1:2852074/1‑67 (MQ=255)
                                           ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                 >  1:3103408/1‑67 (MQ=255)
                                           ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                             >  1:1548159/1‑71 (MQ=255)
                                               aGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTccctaa                         >  1:2973426/1‑68 (MQ=12)
                                                gCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAg                                                   <  1:2314351/44‑1 (MQ=255)
                                                  aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCaa                      <  1:1970270/71‑1 (MQ=255)
                                                  aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCaa                      <  1:3817477/71‑1 (MQ=255)
                                                   ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCa                       <  1:2346159/69‑1 (MQ=255)
                                                    cATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAAt                     <  1:3493206/70‑1 (MQ=255)
                                                     aTCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                                  >  1:1158617/1‑40 (MQ=255)
                                                     aTCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                                  >  1:2067735/1‑40 (MQ=255)
                                                     aTCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                        <  1:3587175/66‑1 (MQ=255)
                                                      tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                         >  1:573019/1‑48 (MQ=255)
                                                      tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                   <  1:1937508/54‑1 (MQ=255)
                                                      tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                    <  1:3905647/53‑1 (MQ=255)
                                                      tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAgg                                <  1:2078536/57‑1 (MQ=255)
                                                      tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATac                   >  1:697547/1‑70 (MQ=255)
                                                       cccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATacca                 >  1:3915903/1‑71 (MQ=255)
                                                        ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCaa                                            >  1:1206526/1‑43 (MQ=255)
                                                        ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                         >  1:1870261/1‑46 (MQ=255)
                                                        ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCaa                      >  1:3004284/1‑65 (MQ=255)
                                                           aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCa                                             <  1:2427382/39‑1 (MQ=255)
                                                           aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCa                                             <  1:4087896/39‑1 (MQ=255)
                                                            aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg             <  1:3528148/70‑1 (MQ=255)
                                                            aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg             >  1:2811118/1‑70 (MQ=255)
                                                            aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGt            <  1:2309965/71‑1 (MQ=255)
                                                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                           <  1:4055097/39‑1 (MQ=255)
                                                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCaa                      >  1:1072439/1‑60 (MQ=255)
                                                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAAt                     >  1:798581/1‑61 (MQ=255)
                                                             cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTa           >  1:2668064/1‑71 (MQ=255)
                                                              aTTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTAc          >  1:1860586/1‑71 (MQ=255)
                                                               ttGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACg         <  1:1957162/71‑1 (MQ=255)
                                                               ttGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACg         <  1:716954/71‑1 (MQ=255)
                                                                 gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACg                                       >  1:3439289/1‑39 (MQ=255)
                                                                 gTTAGCGTGTTTATCGAGAAGCAAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACgg        >  1:3266442/1‑70 (MQ=255)
                                                                   tAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCaa     >  1:3605295/1‑71 (MQ=255)
                                                                    aGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAAt    >  1:1611211/1‑71 (MQ=255)
                                                                    aGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAAt    >  1:3411488/1‑71 (MQ=255)
                                                                     gCGTGTTAATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTAc          <  1:1317124/64‑1 (MQ=255)
                                                                      cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCAACGGTACGGCAATgg  <  1:885078/71‑1 (MQ=255)
                                                                      |                                                                      
CGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGG  >  minE/311627‑311767

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: