Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,584,580 C→A 6.9% V899V (GTG→GTT mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,584,5800CA 6.9% 335.8 / 5.7 101V899V (GTG→GTTmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base C (49/45);  new base A (5/2);  total (54/47)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 102 total aligned reads displayed.

GTCATTACTTAAGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAAT  >  minE/2584511‑2584648
                                                                     |                                                                    
gTCATTACTTAAGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATc                                                                      <  1:1724757/70‑1 (MQ=255)
   aTTACTTAAGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCa                                                                     >  1:278888/1‑68 (MQ=255)
   aTTACTTAAGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCAcc                                                                   >  1:1201786/1‑70 (MQ=255)
   aTTACTTAAGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCAcc                                                                   >  1:915757/1‑70 (MQ=255)
   aTTACTTAAGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCAcc                                                                   >  1:4083558/1‑70 (MQ=255)
      aCTTAAGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCgag                                                                >  1:2465948/1‑70 (MQ=255)
      aCTTAAGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCgag                                                                >  1:4107570/1‑70 (MQ=255)
      aCTTAAGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCgag                                                                >  1:2626370/1‑70 (MQ=255)
          aaGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACg                                                            >  1:3816358/1‑70 (MQ=255)
               gcgcAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTg                                                      <  1:2309241/71‑1 (MQ=255)
               gcgcAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTg                                                      <  1:3114887/71‑1 (MQ=255)
                cgcAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGa                                                     <  1:3031139/71‑1 (MQ=255)
                 gcAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGAc                                                    >  1:746796/1‑71 (MQ=255)
                 gcAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGAc                                                    <  1:3844953/71‑1 (MQ=255)
                 gcAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGAc                                                    >  1:49976/1‑71 (MQ=255)
                    gATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGa                                                     <  1:2522270/67‑1 (MQ=255)
                     aTCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGc                                                <  1:2987539/71‑1 (MQ=255)
                      tCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCgaga                                                               >  1:185445/1‑55 (MQ=255)
                       cGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctc                                              >  1:3018857/1‑71 (MQ=255)
                       cGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctc                                              >  1:1970207/1‑71 (MQ=255)
                        ggTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGaa                                                         >  1:1437207/1‑59 (MQ=255)
                        ggTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctct                                             <  1:1278301/71‑1 (MQ=255)
                        ggTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctct                                             <  1:1498835/71‑1 (MQ=255)
                        ggTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctc                                              <  1:2496733/70‑1 (MQ=255)
                        ggTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctc                                              <  1:3950942/70‑1 (MQ=255)
                         gTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGctct                                             <  1:2787540/70‑1 (MQ=255)
                         gTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCtctc                                            <  1:2232431/71‑1 (MQ=255)
                          tGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCa                                           <  1:2627101/71‑1 (MQ=255)
                          tGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCa                                           <  1:4005309/71‑1 (MQ=255)
                            gCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCAc                                                                    >  1:1149581/1‑44 (MQ=255)
                            gCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCga                                                                 >  1:46931/1‑47 (MQ=255)
                             ccAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATag                                        <  1:1324483/71‑1 (MQ=255)
                             ccAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATag                                        <  1:4033511/71‑1 (MQ=255)
                              cAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATa                                         <  1:1739454/69‑1 (MQ=255)
                               aGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGATTTGACCAGCTCTCATa                                         <  1:1353589/68‑1 (MQ=255)
                                 tAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCa                                           >  1:3209795/1‑64 (MQ=255)
                                 tAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCa                                           >  1:2635989/1‑64 (MQ=255)
                                 tAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGt                                     >  1:3538628/1‑70 (MQ=255)
                                  aaTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGaa                                                         >  1:1815115/1‑49 (MQ=255)
                                  aaTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgc                                   >  1:1608125/1‑71 (MQ=255)
                                  aaTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTg                                    >  1:3692766/1‑70 (MQ=255)
                                   aTGCGCCAACGACGCCTAACGGAACAACCAACATCAcc                                                                   >  1:2579359/1‑38 (MQ=255)
                                   aTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATag                                        >  1:2891432/1‑65 (MQ=255)
                                   aTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgcg                                  <  1:1881534/71‑1 (MQ=255)
                                   aTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgcg                                  <  1:677872/71‑1 (MQ=255)
                                   aTGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgcg                                  >  1:3498591/1‑71 (MQ=255)
                                    tGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCAt                                          <  1:3256360/62‑1 (MQ=255)
                                     gcgcCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATaga                                       >  1:141953/1‑64 (MQ=255)
                                      cgcCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgcggc                                <  1:3984823/70‑1 (MQ=255)
                                        ccAACGACGCCTAACGGAACAACCAACATCACCgag                                                                <  1:2688693/36‑1 (MQ=255)
                                        ccAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGAc                                                    <  1:2646792/48‑1 (MQ=255)
                                        ccAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGct                                               >  1:2485047/1‑53 (MQ=255)
                                         cAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCggcgag                             <  1:839371/70‑1 (MQ=255)
                                         cAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCggcgag                             <  1:807771/70‑1 (MQ=255)
                                          aaCGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCggcgaggg                           >  1:2760393/1‑71 (MQ=255)
                                           acgacgCCTAACGGAACAACCAACATCACCGAGAACgg                                                           <  1:2787243/38‑1 (MQ=255)
                                           acgacgCCTAACGGAACAACCAACATCACCGAGAACgg                                                           <  1:408824/38‑1 (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)
                                                     aCGGAACAACCAACATCACCGAGAACGGAATTGAc                                                    >  1:3373599/1‑35 (MQ=255)
                                                     aCGGAACAACCAACATCACCGAGAACGGAATTGAc                                                    >  1:3378448/1‑35 (MQ=255)
                                                     aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATagag                                      >  1:688725/1‑49 (MQ=255)
                                                     aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTgcgg                                 >  1:902489/1‑54 (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:1266403/1‑70 (MQ=255)
                                                     aCGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAacac                 >  1:273260/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:3097720/3‑40 (MQ=16)
                                                           ctaccaacATAACCGAGACAGGAATTGACCAGCTTTCATa                                         >  1:3534615/3‑40 (MQ=16)
                                                           ctaccaacATAACCGAGACAGGAATTGACCAGCTTTCATa                                         >  1:3486312/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:2710466/71‑1 (MQ=255)
                                                                aacATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGa     <  1:3306459/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)
                                                                     |                                                                    
GTCATTACTTAAGCCGCGCAGATCGGTGGCCAGTAATGCGCCAACGACGCCTAACGGAACAACCAACATCACCGAGAACGGAATTGACCAGCTCTCATAGAGTGCGGCGAGGGCGAGGAACACCACGACCAGTGAAAT  >  minE/2584511‑2584648

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: