Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,747,876 G→A 28.6% S44L (TCA→TTA)  yfiE ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,747,8760GA28.6% 77.5 / 34.9 49S44L (TCA→TTA) yfiEpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (8/27);  new base A (6/8);  total (14/35)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTA  >  minE/1747809‑1747926
                                                                   |                                                  
cgtcagcaACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGaa                                                  <  1:4029465/67‑1 (MQ=255)
   cagcaACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGaa                                                  <  1:385568/67‑1 (MQ=255)
   cagcaACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACt                                              <  1:1074357/71‑1 (MQ=255)
          ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGaa                                                  <  1:880906/60‑1 (MQ=255)
          ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGaa                                                  <  1:546619/60‑1 (MQ=255)
          ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGaa                                                  <  1:3810701/60‑1 (MQ=255)
          ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGaa                                                  <  1:3102781/60‑1 (MQ=255)
                         gTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGctgctg                              <  1:3084331/65‑1 (MQ=255)
                         gTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGctgctg                              <  1:580049/65‑1 (MQ=255)
                            aGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGaaa                     <  1:2808832/71‑1 (MQ=255)
                            aGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGaa                      <  1:231899/70‑1 (MQ=255)
                            aGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAAATCCTGCTCAAGCTGCTGAATATGaaa                     <  1:1942029/71‑1 (MQ=255)
                                acaTCCGTCGGCCAATTTTCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:1882727/71‑1 (MQ=255)
                                acaTCCGTCGGCCAATTTTCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:2777339/71‑1 (MQ=255)
                                acaTCCGTCGGCCAATTTTCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:817752/71‑1 (MQ=255)
                                  aTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:3108976/69‑1 (MQ=255)
                                  aTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:1289111/69‑1 (MQ=255)
                                  aTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:1965817/69‑1 (MQ=255)
                                  aTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:3970946/69‑1 (MQ=255)
                                         ggCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGa         <  1:3955729/70‑1 (MQ=255)
                                         ggCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGa         <  1:1855066/70‑1 (MQ=255)
                                         ggCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGa         <  1:1416738/70‑1 (MQ=255)
                                         ggCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGa         <  1:3239390/70‑1 (MQ=255)
                                         ggCCAATTTACTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGa         <  1:952236/70‑1 (MQ=255)
                                          gCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGa         <  1:3709257/69‑1 (MQ=255)
                                           ccAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATtg      <  1:2457986/71‑1 (MQ=255)
                                            cAATTTTCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATg                        <  1:2650748/52‑1 (MQ=255)
                                            cAATTTTCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATg                        <  1:1487964/52‑1 (MQ=255)
                                              aTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATtgtg    >  1:430156/1‑70 (MQ=255)
                                              aTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATtgtg    >  1:1623502/1‑70 (MQ=255)
                                              aTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATtgtg    >  1:3108859/1‑70 (MQ=255)
                                              aTTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATtgtg    >  1:624473/1‑70 (MQ=255)
                                                tttCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTa  >  1:1034639/1‑70 (MQ=255)
                                                tttCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTa  >  1:1005428/1‑70 (MQ=255)
                                                 ttCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACTGTCGAt        >  1:2262723/1‑63 (MQ=255)
                                                 ttCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGAt        >  1:3890850/1‑63 (MQ=255)
                                                 ttCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGAt        >  1:2913126/1‑63 (MQ=255)
                                                 ttCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGAt        >  1:207916/1‑63 (MQ=255)
                                                 ttCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGAt        >  1:1567877/1‑63 (MQ=255)
                                                 ttCTCAAATAACTGGACTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGAt        >  1:1082828/1‑63 (MQ=255)
                                                         tAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCAc               >  1:304682/1‑48 (MQ=255)
                                                         tAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCAc               >  1:1212474/1‑48 (MQ=255)
                                                                aCTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATtgtg    <  1:3454320/52‑1 (MQ=255)
                                                                aCTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATtgtg    <  1:3010383/52‑1 (MQ=255)
                                                                aCTAAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATtgtg    <  1:1151017/52‑1 (MQ=255)
                                                                  tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:3972837/37‑1 (MQ=255)
                                                                  tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:2186239/37‑1 (MQ=255)
                                                                  tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:4058469/37‑1 (MQ=255)
                                                                  tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTc                 <  1:2181401/37‑1 (MQ=255)
                                                                   |                                                  
CGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTA  >  minE/1747809‑1747926

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: