vcf-compare and making Venn diagram
0
0
Entering edit mode
22 months ago
arshad1292 ▴ 70

Hello,

I want to compare two vcf files using vcf tools, and draw the Venn diagrams to visualize the comparison.

I ran the following command.

vcf-compare file1_sniffles_filtered.vcf file2_Chr.vcf --ignore-indels -r -p

Here is the output:

VN 'Venn-Diagram Numbers'. Use grep ^VN | cut -f 2- to extract this part.

VN The columns are:

VN 1 .. number of sites unique to this particular combination of files

VN 2- .. combination of files and space-separated number, a fraction of sites in the file

VN 1 file2_Chr.vcf (0.0%) file1_sniffles_filtered.vcf (0.0%) VN 2433 file1_sniffles_filtered.vcf (100.0%) VN 708961 file2_Chr.vcf (100.0%)

SN Summary Numbers. Use grep ^SN | cut -f 2- to extract this part.

SN Number of REF matches: 0 SN Number of ALT matches: 0 SN Number of REF mismatches: 1 SN Number of ALT mismatches: 0 SN Number of samples in GT comparison: 0

Number of sites lost due to grouping (e.g. duplicate sites): lost, %lost, read, reported, file

SN Number of lost sites: 268 9.9% 2702 2434 file1_sniffles_filtered.vcf SN Number of lost sites: 449 0.1% 709411 708962 file2_Chr.vcf

Question 1: Is that the right output that I am supposed to get? Question 2: How can I create Venn diagram using these results?

I tried looking up online and found some 5-6 years old thread but that method doesn't work anymore, I guess.

Please help this newbie. Many thanks in advance.

SNP • 493 views
ADD COMMENT

Login before adding your answer.

Traffic: 1733 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