SURVIVOR merge function usage
0
0
Entering edit mode
3.2 years ago
dlekrud456 • 0

Hello,

I have issues with SURVIVOR merge function that only shows coincidental cases when I draw VenDiagram. Can anyone please help me with this issue?

used survivor merge and genComp function to compare two vcfs but it only shows coincidental cases. The below is my code.

In bash,

SURIVOR merge sample_files 1000 2 1 1 0 30 sample_merged.vcf
genComp sample_merged.vcf sample_merged.mat.txt
perl -ne 'print "$1\n" if /SUPP_VEC=([^,;]+)/' sample_merged.vcf | sed -e 's/\(.\)/\1 /g' > sample_merged_overlapp.txt

In R,

t=read.table("sample_merged_overlapp.txt",header=F)

library(VennDiagram)

venn.diagram(list(N15=which(t[,1]==1), T15=which(t[,2]==1)), fill = c("gray", "orange") , alpha = c(0.5, 0.5), cex = 2, lty =2, filename = "my_sample")

And the result shows following,

enter image description here

I'm new to programming and not sure how to solve this issue, your help is very much appreciated!

NGS genome genomics • 1.1k views
ADD COMMENT

Login before adding your answer.

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