Hi,

the commands are obsolete, please see here instructions for the recommended workflow:
http://samtools.github.io/bcftools/howtos/variant-calling.html

Best wishes,
Petr

On 05/07/2020 20:47, White, Noor (NIH/NEI) [F] via Samtools-help wrote:

Hello,

I am repeatedly getting the below error running samtools (v1.10) and python (v3.7, though I get the same error when I re-run with python v2.7). I do not believe the error has to do with my SAM-format output, original sample or query reference, as I have switched all out. I found some similar error messages online, but no resolution. Thanks in advance for any possible help.

Command:

samtools mpileup -Q 20 -q 5 -d 5000 -uf Ggal_exon_sequences.fasta Mapped_AG_WGS_Ggal_exon.bam | bcftools call -c | vcfutils.pl vcf2fq -d 5 -Q 20 -l 1 >Consensus_AG_WGS_Ggal_exon.fq

Error:

Use of uninitialized value $l in numeric lt (<) at /usr/local/apps/samtools/1.10/bin/vcfutils.pl line 566, <> line 412.

--

Noor D. White, PhD

/Postdoctoral Fellow/, National Science Foundation

/Postdoctoral Fellow/, National Eye Institute, National Institutes of Health

/Research Collaborator/, National Museum of Natural History, Smithsonian Institution

NNRL [nei.nih.gov] <https://urldefense.proofpoint.com/v2/url?u=https-3A__nei.nih.gov_research_research-2Dlabs-2Dand-2Dbranches_neurobiology-2Dneurodegeneration-2Drepair-2Dlaboratory&d=DwMGaQ&c=D7ByGjS34AllFgecYw0iC6Zq7qlm8uclZFI0SqQnqBo&r=xdvdTaAZDWitAtUqWIZL0A&m=qgRWAvs_jKP4lIS9VzBJAgoJHo08by44S-lL86l66fA&s=JndXwLAWFmql0Kr_yBqD4a2EYboL7xhaXEvE9qNG6bU&e=>, National Eye Institute [nei.nih.gov] <https://urldefense.proofpoint.com/v2/url?u=https-3A__nei.nih.gov_&d=DwMGaQ&c=D7ByGjS34AllFgecYw0iC6Zq7qlm8uclZFI0SqQnqBo&r=xdvdTaAZDWitAtUqWIZL0A&m=qgRWAvs_jKP4lIS9VzBJAgoJHo08by44S-lL86l66fA&s=CiIsXtp7tv5VQKf9ggJXbzoevgsw-dhuNTwbbG7bhYc&e=>, NIH [nih.gov] <https://urldefense.proofpoint.com/v2/url?u=https-3A__www.nih.gov_&d=DwMGaQ&c=D7ByGjS34AllFgecYw0iC6Zq7qlm8uclZFI0SqQnqBo&r=xdvdTaAZDWitAtUqWIZL0A&m=qgRWAvs_jKP4lIS9VzBJAgoJHo08by44S-lL86l66fA&s=nY2wuoJEluVP6RNxh444Z2Vln53lUxkAwVLmjWOax-o&e=>

6 Center Dr.

Bethesda, MD 20892

(301) 402-5734

www.noorwhite.com [noorwhite.com] <https://urldefense.proofpoint.com/v2/url?u=http-3A__www.noorwhite.com_&d=DwMGaQ&c=D7ByGjS34AllFgecYw0iC6Zq7qlm8uclZFI0SqQnqBo&r=xdvdTaAZDWitAtUqWIZL0A&m=qgRWAvs_jKP4lIS9VzBJAgoJHo08by44S-lL86l66fA&s=b1ZilueILbTAVNfGoopYrZHUzNwBqr1sBj-w1idYmCk&e=>



_______________________________________________
Samtools-help mailing list
Samtools-help@lists.sourceforge.net
https://urldefense.proofpoint.com/v2/url?u=https-3A__lists.sourceforge.net_lists_listinfo_samtools-2Dhelp&d=DwICAg&c=D7ByGjS34AllFgecYw0iC6Zq7qlm8uclZFI0SqQnqBo&r=xdvdTaAZDWitAtUqWIZL0A&m=qgRWAvs_jKP4lIS9VzBJAgoJHo08by44S-lL86l66fA&s=6c1DiAcPCS852_j1kszrjTGW87PrNpir4Wp6et8xac8&e=




--
The Wellcome Sanger Institute is operated by Genome Research Limited, a charity registered in England with number 1021457 and a company registered in England with number 2742969, whose registered office is 215 Euston Road, London, NW1 2BE.
_______________________________________________
Samtools-help mailing list
Samtools-help@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/samtools-help

Reply via email to