Question: RRBS: import Bismark output to RnBeads
0
gravatar for cvu
2.1 years ago by
cvu150
India
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

bismark2report

bismark2summary

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.

Help
Access

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