Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,673 G→A 9.6% L400L (CTC→CTT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,6730GA 9.6% 228.5 / 10.8 73L400L (CTC→CTTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base A (2/5);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

TCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCC  >  minE/311605‑311740
                                                                    |                                                                   
tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                                                   >  1:286891/1‑71 (MQ=255)
tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                                                   >  1:2272152/1‑71 (MQ=255)
 cttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                                                   >  1:2593503/1‑70 (MQ=255)
 cttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                                                   <  1:2022659/70‑1 (MQ=255)
    ccGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGc                                                                  <  1:158362/68‑1 (MQ=255)
     cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCAcc                                                               <  1:2438285/70‑1 (MQ=255)
          tAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCg                                                                     <  1:2766905/59‑1 (MQ=255)
           aaCACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATc                                                            >  1:495542/1‑67 (MQ=255)
            acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGa                                                                    >  1:999226/1‑58 (MQ=255)
            acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                                       >  1:2754868/1‑71 (MQ=255)
            acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                                       >  1:4083210/1‑71 (MQ=255)
              acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                                     <  1:277212/71‑1 (MQ=255)
              acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                                                     <  1:1344686/71‑1 (MQ=255)
               cGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCg                                                                     >  1:1691291/1‑54 (MQ=255)
               cGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCg                                                         <  1:1523054/66‑1 (MQ=255)
                gCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAtt                                                   >  1:1042670/1‑71 (MQ=255)
                  tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCa                                                                 >  1:2775837/1‑55 (MQ=255)
                  tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGt                                                 >  1:249664/1‑71 (MQ=255)
                      cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCa                                                              <  1:3409455/54‑1 (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:3510059/1‑44 (MQ=255)
                                               aTCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTa                                               >  1:304189/1‑44 (MQ=255)
                                                 ccaccaACAGGCCGATGGCGAGCACCATCCCGAACa                                                     <  1:1693656/36‑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)
                                                      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:190568/1‑70 (MQ=255)
                                                          ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc          >  1:3299922/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:2469561/71‑1 (MQ=255)
                                                           gCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa        <  1:1182417/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:2602027/69‑1 (MQ=255)
                                                            ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa         >  1:2581942/1‑69 (MQ=255)
                                                             cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa         >  1:2716595/1‑68 (MQ=255)
                                                             cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg      <  1:2415528/71‑1 (MQ=255)
                                                             cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg      >  1:774543/1‑71 (MQ=255)
                                                             cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg      <  1:934981/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)
                                                                    |                                                                   
TCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCC  >  minE/311605‑311740

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: