Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,860,420 | G→T | 42.1% | L449L (CTG→CTT) | mutL → | methyl‑directed mismatch repair protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,860,420 | 0 | G | T | 42.1% | 49.7 / 74.1 | 83 | L449L (CTG→CTT) | mutL | methyl‑directed mismatch repair protein |
Reads supporting (aligned to +/- strand): ref base G (3/45); new base T (23/12); total (26/57) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.18e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 |
GGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCAC > minE/2860360‑2860484 | ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAgcgc < 1:825617/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAgcgc < 1:1878135/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAgcgc < 1:1867123/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAgcgc < 1:502237/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:317424/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:1791249/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:1809246/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:1494761/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:1420487/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:1402120/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:1330371/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:3208810/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:3434730/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:17328/67‑1 (MQ=255) ggcgAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGGTGCTGGAgcgc < 1:2065550/67‑1 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAgc > 1:2998873/1‑48 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGc > 1:1947663/1‑56 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACa > 1:2128535/1‑60 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACa > 1:1620093/1‑60 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:1916601/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:2067972/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:1920778/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:1077786/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:2099653/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:2303976/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:2597839/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:2920261/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:3127185/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:3138428/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:3316540/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:545512/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:617553/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:858870/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:1196012/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:1611910/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:1657237/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:178659/1‑61 (MQ=255) ttGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACAt > 1:1836082/1‑61 (MQ=255) gTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACtttt < 1:1698779/67‑1 (MQ=255) tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt < 1:1249398/68‑1 (MQ=255) tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt < 1:3224927/68‑1 (MQ=255) tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt < 1:2209277/68‑1 (MQ=255) tCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTAt < 1:658208/68‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc < 1:3129408/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc < 1:2294070/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc < 1:1593080/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc < 1:2351107/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc < 1:1099131/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc < 1:1560742/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc < 1:3327439/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc < 1:538611/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc < 1:640577/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc < 1:682065/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc < 1:206960/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc < 1:484351/69‑1 (MQ=255) tACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc < 1:1849690/69‑1 (MQ=255) aCTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGcc > 1:2159661/1‑69 (MQ=255) tGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc < 1:661426/38‑1 (MQ=255) ccATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAAc < 1:1454079/68‑1 (MQ=255) ttCCGACTGTCCGTTCTTGGAGCGCGACGGCAACAATTCACTTTTATCCTTGc < 1:1558319/53‑1 (MQ=255) tCCGACTGTGCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc < 1:865534/52‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:537392/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:1202238/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:70436/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:598319/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:928653/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:1023358/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:1900866/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:494510/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:493841/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:120397/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:3309405/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:2980272/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:1571837/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:2479329/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:2472657/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:1937381/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:1895417/46‑1 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTgg > 1:3048764/1‑68 (MQ=255) ccGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGc > 1:2009470/1‑69 (MQ=255) cGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc < 1:1168025/45‑1 (MQ=255) gCGTTGCTTGAGCGCGATGGCAACATTTCACTTTTATCCTTGc < 1:1237948/43‑1 (MQ=38) gTTGCTGGAGCGCGACGTCAACATTTCACTTTTATCCTTGCCAGTGGCAGAAc < 1:1558147/53‑1 (MQ=255) ggAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGcac < 1:1203891/65‑1 (MQ=255) | GGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCAC > minE/2860360‑2860484 |
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: ‑ |