Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 346,941 C→G 47.4% W224S (TGG→TCG)  ybbO ← predicted oxidoreductase with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE346,9410CG47.4% 1.9 / 66.3 57W224S (TGG→TCG) ybbOpredicted oxidoreductase with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base C (12/18);  new base G (7/20);  total (19/38)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

TTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATC  >  minE/346903‑347003
                                      |                                                              
ttCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGt                                                 <  1:1606600/54‑1 (MQ=255)
ttCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGt                                                 <  1:1600141/54‑1 (MQ=255)
ttCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAgcgg                                  >  1:1318239/1‑69 (MQ=255)
ttCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAgcgg                                  >  1:409230/1‑69 (MQ=255)
ttCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAgcgg                                  >  1:413252/1‑69 (MQ=255)
ttCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAgcgg                                  >  1:100794/1‑69 (MQ=255)
ttCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAgcgg                                  >  1:2234890/1‑69 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:616847/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:3076560/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:1247461/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:691241/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:1154426/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:2118639/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:2962193/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:2323055/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:1966882/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:1852809/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGggcggc                             <  1:2493101/69‑1 (MQ=255)
     cTTCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCGGTCCCAACGTAAAGCGGgcggcg                            <  1:530892/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGTTCCCAACGTAAAGCGGgcggcg                            <  1:1838553/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:2609368/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:2869927/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:2969644/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            >  1:2473777/1‑69 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            >  1:3041730/1‑69 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            >  1:44408/1‑69 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:476639/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:604670/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            >  1:672423/1‑69 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:792759/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:965883/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:187028/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:1115438/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:1177022/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:1205063/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:1244777/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            >  1:1613798/1‑69 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:1733536/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:1926541/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            >  1:214049/1‑69 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:243324/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:2257716/69‑1 (MQ=255)
      ttCTCGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            >  1:2142402/1‑69 (MQ=255)
       tctcGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:3262494/68‑1 (MQ=255)
       tctcGCTAATAAAAGCATGGCGCACTTTGTCGACCACCGCTTCCGGTCCCAACGTAAAGCGGgcggcg                            <  1:652614/68‑1 (MQ=255)
          cGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGc                        >  1:2814209/1‑69 (MQ=255)
          cGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGc                        >  1:2340329/1‑69 (MQ=255)
          cGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGc                        >  1:3066401/1‑69 (MQ=255)
                 aaaaGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGAtt                 >  1:515492/1‑69 (MQ=255)
                 aaaaGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGAtt                 >  1:1801472/1‑69 (MQ=255)
                 aaaaGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGAtt                 >  1:2878208/1‑69 (MQ=255)
                 aaaaGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGAtt                 >  1:1192052/1‑69 (MQ=255)
                             tACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGttt     <  1:2536551/69‑1 (MQ=255)
                             tACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGttt     <  1:780959/69‑1 (MQ=255)
                             tACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGttt     <  1:3162431/69‑1 (MQ=255)
                              aCTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTTGttt     <  1:3328116/68‑1 (MQ=255)
                                 ttGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATc  <  1:1133827/68‑1 (MQ=255)
                                      |                                                              
TTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTCCGGTCCCAACGTAAAGCGGGCGGCGATGCCGGGATTTTCGACTGGTTTATC  >  minE/346903‑347003

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: