BCFTOOLS VIEW parameter?
I am trying to split my multiple samples vcf files to single one with Bcftools view, and -O z which can generate the compressed VCF.
bcftools view -s "$id" -g het -v snps -o "$id/${id}_${vcf%.vcf.gz}_hetsnps.vcf" -O z "$vcf" --threads 50
However, the output files are uncompressed?
• 1,281 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
charlie •
•
written
by
charlie •
1 answer
you should add a '.gz' to the output file name.
-o "$id/${id}_${vcf%.vcf.gz}_hetsnps.vcf.gz" -O z
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
after gatk VariantAnnotator -V *_com_norm.vcf -A AlleleFraction -O *_norm_AB.vcf There "nan,nan"…
written by zhuo •
-
Does a merged VCF file (using bcftools merge) contain unique SNPs or INDELs from all merged VCF fi…
written by mohsamir2016 •
-
BCFtools consensus
written by Peerzada •
-
bcftools merge of over 9000+ vcf files
written by Vanish007 •
-
for loop in a bash script
written by Erika •
-
PCA for combined set
written by raalsuwaidi •
-
Only output SNPs in bcftools
written by NewtoBioinfo •
-
non-SNP variants in a vcf merge file
written by evelyn
-
bcftools compressing and indexing vcf files
written by Inquisitive8995
-
bcftools merge, error "Could not parse the region(s)"
written by agathejouet •
Please do not paste screenshots of plain text content, it is counterproductive. You can copy paste the content directly here (using the code formatting option shown below), or use a GitHub Gist if the content volume exceeds allowed length here.