| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,014,491 | A→C | 26.7% | L526L (CTT→CTG) | 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,491 | 0 | A | C | 26.7% | 29.2 / 20.9 | 30 | L526L (CTT→CTG) | fumA | aerobic Class I fumarate hydratase |
| Reads supporting (aligned to +/- strand): ref base A (9/13); new base C (3/5); total (12/18) | |||||||||||
| 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.04e-01 | |||||||||||
TATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAA > minE/1014423‑1014561 | tATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAg > 1:802756/1‑71 (MQ=255) ttCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACaa < 1:574330/67‑1 (MQ=255) cacacaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa < 1:1409133/70‑1 (MQ=255) cacaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAAc < 1:2397066/71‑1 (MQ=255) acaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACg > 1:766524/1‑71 (MQ=255) aTTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCt < 1:879644/71‑1 (MQ=255) gtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTcc < 1:3225121/69‑1 (MQ=255) gAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACg > 1:520733/1‑53 (MQ=255) gAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACg > 1:3318998/1‑53 (MQ=255) gTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAt < 1:786651/70‑1 (MQ=255) tATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAttt < 1:3271043/71‑1 (MQ=255) gctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCa < 1:1875849/70‑1 (MQ=255) gctgGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGaa < 1:1684018/49‑1 (MQ=255) ctgGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCCGGAAATCTTCCACTTCAATTTTCCAg < 1:1928436/70‑1 (MQ=255) gaagaaGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGAt < 1:1659224/69‑1 (MQ=255) aagaagTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGc < 1:2995623/71‑1 (MQ=255) agaagTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCt < 1:1022793/71‑1 (MQ=255) agaagTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCt < 1:3026090/71‑1 (MQ=255) aagTCATTTCCTTTATCATCCACAAGGATAAACGc > 1:879320/1‑35 (MQ=255) gTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCt < 1:2939095/43‑1 (MQ=255) tCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCa < 1:1330090/70‑1 (MQ=255) cATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCAtt > 1:22054/1‑71 (MQ=255) aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTc > 1:3510638/1‑71 (MQ=255) aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTc < 1:1984600/71‑1 (MQ=255) ttCCTTTATCATCCACAAGGATAAACGCCGGGAGATCTTCCACTTCaa > 1:2030932/1‑48 (MQ=255) accTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATAcc > 1:3306792/2‑69 (MQ=255) tttGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAgct > 1:199777/1‑69 (MQ=255) ttgtcatcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTc > 1:122209/4‑71 (MQ=255) atcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGat < 1:2037959/71‑1 (MQ=255) atcCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGa < 1:3450025/70‑1 (MQ=255) atcCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGa > 1:412057/1‑70 (MQ=255) acaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc < 1:3181219/71‑1 (MQ=255) aaGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCaa > 1:2741326/1‑71 (MQ=255) aaGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCaa < 1:2857610/71‑1 (MQ=255) aaGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCaa > 1:843968/1‑71 (MQ=255) | TATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAA > minE/1014423‑1014561 |
| 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: ‑ |