Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 452,584 | T→G | 26.8% | P183P (CCT→CCG) | phr → | deoxyribodipyrimidine photolyase, FAD‑binding |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 452,584 | 0 | T | G | 26.8% | 64.2 / 31.1 | 56 | P183P (CCT→CCG) | phr | deoxyribodipyrimidine photolyase, FAD‑binding |
Reads supporting (aligned to +/- strand): ref base T (16/25); new base G (8/7); total (24/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.75e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01 |
GATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAA > minE/452529‑452651 | gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat > 1:521613/1‑66 (MQ=255) gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat > 1:2542580/1‑66 (MQ=255) gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat > 1:2519245/1‑66 (MQ=255) gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat > 1:2660126/1‑66 (MQ=255) gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat > 1:612876/1‑66 (MQ=255) gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat > 1:2031065/1‑66 (MQ=255) gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat > 1:2737384/1‑66 (MQ=255) ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg > 1:316003/1‑69 (MQ=255) ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg > 1:2482986/1‑69 (MQ=255) ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg > 1:798798/1‑69 (MQ=255) ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg > 1:2247162/1‑69 (MQ=255) ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg > 1:2243116/1‑69 (MQ=255) ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg > 1:2131850/1‑69 (MQ=255) ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg > 1:1023626/1‑69 (MQ=255) ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg > 1:171342/1‑69 (MQ=255) gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGATCCCGCGCCATCCATTACGCTGAATTATc < 1:347568/69‑1 (MQ=255) gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc < 1:861349/69‑1 (MQ=255) gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc < 1:3322340/69‑1 (MQ=255) gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc < 1:1967639/69‑1 (MQ=255) gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc < 1:3016770/69‑1 (MQ=255) gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc < 1:2411340/69‑1 (MQ=255) gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc < 1:1356270/69‑1 (MQ=255) cGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATcc > 1:2829707/1‑69 (MQ=255) gcgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg < 1:3282142/69‑1 (MQ=255) gcgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg < 1:1255017/69‑1 (MQ=255) gcgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg < 1:1499769/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt < 1:2214036/67‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt < 1:3371760/67‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt < 1:2937722/67‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt < 1:206218/67‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt < 1:1754501/67‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt < 1:3172490/67‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:584911/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:596540/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:869025/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:3330189/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:3240963/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:912557/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:1749968/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:2110332/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:2207286/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:2327249/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:2655848/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:2431938/69‑1 (MQ=255) ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCT‑‑CCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:3069582/67‑1 (MQ=255) aaaaGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt < 1:743626/67‑1 (MQ=255) tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCAtttcc > 1:2363561/1‑65 (MQ=255) tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTccc > 1:2850948/1‑68 (MQ=255) tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg > 1:3107005/1‑69 (MQ=255) tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg > 1:2632027/1‑69 (MQ=255) tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg > 1:2788698/1‑69 (MQ=255) tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg > 1:332464/1‑69 (MQ=255) tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg > 1:294648/1‑69 (MQ=255) tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg > 1:3109170/1‑69 (MQ=255) agCCCTCGCCATACATTACGCTGAATTATCCTCGTCAGTCt < 1:2529100/41‑1 (MQ=255) cTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGaaaa < 1:820525/69‑1 (MQ=255) | GATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAA > minE/452529‑452651 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |