vcf to venn diagram

Hi, I tried to reproduce this article but this command failed:

$ bcftools isec -n +2 calls-16c-gt-q20.vcf  calls-SA-gt-q20.vcf > isec_file1-v-2_out.vcf
[main] Unrecognized command.

I am using samtools 1.9.

What did I miss?

Thank you in advance,

snp

• 1,252 views • link

written 7.6 years ago by Ric

Can you try running which bcftools followed by just bcftools and paste both outputs here? This might be a version mismatch.

Also, just to eliminate a possibility, you're not typing the $ as part of the command, right?

• 0 views • link

written 7.6 years ago by Ram

0 answers

No answers yet.

More posts like this

gatk LiftoverVCF finishes with no output VCF nor error.
written 3.6 years ago by jpablo •

Hi, I'm trying to liftover a VCF from the ALFA project in order to get the frequencies in HG19 instead of HG38, but LiftoverVCF (latest …
How to handle VCFs from the same sample but using different aligners and variant callers?
written 4.7 years ago by aldhairmedico •

Hi, I'm using whole-exome sequencing (WES) for somatic variant calling. During the process, I tried to follow the approach described here: https://pubmed.ncbi.nlm.nih.gov/28420412/ **Basically my workflow …
Bcftools isec output only contains one sample name - how to identify which samples share the same v…
written 5.1 years ago by Linda •

I have read this post: https://www.biostars.org/p/136937 to further understand the output generated by bcftools isec. I'm trying to find variants that are common by location …
GATK4 ASEReadCounter returns nothing
written 6.6 years ago by alfonsosaera •

Dear BioStar Users, I have a problem with `GAT4 ASEReadCounter` and could not find help googling or in GATK4 forum. My code: gatk ASEReadCounter -I …
bcftools filtering INFO column issues (empty output)
written 7.7 years ago by Molly_K

How can I filter out variants with R2!="." in the INFO column? > bcftools view -i 'INFO/R2!="."' -Oz $FILE -o $OUTFILTERED The above generated a …
Tutorial: Produce PCA bi-plot for 1000 Genomes Phase III - Version 2
written 7.9 years ago by Kevin Blighe

**Note1 - Previous version: https://www.biostars.org/p/271694/** **Note2 - this data is for hg19 / GRCh37** **Note3 - GRCh38 data is available [HERE][1]** ------------------------------- The tutorial has …
BAN - Best alignment normalisation error
written 8.5 years ago by pinn

Hi While performing vcf normalisation I go this error using BAN software. Is it mummer showing the error ? Is their any better vcf normalisation …
Extract Allele At Particular Positions From Bam Files
written 12.4 years ago by guns.guru •

<p>Is there an easy way to extract the allele (monomorphic or polymorphic wrt the reference) for a sample across all reads at specific locations, given …
Tool: Gemini: Integrative Exploration Of Genetic Variation And Genome Annotations
written 13.3 years ago by Aaronquinlan

We have developed a new tool called GEMINI (GEnome MINIng) to facilitate the exploration of genetic variation in the context of a wide range of …
Parse Blast Output ..
written 14.4 years ago by Suraj •

<p>i want to parse this blast output but compilers giving NULLPOINTEREXCEPTION .. input file and source code of java program is given below ... plz …