Compare 4 vcf files with vcftools
0
0
Entering edit mode
8.4 years ago

mmm I tried with vcftools and --diff to see which snps, indels... share 4 individuals from a family. But with this code seems just compare two individuals

vcftools --vcf 2.vcf --diff  2.vcf --diff 3.vcf --diff 4.vcf --out compareall.xls

so at end I did like someone posted in this forum

cat 2.vcf 4.vcf 10.vcf > 2410.vcf

vcftools --vcf 2410.vcf --recode --recode-INFO-all --diff padre.vcf --out compareall.xls

I do right?

Isn't another way to do like on first stage?

Thanks :)

vcf • 4.4k views
ADD COMMENT
1
Entering edit mode

see if this tool will work.

ADD REPLY
0
Entering edit mode

well...I dont know but always I get some trouble to insatll Python packages

I do like say

pip install https://github.com/Clinical-Genomics/taboo/zipball/master

I sort my samples, but when I run this code in terminal where I have my files

taboo compare sample*.sorted.vcf > results.txt

says taboo command not found ;_;

ADD REPLY
0
Entering edit mode

well no idea how...but I managed to install taboo

now i get this

taboo compare sample2.sorted.vcf sample10.sorted.vcf > results.txt
taboo compare sample*.sorted.vcf > results.txt


both codes Aborted!
ADD REPLY

Login before adding your answer.

Traffic: 2387 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6