I'm trying to find a way to convert a GEO WGBS data set into a methylKit-workable form. I've used methylKit for ERRBS data in the past but this data set has very different formatting. This data was processed by the authors with Bismark methylation extractor and is in a txt file. It looks like one of the files is sorted and the other one isn't but it's hard to understand exactly what processing they did. From doing some reading it looks like the data needs to be in a sorted SAM or BAM file format in order to use processBismarkAln in methylKit. I'm not familiar with Bismark and I haven't used Perl very much. Does anyone have any advice as to how to convert it based on it's current form? Does anyone have an example of the SAM file format that will work with with processBismarkAln?
This is the first data set:
The second file:
0 answers
No answers yet.
Log in to answer this question.
-
Error: free(): invalid next size (fast) and Aborted (core dumped) while using R in Linux server
written by Sowmya Pulapet •
-
Use bismark files to run methylKit
written by anasjamshed
-
Remove Chr from methylkit analysis
written by janayfox95 •
-
Methylkit (differential methylation analysis) error in data frame
written by harshraje19 •
-
Bismark to MethylKit
written by Arindam Ghosh
-
From Bismark to methylkit, which file is better to use?
written by iraia.munoa
-
Understanding Bismark and methylkit output
written by pooja_shah •
-
error in reading bam file in methylKit
written by chaudharyc61 •
-
Interpreting Bismark Methylation Extractor output
written by Yasin Uzun •
-
apply cufflinks to ENCODE RNA-Seq data
written by liu4gre
Please do not paste screenshots of text. You can use the code formatting option to showcase text better, and even format tabular text so it's easy on the eyes.