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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,6910TC 8.2% 235.4 / 10.2 85T394T (ACA→ACGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (41/37);  new base C (2/5);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACG  >  minE/311621‑311760
                                                                      |                                                                     
gCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAtt                                                                       >  1:1042670/1‑71 (MQ=255)
  tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                                                                     >  1:249664/1‑71 (MQ=255)
      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)
               ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTa                                                                   <  1:1212382/60‑1 (MQ=255)
               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:988806/1‑69 (MQ=255)
                    aCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag                                                     >  1:2723166/1‑69 (MQ=255)
                      gATGGCGTCATCCACCATCAGGCCGATGGCGAGCTCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAg                                                  <  1:4105831/70‑1 (MQ=255)
                          gcatCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                              <  1:2358204/67‑1 (MQ=9)
                          gcatCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                              >  1:3728942/4‑70 (MQ=9)
                             tCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGtt                                                                    <  1:1742676/45‑1 (MQ=12)
                              cATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                                                          <  1:2515136/54‑1 (MQ=255)
                               aTCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTa                                                                   >  1:304189/1‑44 (MQ=255)
                               aTCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTa                                                                   >  1:3510059/1‑44 (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)
                                      aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAg                                                                  >  1:389126/1‑38 (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:3879823/1‑56 (MQ=255)
                                           gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                           >  1:93436/1‑56 (MQ=255)
                                           gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                            <  1:1963441/71‑1 (MQ=255)
                                           gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                            <  1:2469561/71‑1 (MQ=255)
                                           gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                            <  1:2540081/71‑1 (MQ=255)
                                           gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                            <  1:1182417/71‑1 (MQ=255)
                                           cccGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaaa                                             <  1:3136104/53‑1 (MQ=12)
                                            ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                             <  1:2602027/69‑1 (MQ=255)
                                            ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                             >  1:2581942/1‑69 (MQ=255)
                                             cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                             >  1:2716595/1‑68 (MQ=255)
                                             cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                          <  1:934981/71‑1 (MQ=255)
                                             cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                          >  1:774543/1‑71 (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:3103408/1‑67 (MQ=255)
                                                 ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                          >  1:2852074/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:2067735/1‑40 (MQ=255)
                                                           aTCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                           >  1:1158617/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:4087896/39‑1 (MQ=255)
                                                                 aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCa                                      <  1:2427382/39‑1 (MQ=255)
                                                                  aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg      >  1:2811118/1‑70 (MQ=255)
                                                                  aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg      <  1:3528148/70‑1 (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:716954/71‑1 (MQ=255)
                                                                     ttGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACg  <  1:1957162/71‑1 (MQ=255)
                                                                      |                                                                     
GCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACG  >  minE/311621‑311760

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: