| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,014,521 | A→G | 27.3% | I516I (ATT→ATC) | 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,521 | 0 | A | G | 27.3% | 46.4 / 19.6 | 33 | I516I (ATT→ATC) | fumA | aerobic Class I fumarate hydratase |
| Reads supporting (aligned to +/- strand): ref base A (7/17); new base G (4/5); total (11/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.38e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01 | |||||||||||
GTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGC > minE/1014454‑1014590 | gTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAt < 1:786651/70‑1 (MQ=255) tATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAttt < 1:3271043/71‑1 (MQ=255) gctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCa < 1:1875849/70‑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) tCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCa < 1:1330090/70‑1 (MQ=255) cATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCAtt > 1:22054/1‑71 (MQ=255) aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTc < 1:1984600/71‑1 (MQ=255) aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTc > 1:3510638/1‑71 (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:412057/1‑70 (MQ=255) atcCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGa < 1:3450025/70‑1 (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) aTAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGata < 1:2356006/62‑1 (MQ=255) tAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc < 1:2867404/64‑1 (MQ=255) tAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACAt > 1:1709500/1‑71 (MQ=255) aaCGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc < 1:2422416/62‑1 (MQ=255) aCGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGat > 1:1281622/1‑57 (MQ=255) cGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa < 1:2376214/71‑1 (MQ=255) cGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGata < 1:3744670/54‑1 (MQ=255) cGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGc < 1:857845/71‑1 (MQ=255) gggAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTc < 1:3738673/65‑1 (MQ=255) aaTCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCgg > 1:1256131/1‑46 (MQ=255) ccACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATAc < 1:1985466/69‑1 (MQ=255) aaTTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTgc < 1:1216150/70‑1 (MQ=255) | GTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGC > minE/1014454‑1014590 |
| 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: ‑ |