Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,717,162 | G→C | 32.1% | V344V (GTC→GTG) | glyA ← | serine hydroxymethyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,717,162 | 0 | G | C | 32.1% | 58.1 / 38.8 | 56 | V344V (GTC→GTG) | glyA | serine hydroxymethyltransferase |
Reads supporting (aligned to +/- strand): ref base G (23/15); new base C (10/8); total (33/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.77e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 |
CCTACACGAATACCGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTT > minE/1717100‑1717218 | ccTACACGAATACCGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTg > 1:1008408/1‑69 (MQ=255) ccTACACGAATACCGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTg > 1:3225685/1‑69 (MQ=255) cGAATACCGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGt > 1:3414892/1‑62 (MQ=255) ccGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACgg > 1:3193094/1‑69 (MQ=255) ccGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACgg > 1:2171560/1‑69 (MQ=255) ccGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACgg > 1:2421611/1‑69 (MQ=255) ccGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACgg > 1:1856404/1‑69 (MQ=255) ccGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACgg > 1:2598089/1‑69 (MQ=255) ccGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACgg > 1:2694814/1‑69 (MQ=255) ccGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACgg > 1:2628992/1‑69 (MQ=255) ccGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGATAGCACgg > 1:2451306/1‑69 (MQ=255) cGGAAGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACgg > 1:2862560/1‑68 (MQ=255) cGGAAGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACgg > 1:2256369/1‑68 (MQ=255) cGGAAGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACgg > 1:2298169/1‑68 (MQ=255) cGGAAGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACgg > 1:1965979/1‑68 (MQ=255) cGGAAGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACgg > 1:2705520/1‑68 (MQ=255) cGGAAGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACgg > 1:2911322/1‑68 (MQ=255) cGGAAGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACgg > 1:155033/1‑68 (MQ=255) cGGAAGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACgg > 1:3014683/1‑68 (MQ=255) cGGAAGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACgg > 1:536300/1‑68 (MQ=255) cGGAAGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACgg > 1:1016895/1‑68 (MQ=255) cacaAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAgc > 1:926408/1‑69 (MQ=255) cacaAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAgc > 1:2940335/1‑69 (MQ=255) acaAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAgcg > 1:1147538/1‑69 (MQ=255) acaAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAgcg > 1:1640017/1‑69 (MQ=255) aaaCGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAgcggc < 1:3178509/69‑1 (MQ=255) ctTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGcttc < 1:116412/69‑1 (MQ=255) ctTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGcttc < 1:2480521/69‑1 (MQ=255) cGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAgcggc > 1:3251469/1‑57 (MQ=255) cGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAgcg > 1:241721/1‑55 (MQ=255) cGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGt > 1:35184/1‑59 (MQ=255) cGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCttt > 1:1443585/1‑69 (MQ=255) aTCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTAcc < 1:1191983/69‑1 (MQ=255) aTCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTAcc < 1:1391374/69‑1 (MQ=255) aTCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTAcc < 1:1593955/69‑1 (MQ=255) aTCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTAcc < 1:3084066/69‑1 (MQ=255) aTCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTAcc < 1:2962308/69‑1 (MQ=255) aTCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTAcc < 1:300092/69‑1 (MQ=255) tCGTTCGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTAcc < 1:671070/68‑1 (MQ=255) cGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAgcg > 1:2594118/1‑46 (MQ=255) cGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGt > 1:2657554/1‑50 (MQ=255) cGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAg < 1:1496809/69‑1 (MQ=255) cGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAg > 1:2927510/1‑69 (MQ=255) cGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAg < 1:1418598/69‑1 (MQ=255) cGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAg < 1:1816187/69‑1 (MQ=255) cGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAg > 1:163935/1‑69 (MQ=255) cGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAg < 1:1980043/69‑1 (MQ=255) cGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAg < 1:615667/69‑1 (MQ=255) cGGTACGCTGTTTTTGTTCACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTAcc < 1:1625680/63‑1 (MQ=255) cTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGTCGTCTGCTTCTTTACCGGTCAGGttttt < 1:917123/68‑1 (MQ=255) cTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGttttt < 1:3046363/68‑1 (MQ=255) cTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGttttt < 1:3080202/68‑1 (MQ=255) cTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGttttt < 1:3263162/68‑1 (MQ=255) cTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGttttt < 1:1744016/68‑1 (MQ=255) cTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGttttt < 1:847329/68‑1 (MQ=255) cTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGttttt < 1:2251001/68‑1 (MQ=255) gttgACGGTGATGTTAGCACGGCCCAGAGCGGCGt > 1:1665577/1‑35 (MQ=255) | CCTACACGAATACCGGAGGTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTT > minE/1717100‑1717218 |
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: ‑ |