Question: RRBS: import Bismark output to RnBeads
gravatar for cvu
2.1 years ago by
cvu150 wrote:

Hello All,

I am running Bisulfide sequencing samples. I've directional libraries with four single-ended fastq files for each sample and completed with bismark and got four bam files which i've merged and finally got one bam file for each sample.

These are the commands I ran:

bismark --genome_folder /home/xyz/mm10 read1.fastq,read2.fastq,read3.fastq,read4.fastq -o /home/out/ --parallel 3 -p 10 --temp_dir /home/tmp_dir/ 

samtools merge out.bam read1.bam read2.bam read3.bam read4.bam

bismark_methylation_extractor --bedGraph -o /home/out/ out.bam



bam2nuc --genome_folder  /home/xyz/mm10 out.bam --dir /home/out/

I have got around 26 output files. now I want to analyse outputs with RnBeads R package, but I'm not able to import bismark result to RnBeads. Please suggest which bismark output file I should use for RnBeads.

Thanks in advance!! Any help appreciated!!

ADD COMMENTlink modified 2.1 years ago by Biostar ♦♦ 20 • written 2.1 years ago by cvu150
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2115 users visited in the last hour