Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | 40 | 30,807 | 2 bp→GA | intergenic (–/–) | – / – | –/– |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | 40 | 30,807 | 0 | C | G | 100.0% | 115.6 / NA | 37 | intergenic (–/–) | –/– | –/– |
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (20/17); total (20/17) | |||||||||||
* | 40 | 30,808 | 0 | C | A | 97.3% | 113.2 / NA | 37 | intergenic (–/–) | –/– | –/– |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base A (20/16); total (20/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.59e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 |
TCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGCGCGTGTCATACAAGGTTTCGGCGGGGGCGGCCCCCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAGGTGAACGGCAAAAC > 40/30712‑30908 || gaataccccgaagataatcgggttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGAccgc > 1:1076218‑M2/91‑101 (MQ=255) gaataccccgaagataatcgggttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGAccgc > 2:1590623‑M2/91‑101 (MQ=255) aataccccgaagataatcgggttggaaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGAccgcc < 2:199144‑M2/12‑1 (MQ=255) ataccccgaagataatcgggttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGAccgccg > 2:1590810‑M2/89‑101 (MQ=255) ataccccgaagataatcgggttgggaagcacgagccgcacgattttagtttgacccgtcgagatctgattaaagtgagcgccgcaacgGCGGCGAccgccg > 2:1709757‑M2/89‑101 (MQ=255) agataatcgggttgggaagcacgagccgcactattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATc < 1:1078725‑M2/22‑1 (MQ=255) agataatcgggttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgTCGGCGACCGCCGTGGTTTATc > 1:465430‑M2/80‑101 (MQ=255) atcgggttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCAt > 1:1689069‑M2/75‑101 (MQ=255) tggggtgggaaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaccgGCGGCGACCGCCGTGGTTTATCCTCAtt < 1:607838‑M2/28‑1 (MQ=255) gttgggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTAc > 1:991082‑M2/70‑101 (MQ=255) ggaagcacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGCCCGCCGTGGTTTATCCTCATTCTACGCtt < 2:1689063‑M2/36‑1 (MQ=255) cacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGc > 2:1177343‑M2/61‑101 (MQ=255) cacgagccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGc > 1:1088065‑M2/61‑101 (MQ=255) ccgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCNTTCTANGCTTGCGGCAAGCGt > 2:1069908‑M2/55‑101 (MQ=255) cgcacgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGtt > 2:1513357‑M2/54‑101 (MQ=255) cacgattttagtttgccccgtcgagatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTcc < 1:1709763‑M2/50‑1 (MQ=255) acgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCg < 1:1710387‑M2/51‑1 (MQ=255) acgattttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCg < 1:1106823‑M2/51‑1 (MQ=255) tttcgtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGCCGACCGCCGTGGTATATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGc > 2:462125‑M2/46‑101 (MQ=255) tttagtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGc > 2:1626794‑M2/46‑101 (MQ=255) gtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTaca < 2:1459561‑M2/60‑1 (MQ=255) gtttgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGCTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTaca > 2:747359‑M2/42‑101 (MQ=255) tgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACAccc > 1:918702‑M2/39‑101 (MQ=255) tgacccgtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACAccc < 2:245906‑M2/63‑1 (MQ=255) gtcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCa > 1:946472‑M2/33‑101 (MQ=255) tcgcgatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTCCCGGCAGCTACACCCGCGCCag > 1:231755‑M2/32‑101 (MQ=255) gcgatctgattaaagtgaccgccgcagcgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCagag < 2:232834‑M2/72‑1 (MQ=255) gatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGAt > 1:1715335‑M2/28‑101 (MQ=255) gatctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGAt > 1:1371656‑M2/28‑101 (MQ=255) tctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATaa < 2:1353903‑M2/76‑1 (MQ=255) tctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATaa < 2:1354578‑M2/76‑1 (MQ=255) ctgattaaagtgagcgccgcaacgTCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAAt < 2:465425‑M2/77‑1 (MQ=255) ctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAAt < 2:1076213‑M2/77‑1 (MQ=255) ctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAAt < 1:1349251‑M2/77‑1 (MQ=255) ctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAAt > 1:1296569‑M2/25‑101 (MQ=255) ctgattaaagtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAAt < 1:1702168‑M2/77‑1 (MQ=255) aagtgagngcngcaangGCGGNGACCGCCGTGGTTTATNCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCTCCAGAGATAATGCCCCTg < 2:1203060‑M2/84‑1 (MQ=255) gtgagcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGac > 2:973106‑M2/16‑101 (MQ=255) agcgccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACt < 2:991077‑M2/89‑1 (MQ=255) ccgcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCAGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAg > 2:769146‑M2/9‑101 (MQ=255) gcaacgGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAGGt > 1:763240‑M2/7‑101 (MQ=255) gGCGGCGACCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTCCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAGGTGAACg < 2:231755‑M2/100‑1 (MQ=255) cccGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAGGTGAACGGCAAAAc < 1:519009/101‑1 (MQ=255) || TCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGCGCGTGTCATACAAGGTTTCGGCGGGGGCGGCCCCCGCCGTGGTTTATCCTCATTCTACGCTTGCGGCAAGCGTTCCGGCAGCTACACCCGCGCCAGAGATAATGCCCCTGACACTTAAGGTGAACGGCAAAAC > 40/30712‑30908 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |