SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
best-practices for assigning read-groups for snp-caller brentp Bioinformatics 2 12-20-2011 07:34 AM
GATK unified SNP caller error dakl Bioinformatics 4 11-23-2010 09:27 PM
Aligner/SNP caller for Illumina data sklages Bioinformatics 2 10-06-2010 04:48 AM
SNP Caller for BAM Files AnamikaDarwin Bioinformatics 0 01-22-2010 05:53 AM
SNP caller for RNA-Seq data shurjo Bioinformatics 0 12-08-2009 02:58 PM

Reply
 
Thread Tools
Old 10-20-2010, 06:19 PM   #1
BertieWooster
Member
 
Location: Sydney, Australia

Join Date: Jul 2009
Posts: 14
Question Samtools SNP Caller versus MAQ SNP caller

I have been using MAQ's variant caller to call SNPs. I now want to use Samtools.pl varfilter and want to make sure that I maintain the same conditions I had for MAQ. My question is which of the parameters are common to MAQ and Samtools' variant caller?


Thanks,
BW
BertieWooster is offline   Reply With Quote
Old 10-21-2010, 05:25 AM   #2
drio
Senior Member
 
Location: 4117'49"N / 24'42"E

Join Date: Oct 2008
Posts: 323
Default

Take a look to this .
__________________
-drd
drio is offline   Reply With Quote
Old 05-18-2011, 07:42 AM   #3
egatti
Junior Member
 
Location: Milan

Join Date: May 2011
Posts: 3
Default MAQ SNP caller

Dear all,

In need to know if MAQ SNP caller takes .bam files in input.

I have aligned reads by Novoalign and I want to use MAQ only as a SNP caller, is it possible or the alignment with MAQ is mandatory?

Waiting for your feedback,
Thanks a lot

Elena
egatti is offline   Reply With Quote
Old 05-18-2011, 09:00 AM   #4
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

No, MAQ does not take BAM files as input. Use "samtools pileup -A", however "pileup" will be removed from samtools in the next release, so you might as well move over to "samtools mpileup".
nilshomer is offline   Reply With Quote
Old 05-18-2011, 11:19 PM   #5
egatti
Junior Member
 
Location: Milan

Join Date: May 2011
Posts: 3
Default

Thanks a lot,
I used samtools pileup and I had no problem, but I would like to try with an other SNP caller, for example MAQ.
So, as far as I understood, It is not possible to call variants with MAQ starting with reads aligned by Novoalign. If I want to call variants with MAQ I had to re-align my reads with MAQ. Is it correct?

Thank you in advance for your suggestions.

Elena
egatti is offline   Reply With Quote
Old 05-18-2012, 01:41 PM   #6
bgulko
Member
 
Location: New York

Join Date: Nov 2009
Posts: 13
Default

This may be a dead thread, but just in case, I believe the old version of samtools pileup (the new version is samtools mpileup) can be forced to produce the MAQ consensus call (it requires a special option, -c). Also, if you are interested in another caller, the pileup can be used as input to the BSNP genotype caller at http://compgen.bscb.cornell.edu/GPhoCS/BSNP/ .
--Brad
bgulko is offline   Reply With Quote
Reply

Tags
samtools.pl varfilger

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 07:15 AM.


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