metrics on BAM file

Moip · 09-26-2018, 11:13 AM

Hi everyone,
I am doing exome sequencing, paired end, but not used to analyze reads.
After alignement, how do i proceed the BAM file?
Thank you

Richard Finney · 09-26-2018, 05:19 PM

How many samples do you have?
How did you align the reads?
Are they tumor and normal samples?
Are looking for germline characteristics or somatic changes?

Moip · 09-27-2018, 12:54 AM

I have three normal samples and looking for germline characteristics.
Mapping with BWA-MEM.
Thank you for you help Richard

Richard Finney · 09-27-2018, 07:20 AM

You will probably want to do some variant calling on the bam files. This may include SNPS, fusions, amplification or deletions.

There are many tools to do this. Searching for these tools is straightforward:

https://www.google.com/search?q=germ...g+on+bam+files

You may want to use favorite, known tools already used at your organization, so ask around. A new tool may require some non-trivial installation.

Moip · 09-27-2018, 07:48 AM

Thank you for you reply.
i need to understand some details, once i have a BAM file, before variant calling i have to remove duplicate, recalibrate BAM? run other metrics ?
sorry i am completely new to this.

Richard Finney · 09-27-2018, 05:01 PM

You could do some quality checks and remove duplicates and recalibrate; but if you are just starting off, just go ahead and assume the data is good and try and see if you can get a report on variations using a commonly used tool. After you get that working, you can check out the other stuff and re-run a new workflow. Order of importance would be : QC, duplicate removal or marking , then re calibration. Many here at work say recalibration is not worth the effort. Theoretically marking or removing duplicates may be important upstream to a variant caller, and again, theoretically a variant caller could do that work for you so read the docs for your tools and follow their advice.

Do try and run a QC tool on your fastqs early.
Absolute easiest QC on bam files is "samtools flagstat" (google it).
Study the results and see if you can make sense of them.

Moip · 09-28-2018, 12:12 AM

Thank you, i ll give a try. Do i ignore UMi for this analysis?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Cuffdiff error [bam_header_read] invalid BAM binary header (this is not a BAM file).	adrian	Bioinformatics	0	12-29-2016 02:23 PM
Write the subset of reads from BAM file into new SAM/BAM file, using R tools.	Old Pioneer	Bioinformatics	0	01-27-2016 05:41 AM
Convert merged BAM back to per lane BAM or FASTQ file	danielsbrewer	Bioinformatics	6	10-03-2013 08:29 AM
Problems with Picard metrics file	cristinavalente	Bioinformatics	0	05-11-2013 03:19 AM

09-26-2018, 11:13 AM	#1
Moip Junior Member Location: canada Join Date: Sep 2018 Posts: 5	metrics on BAM file Hi everyone, I am doing exome sequencing, paired end, but not used to analyze reads. After alignement, how do i proceed the BAM file? Thank you

09-28-2018, 12:12 AM	#7
Moip Junior Member Location: canada Join Date: Sep 2018 Posts: 5	Thank you, i ll give a try. Do i ignore UMi for this analysis? Last edited by Moip; 09-28-2018 at 11:34 PM.

09-26-2018, 05:19 PM	#2
Richard Finney Senior Member Location: bethesda Join Date: Feb 2009 Posts: 700	How many samples do you have? How did you align the reads? Are they tumor and normal samples? Are looking for germline characteristics or somatic changes?

09-27-2018, 12:54 AM	#3
Moip Junior Member Location: canada Join Date: Sep 2018 Posts: 5	I have three normal samples and looking for germline characteristics. Mapping with BWA-MEM. Thank you for you help Richard

09-27-2018, 07:20 AM	#4
Richard Finney Senior Member Location: bethesda Join Date: Feb 2009 Posts: 700	You will probably want to do some variant calling on the bam files. This may include SNPS, fusions, amplification or deletions. There are many tools to do this. Searching for these tools is straightforward: https://www.google.com/search?q=germ...g+on+bam+files You may want to use favorite, known tools already used at your organization, so ask around. A new tool may require some non-trivial installation.

09-27-2018, 07:48 AM	#5
Moip Junior Member Location: canada Join Date: Sep 2018 Posts: 5	Thank you for you reply. i need to understand some details, once i have a BAM file, before variant calling i have to remove duplicate, recalibrate BAM? run other metrics ? sorry i am completely new to this.

09-27-2018, 05:01 PM	#6
Richard Finney Senior Member Location: bethesda Join Date: Feb 2009 Posts: 700	You could do some quality checks and remove duplicates and recalibrate; but if you are just starting off, just go ahead and assume the data is good and try and see if you can get a report on variations using a commonly used tool. After you get that working, you can check out the other stuff and re-run a new workflow. Order of importance would be : QC, duplicate removal or marking , then re calibration. Many here at work say recalibration is not worth the effort. Theoretically marking or removing duplicates may be important upstream to a variant caller, and again, theoretically a variant caller could do that work for you so read the docs for your tools and follow their advice. Do try and run a QC tool on your fastqs early. Absolute easiest QC on bam files is "samtools flagstat" (google it). Study the results and see if you can make sense of them.