Converting Bisulfite-seq BAM files to fastq
Entering edit mode
3.5 years ago
divyswar01 • 0

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.

Methylation Bisulfite-Seq TCGA BAM • 897 views
Entering edit mode

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.


Login before adding your answer.

Traffic: 2007 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6