This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SnpSift match overlapping mutations with a given bed file

Hi,

I know how to get overlapping mutations in my vcf with a given bed file. But is there a way to get a simple number from a log file ? like in the vcftools ?

Best,

Tunc.

snpsift vcftools
• 2,323 views • link
written by morovatunc morovatunc

can you give example of what exactly you want

• 1 views • link
written by Medhat Medhat

On one hand I have bed file. On the other hand I have my vcf file.

I want to extract mutations on that vcf file within that bed file locations as a VCF + I want to know how many of them are there.

Thank you for your response.

Tunc.

• 1 views • link
written by morovatunc morovatunc

Did you tried bedtools intersect -a yourVcfFile.vcf -b yourBed.bed > result.vcf

and for how many "number of lines" wc -l result.vcf

If I understood right your question.

• 1 views • link
written by Medhat Medhat

With this method I got that number that I want, but I will try to keep it as a last option because snpsift itself can do that isolation as well, I just wanted to know if it has a option like vcftools's log file.

Thank you very much for your useful solution.

• 1 views • link
written by morovatunc morovatunc

I think yes, you will need to use one of those according to your needs

  • countHom() Count number of homozygous genotypes No arguments (countHom() > 0)
  • countHet() Count number of heterozygous genotypes No arguments (countHet() > 2)
  • countVariant() Count number of genotypes that are variants (i.e. not reference 0/0) No arguments (countVariant() > 5)
  • countRef() Count number of genotypes that are NOT variants (i.e. reference 0/0)
• 1 views • link
written by Medhat Medhat

0 answers

No answers yet.

Log in to answer this question.

More posts like this