Comments (2)
The reason for entirely removing C/T SNP positions in a bisulfite experiment is obviously that you cannot trust the result at all.
Let's assume genome 1 would normally be a C, but due to its methylation state you might see it as either a C (methylated) or T (unmethylated) after bisulfite treatment. Genome 2 would always have a T at that position in any case.
While a C in a read would have to be a (methylated) genome 1 read, there is no way to discriminate between a genome 2 or unmethylated genome 1 read, unless the read contains additional informative SNPs at the same time. In other words, the methylation state might not be a real methylation state but could be confounding signal coming simply from the underlying genomic sequence (that just looks like non-methylation).
If you still wanted to include reads covering those C/T positions (I am not sure you should be wanting that...) and I am not mistaken, you should be able to pre-filter your list of SNPs to exclude C/T SNPs, and then run the N-masking and indexing again. This way you wouldn't N-mask and thus include all C/T positions, and therefore also get 'methylation calls'. You should not observe any mapping bias in this case because Cs are in-silico C-to-T converted to Ts for the mapping procedure anyway. Since those positions were not N-masked they would also not participate in allele (or haplotype) assignment, but please be aware that the 'methylation' you observe at these positions might be confounded by genetics, i.e. the T SNP.
Do I make sense?
from snpsplit.
Yes, I do see your point about the confound of looking at methylation at C/T SNPs. The good thing is that I'm studying F1 hybrids, and we have a very good list of non-CG fixed SNPs to assign reads to haplotype, so dissecting out methylation vs. SNP is possible with my particular system. The approach you suggested, excluding C/T SNPs, worked for me. Thank you!
from snpsplit.
Related Issues (20)
- Two much variation in the Two strain specific reads percentage HOT 1
- Using SNPsplit to analysis Bisulfite-Seq data in plant HOT 2
- SNPsplit for mouse RNA-seq data HOT 8
- 0-indexed or 1-indexed for the SNP output? HOT 2
- single-end mode does not work HOT 2
- N-masking human genome HOT 2
- How to phase Hi-C data using SNPs HOT 2
- Why I can't split bisulfite data ? HOT 3
- Mask reference genome for human genome build 37 using bisulphite fasta file HOT 6
- single end mode not working HOT 2
- high reads did not contain one of the expected bases at known SNP positions HOT 7
- SNP level counting showing bi and mono allele count for same gene HOT 2
- Using SNPsplit for single cell data HOT 4
- minimap2 HOT 4
- Genome preparation:129S1 (paternal) mice × B6D2F1/J (F1 of C57BL6NJ × DBA2J, maternal) mice HOT 7
- Update SNPsplit to v8 release of the Mouse Genomes Project HOT 2
- Tutorial for SNPsplit HOT 2
- Getting 0 reads processed HOT 11
- Bowtie1 as the aligner of Bismark HOT 1
- Zero segregation for scRNAseq HOT 5
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from snpsplit.