Biostar Beta. Not for public use.
Question: To Avoid The 'Argument "F" Isn'T Numeric In Numeric
0
Entering edit mode

Hello, I am running the command like:

samtools mpileup -uf hg18.fa my.bam | bcftools view -cg - | vcfutils.pl vcf2fq > cns.fq

However I got many many warnings on the screen which are like

Argument "chr1^I22217285^I.^IT^I.^I28.2^I.^IDP=1;;CI95=1.5,0;FQ=-3..." isn't numeric in numeric lt (<) at /mypath/samtools/bcftools/vcfutils.pl line 503, <> line 21409998.
Argument "chr1^I22250317^I.^IT^I.^I28.2^I.^IDP=9;;CI95=1.5,0;FQ=-3..." isn't numeric in numeric lt (<) at /Mypath/samtools/bcftools/vcfutils.pl line 503, <> line 21443057.

What is wrong and how to avoid it? My reference fa file likes

>chr10

NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN

Thanks.

ADD COMMENTlink 8.3 years ago Zhshqzyc • 490 • updated 15 months ago Biostar 20
Entering edit mode
0

It works on a BAM file for me - how did you make your BAM file. ^I is often a tab - or "t" in perl.

ADD REPLYlink 8.3 years ago
Chris Penkett
• 480
1
Entering edit mode

Come on. You have to try a little. If you are having trouble with a command that pipes a whole bunch of commands together, you have to split them apart, to figure out which command is the problem. Given that the problem looks like mpileup is balking at your .bam file, you really ought to convert some of it to .sam, and put up the first few offending lines of the .sam.

ADD COMMENTlink 8.2 years ago Swbarnes2 ♦ 1.5k
0
Entering edit mode

(not tested) samtools mpileup is missing the option -g ("generate BCF output ") ?

ADD COMMENTlink 8.3 years ago Pierre Lindenbaum 120k
Entering edit mode
0

I tried the command samtools mpileup -uf hg18.fa -g my.bam | bcftools view -cg - | vcfutils.pl vcf2fq > cns.fq But it still failed. Can you list your command?

ADD REPLYlink 8.3 years ago
Zhshqzyc
• 490

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0