I have bisulfite sequencing data that was aligned and provided in bigwig format. I've since converted the files to bedgraph format, and based off the naming of the files, the BS-seq data was split by +/- strand (i.e., sample1 has two bedgraph files, one bedgraph for the reads on the plus strand, and another bedgraph for reads on the negative strand). I noticed that the start and end coordinates on the plus strand sample to be 1 position lower than the negative strand samples. I appreciate any insight in explaining this discrepancy.
•
written
by
biohack92
1 answer
If the methylation call or mark or event is generated on the 5' side of the cytosine nucleotide on the 5'-CG-3' on the forward/sense strand, the position of the mark or event on the reverse/antisense strand — to make it on the 5' side of the cytosine on the 3'-GC-5' on the reverse stand — would probably need to be one base "forward" as measured with respect to the methylation event on the forward/sense strand.
Log in to answer this question.
-
bedGraphToBigWig: Missing Genome Coordinates
written by RelativeOptics •
-
STAR split reads into two strands
written by Alexander •
-
Output bedGraph from genomeCoverageBed is missing "chr"
written by m98
-
Splitting paired-end RNA-seq data in bedgraph format according to strand
written by virlana.shchuka •
-
BigWig file of strand-specific RNA-seq data
written by yuabrahamliu •
-
Is there any reason to analyze plus/minus strands separately in bisulfite seq/methylation data?
written by biohack92
-
How can you distinguish between whole genome and targeted bisulfite sequencing data?
written by biohack92
-
Multiple files for each sample in bisulfite sequencing dataset
written by biohack92
-
How to find differentially methylated regions from bigwig/wig/bedgraph files?
written by biohack92
-
Negative Strand Coordinates in Fasta Files?
written by bn •
