Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,584,325 C→T 13.6% V984V (GTG→GTA mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,584,3250CT13.6% 248.6 / 22.7 88V984V (GTG→GTAmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base C (38/38);  new base T (5/7);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

ATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATT  >  minE/2584260‑2584394
                                                                 |                                                                     
aTCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                                                  <  1:606367/71‑1 (MQ=255)
   aCGCCGTTACGTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                                                  >  1:8647/1‑68 (MQ=255)
   aCGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                                                  >  1:1972868/1‑68 (MQ=255)
   aCGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                                                  >  1:761617/1‑68 (MQ=255)
   aCGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                                                  >  1:1453642/1‑68 (MQ=255)
      ccGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAAt                                                            <  1:2638641/71‑1 (MQ=255)
       cGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAAt                                                            <  1:389211/70‑1 (MQ=255)
        gatacccacaGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATaa                                                          <  1:3682627/69‑1 (MQ=9)
         gTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                                           >  1:1324176/1‑69 (MQ=255)
         gTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                                           >  1:2777140/1‑69 (MQ=255)
           acctaccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                                                          >  1:1176241/1‑68 (MQ=255)
           acctaccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                                                          >  1:548358/1‑68 (MQ=255)
             ctaccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                                           <  1:2550853/65‑1 (MQ=255)
             ctaccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGc                                                      >  1:3605113/1‑70 (MQ=255)
             ctaccGCGTTTTGCCCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGcc                                                     <  1:3879304/71‑1 (MQ=255)
               acagcgTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGag                                                   <  1:3860792/68‑1 (MQ=9)
                ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGgcagca                                                                      >  1:72767/1‑51 (MQ=255)
                ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCag                                                   <  1:78835/70‑1 (MQ=255)
                 cgcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagag                                                 <  1:2617272/71‑1 (MQ=255)
                 cgcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                                                    >  1:825207/1‑68 (MQ=255)
                 cgcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                                                    >  1:1033604/1‑68 (MQ=255)
                    gTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCAc                                                                     >  1:458692/1‑48 (MQ=255)
                    gTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAgg                                              >  1:1638040/1‑71 (MQ=255)
                    gTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAgg                                              >  1:798087/1‑71 (MQ=255)
                       ttGCGCGCCAGAACCGGCACCATGCCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGCGGtca                                           >  1:683703/1‑71 (MQ=255)
                       ttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAgg                                                       >  1:2058887/1‑59 (MQ=255)
                       ttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAgg                                                       >  1:2755116/1‑59 (MQ=255)
                       ttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                  >  1:3515077/1‑64 (MQ=255)
                       ctGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtca                                           >  1:3009935/2‑71 (MQ=255)
                        tGCGCGCCAGAACCGGCACCATGAGTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtca                                           <  1:449981/70‑1 (MQ=255)
                         gcgcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtcat                                          <  1:1818497/70‑1 (MQ=255)
                         gcgcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtcat                                          <  1:2212633/70‑1 (MQ=255)
                          cgcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagagag                                               >  1:2061067/1‑64 (MQ=255)
                            cgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                                      <  1:1064270/71‑1 (MQ=255)
                            cgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAg                                       >  1:778072/1‑70 (MQ=255)
                            cgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAg                                       <  1:3726042/70‑1 (MQ=255)
                             gcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagag                                                 <  1:995071/59‑1 (MQ=255)
                             gcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagag                                                 <  1:3366805/59‑1 (MQ=255)
                             gcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                                     <  1:2138555/71‑1 (MQ=255)
                              ccAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAgg                                              <  1:1522864/61‑1 (MQ=255)
                              ccAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt                                    <  1:2543743/71‑1 (MQ=255)
                               cAGAACCGGCACCATGACTGATAACCAGCGgcagca                                                                      <  1:2425462/36‑1 (MQ=255)
                               cAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                                                         <  1:1447311/49‑1 (MQ=255)
                                aGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                                                          >  1:2330768/1‑47 (MQ=255)
                                aGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTg                                  >  1:1199883/1‑71 (MQ=255)
                                aGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTg                                  >  1:2286118/1‑71 (MQ=255)
                                 gAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTgg                                 >  1:2328731/1‑71 (MQ=255)
                                 gAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTgg                                 >  1:806792/1‑71 (MQ=255)
                                 gAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTgg                                 >  1:3322800/1‑71 (MQ=255)
                                  aaCCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTg                                  <  1:1525240/69‑1 (MQ=255)
                                   gccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCgagaga                                                >  1:880052/2‑54 (MQ=255)
                                   gccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGAta                                   >  1:243958/2‑66 (MQ=12)
                                   aCCGGCACCATGACTGATAACCAGCGGCAGCACGCCga                                                                >  1:2378583/1‑38 (MQ=255)
                                     cGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGt                             <  1:1801768/71‑1 (MQ=255)
                                     cGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGt                             <  1:3420762/71‑1 (MQ=255)
                                      ggCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACgca                            <  1:3777631/71‑3 (MQ=9)
                                      ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa                            <  1:3368213/71‑1 (MQ=255)
                                       gCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACgca                            >  1:1112377/1‑68 (MQ=9)
                                       gCACCATGACTGATAACCAGCGGCAGCACGCCgaga                                                              >  1:1350824/1‑36 (MQ=255)
                                       gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAtt                                   <  1:2733812/63‑1 (MQ=255)
                                       gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaa                           <  1:1463295/71‑1 (MQ=255)
                                       gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaa                           <  1:1981888/71‑1 (MQ=255)
                                       gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaa                           <  1:76557/71‑1 (MQ=255)
                                       gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa                            >  1:1854631/1‑70 (MQ=255)
                                       gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa                            >  1:1447028/1‑70 (MQ=255)
                                         aCCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAAc                         >  1:3143752/1‑71 (MQ=255)
                                         aCCATGACTGATAACCAGCGGCAGCAAGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                          >  1:478243/1‑70 (MQ=255)
                                           cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                       <  1:3305821/71‑1 (MQ=255)
                                            aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                      >  1:1044603/1‑71 (MQ=255)
                                             ttaCTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGgacgcagacgca                      >  1:318409/3‑70 (MQ=2)
                                               aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATcc                   <  1:2357308/71‑1 (MQ=255)
                                               aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATcc                   <  1:2894015/71‑1 (MQ=255)
                                                cTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAggac                               <  1:263185/58‑3 (MQ=9)
                                                cTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAggac                               <  1:3151316/58‑3 (MQ=9)
                                                cTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAggac                               <  1:3012318/58‑3 (MQ=9)
                                                cTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACGCAtac                   <  1:627934/70‑3 (MQ=9)
                                                cTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATcc                   <  1:35348/70‑1 (MQ=255)
                                                 tGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTgg                                 <  1:3715236/55‑1 (MQ=255)
                                                  gATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAtt                                   <  1:2947039/52‑1 (MQ=255)
                                                   aTAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                          >  1:2626189/1‑60 (MQ=255)
                                                   aTAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCAt                     >  1:1428238/1‑65 (MQ=255)
                                                   aTAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCcgcg                >  1:2795092/1‑70 (MQ=255)
                                                    tAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCgcgc               <  1:634815/70‑1 (MQ=255)
                                                    tAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgcc              <  1:264338/71‑1 (MQ=255)
                                                      aCCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgccgc            <  1:3082112/71‑1 (MQ=255)
                                                        cAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCtt          <  1:2234995/71‑1 (MQ=255)
                                                         aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTc         >  1:46695/1‑71 (MQ=255)
                                                          gCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTc         >  1:2027311/1‑70 (MQ=255)
                                                          gCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCg        <  1:3706929/71‑1 (MQ=255)
                                                           cGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCg        <  1:1936266/70‑1 (MQ=255)
                                                             gcagcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                                     >  1:2192197/1‑39 (MQ=255)
                                                              cagcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAAc                         <  1:320007/50‑1 (MQ=255)
                                                                gcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGAtt  <  1:3651375/71‑1 (MQ=255)
                                                                 |                                                                     
ATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATT  >  minE/2584260‑2584394

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: