Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,014,650 | G→C | 12.5% | A473A (GCC→GCG) | 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,650 | 0 | G | C | 12.5% | 83.7 / 4.3 | 32 | A473A (GCC→GCG) | fumA | aerobic Class I fumarate hydratase |
Reads supporting (aligned to +/- strand): ref base G (13/15); new base C (2/2); total (15/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01 |
TACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGC > minE/1014580‑1014716 | tACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAg > 1:910951/1‑71 (MQ=255) ccTGCGCCAATACAGCGGCCGGCCCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATcc > 1:343947/1‑71 (MQ=255) gcgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAg > 1:621400/1‑62 (MQ=255) gtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgc > 1:2303198/2‑71 (MQ=255) tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc < 1:1881849/64‑1 (MQ=255) aCAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgctg > 1:4007781/1‑71 (MQ=255) aCAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgctg > 1:2401664/1‑71 (MQ=255) gCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctggct > 1:1782672/1‑71 (MQ=255) cGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTTCAGGCATCCGTCACCTgctggct > 1:1105988/1‑70 (MQ=255) ggCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGc > 1:3350080/1‑71 (MQ=255) ccGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATc > 1:3629039/1‑52 (MQ=255) ccaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGc > 1:4037026/1‑44 (MQ=255) ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc > 1:446543/1‑44 (MQ=12) cGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGt < 1:708229/46‑1 (MQ=255) aTACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGc < 1:1980606/71‑1 (MQ=255) cTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgc > 1:2080275/1‑48 (MQ=255) ccgaGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCt < 1:2282755/57‑1 (MQ=255) aaGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGa > 1:586954/1‑71 (MQ=255) aGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGTCTGCGGTTGCCTTTCGCCAGCATGAt < 1:2841483/71‑1 (MQ=255) aGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGc < 1:3700317/66‑1 (MQ=255) aGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGc < 1:2502461/66‑1 (MQ=255) ggTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgct < 1:2800780/42‑1 (MQ=255) gAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATa < 1:562015/70‑1 (MQ=255) gccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGc > 1:3228494/1‑71 (MQ=255) gtgtTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGAt < 1:227388/57‑1 (MQ=255) ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGa < 1:1671168/53‑1 (MQ=255) ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcg > 1:4009513/1‑71 (MQ=255) ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgc < 1:1063912/70‑1 (MQ=255) ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATc > 1:2770714/1‑53 (MQ=255) ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcgc < 1:4126777/70‑1 (MQ=255) tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg < 1:1973623/70‑1 (MQ=255) tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg < 1:583684/70‑1 (MQ=255) tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg < 1:2135436/70‑1 (MQ=255) aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc < 1:426451/70‑1 (MQ=255) | TACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGC > minE/1014580‑1014716 |
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: ‑ |