Question: Checking bam quality
0
gravatar for banerjeeshayantan
2.8 years ago by
banerjeeshayantan190 wrote:

I have multiple bam files from cancer datasets and I want to check their quality. What I did was, I converted the bam to fastq using bedtools. I got two fastq read files (paired end) per bam file. Now I checked the quality of the fastq files generated using fastqc tool. This is how I am trying the find a good dataset for my analysis. Is this approach ok?

alignment sequence • 3.5k views
ADD COMMENTlink modified 2.8 years ago by venu6.7k • written 2.8 years ago by banerjeeshayantan190
0
gravatar for venu
2.8 years ago by
venu6.7k
Germany
venu6.7k wrote:

No. It's not ok.

You are checking the quality of fastq file but what you want is BAM quality. There is a difference between two. Check Qualimap to assess the quality of aligned BAM.

If it is ChIP-seq data, explore deepTools. It has very nice functions.

ADD COMMENTlink modified 2.8 years ago • written 2.8 years ago by venu6.7k

Thanks for your reply. I get your point. Let me rephrase. So now I have only bam files and I want to check the fastq quality of the files that produced these bam files, then is the above method ok?

ADD REPLYlink written 2.8 years ago by banerjeeshayantan190
1

You can directly run fastqc on bam files to get fastq quality. Skip time consuming BAM -> fastq.

ADD REPLYlink written 2.8 years ago by venu6.7k

But my reads are paired end. So to check quality of both reads, I need to make read 1 and read 2 for each bam file and then check it using fastqc tool. Isn't it?

ADD REPLYlink written 2.8 years ago by banerjeeshayantan190

So to check quality of both reads

Please define what you mean by "quality". There are so many things one can look for and it maybe depends on the goal of your analysis.

fin swimmer

ADD REPLYlink written 2.8 years ago by finswimmer14k

Per base sequence quality

ADD REPLYlink written 2.8 years ago by banerjeeshayantan190
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: 1343 users visited in the last hour