Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 910,465 G→A 71.5% A126A (GCG→GCA ycjU → predicted beta‑phosphoglucomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE910,4650GA71.5% 23.2 / 8.8 28A126A (GCG→GCAycjUpredicted beta‑phosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base A (18/2);  total (21/7)
Fisher's exact test for biased strand distribution p-value = 9.47e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACAGCAGATCTCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCGGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGATGCTTCCCAACTT  >  minE/910406‑910516
                                                           |                                                   
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:2036879/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:37557/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:3317679/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:3286040/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:3253909/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:323445/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:3167045/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:3047790/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:2996735/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:38560/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:1819036/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:1748484/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:1744622/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:1693819/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:157201/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:1221360/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:1175696/1‑69 (MQ=255)
gCACAGCAGATTCCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCAGCGCTGGAg                                            >  1:782169/1‑69 (MQ=255)
                               tctcCCTGAATGCGCCGACGATTTTAGCGGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGa               <  1:1219822/67‑1 (MQ=255)
                                       aaTGCGCCGACGATTTTAGCGGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGATGCTTCCCaa     <  1:3160213/69‑1 (MQ=255)
                                         tGCTCCGACGATTTTAGCAGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGACGCTTCCCAAc    <  1:2012022/68‑1 (MQ=255)
                                         tGCGCCGACGATTTTAGCGGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGATGCTTCCCaa     <  1:101524/67‑1 (MQ=255)
                                         tGCGCCGACGATTTTAGCGGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGATGCTTCCCaa     <  1:3047644/67‑1 (MQ=255)
                                         tGCGCCGACGATTTTAGCGGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGATGCTTCCCaa     <  1:159540/67‑1 (MQ=255)
                                         tGCGCCGACGATTTTAGCAGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGACGCTTCCCAAc    <  1:1913131/68‑1 (MQ=255)
                                           cgcCGACGATTTTAGCGGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGATGCTTCCCAACtt  >  1:224940/1‑68 (MQ=255)
                                           cgcCGACGATTTTAGCGGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGATGCTTCCCAACtt  >  1:2641166/1‑68 (MQ=255)
                                           cgcCGACGATTTTAGCGGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGATGCTTCCCAACtt  >  1:1097591/1‑68 (MQ=255)
                                                           |                                                   
GCACAGCAGATCTCGGTTGGGCTGGCTTCTGTCTCCCTGAATGCGCCGACGATTTTAGCGGCGCTGGAGCTGCGCGAGTTTTTCACCTTCTGCGCGGATGCTTCCCAACTT  >  minE/910406‑910516

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: