Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,002 T→G 6.3% L425L (CTA→CTC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0020TG 6.3% 164.8 / 4.8 64L425L (CTA→CTCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base T (27/33);  new base G (2/2);  total (29/35)
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.81e-01

ACGGCAATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTA  >  minE/2585937‑2586072
                                                                 |                                                                      
acGGCAATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCtt                                                                   <  1:807007/71‑1 (MQ=255)
 cGGCAATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTa                                                                  >  1:2298244/1‑71 (MQ=255)
  ggCAATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTAt                                                                 <  1:1850914/71‑1 (MQ=255)
   gCAATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCtt                                                                   >  1:3621996/1‑68 (MQ=255)
   gCAATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATc                                                                >  1:2736605/1‑71 (MQ=255)
    cAATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCt                                                               <  1:562269/71‑1 (MQ=255)
    cAATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCt                                                               <  1:3597544/71‑1 (MQ=255)
       tACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTcc                                                            >  1:192069/1‑71 (MQ=255)
       caCCGACCAGCGCACGTTGGATCTGCCCCATCTATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTcc                                                            <  1:3480439/70‑1 (MQ=255)
          cGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCa                                                         <  1:3548253/71‑1 (MQ=255)
            aCCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCtt                                                                   >  1:2969976/1‑59 (MQ=255)
                 cgcACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGa                                                     >  1:2228531/1‑68 (MQ=255)
                     cGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGc                                                          >  1:726840/1‑59 (MQ=255)
                      gTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTc                                              >  1:3465481/1‑70 (MQ=255)
                      gTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTc                                              >  1:2087638/1‑70 (MQ=255)
                      gTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCg                                             <  1:2147322/71‑1 (MQ=255)
                        tGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGAc                                           >  1:3050303/1‑71 (MQ=255)
                        tGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGATCGa                                            >  1:737124/1‑70 (MQ=255)
                                ccccATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTtctc                                     <  1:1335659/69‑1 (MQ=255)
                                ccccATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTcca                                   <  1:1763130/71‑1 (MQ=255)
                                  ccATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGa                                                     >  1:1062587/1‑51 (MQ=255)
                                  ccATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccacc                                 <  1:3363532/71‑1 (MQ=255)
                                   cATCGATTTTTCCGTCGCTTCTTTTGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTccacca                                >  1:1872414/1‑71 (MQ=255)
                                    aTCGATTTTTCCGTCGCTTCTTTTGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTcc                                    <  1:3054131/66‑1 (MQ=255)
                                    aTCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                               >  1:3949271/1‑71 (MQ=255)
                                    aTCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                               >  1:1533686/1‑71 (MQ=255)
                                    aTCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                               >  1:20320/1‑71 (MQ=255)
                                     tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                               >  1:3540943/1‑70 (MQ=255)
                                     tCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccaccac                               >  1:1977280/1‑70 (MQ=255)
                                       gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGac                                                    >  1:3065355/1‑47 (MQ=255)
                                       gATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                            >  1:2219701/1‑71 (MQ=255)
                                        aTTTTTCCGTCGCTTCTTTTGGCGGGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTccaccacca                             >  1:1671793/1‑69 (MQ=255)
                                        aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                            >  1:3886668/1‑70 (MQ=255)
                                        aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACg                           <  1:2485590/71‑1 (MQ=255)
                                        aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACg                           <  1:559040/71‑1 (MQ=255)
                                        aTTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACg                           <  1:3494392/71‑1 (MQ=255)
                                         tttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCg                                             <  1:437012/52‑1 (MQ=255)
                                         tttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                          <  1:2510934/71‑1 (MQ=255)
                                         tttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                          <  1:2038118/71‑1 (MQ=255)
                                          ttATGCGTCGCTTCCTTCGGTGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGa                          >  1:225237/1‑70 (MQ=255)
                                          ttATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGAt                         >  1:1175403/1‑71 (MQ=255)
                                           tATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGAt                         <  1:407216/70‑1 (MQ=255)
                                            aTGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATg                        <  1:3945938/70‑1 (MQ=255)
                                            aTGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATg                        <  1:1426324/70‑1 (MQ=255)
                                             tGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATgg                       <  1:274622/70‑1 (MQ=255)
                                             tGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGc                      <  1:2282291/71‑1 (MQ=255)
                                              gCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCa                                                         <  1:2126999/35‑1 (MQ=255)
                                                 tCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACg                                          <  1:689279/47‑1 (MQ=255)
                                                 tCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTccac                                  <  1:184344/55‑1 (MQ=255)
                                                    cttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcc               <  1:2360006/71‑1 (MQ=255)
                                                    cttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcc               <  1:576957/71‑1 (MQ=255)
                                                      tccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac             >  1:1285240/1‑71 (MQ=255)
                                                      tccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac             >  1:625026/1‑71 (MQ=255)
                                                       ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTTATccac             <  1:511122/70‑1 (MQ=255)
                                                       ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac             <  1:2412808/70‑1 (MQ=255)
                                                          ttggcggGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAg          <  1:2026811/68‑1 (MQ=12)
                                                          tcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCcaccaccac                            <  1:1894718/52‑1 (MQ=255)
                                                           cggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCGCCAGt         >  1:1686169/1‑70 (MQ=255)
                                                            ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa       >  1:3302608/1‑71 (MQ=255)
                                                            ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa        <  1:1024032/70‑1 (MQ=255)
                                                            ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa        <  1:1783259/70‑1 (MQ=255)
                                                            ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa        <  1:3740850/70‑1 (MQ=255)
                                                             gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc      >  1:3736467/1‑71 (MQ=255)
                                                               ggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca           <  1:1971060/64‑1 (MQ=255)
                                                                 tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa  <  1:290303/71‑1 (MQ=255)
                                                                 tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa  <  1:3595756/71‑1 (MQ=255)
                                                                 |                                                                      
ACGGCAATACCGACCAGCGCACGTTGGATCTGCCCCATCGATTTATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTA  >  minE/2585937‑2586072

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: