Question: Does whole exome bam files in TCGA have unpaired reads?
gravatar for akhattri
2.9 years ago by
akhattri0 wrote:

I am trying to slice some whole exome bam files from TCGA and for my pipeline I need to extract read pairs. I used samtools sort -n to sort the bam file and using bedtools bamtofastq to extract paired reads but it seems there are lot of reads for which there is no mate in the bam file. I thought I could extract the fastq files from the bam files which will be equivalent to raw sequencing data but it does not seems to be the case. Please share your experience or if you have better method to extract reads. I am going to try Sambamba which some folks suggest to be much faster than samtools.

Thanks, Arun

bam tcga ngs exome • 1.0k views
ADD COMMENTlink modified 2.9 years ago • written 2.9 years ago by akhattri0
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: 1717 users visited in the last hour