Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,587,282 T→C 10.9% intergenic (‑6/+19) mdtF ← / ← mdtE multidrug transporter, RpoS‑dependent/multidrug resistance efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,587,2820TC10.9% 119.8 / 6.3 46intergenic (‑6/+19)mdtF/mdtEmultidrug transporter, RpoS‑dependent/multidrug resistance efflux transporter
Reads supporting (aligned to +/- strand):  ref base T (21/20);  new base C (3/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCCTGGCTGGAGGAAATTGCTCG  >  minE/2587212‑2587351
                                                                      |                                                                     
gCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGt                                                                       >  1:1146847/1‑71 (MQ=255)
      atcatAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCtt                                                                 <  1:599994/71‑1 (MQ=255)
         ataataATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAg                                                              <  1:1925123/71‑1 (MQ=255)
           aataatGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGcc                                                            >  1:2504207/1‑71 (MQ=255)
            ataatGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCt                                                           >  1:1863013/1‑71 (MQ=255)
               atGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGc                                                         <  1:2116179/70‑1 (MQ=255)
                tGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCaa                                                       <  1:3002966/71‑1 (MQ=255)
                    aaGTACCCAGGCAAAAACCGGGCCATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGt                                                    >  1:718176/1‑70 (MQ=255)
                         cccAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGt                                              <  1:684563/71‑1 (MQ=255)
                          ccAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGt                                              <  1:75064/70‑1 (MQ=255)
                           cAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGttt                                            >  1:3605232/1‑71 (MQ=255)
                              gCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGa                                         >  1:3018641/1‑71 (MQ=255)
                                aaaaaCCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTc                                           <  1:1840941/67‑1 (MQ=255)
                                 aaaaCCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGttt                                            <  1:2577515/65‑1 (MQ=255)
                                 aaaaCCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCg                                          <  1:1441160/67‑1 (MQ=255)
                                  aaaCCGGGCGATCAATAAAATAGTTAGCCATGAAAGTccc                                                                    >  1:608129/1‑40 (MQ=255)
                                  aaaCCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGc                                                             <  1:3449815/47‑1 (MQ=255)
                                  aaaCCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCg                                     <  1:3626268/71‑1 (MQ=255)
                                   aaCCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCgg                                    >  1:757200/1‑71 (MQ=255)
                                     ccGGGCGATCAATAAAATAGTTAGCCATGAAAGCCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGt                                   >  1:3677879/1‑70 (MQ=255)
                                       gggCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAAc                                                      >  1:3254242/1‑49 (MQ=255)
                                       gggCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCt                                >  1:3626221/1‑71 (MQ=255)
                                       gggCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCt                                >  1:2721640/1‑71 (MQ=255)
                                        ggCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTg                               >  1:45585/1‑71 (MQ=255)
                                         gCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCg                                     >  1:1684710/1‑64 (MQ=255)
                                         gCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTgg                              <  1:55577/71‑1 (MQ=255)
                                             tCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGAtt                                       <  1:2492860/58‑1 (MQ=255)
                                             tCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGtt                          >  1:3419888/1‑71 (MQ=255)
                                               aaTAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGtttt                        <  1:314787/71‑1 (MQ=255)
                                               aaTAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGttt                         <  1:1496404/70‑1 (MQ=255)
                                                  aaaaTAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGc                                 >  1:200330/1‑59 (MQ=255)
                                                  aaaaTAGTTAGCCATGAAAGCCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCct                     >  1:3443980/1‑71 (MQ=255)
                                                   aaaTAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCctg                    >  1:1831889/1‑71 (MQ=255)
                                                   aaaTAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCct                     <  1:2585674/70‑1 (MQ=255)
                                                    aaTAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTa                                                  >  1:829744/1‑40 (MQ=255)
                                                    aaTAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTa                                                  >  1:1714123/1‑40 (MQ=255)
                                                      tAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCctggct                 <  1:897587/71‑1 (MQ=255)
                                                       aGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTc                       >  1:808967/1‑64 (MQ=255)
                                                               aTGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCCTGGCTGGAGGAAAt        >  1:1886671/1‑71 (MQ=255)
                                                               aTGAAAGCCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCCTGGCTGGAGGAAAt        >  1:2136880/1‑71 (MQ=255)
                                                                tGAAAGCCCCCTTAAGCCTGCAACGTTATTGTTTCGATTc                                      <  1:1459178/40‑1 (MQ=255)
                                                                tGAAAGCCCCCTTAAGCCTGCAACGTTATTGTTTCGATTc                                      <  1:3291854/40‑1 (MQ=255)
                                                                  aaaGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTcc                      >  1:2178674/1‑54 (MQ=255)
                                                                  aaaGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCCTGGCTGGAGGAAATTGc     >  1:3850437/1‑71 (MQ=255)
                                                                  aaaGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCCTGGCTGGAGGAAATTGc     <  1:3350851/71‑1 (MQ=255)
                                                                     gTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCCTGGCTGGAGGaa          <  1:3214516/63‑1 (MQ=255)
                                                                     gTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCCTGGCTGGAGGAAATTGCTCg  <  1:2140962/71‑1 (MQ=255)
                                                                      tCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCCTGGCTggagga           <  1:2299016/61‑1 (MQ=255)
                                                                      |                                                                     
GCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATAGTTAGCCATGAAAGTCCCCTTAAGCCTGCAACGTTATTGTTTCGATTCGGTGCTGGCGTTTTCCTGGCTGGAGGAAATTGCTCG  >  minE/2587212‑2587351

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: