Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,647,952 C→T 55.8% R120R (CGC→CGT glgB → 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,647,9520CT55.8% 6.8 / 39.6 52R120R (CGC→CGTglgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base C (11/12);  new base T (12/17);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

CTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCAGATACTATGGATGG  >  minE/2647902‑2647993
                                                  |                                         
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:874587/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:2701965/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:2723883/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:2788348/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:2473800/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:2367068/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:233051/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:2318427/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:2085944/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:1773504/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:1765614/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:3080091/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:144670/69‑1 (MQ=255)
cTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:607829/69‑1 (MQ=255)
 tAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:1422830/68‑1 (MQ=255)
 tAATCCAGGAAATGGATGCATGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:2447357/68‑1 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTa                                      >  1:670570/1‑54 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgcg                        >  1:3026503/1‑68 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgcg                        >  1:2896390/1‑68 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgcg                        >  1:2832014/1‑68 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgcg                        >  1:774941/1‑68 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgcg                        >  1:2713325/1‑68 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgcg                        >  1:995189/1‑68 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgcg                        >  1:2595551/1‑68 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgcg                        >  1:1418561/1‑68 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgcg                        >  1:1964427/1‑68 (MQ=255)
  aaTCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         >  1:768225/1‑67 (MQ=255)
  aaTCCAAGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgcg                        >  1:439455/1‑68 (MQ=255)
      cAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCa                    >  1:2606368/1‑68 (MQ=255)
      cAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCa                    >  1:1395116/1‑68 (MQ=255)
      cAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCa                    >  1:3097030/1‑68 (MQ=255)
      cAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCa                    >  1:1419273/1‑68 (MQ=255)
      cAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCa                    >  1:3432582/1‑68 (MQ=255)
      cAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTAACCTGCGCCAGTATGaaa                                 >  1:637500/1‑55 (MQ=255)
         gAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCa                <  1:2543268/69‑1 (MQ=255)
         gAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCa                <  1:2892344/69‑1 (MQ=255)
         gAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCa                <  1:2045338/69‑1 (MQ=255)
         gAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCa                <  1:1706663/69‑1 (MQ=255)
              ggATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCAGATAc           <  1:3236288/69‑1 (MQ=255)
              ggATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCAGATAc           <  1:2644119/69‑1 (MQ=255)
              ggATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCAGATAc           <  1:253830/69‑1 (MQ=255)
                     tGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGa                                   >  1:2944204/1‑38 (MQ=255)
                     tGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGgcgcg                      >  1:889498/1‑51 (MQ=255)
                     tGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGgcgcg                      >  1:1075233/1‑51 (MQ=255)
                       gCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCAGATACTatggatgg  >  1:5806/1‑69 (MQ=255)
                       gCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGAGCGCATGCAGATACTatggatgg  >  1:721558/1‑69 (MQ=255)
                             atCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCa                <  1:3086265/49‑1 (MQ=255)
                             atCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCa                <  1:751310/49‑1 (MQ=255)
                             atCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCa                <  1:2838652/49‑1 (MQ=255)
                             atCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCa                <  1:1329120/49‑1 (MQ=255)
                                tGAAGGTACTCACCTACGTCCGTATGAAACCTTAGgc                         <  1:955417/37‑1 (MQ=255)
                                     gTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCAGATAc           <  1:3107282/46‑1 (MQ=255)
                                                  |                                         
CTAATCCAGGAAATGGATGCCTGGCTATTATCTGAAGGTACTCACCTGCGCCCGTATGAAACCTTAGGCGCGCATGCAGATACTATGGATGG  >  minE/2647902‑2647993

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: