-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmethylation workflow clinvar.dot
42 lines (22 loc) · 1.16 KB
/
methylation workflow clinvar.dot
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
digraph workflow {
"ClinVar VCF" [shape=hexagon];
"GRCh38 -> T2Tv1.0 chain file" [shape=hexagon];
"T2Tv1.0 -> T2Tv1.1 chain file" [shape=hexagon];
"GENCODE GFF" [shape=hexagon];
"ClinVar VCF" -> "ClinVar BED" [label="vcf2bed"];
"GENCODE GFF" -> "GENCODE BED" [label="gff2bed"];
"ClinVar GRCh38 regions BED" [shape=parallelogram];
"ClinVar T2Tv1.1 regions BED" [shape=parallelogram];
"CrossMap1" [shape=diamond, color=gray85, style=filled, label="CrossMap.py bed <chain file> <BED file>"];
"CrossMap2" [shape=diamond, color=gray85, style=filled, label="CrossMap.py bed <chain file> <BED file>"];
"GRCh38 -> T2Tv1.0 chain file" -> "CrossMap1";
"bedtools intersect -u" -> "ClinVar GRCh38 regions BED";
"bedtools intersect -u" [shape=diamond, color=gray85, style=filled, label="bedtools intersect -u \n-a <GENCODE> -b <ClinVar>"];
"GENCODE BED" -> "bedtools intersect -u";
"ClinVar BED" -> "bedtools intersect -u";
"ClinVar GRCh38 regions BED" -> "CrossMap1";
"CrossMap1" -> "ClinVar T2Tv1.0 regions BED";
"ClinVar T2Tv1.0 regions BED" -> "CrossMap2";
"T2Tv1.0 -> T2Tv1.1 chain file" -> "CrossMap2";
"CrossMap2" -> "ClinVar T2Tv1.1 regions BED";
}