Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,781 C→T 10.9% A364A (GCG→GCA acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,7810CT10.9% 342.9 / 20.3 120A364A (GCG→GCAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (42/62);  new base T (7/6);  total (52/68)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

CAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTA  >  minE/311714‑311850
                                                                   |                                                                     
   aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTcgcgcg                                                                  >  1:4846186/1‑70 (MQ=255)
   aGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCgg                                                                 <  1:3663923/71‑1 (MQ=255)
    ggcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGa                                                                <  1:2655864/71‑1 (MQ=255)
     gcggcAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCACCGTCGCGCgg                                                                 >  1:3497152/1‑69 (MQ=255)
          aaGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCt                                                          <  1:1388741/71‑1 (MQ=255)
           aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                         >  1:2082518/1‑71 (MQ=255)
           aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                         <  1:1954891/71‑1 (MQ=255)
           aGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTg                                                         >  1:2930044/1‑71 (MQ=255)
                ggCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGa                                                    <  1:3631523/71‑1 (MQ=255)
                    aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                                >  1:5215938/1‑71 (MQ=11)
                    aaGGTCCCCAGCAGCACCACCGGTACGGCGATGGTCGGGATCAGCGTGGCGCGGAAGTTCTGCAGGAACAg                                                >  1:5237249/1‑71 (MQ=11)
                     aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                              <  1:2496179/56‑1 (MQ=255)
                     aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                              <  1:13541/56‑1 (MQ=255)
                     aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                              <  1:2566969/56‑1 (MQ=255)
                     aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                              <  1:4773430/56‑1 (MQ=255)
                     aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                              <  1:2961535/56‑1 (MQ=255)
                     aGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAg                                                              <  1:4473877/56‑1 (MQ=255)
                          ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                               >  1:4441074/1‑50 (MQ=255)
                          ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                                                 <  1:122985/64‑1 (MQ=255)
                          ccGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACa                                           <  1:3948426/70‑1 (MQ=255)
                           cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAcatacata                                         <  1:128063/71‑1 (MQ=255)
                           cGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                                         <  1:2566942/71‑1 (MQ=255)
                             aacAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                       >  1:4922309/3‑71 (MQ=2)
                             aacAACACCACGGGTACCGCAATGGTGGGGATCAGAGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                       >  1:5222363/3‑71 (MQ=2)
                             aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                   >  1:2727410/1‑59 (MQ=255)
                             aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                   >  1:2609659/1‑59 (MQ=255)
                             aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGaa                                                   >  1:2780819/1‑59 (MQ=255)
                             aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                                       <  1:966593/71‑1 (MQ=255)
                             aGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAc                                       <  1:2970212/71‑1 (MQ=255)
                              gCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAAcc                                      <  1:1008852/71‑1 (MQ=255)
                              acAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                       >  1:1975921/2‑70 (MQ=2)
                              acAACACCACGGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAc                                       >  1:1661494/2‑70 (MQ=2)
                                aacaccacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAg                                    >  1:4596896/4‑71 (MQ=2)
                                aaTACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                          <  1:511256/65‑1 (MQ=255)
                                   accacgGGTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAgg                                   <  1:1106224/69‑1 (MQ=9)
                                   accaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                               >  1:822193/1‑41 (MQ=255)
                                   accaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                               >  1:295208/1‑41 (MQ=255)
                                   accaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                               >  1:2786853/1‑41 (MQ=255)
                                    ccaccGGTACGGCAATGGTCGGATTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                 >  1:4568543/1‑70 (MQ=255)
                                    ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                               >  1:1083101/1‑40 (MQ=255)
                                    ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGaa                                                               >  1:1165736/1‑40 (MQ=255)
                                    ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                          >  1:4393492/1‑61 (MQ=255)
                                    ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGaa                                 >  1:4377344/1‑70 (MQ=255)
                                    ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAac                                <  1:403968/71‑1 (MQ=255)
                                    ccaccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGAAGGAACAg                                                >  1:2907685/1‑55 (MQ=255)
                                     caccGGTACGGCAATGGTCGGAATCAACGTCGCGCgg                                                                 >  1:4014668/1‑37 (MQ=255)
                                     caccGGTACGGCAATGGTCGGAATCAACGTCGCGCgg                                                                 >  1:5433506/1‑37 (MQ=255)
                                     caccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                                  <  1:5398554/52‑1 (MQ=255)
                                      acgggTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAg                                    >  1:3124164/4‑65 (MQ=2)
                                      acgggTACCGCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAacac                              <  1:5493422/68‑1 (MQ=9)
                                      accGGTACGGCAATGGTCGGAATCAACTTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                              <  1:4566776/71‑1 (MQ=255)
                                      accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                   <  1:4882378/66‑1 (MQ=255)
                                      accGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAacac                              <  1:3584368/71‑1 (MQ=255)
                                      accGGCACCGCGATCGTCGGGATCAGCGTCGCCCGCCAATTCTGCAGGAACAGGAACATCACCAGGAagac                              <  1:4785712/71‑1 (MQ=2)
                                       ccGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAg                                                >  1:266931/1‑52 (MQ=255)
                                       ccGGTACGGCAATGGTAGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACa                                                 >  1:4796944/1‑51 (MQ=255)
                                         ggTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                           >  1:1448485/1‑71 (MQ=255)
                                           tACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGga                         >  1:2878751/1‑71 (MQ=255)
                                            aCGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgat                        >  1:4248622/1‑71 (MQ=255)
                                             ccgCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGcat                        <  1:786931/68‑1 (MQ=9)
                                             ccgCAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGcat                        <  1:3775086/68‑1 (MQ=9)
                                             cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                       <  1:109107/71‑1 (MQ=255)
                                             cGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                       <  1:4147432/71‑1 (MQ=255)
                                               gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                                  >  1:2497879/1‑42 (MQ=255)
                                               gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAAc                                                  >  1:1191143/1‑42 (MQ=255)
                                               gCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatgat                     <  1:5089969/71‑1 (MQ=255)
                                                cAATGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCAcc                                      <  1:1768456/53‑1 (MQ=2)
                                                 aaTGGTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGACCACCAGcataatg                    >  1:5560352/1‑69 (MQ=2)
                                                  aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGa                                               >  1:1807280/1‑42 (MQ=255)
                                                  aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                          <  1:1387417/47‑1 (MQ=255)
                                                  aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACAt                                          <  1:1821335/47‑1 (MQ=255)
                                                  aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATa                                         <  1:2646836/48‑1 (MQ=255)
                                                  aTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                             <  1:2180989/60‑1 (MQ=255)
                                                   tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                   <  1:1423131/53‑1 (MQ=255)
                                                   tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                           <  1:2015780/61‑1 (MQ=255)
                                                   tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                           <  1:3730336/61‑1 (MQ=255)
                                                   tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                           <  1:5267957/61‑1 (MQ=255)
                                                   tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                           <  1:1524122/61‑1 (MQ=255)
                                                   tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                           <  1:3154050/61‑1 (MQ=255)
                                                   tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                           <  1:3283076/61‑1 (MQ=255)
                                                    ggTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                           >  1:4168476/1‑60 (MQ=255)
                                                    ggTCGGAATCAACGTCGCGAGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCtt                >  1:5645169/1‑71 (MQ=255)
                                                     gTCGGATTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAg                           <  1:4376172/59‑1 (MQ=255)
                                                       cGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa             >  1:1014996/1‑71 (MQ=255)
                                                          aaTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa          <  1:5031656/71‑1 (MQ=255)
                                                           aTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAg         >  1:1071591/1‑71 (MQ=255)
                                                            tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACg                             <  1:1713518/50‑1 (MQ=255)
                                                            tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                       <  1:2094437/56‑1 (MQ=255)
                                                            tCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGgatg                       <  1:1918509/56‑1 (MQ=255)
                                                              aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg       <  1:3908644/70‑1 (MQ=255)
                                                                cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAgg                                   <  1:4149973/40‑1 (MQ=255)
                                                                cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt    >  1:2754883/1‑71 (MQ=255)
                                                                 gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt   <  1:3234591/71‑1 (MQ=255)
                                                                 gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt   <  1:3121986/71‑1 (MQ=255)
                                                                 gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt   <  1:1484380/71‑1 (MQ=255)
                                                                 gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt   <  1:950230/71‑1 (MQ=255)
                                                                  tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg       <  1:1937636/66‑1 (MQ=255)
                                                                  tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa  <  1:3918965/71‑1 (MQ=255)
                                                                  tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa  <  1:1121907/71‑1 (MQ=255)
                                                                  tCGCGCGGAAGTTCTGAAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg       <  1:2248828/66‑1 (MQ=255)
                                                                   |                                                                     
CAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTA  >  minE/311714‑311850

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: