Missing -d option in new BCFTools?

In the old version of BCFTools (0.1.19), there was an option "-d skip loci where less than FLOAT fraction of samples covered [0]"

It seems to have disappeared in the new version (1.1). Is there an alternative that does the same thing?

Thanks! :-)

sequence snp

• 1,620 views • link

updated 11.6 years ago by Ram

• written 11.6 years ago by Tails

1 answer

From a quick look at the manual and your question, I'm guessing you could use the -e option to exclude sites where number of samples (N_SAMPLES) is less than fraction * total number of samples (a calculated constant)

• 0 views • link

written 11.6 years ago by Ram

More posts like this

bcftools mpileup returns vcf files with only headers
written 1.5 years ago by Brenden •

I've been working on a project the past few weeks where I'm analyzing SAM files for specific point mutations. I'm aware that bcftools has the …
bcftools consensus still returns "Could not parse the header" error
written 5.0 years ago by shpak.max •

I attempted to create a consensus fasta file using bcftools, i.e. bgzip -c All_SRR_SNP_Clean.vcf > All_SRR_SNP_Clean.vcf.gz tabix All_SRR_SNP_Clean.vcf.gz cat $ref| bcftools consensus $vcf_dir/All_SRR_SNP_Clean.vcf.gz > consensus.fasta …
bcftools call variant calling model
written 5.6 years ago by igor

This seems like a question that should have been addressed before, but somehow I cannot find anything. There is a [post][1] on samtools SourceForge page …
Does STAR identify new splicing?
written 9.0 years ago by zhangzhenyusz •

At beginning I was using Tophat-Cufflinks and had processed several sets of data. But it was way too slow so I switched to STAR. However, …
SNP calls all seem to have quality of zero
written 10.1 years ago by novice

I used the following pipeline to call SNPs on an alignment file: $ samtools mpileup -uf s_cerevisiae.fasta aln.sorted.bam | bcftools view -O b - > …
samtools/bcftools strange (error ?) lack of call
written 12.0 years ago by mbourgey •

Hi all, I got a strange lack of call when I used samtolls mpileup + bcftools to call my SNP (version 0.19 or 1.0 give …
A Quick Way To Annotate Entire Bacterial Genomes To Go?
written 12.5 years ago by mgalactus

Is there any fast and reliable way to annotate protein sequences from bacterial genomes to the GO terms? The best option so far seems to …
When Would You Use Sage Over Rna-Seq
written 12.7 years ago by Ying W

Is there ever a situation where you would use SAGE (or related protocols) over RNA-seq in organism where all transcripts are known? The only situation …
Meaning Of -A In Vcftools Vcfutils.Pl Varfilter?
written 13.4 years ago by Dan

The command line says: <pre><code>Usage: vcfutils.pl varFilter [options] <in.vcf> Options: -Q INT minimum RMS mapping quality for SNPs [10] -d INT minimum read depth [2] …
Samtools Mpileup And Bcftools Parameter
written 13.8 years ago by Ashutosh Pandey

Hi, I am currently using new samtools 0.1.17. In older versions of samtools there was a -r parameter explained below: <pre><code>-r FLOAT prior of a …