Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,592,861 C→G 5.1% M49I (ATG→ATC) ‡ yhiF ← predicted DNA‑binding ranscriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,592,8610CG 5.1% 301.9 / 2.9 96M49I (ATG→ATC) ‡yhiFpredicted DNA‑binding ranscriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (42/49);  new base G (2/3);  total (44/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAATAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTGTATAT  >  minE/2592792‑2592931
                                                                     |                                                                      
gaAAATAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCa                                                                       >  1:43629/1‑71 (MQ=255)
gaAAATAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCa                                                                       >  1:3216908/1‑71 (MQ=255)
gaAAATAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCa                                                                       >  1:3324637/1‑71 (MQ=255)
  aaaTAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCAt                                                                      >  1:3761343/1‑70 (MQ=255)
  aaaTAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCAt                                                                      >  1:303299/1‑70 (MQ=255)
  aaaTAATGACATGAACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCAt                                                                      >  1:1815246/1‑70 (MQ=255)
        tGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTa                                                               <  1:3236935/71‑1 (MQ=255)
        tGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTa                                                               <  1:807356/71‑1 (MQ=255)
          aCATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTAt                                                              >  1:3611806/1‑70 (MQ=255)
            aTGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCaa                                                           >  1:28275/1‑71 (MQ=255)
            aTGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCaa                                                           >  1:2289690/1‑71 (MQ=255)
              ggACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCaat                                                          >  1:3344293/1‑70 (MQ=255)
               gACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCaataa                                                        >  1:1926364/1‑71 (MQ=255)
               gACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCaata                                                         <  1:1608663/70‑1 (MQ=255)
               gACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCaata                                                         <  1:2262374/70‑1 (MQ=255)
                aCAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCaataat                                                       >  1:1278632/1‑71 (MQ=255)
                aCAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCaataat                                                       <  1:38171/71‑1 (MQ=255)
                 cAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATg                                                      <  1:1672743/71‑1 (MQ=255)
                    gTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACg                                                   <  1:2557239/71‑1 (MQ=255)
                    gTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACg                                                   <  1:1647664/71‑1 (MQ=255)
                     ttttAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGa                                                  <  1:806003/71‑1 (MQ=255)
                     ttttAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGa                                                  <  1:2344513/71‑1 (MQ=255)
                      tttAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTGATGTATCAATAATGACg                                                   >  1:2738810/1‑69 (MQ=255)
                      tttAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACg                                                   >  1:806147/1‑69 (MQ=255)
                      tttAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACg                                                   >  1:1979347/1‑69 (MQ=255)
                      tttAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGa                                                  <  1:673487/70‑1 (MQ=255)
                        tAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAAt                                                                  >  1:215471/1‑52 (MQ=255)
                           tCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCaat                                                          <  1:2337212/57‑1 (MQ=255)
                           tCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCa                                             <  1:2320471/70‑1 (MQ=255)
                              ttCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATg                                                      <  1:3686424/58‑1 (MQ=255)
                               tCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCt                                         >  1:117183/1‑70 (MQ=255)
                               tCAATTTGAGTGAGAAAATTAGAATGAAATACATTAATGAATAATGTAtt                                                             >  1:424483/1‑49 (MQ=25)
                               tCAATTTGAGTGAGAAAATTAGAATGAAATACATTAATGAATAATGTAtt                                                             >  1:2445290/1‑49 (MQ=25)
                                cAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGa                                       >  1:1848234/1‑71 (MQ=255)
                                 aaTTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGa                                       <  1:2599165/70‑1 (MQ=255)
                                 aaTTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGAt                                      >  1:2857707/1‑71 (MQ=255)
                                  aTTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTa                                                               <  1:2002778/45‑1 (MQ=255)
                                  aTTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATc                                                             <  1:514463/47‑1 (MQ=255)
                                     tgagtgagAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCaa                                  >  1:262060/1‑71 (MQ=255)
                                     tgagtgagAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCaa                                  >  1:1615748/1‑71 (MQ=255)
                                      gagtgagAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCaa                                  <  1:2563118/70‑1 (MQ=255)
                                       agtgagAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACa                                <  1:470235/71‑1 (MQ=255)
                                        gtgagAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACAt                               <  1:264206/71‑1 (MQ=255)
                                         tgagAAAATTAGAATGAAATACATTATTCATTAATGTATCaat                                                          <  1:2032589/43‑1 (MQ=255)
                                          gagaAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTc                                          <  1:3951163/58‑1 (MQ=255)
                                           agaAAATTAGAATGATATACATTATTCATTAATGTATCtataa                                                        <  1:2653233/43‑1 (MQ=255)
                                           agaAAATTAGAATGAAATAGATTATTCATTAATGTATCAATAATGACGAAGGCAt                                            <  1:3096229/55‑1 (MQ=255)
                                           agaAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCAt                                            <  1:2485101/55‑1 (MQ=255)
                                           agaAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTa                             >  1:3281355/1‑70 (MQ=255)
                                           agaAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTAc                            >  1:1588316/1‑71 (MQ=255)
                                           agaAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTCATg                                     >  1:2983063/1‑62 (MQ=255)
                                            gaAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACg                           <  1:1038073/71‑1 (MQ=255)
                                             aaaaTTAGAATGAAATACATTATTCATTGATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGt                          <  1:668691/71‑1 (MQ=255)
                                             aaaaTTAGAATGAAATACATTATTCATTAATGTATCAATAATGACg                                                   <  1:154280/46‑1 (MQ=255)
                                             aaaaTTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAgg                                               <  1:3738218/50‑1 (MQ=255)
                                             aaaaTTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAgg                                               <  1:2797535/50‑1 (MQ=255)
                                             aaaaTTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAgg                                               <  1:3455950/50‑1 (MQ=255)
                                              aaaTTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTc                         <  1:1660064/71‑1 (MQ=255)
                                                aTTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAgg                                               <  1:38115/47‑1 (MQ=255)
                                                aTTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCg                        >  1:3191171/1‑70 (MQ=255)
                                                aTTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGa                       >  1:517765/1‑71 (MQ=255)
                                                  tAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACg                           >  1:1449385/1‑65 (MQ=255)
                                                  tAGAATGAAATACATTAGTCATTAATGTATCaataat                                                       >  1:19590/1‑37 (MQ=255)
                                                     aaTGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATctct                  >  1:1201528/1‑71 (MQ=255)
                                                      aTGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATctct                  <  1:189994/70‑1 (MQ=255)
                                                       tGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttc                >  1:3774142/1‑71 (MQ=255)
                                                       tGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttc                >  1:3796929/1‑71 (MQ=255)
                                                         aaaTACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttct               <  1:2930656/70‑1 (MQ=255)
                                                         aaaTACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttct               <  1:2263095/70‑1 (MQ=255)
                                                         aaaTACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttct               <  1:27849/70‑1 (MQ=255)
                                                         aaaTACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttct               >  1:2357981/1‑70 (MQ=255)
                                                          gaTACATTAATGAATAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTAc                            <  1:42363/55‑1 (MQ=37)
                                                          gaTACATTAATGAATAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGAtc                     <  1:1835336/62‑1 (MQ=37)
                                                          gaTACATTAATGAATAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGAtc                     <  1:128008/62‑1 (MQ=37)
                                                          aaTACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCttcttc             >  1:2046257/1‑71 (MQ=255)
                                                           aTACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCttcttc             >  1:2136852/1‑70 (MQ=255)
                                                           aTACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCg            >  1:2415162/1‑71 (MQ=255)
                                                            tACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCg                        <  1:2289758/58‑1 (MQ=255)
                                                            tACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttct               <  1:1163663/67‑1 (MQ=255)
                                                            tACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGt           <  1:1396905/71‑1 (MQ=255)
                                                            tACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGt           <  1:2416219/71‑1 (MQ=255)
                                                            tACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGt           >  1:2208794/1‑71 (MQ=255)
                                                             aCATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCg                        <  1:1518898/57‑1 (MQ=255)
                                                              cATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTc         <  1:1766024/71‑1 (MQ=255)
                                                               attattCATTAATGTATCAATAATGACGAAGGCATTCTGa                                       >  1:3543313/1‑40 (MQ=255)
                                                                ttattCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTg       <  1:2616465/71‑1 (MQ=255)
                                                                ttattCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTg       >  1:1679885/1‑71 (MQ=255)
                                                                 tattCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttct               >  1:2955944/1‑62 (MQ=255)
                                                                 tattCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTg       >  1:63775/1‑70 (MQ=255)
                                                                  attCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTc                         >  1:925284/1‑51 (MQ=255)
                                                                  attCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCttcttc             <  1:2205570/63‑1 (MQ=255)
                                                                  attCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCttcttc             <  1:3510403/63‑1 (MQ=255)
                                                                  attCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTGt      <  1:2317219/70‑1 (MQ=255)
                                                                  attCATTAATGTATCAATAATGAAGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTGta     >  1:2249072/1‑71 (MQ=255)
                                                                   ttCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTGtat    <  1:41333/71‑1 (MQ=255)
                                                                   ttCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTGtat    <  1:1338920/71‑1 (MQ=255)
                                                                     cATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTGTatat  >  1:2117104/1‑71 (MQ=255)
                                                                     |                                                                      
GAAAATAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTGTATAT  >  minE/2592792‑2592931

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: