samtools mpileup: bam file to consensus sequence generation error?
0
0
Entering edit mode
5.3 years ago
Kumar ▴ 120

I have been doing reference based assembly of illumina based paired end reads. For which, I indexed my reference genome and generated sam, bam and sorted.bam files, further proceeded for consensus sequence generation using the following command,

samtools mpileup -uf reference.fasta query.sorted.bam | bcftools view -cg - | perl vcfutils.pl vcf2fq | seqtk seq -A -l 70 - | sed -r "s/(>)(.*)/\1 query.bam.consensus/g" > consensus.fasta

and end up with error as follows,

Error: Could not parse --min-ac g
Use of uninitialized value $l in numeric lt (<) at vcfutils.pl line 566.
Use of uninitialized value $l in numeric lt (<) at vcfutils.pl line 566.
[mpileup] 1 samples in 1 input files
<mpileup> Set max per-file depth to 8000

The strange thing is that, the same command line is working well (without showing error) in ubuntu 14 and 16 LTS OS version, but showing above mentioned error in ubuntu 18 LTS OS version. Please help me to fix this issue.

Assembly genome ngs samtools bcftools • 1.6k views
ADD COMMENT
0
Entering edit mode

There is an issue on github related to his https://github.com/samtools/bcftools/issues/391

ADD REPLY

Login before adding your answer.

Traffic: 1647 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