Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools pileup vs mpileup for computing per-base coverage SamH Bioinformatics 5 08-10-2011 03:14 PM
quantitative measure of coverage and reference genome issue sara_ General 1 03-14-2011 08:37 PM
Why are there all these As in my consensus sequence: An issue with samtools pileup dagarfield Bioinformatics 3 02-03-2011 11:14 AM
samtools pileup coverage guavajuice Bioinformatics 4 05-27-2010 03:30 PM
Coverage issue hersh Bioinformatics 0 03-18-2009 11:17 PM

Thread Tools
Old 06-11-2010, 05:46 AM   #1
Senior Member
Location: Germany

Join Date: Oct 2008
Posts: 415
Default Samtools pileup coverage issue

Hi everyone,

I am trying to assess the alignment quality of two aligners.

It seems the samtools pileup command does not give data for every position in the genome.

Specifically, it does not print low coverage regions - coverage <3 from what I can see.

Is there a way to modify this behaviour ?

I see the -l FILE
may be helpful?

I want to include _all_ positions in the pileup.

colindaven is offline   Reply With Quote
Old 06-11-2010, 07:31 AM   #2
Senior Member
Location: Boston

Join Date: Feb 2008
Posts: 693

No. You have to write a script to do this by yourself.
lh3 is offline   Reply With Quote
Old 06-11-2010, 08:49 AM   #3
Location: wisconsin

Join Date: Jan 2010
Posts: 12

there's a good perl API for samtools if you want to go that route. google Bio:B::Samtools or something along those lines and you'll find it.
bbimber is offline   Reply With Quote

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 05:50 PM.

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