Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Yes .. BBMap can do that! GenoMax Bioinformatics 211 03-21-2018 04:59 PM
Several questions regarding BBMap/BBSplit MSchm Bioinformatics 2 12-18-2017 10:01 AM
2 pileup format questions krobison Bioinformatics 1 11-03-2010 10:52 AM
Samtools pileup SNPs questions Melissa Bioinformatics 2 03-24-2010 10:15 AM

Thread Tools
Old 12-20-2017, 08:21 AM   #1
Junior Member
Location: Sweden

Join Date: May 2016
Posts: 7
Smile BBMap pileup questions

I am interested in metrics that summarise coverage from my BAM file. I came across BBMap pileup which does a decent job. But I have a couple of questions.

1. It seems like if a window is provided to pileup, it returns the sum over the given window. Can other metrics such as sum/mean/median/min/max/sd/mad be calculated? I am particularly interested in the median.

2. Does pileup do only fixed windows or can it do sliding window metrics?

3. BBMap pileup returns a summary over the whole file to the console like below:

Average coverage: 71.19
Standard deviation: 198.58
Percent scaffolds with any coverage: 100.00
Percent of reference bases covered: 97.66

Can this be saved to a file? I tried '> summary.txt' at the end of my script but that didn't work.

rmf is offline   Reply With Quote
Old 12-20-2017, 10:10 AM   #2
Senior Member
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,731

Have you tried using "out=" to save the results to a file.
GenoMax is offline   Reply With Quote

bbmap, coverage, mapping coverage, pileup

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 10:07 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO