Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,584,604 C→T 8.5% E891E (GAG→GAA mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,584,6040CT 8.5% 263.4 / 10.4 82E891E (GAG→GAAmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base C (41/34);  new base T (5/2);  total (46/36)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

CGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCAGGAGCCTGATTT  >  minE/2584534‑2584673
                                                                      |                                                                     
cGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctc                                                                       >  1:1970207/1‑71 (MQ=255)
cGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctc                                                                       >  1:3018857/1‑71 (MQ=255)
 ggTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctct                                                                      <  1:1278301/71‑1 (MQ=255)
 ggTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctct                                                                      <  1:1498835/71‑1 (MQ=255)
 ggTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctc                                                                       <  1:3950942/70‑1 (MQ=255)
 ggTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctc                                                                       <  1:2496733/70‑1 (MQ=255)
  gTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctct                                                                      <  1:2787540/70‑1 (MQ=255)
  gTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCtctc                                                                     <  1:2232431/71‑1 (MQ=255)
   tGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCa                                                                    <  1:4005309/71‑1 (MQ=255)
   tGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCa                                                                    <  1:2627101/71‑1 (MQ=255)
      ccAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATag                                                                 <  1:4033511/71‑1 (MQ=255)
      ccAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATag                                                                 <  1:1324483/71‑1 (MQ=255)
       cAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATa                                                                  <  1:1739454/69‑1 (MQ=255)
        aGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGATTTGACCAGCTCTCATa                                                                  <  1:1353589/68‑1 (MQ=255)
          tAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCa                                                                    >  1:2635989/1‑64 (MQ=255)
          tAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCa                                                                    >  1:3209795/1‑64 (MQ=255)
          tAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGt                                                              >  1:3538628/1‑70 (MQ=255)
           aaTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgc                                                            >  1:1608125/1‑71 (MQ=255)
           aaTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTg                                                             >  1:3692766/1‑70 (MQ=255)
            aTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATag                                                                 >  1:2891432/1‑65 (MQ=255)
            aTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgcg                                                           <  1:677872/71‑1 (MQ=255)
            aTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgcg                                                           >  1:3498591/1‑71 (MQ=255)
            aTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgcg                                                           <  1:1881534/71‑1 (MQ=255)
             tGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCAt                                                                   <  1:3256360/62‑1 (MQ=255)
              gcgcCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATaga                                                                >  1:141953/1‑64 (MQ=255)
               cgcCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgcggc                                                         <  1:3984823/70‑1 (MQ=255)
                  cAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCggcgag                                                      <  1:807771/70‑1 (MQ=255)
                  cAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCggcgag                                                      <  1:839371/70‑1 (MQ=255)
                   aaCGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCggcgaggg                                                    >  1:2760393/1‑71 (MQ=255)
                     cgacgCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCa                                                                    <  1:217846/53‑1 (MQ=255)
                      gacgCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATaga                                                                >  1:2307951/1‑56 (MQ=255)
                       acgCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCggcgagggcgag                                                >  1:3807854/1‑71 (MQ=255)
                         gCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGgcgagggcgagg                                               <  1:2633125/70‑1 (MQ=255)
                          ccTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGaa                                             <  1:1289382/71‑1 (MQ=255)
                           cTAACGGCACTACCAACATAACCGAGACAGGAATTGACCAGCTTTCATag                                                                 <  1:2395755/50‑1 (MQ=16)
                           cTAACGGAACAACCAACATCTCCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGaa                                             >  1:59456/1‑70 (MQ=255)
                           cTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctct                                                                      <  1:1328816/45‑1 (MQ=255)
                            tAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCa                                                                    <  1:1286912/46‑1 (MQ=255)
                            tAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGgcgagggcgagg                                               <  1:3060834/67‑1 (MQ=255)
                            tAACGGAACAACCAACATCACCGAGAACGGAATTGACAGCTCTCATagag                                                               >  1:3056635/1‑50 (MQ=38)
                            tAACGGAACAACCAACATCACCGACAACGGAATTGACCAGCTCTCATa                                                                  <  1:2978922/48‑1 (MQ=255)
                             aaCGGAACAACCAACATCACCGAGGACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAacac                                          >  1:3289650/1‑71 (MQ=255)
                             aaCGGAACAACCAACATCACCGAGAACGGAATTGACCAGctct                                                                      <  1:1780761/43‑1 (MQ=255)
                              aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATagag                                                               >  1:688725/1‑49 (MQ=255)
                              aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgcgg                                                          >  1:902489/1‑54 (MQ=255)
                              aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTg                                                             >  1:961289/1‑51 (MQ=255)
                              aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTg                                                             >  1:3219363/1‑51 (MQ=255)
                              aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCggcgagggcga                                                 >  1:331539/1‑63 (MQ=255)
                              aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCggcgagggcga                                                 >  1:1739379/1‑63 (MQ=255)
                              aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAacac                                          >  1:273260/1‑70 (MQ=255)
                              aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAacac                                          >  1:1266403/1‑70 (MQ=255)
                              aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAacac                                          >  1:3177347/1‑70 (MQ=255)
                                 gAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGt                                                              <  1:481258/47‑1 (MQ=255)
                                   acaacCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAacac                                          <  1:1918614/65‑1 (MQ=255)
                                    ctaccaacATAACCGAGACAGGAATTGACCAGCTTTCATa                                                                  >  1:3534615/3‑40 (MQ=16)
                                    ctaccaacATAACCGAGACAGGAATTGACCAGCTTTCATa                                                                  >  1:3486312/3‑40 (MQ=16)
                                    ctaccaacATAACCGAGACAGGAATTGACCAGCTTTCATa                                                                  >  1:3097720/3‑40 (MQ=16)
                                     taccaacATAACCGAGACAGGAATTGACCAGCTTTCATa                                                                  <  1:4083989/38‑1 (MQ=18)
                                     taccaacATAACCGAGACAGGAATTGACCAGCTTTCATAGAGTGCagcaag                                                      >  1:222464/2‑51 (MQ=19)
                                     aaccaacATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGgcgagggcgagg                                               <  1:3710034/58‑1 (MQ=255)
                                     aaccaacATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGa                                              <  1:596119/59‑1 (MQ=255)
                                       ccaacATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGt                                >  1:2438232/1‑71 (MQ=255)
                                        caacATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGaa                                             <  1:2762973/57‑1 (MQ=255)
                                         aacATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGa                              <  1:3306459/71‑1 (MQ=255)
                                         aacATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGa                              <  1:2710466/71‑1 (MQ=255)
                                         aacATAACCGAGACAGGAATTGACCAGCTTTCATaga                                                                >  1:1599954/1‑37 (MQ=21)
                                          acATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGaa                             <  1:1374061/71‑1 (MQ=255)
                                          acATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGa                              <  1:1214633/70‑1 (MQ=255)
                                          acATCACCGAGAACGGAATTGACCACCTCTCATaga                                                                >  1:314943/1‑36 (MQ=255)
                                           cATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGaaa                            >  1:3094205/1‑71 (MQ=255)
                                             tCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAAt                           >  1:2990247/1‑70 (MQ=255)
                                               aCCGAGAACGGAATTGACCAGCTCTCATAGAGTgc                                                            >  1:12477/1‑35 (MQ=255)
                                               aCCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGa                                     >  1:3416754/1‑58 (MQ=255)
                                               aCCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATc                          >  1:2364050/1‑69 (MQ=255)
                                                ccGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCa                       >  1:350571/1‑71 (MQ=255)
                                                ccGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCa                       >  1:3219231/1‑71 (MQ=255)
                                                ccGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCa                       >  1:69727/1‑71 (MQ=255)
                                                ccGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCa                       >  1:2438345/1‑71 (MQ=255)
                                                 cGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGa                              <  1:39606/63‑1 (MQ=255)
                                                    gaACGTAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAAcacca                                        >  1:1008266/1‑50 (MQ=255)
                                                      aCGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAAcacc                                         <  1:905322/47‑1 (MQ=255)
                                                       cGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATc                          <  1:457924/61‑1 (MQ=255)
                                                         gAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCa              <  1:2209047/71‑1 (MQ=255)
                                                          aaTTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAgcgc               <  1:2157949/69‑1 (MQ=255)
                                                          aaTTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCAg             >  1:3930261/1‑71 (MQ=255)
                                                            ttGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCg                         <  1:2549544/57‑1 (MQ=255)
                                                              gACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCAGGAGc         >  1:4049698/1‑71 (MQ=255)
                                                               aCCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCAGGAGc         <  1:3311913/70‑1 (MQ=255)
                                                                  aGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCAGGAGCCTGa     >  1:3133806/1‑71 (MQ=255)
                                                                   gCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCAGGAGCCTGAt    >  1:2401641/1‑71 (MQ=255)
                                                                   gCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCAGGAGCCTGAt    >  1:3050954/1‑71 (MQ=255)
                                                                   gCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCAGGAGCCTGAt    >  1:1753335/1‑71 (MQ=255)
                                                                     tctcATAGTGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCAGGAGCCTGAttt  <  1:2312535/71‑1 (MQ=255)
                                                                      |                                                                     
CGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAATCGCATACAGCGCAGGAGCCTGATTT  >  minE/2584534‑2584673

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: