Question: Converting Bisulfite-seq BAM files to fastq
gravatar for divyswar01
2.1 years ago by
divyswar010 wrote:

TCGA has some DNA Methylation files here But they only have BAM files here. I need fastq files for these BAM files. I have contacted GDC and they said they only have BAMs available. Since fastqs in bisulfite context are aligned with Bisulfite-aware aligner, I don't think converting BAM to fastq is straightforward. I tried using samtools bam2fq function and it produces fastq files containing reads with very low base quality and excessive amount of repeats. When I aligned these converted fastqs I got almost 0 alignment rate.

I am wondering if there is a way to convert bisulfite sequencing BAMs back to fastqs. It would be more convenient if TCGA had fastqs available but since that is not an option I am not sure how to regenerate the fastqs.

ADD COMMENTlink modified 23 months ago by Biostar ♦♦ 20 • written 2.1 years ago by divyswar010

Hi, is it an old sequencing ? Do you know the technology used (maybe the aligner consider an older phred score ? ) ? I'm a bit surprise because the BAM files is generated in a way to conserved the original quality of sequencing bases so i couldn't imagine an original fastq with better bases qualities.

ADD REPLYlink written 23 months ago by Titus910
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: 765 users visited in the last hour