-
Notifications
You must be signed in to change notification settings - Fork 260
Description
Greetings,
The bcftools consensus -S setting is not working for me.
I am running:
samtools faidx reference.fa -r regions.txt | bcftools consensus samples.vcf.gz -S subsamples.txt -o consensus.fasta
my samples.vcf.gz file has 347 samples
my subsamples.txt has 2 samples formatted as follows:
sample1
sample2
However the stream says, "Note: applying IUPAC codes based on FORMAT/GT in 347 samples
Applied 37 variants" instead of two samples
If I write out the samples as in
samtools faidx reference.fa -r regions.txt | bcftools consensus samples.vcf.gz -s sample1,sample2 -o consensus2.fasta
the stream correctly says, "Note: applying IUPAC codes based on FORMAT/GT in 2 samples
Applied 30 variants."
Am I missing something about using -S or the formatting of the subsamples.txt? Or is this an issue with bcftools?
Thanks!