SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Why MAQ consensus seq better than SAMtools consensus ?? av_d Genomic Resequencing 5 10-18-2015 04:44 AM
Consensus from mpileup for haploid sequences (forcing base calls - no ambiguities) ericarcher Bioinformatics 6 01-24-2014 09:52 AM
Generate consensus sequence with Samtools mpileup? Heisman Bioinformatics 4 10-28-2013 05:57 PM
mpileup consensus behaviour rwness Bioinformatics 5 07-24-2012 08:12 AM
Problem calling consensus sequence with Samtools mpileup? Lspoor Bioinformatics 3 10-24-2011 08:43 AM

Reply
 
Thread Tools
Old 07-14-2011, 09:07 AM   #1
prm36
Member
 
Location: Edinburgh

Join Date: Mar 2010
Posts: 16
Default consensus from mpileup

Hi,

I am trying to call the consensus sequence from an mpileup file using samtools (specifically the vcf2fq command).

The consensus fastq file is produced, and looks correct, but during the running of vcf2fast I receive the following error(?) numerous times:


Argument "MRSA252^I2001121^I.^IA^I.^I28.2^I.^IDP=1;;CI95=1.5,0;FQ=..." isn't numeric in numeric lt (<) at /samtools-0.1.16/bcftools/vcfutils.pl line 503, <> line 1958656.

Could anyone explain to me what this means, and whether I am doing something wrong?!

Cheers

Paul
prm36 is offline   Reply With Quote
Old 07-14-2011, 09:08 AM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

I don't think anything is wrong. It does that at every point where the DP is 1 or 0. The software will fill in N's at those points. Just let it finish, the look at it, it should be alright.
swbarnes2 is offline   Reply With Quote
Old 07-14-2011, 09:10 AM   #3
prm36
Member
 
Location: Edinburgh

Join Date: Mar 2010
Posts: 16
Default

Thank you!
prm36 is offline   Reply With Quote
Old 07-14-2011, 09:11 AM   #4
prm36
Member
 
Location: Edinburgh

Join Date: Mar 2010
Posts: 16
Default

Thank you!!
prm36 is offline   Reply With Quote
Reply

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:51 AM.


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