Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,014,734 | A→G | 13.6% | D445D (GAT→GAC) | fumA ← | aerobic Class I fumarate hydratase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,014,734 | 0 | A | G | 13.6% | 105.9 / 9.3 | 44 | D445D (GAT→GAC) | fumA | aerobic Class I fumarate hydratase |
Reads supporting (aligned to +/- strand): ref base A (20/18); new base G (2/4); total (22/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.64e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01 |
GCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCG > minE/1014664‑1014804 | gctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGaa > 1:345416/1‑71 (MQ=255) gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatcc < 1:664677/70‑1 (MQ=255) actgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtcca < 1:3597046/70‑1 (MQ=255) ctgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtccat < 1:1565129/71‑1 (MQ=255) tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccgtc > 1:453332/1‑71 (MQ=255) tgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtcc < 1:1060785/68‑1 (MQ=255) cGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccg > 1:2532283/1‑71 (MQ=255) ggTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgt < 1:816599/71‑1 (MQ=255) gCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatc < 1:2955553/60‑1 (MQ=255) accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcctg > 1:723649/2‑69 (MQ=255) accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcct < 1:3140327/69‑2 (MQ=255) ccTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcggc > 1:1770302/1‑71 (MQ=255) cTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtccat < 1:2601338/61‑1 (MQ=255) cTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtc > 1:1039793/1‑66 (MQ=255) tttCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCctgc < 1:586561/69‑1 (MQ=255) ttCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtccatac > 1:118187/1‑59 (MQ=255) ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt > 1:685522/1‑71 (MQ=255) ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt > 1:308783/1‑71 (MQ=255) ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt > 1:2013825/1‑71 (MQ=255) tCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGtg > 1:3941772/1‑70 (MQ=255) gCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatcc > 1:654576/1‑54 (MQ=255) ccAGCATGATCATGCTGCCGCCGTG‑GATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGtt < 1:852026/70‑1 (MQ=255) ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGtt > 1:2745306/1‑71 (MQ=255) ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGtt < 1:3893653/71‑1 (MQ=255) cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg < 1:2163215/71‑1 (MQ=255) cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg < 1:2047111/71‑1 (MQ=255) aGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg > 1:1354929/1‑70 (MQ=255) gCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcggc > 1:62955/1‑61 (MQ=255) aTGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc < 1:1716779/70‑1 (MQ=255) tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTgg > 1:3926520/1‑67 (MQ=255) cATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCaa > 1:643730/1‑68 (MQ=255) tCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCAGGCGGTCGTTGGGCCAAGAGAAcc > 1:474120/1‑69 (MQ=255) ccGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg > 1:960034/1‑69 (MQ=255) ccGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg > 1:1328016/1‑69 (MQ=255) gCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAgg < 1:3066225/71‑1 (MQ=255) ccTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAg > 1:1075183/1‑68 (MQ=255) tGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcggcg < 1:3785807/42‑1 (MQ=255) tGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcggcg < 1:2062144/42‑1 (MQ=255) tGCCCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCAt > 1:463026/1‑70 (MQ=255) cgcTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATa < 1:1936669/69‑1 (MQ=255) gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCa > 1:684540/1‑62 (MQ=255) gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCtt < 1:2913606/69‑1 (MQ=255) gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTc < 1:2428859/70‑1 (MQ=255) cAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCg < 1:3118911/71‑1 (MQ=255) gTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGc > 1:1241017/1‑71 (MQ=255) gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTgg < 1:1774378/37‑1 (MQ=255) gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc > 1:3009677/1‑40 (MQ=255) aTCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGttt < 1:2785296/71‑1 (MQ=255) aTCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGttt < 1:3748408/71‑1 (MQ=255) aTCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAag < 1:1429219/42‑1 (MQ=255) tCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCt < 1:467292/60‑1 (MQ=255) gACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTg < 1:2334704/70‑1 (MQ=255) aTAAGAATCCATCCGTCCGGCGGTCGTTGGGCCaa > 1:1737134/1‑35 (MQ=255) tAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCgg < 1:3038719/71‑1 (MQ=255) atccatccGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCcg < 1:2553805/71‑1 (MQ=255) | GCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCG > minE/1014664‑1014804 |
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: ‑ |