Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,610,889 G→A 46.5% C283C (TGC→TGT eutB ← ethanolamine ammonia‑lyase, large subunit, heavy chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,610,8890GA46.5% 2.9 / 33.1 43C283C (TGC→TGTeutBethanolamine ammonia‑lyase, large subunit, heavy chain
Reads supporting (aligned to +/- strand):  ref base G (14/9);  new base A (7/13);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

AGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCT  >  minE/1610844‑1610937
                                             |                                                
aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt                            >  1:2779911/1‑68 (MQ=255)
aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt                            >  1:1353132/1‑68 (MQ=255)
aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt                            >  1:1355978/1‑68 (MQ=255)
aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt                            >  1:3273839/1‑68 (MQ=255)
aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt                            >  1:2245818/1‑68 (MQ=255)
aGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGt                            >  1:2203472/1‑68 (MQ=255)
            tAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGcc               <  1:2002504/69‑1 (MQ=255)
            tAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGcc               <  1:1998471/69‑1 (MQ=255)
               cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc            >  1:1269724/1‑69 (MQ=255)
               cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc            >  1:3192928/1‑69 (MQ=255)
               cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc            >  1:2768371/1‑69 (MQ=255)
               cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc            >  1:2537355/1‑69 (MQ=255)
               cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc            >  1:82689/1‑69 (MQ=255)
               cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc            >  1:1569268/1‑69 (MQ=255)
               cgcCGAGCCTTGTCCGGTTTCGAAGTAGAGACAGGTTTCCCCGGCGATACGGTTGAACTCCGCGCCCAc            >  1:3220171/1‑69 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               >  1:2951729/1‑44 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               >  1:1344329/1‑44 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               <  1:1039022/44‑1 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               >  1:2732979/1‑44 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               <  1:1481749/44‑1 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               <  1:2696627/44‑1 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               >  1:255887/1‑44 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               <  1:1877697/44‑1 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               <  1:3301059/44‑1 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               >  1:2472639/1‑44 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               >  1:2356273/1‑44 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               <  1:827701/44‑1 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               >  1:2242668/1‑44 (MQ=255)
                     gCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               >  1:3404716/1‑44 (MQ=255)
                      ccTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATAc                               <  1:2834453/43‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:828916/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:996661/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:568803/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:2860712/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:2803665/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:2711342/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:2477371/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:2059545/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:1695411/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:159905/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:1539465/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:1504128/69‑1 (MQ=255)
                         tGTCCGGTTTCGAAGTAGAGACAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGTGCt  <  1:1469656/69‑1 (MQ=255)
                                             |                                                
AGCGCCAGCGGATAGCGCAGAGCCTTGTCCGGTTTCGAAGTAGAGGCAGTTTTCCCCGGCGATACGGTTGAACTCCGCGCCCACTGCGCGCGCT  >  minE/1610844‑1610937

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: