To Avoid The 'Argument "F" Isn'T Numeric In Numeric
2
0
Entering edit mode
9.4 years ago
Zhshqzyc ▴ 490

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.

samtools perl • 3.5k views
ADD COMMENT
0
Entering edit mode

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 REPLY
1
Entering edit mode
9.4 years ago
Swbarnes2 ★ 1.5k

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 COMMENT
0
Entering edit mode
9.4 years ago

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

ADD COMMENT
0
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 2067 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6