Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,592,859 T→A 6.0% N50I (AAT→ATT)  yhiF ← predicted DNA‑binding ranscriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,592,8590TA 6.0% 231.2 / 5.8 84N50I (AAT→ATT) yhiFpredicted DNA‑binding ranscriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (35/43);  new base A (2/3);  total (38/46)
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.76e-01

GGAGAAAATAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTGTAT  >  minE/2592789‑2592929
                                                                      |                                                                      
ggAGAAAATAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACattat                                                                        <  1:779615/71‑1 (MQ=255)
ggAGAAAATAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACattat                                                                        <  1:3493711/71‑1 (MQ=255)
  agaAAATAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACattatt                                                                       >  1:2014740/1‑70 (MQ=255)
   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:2445290/1‑49 (MQ=25)
                                  tCAATTTGAGTGAGAAAATTAGAATGAAATACATTAATGAATAATGTAtt                                                           >  1:424483/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:2797535/50‑1 (MQ=255)
                                                aaaaTTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAgg                                             <  1:3738218/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:2263095/70‑1 (MQ=255)
                                                            aaaTACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttct             <  1:2930656/70‑1 (MQ=255)
                                                            aaaTACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttct             >  1:2357981/1‑70 (MQ=255)
                                                            aaaTACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCtcttct             <  1:27849/70‑1 (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)
                                                                      |                                                                      
GGAGAAAATAATGACATGGACAGGTTTTAATCGTTCAATTTGAGTGAGAAAATTAGAATGAAATACATTATTCATTAATGTATCAATAATGACGAAGGCATTCTGATGCAACATTACGTCGATCTCTTCTTCGTTCTGTAT  >  minE/2592789‑2592929

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: