Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Variant Calling with mpileup asebastian Bioinformatics 0 03-31-2013 08:53 PM
Resequencing & variant calling issue Jluis Bioinformatics 5 02-05-2013 02:12 AM
BFAST and Variant Calling nexgengirl Bioinformatics 8 01-02-2013 05:03 AM
variant calling with SAMtools krueml Bioinformatics 2 09-28-2012 05:36 AM
variant calling kjaja Bioinformatics 1 11-04-2011 07:16 AM

Thread Tools
Old 04-01-2013, 12:07 PM   #1
Junior Member
Location: USA

Join Date: Jul 2012
Posts: 4
Default BioPerl & Variant calling

Hi ,

I am trying to do variant calling with mpileup using BioPerl.
I am using Bio::Tools::Run::Samtools but it says mpileup not registered.
Does this mean BioPerl doesn't support mpileup ?
Can anyone tell me how to run it(or any other variant calling program) through a BioPerl module ?

Thank you
asebastian is offline   Reply With Quote
Old 04-02-2013, 02:37 AM   #2
Location: Oxford, UK

Join Date: Nov 2011
Posts: 17


After looking at CPAN and then at the module code directly it looks as if mpileup is not supported as there is no '-command' option in the code.

This may be because the module was written for an old version of samtools and so mpileup was not implemented.

you could email them [BioPerl] and see if they will write in support as I am guessing it would be useful for other people and relatively simple to do. Otherwise have you considered using system calls in perl to run mpileup?

system("samtools mpileup ......");
davidblaney is offline   Reply With Quote
Old 04-02-2013, 12:57 PM   #3
Junior Member
Location: USA

Join Date: Jul 2012
Posts: 4

Thanks David,
With BioPerl, I am trying to contact them...
As for the system calls, I am trying to avoid direct system commands... Thats why I was looking for a BioPerl wrapper...

asebastian is offline   Reply With Quote

bioinformatics, bioperl, mpileup, variant calling

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 02:23 PM.

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