Question: Checking bam quality
0
gravatar for banerjeeshayantan
13 months ago by
banerjeeshayantan100 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 • 1.3k views
ADD COMMENTlink modified 13 months ago by venu6.1k • written 13 months ago by banerjeeshayantan100
0
gravatar for venu
13 months ago by
venu6.1k
Germany
venu6.1k 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 13 months ago • written 13 months ago by venu6.1k

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 13 months ago by banerjeeshayantan100
1

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

ADD REPLYlink written 13 months ago by venu6.1k

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 13 months ago by banerjeeshayantan100

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 13 months ago by finswimmer11k

Per base sequence quality

ADD REPLYlink written 13 months ago by banerjeeshayantan100
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: 1226 users visited in the last hour