R代碼：

args<-commandArgs(trailingOnly=TRUE)
data<-read.table(args[1],header=F)
library(VennDiagram)
data<-read.table("*.txt",header=T)
grid.newpage()
png("*.png")
for(i in 1:6)
{
vp <- viewport(x = 0.15*i, y = 0.8, width = 0.14, height = 0.18, angle=0)
pushViewport(vp)
grid.rect()
draw.pairwise.venn(area1=data[1,i],area2=data[5,i],cross.area=data[3,i],col='white',fill=c('deepskyblue','brown1'),cat.col=c('deepskyblue','brown1'),category=c(names(data)[i],'GATK_v3.5'),cat.pos = c(0, 180),euler.d = TRUE,sep.dist = 0.03)
upViewport()
vp <- viewport(x = 0.15*i, y = 0.6, width = 0.14, height = 0.18, angle=0)
pushViewport(vp)
grid.rect()
draw.pairwise.venn(area1=data[6,i],area2=data[10,i],cross.area=data[8,i],col='white',fill=c('deepskyblue','brown1'),cat.col=c('deepskyblue','brown1'),category=c(names(data)[i],'GATK_v3.5'),cat.pos = c(0, 180),euler.d = TRUE,sep.dist = 0.03)
upViewport()
vp <- viewport(x = 0.15*i, y = 0.4, width = 0.14, height = 0.18, angle=0)
pushViewport(vp)
grid.rect()
draw.pairwise.venn(area1=data[11,i],area2=data[15,i],cross.area=data[13,i],col='white',fill=c('deepskyblue','brown1'),cat.col=c('deepskyblue','brown1'),category=c(names(data)[i],'GATK_v3.5'),cat.pos = c(0, 180),euler.d = TRUE,sep.dist = 0.03)
upViewport()
vp <- viewport(x = 0.15*i, y = 0.2, width = 0.14, height = 0.18, angle=0)
pushViewport(vp)
grid.rect()
draw.pairwise.venn(area1=data[16,i],area2=data[20,i],cross.area=data[18,i],col='white',fill=c('deepskyblue','brown1'),cat.col=c('deepskyblue','brown1'),category=c(names(data)[i],'GATK_v3.5'),cat.pos = c(0, 180),euler.d = TRUE,sep.dist = 0.03)
upViewport()
}
vp <- viewport(x = 0.5, y = 0.92, width = 0.115, height = 0.18, angle=0)
pushViewport(vp)
grid.text("BWA")
upViewport()
vp <- viewport(x = 0.05, y = 0.8, width = 0.115, height = 0.18, angle=0)
pushViewport(vp)
grid.text("C1")
upViewport()
vp <- viewport(x = 0.05, y = 0.6, width = 0.115, height = 0.18, angle=0)
pushViewport(vp)
grid.text("C2")
upViewport()
vp <- viewport(x = 0.05, y = 0.4, width = 0.115, height = 0.18, angle=0)
pushViewport(vp)
grid.text("T1")
upViewport()
vp <- viewport(x = 0.05, y = 0.2, width = 0.115, height = 0.18, angle=0)
pushViewport(vp)
grid.text("T2")
upViewport()
dev.off()

數(shù)據(jù)：

結(jié)果圖：