analysis_bam.svg

<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"
 "http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
<!-- Generated by graphviz version 2.47.3 (20210619.1520)
 -->
<!-- Title: workflow Pages: 1 -->
<svg width="3522pt" height="1327pt"
 viewBox="0.00 0.00 3522.16 1327.43" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g id="graph0" class="graph" transform="scale(1 1) rotate(0) translate(4 1323.43)">
<title>workflow</title>
<polygon fill="white" stroke="transparent" points="-4,4 -4,-1323.43 3518.16,-1323.43 3518.16,4 -4,4"/>
<!-- grch38 fasta -->
<g id="node1" class="node">
<title>grch38 fasta</title>
<polygon fill="none" stroke="black" points="1268.73,-1214.43 1237.69,-1232.43 1175.61,-1232.43 1144.56,-1214.43 1175.61,-1196.43 1237.69,-1196.43 1268.73,-1214.43"/>
<text text-anchor="middle" x="1206.65" y="-1210.73" font-family="Times,serif" font-size="14.00">grch38 fasta</text>
</g>
<!-- MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2 -->
<g id="node4" class="node">
<title>MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="1507.65,-1145.43 1248.65,-1107.43 1507.65,-1069.43 1766.65,-1107.43 1507.65,-1145.43"/>
<text text-anchor="middle" x="1507.65" y="-1111.23" font-family="Times,serif" font-size="14.00">MethylDackel extract &#45;@ 16 &#45;m &lt;Bismap&gt; </text>
<text text-anchor="middle" x="1507.65" y="-1096.23" font-family="Times,serif" font-size="14.00">&#45;q &lt;cutoff&gt; &#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97</text>
</g>
<!-- grch38 fasta&#45;&gt;MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2 -->
<g id="edge5" class="edge">
<title>grch38 fasta&#45;&gt;MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2</title>
<path fill="none" stroke="black" d="M1244.46,-1200.24C1289.04,-1184.69 1364.46,-1158.38 1423.17,-1137.9"/>
<polygon fill="black" stroke="black" points="1424.52,-1141.14 1432.81,-1134.54 1422.21,-1134.53 1424.52,-1141.14"/>
</g>
<!-- Bismap complete bedGraph -->
<g id="node2" class="node">
<title>Bismap complete bedGraph</title>
<ellipse fill="none" stroke="black" cx="1397.65" cy="-1214.43" rx="110.48" ry="18"/>
<text text-anchor="middle" x="1397.65" y="-1210.73" font-family="Times,serif" font-size="14.00">Bismap complete bedGraph</text>
</g>
<!-- Bismap complete bedGraph&#45;&gt;MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2 -->
<g id="edge2" class="edge">
<title>Bismap complete bedGraph&#45;&gt;MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2</title>
<path fill="none" stroke="black" d="M1415.32,-1196.56C1429.14,-1183.38 1448.9,-1164.51 1466.6,-1147.61"/>
<polygon fill="black" stroke="black" points="1469.22,-1149.95 1474.04,-1140.51 1464.39,-1144.89 1469.22,-1149.95"/>
<text text-anchor="middle" x="1471.15" y="-1167.23" font-family="Times,serif" font-size="14.00">Bismap</text>
</g>
<!-- PREPROC_1 -->
<!-- PREPROC_1&#45;&gt;Bismap complete bedGraph -->
<g id="edge1" class="edge">
<title>PREPROC_1&#45;&gt;Bismap complete bedGraph</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M1397.65,-1283.23C1397.65,-1271.6 1397.65,-1255.98 1397.65,-1242.67"/>
<polygon fill="black" stroke="black" points="1401.15,-1242.61 1397.65,-1232.61 1394.15,-1242.61 1401.15,-1242.61"/>
<text text-anchor="middle" x="1450.65" y="-1254.23" font-family="Times,serif" font-size="14.00">from preprocessing</text>
</g>
<!-- CpG -->
<g id="node11" class="node">
<title>CpG</title>
<ellipse fill="none" stroke="black" cx="1432.65" cy="-1014.43" rx="27.1" ry="18"/>
<text text-anchor="middle" x="1432.65" y="-1010.73" font-family="Times,serif" font-size="14.00">CpG</text>
</g>
<!-- MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2&#45;&gt;CpG -->
<g id="edge13" class="edge">
<title>MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2&#45;&gt;CpG</title>
<path fill="none" stroke="black" d="M1480.36,-1073.32C1470.84,-1061.77 1460.37,-1049.07 1451.67,-1038.52"/>
<polygon fill="black" stroke="black" points="1454.19,-1036.07 1445.13,-1030.58 1448.79,-1040.52 1454.19,-1036.07"/>
</g>
<!-- CHG -->
<g id="node12" class="node">
<title>CHG</title>
<ellipse fill="none" stroke="black" cx="1507.65" cy="-1014.43" rx="29.5" ry="18"/>
<text text-anchor="middle" x="1507.65" y="-1010.73" font-family="Times,serif" font-size="14.00">CHG</text>
</g>
<!-- MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2&#45;&gt;CHG -->
<g id="edge12" class="edge">
<title>MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2&#45;&gt;CHG</title>
<path fill="none" stroke="black" d="M1507.65,-1069.42C1507.65,-1060.42 1507.65,-1050.97 1507.65,-1042.53"/>
<polygon fill="black" stroke="black" points="1511.15,-1042.5 1507.65,-1032.5 1504.15,-1042.5 1511.15,-1042.5"/>
</g>
<!-- CHH -->
<g id="node13" class="node">
<title>CHH</title>
<ellipse fill="none" stroke="black" cx="1583.65" cy="-1014.43" rx="29.5" ry="18"/>
<text text-anchor="middle" x="1583.65" y="-1010.73" font-family="Times,serif" font-size="14.00">CHH</text>
</g>
<!-- MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2&#45;&gt;CHH -->
<g id="edge14" class="edge">
<title>MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2&#45;&gt;CHH</title>
<path fill="none" stroke="black" d="M1535.3,-1073.32C1544.81,-1061.93 1555.26,-1049.43 1564,-1038.96"/>
<polygon fill="black" stroke="black" points="1566.87,-1040.98 1570.59,-1031.06 1561.5,-1036.49 1566.87,-1040.98"/>
</g>
<!-- Methylation sequencing reads -->
<g id="node5" class="node">
<title>Methylation sequencing reads</title>
<polygon fill="none" stroke="black" points="1753.82,-1301.43 1686.23,-1319.43 1551.06,-1319.43 1483.47,-1301.43 1551.06,-1283.43 1686.23,-1283.43 1753.82,-1301.43"/>
<text text-anchor="middle" x="1618.65" y="-1297.73" font-family="Times,serif" font-size="14.00">Methylation sequencing reads</text>
</g>
<!-- sorted alignment BAM -->
<g id="node6" class="node">
<title>sorted alignment BAM</title>
<ellipse fill="none" stroke="black" cx="1618.65" cy="-1214.43" rx="92.88" ry="18"/>
<text text-anchor="middle" x="1618.65" y="-1210.73" font-family="Times,serif" font-size="14.00">sorted alignment BAM</text>
</g>
<!-- Methylation sequencing reads&#45;&gt;sorted alignment BAM -->
<g id="edge3" class="edge">
<title>Methylation sequencing reads&#45;&gt;sorted alignment BAM</title>
<path fill="none" stroke="black" d="M1618.65,-1283.23C1618.65,-1271.6 1618.65,-1255.98 1618.65,-1242.67"/>
<polygon fill="black" stroke="black" points="1622.15,-1242.61 1618.65,-1232.61 1615.15,-1242.61 1622.15,-1242.61"/>
<text text-anchor="middle" x="1637.65" y="-1254.23" font-family="Times,serif" font-size="14.00">aligner</text>
</g>
<!-- sorted alignment BAM&#45;&gt;MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2 -->
<g id="edge11" class="edge">
<title>sorted alignment BAM&#45;&gt;MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2</title>
<path fill="none" stroke="black" d="M1600.81,-1196.56C1586.87,-1183.38 1566.92,-1164.51 1549.06,-1147.61"/>
<polygon fill="black" stroke="black" points="1551.23,-1144.84 1541.56,-1140.51 1546.42,-1149.93 1551.23,-1144.84"/>
</g>
<!-- MapQ cutoff -->
<g id="node7" class="node">
<title>MapQ cutoff</title>
<polygon fill="none" stroke="black" points="1861.73,-1214.43 1828.69,-1232.43 1762.6,-1232.43 1729.56,-1214.43 1762.6,-1196.43 1828.69,-1196.43 1861.73,-1214.43"/>
<text text-anchor="middle" x="1795.65" y="-1210.73" font-family="Times,serif" font-size="14.00">MapQ cutoff</text>
</g>
<!-- MapQ cutoff&#45;&gt;MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2 -->
<g id="edge4" class="edge">
<title>MapQ cutoff&#45;&gt;MethylDackel extract &#45;@ 16 &#45;q &lt;cutoff&gt;\n&#45;&#45;CHH &#45;&#45;CHG &#45;&#45;OT 2,0,0,97_2</title>
<path fill="none" stroke="black" d="M1757.1,-1199.38C1714.52,-1183.86 1644.7,-1158.4 1589.74,-1138.36"/>
<polygon fill="black" stroke="black" points="1590.61,-1134.95 1580.01,-1134.82 1588.21,-1141.53 1590.61,-1134.95"/>
<text text-anchor="middle" x="1710.15" y="-1167.23" font-family="Times,serif" font-size="14.00">cutoff</text>
</g>
<!-- sort &#45;k1,1 &#45;k2,2n_2 -->
<g id="node8" class="node">
<title>sort &#45;k1,1 &#45;k2,2n_2</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="1507.65,-959.43 1408.61,-941.43 1507.65,-923.43 1606.69,-941.43 1507.65,-959.43"/>
<text text-anchor="middle" x="1507.65" y="-937.73" font-family="Times,serif" font-size="14.00">sort &#45;k1,1 &#45;k2,2n</text>
</g>
<!-- sorted methyl calls bedGraph -->
<g id="node9" class="node">
<title>sorted methyl calls bedGraph</title>
<ellipse fill="none" stroke="black" cx="1507.65" cy="-868.43" rx="115.08" ry="18"/>
<text text-anchor="middle" x="1507.65" y="-864.73" font-family="Times,serif" font-size="14.00">sorted methyl calls bedGraph</text>
</g>
<!-- sort &#45;k1,1 &#45;k2,2n_2&#45;&gt;sorted methyl calls bedGraph -->
<g id="edge6" class="edge">
<title>sort &#45;k1,1 &#45;k2,2n_2&#45;&gt;sorted methyl calls bedGraph</title>
<path fill="none" stroke="black" d="M1507.65,-923.25C1507.65,-915.22 1507.65,-905.48 1507.65,-896.5"/>
<polygon fill="black" stroke="black" points="1511.15,-896.46 1507.65,-886.46 1504.15,-896.46 1511.15,-896.46"/>
</g>
<!-- sorted methyl calls bedGraph\n(no headers) -->
<g id="node10" class="node">
<title>sorted methyl calls bedGraph\n(no headers)</title>
<ellipse fill="none" stroke="black" cx="1507.65" cy="-772.56" rx="125.32" ry="26.74"/>
<text text-anchor="middle" x="1507.65" y="-776.36" font-family="Times,serif" font-size="14.00">sorted methyl calls bedGraph</text>
<text text-anchor="middle" x="1507.65" y="-761.36" font-family="Times,serif" font-size="14.00">(no headers)</text>
</g>
<!-- sorted methyl calls bedGraph&#45;&gt;sorted methyl calls bedGraph\n(no headers) -->
<g id="edge7" class="edge">
<title>sorted methyl calls bedGraph&#45;&gt;sorted methyl calls bedGraph\n(no headers)</title>
<path fill="none" stroke="black" d="M1507.65,-850.22C1507.65,-838.91 1507.65,-823.7 1507.65,-809.81"/>
<polygon fill="black" stroke="black" points="1511.15,-809.53 1507.65,-799.53 1504.15,-809.53 1511.15,-809.53"/>
<text text-anchor="middle" x="1535.65" y="-821.23" font-family="Times,serif" font-size="14.00">head &#45;n &#45;3</text>
</g>
<!-- bedtools intersect &#45;u_2 -->
<g id="node14" class="node">
<title>bedtools intersect &#45;u_2</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="912.65,-694.69 626.65,-656.69 912.65,-618.69 1198.65,-656.69 912.65,-694.69"/>
<text text-anchor="middle" x="912.65" y="-660.49" font-family="Times,serif" font-size="14.00">bedtools intersect &#45;u &#45;sorted </text>
<text text-anchor="middle" x="912.65" y="-645.49" font-family="Times,serif" font-size="14.00">&#45;a &lt;methyl calls&gt; &#45;b &lt;low map. clinvar regions&gt;</text>
</g>
<!-- sorted methyl calls bedGraph\n(no headers)&#45;&gt;bedtools intersect &#45;u_2 -->
<g id="edge15" class="edge">
<title>sorted methyl calls bedGraph\n(no headers)&#45;&gt;bedtools intersect &#45;u_2</title>
<path fill="none" stroke="black" d="M1416.31,-754.08C1314.8,-734.66 1150.3,-703.18 1037.08,-681.51"/>
<polygon fill="black" stroke="black" points="1037.45,-678.01 1026.97,-679.57 1036.13,-684.89 1037.45,-678.01"/>
<text text-anchor="middle" x="1297.15" y="-716.49" font-family="Times,serif" font-size="14.00">methyl calls</text>
</g>
<!-- bedtools intersect &#45;u_2b -->
<g id="node15" class="node">
<title>bedtools intersect &#45;u_2b</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="1847.65,-694.69 1602.65,-656.69 1847.65,-618.69 2092.65,-656.69 1847.65,-694.69"/>
<text text-anchor="middle" x="1847.65" y="-660.49" font-family="Times,serif" font-size="14.00">bedtools intersect &#45;u &#45;sorted </text>
<text text-anchor="middle" x="1847.65" y="-645.49" font-family="Times,serif" font-size="14.00">&#45;a &lt;methyl calls&gt; &#45;b &lt;low map. regions&gt;</text>
</g>
<!-- sorted methyl calls bedGraph\n(no headers)&#45;&gt;bedtools intersect &#45;u_2b -->
<g id="edge16" class="edge">
<title>sorted methyl calls bedGraph\n(no headers)&#45;&gt;bedtools intersect &#45;u_2b</title>
<path fill="none" stroke="black" d="M1531.31,-746.08C1543.71,-734.2 1559.8,-720.9 1576.65,-712.69 1613.33,-694.82 1655.52,-682.76 1695.53,-674.61"/>
<polygon fill="black" stroke="black" points="1696.24,-678.04 1705.38,-672.68 1694.89,-671.17 1696.24,-678.04"/>
<text text-anchor="middle" x="1610.15" y="-716.49" font-family="Times,serif" font-size="14.00">methyl calls</text>
</g>
<!-- bedtools intersect &#45;u_2c -->
<g id="node16" class="node">
<title>bedtools intersect &#45;u_2c</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="2766.65,-694.69 2536.65,-656.69 2766.65,-618.69 2996.65,-656.69 2766.65,-694.69"/>
<text text-anchor="middle" x="2766.65" y="-660.49" font-family="Times,serif" font-size="14.00">bedtools intersect &#45;u &#45;sorted </text>
<text text-anchor="middle" x="2766.65" y="-645.49" font-family="Times,serif" font-size="14.00">&#45;a &lt;methyl calls&gt; &#45;b &lt;clinvar regions&gt;</text>
</g>
<!-- sorted methyl calls bedGraph\n(no headers)&#45;&gt;bedtools intersect &#45;u_2c -->
<g id="edge17" class="edge">
<title>sorted methyl calls bedGraph\n(no headers)&#45;&gt;bedtools intersect &#45;u_2c</title>
<path fill="none" stroke="black" d="M1606.17,-755.89C1646.33,-748.64 1693.03,-739.08 1734.65,-727.69 1753.77,-722.46 1757.16,-716.31 1776.65,-712.69 1940.78,-682.23 2361.37,-709.49 2527.65,-694.69 2567.67,-691.13 2611.17,-685.06 2649.95,-678.86"/>
<polygon fill="black" stroke="black" points="2650.92,-682.25 2660.23,-677.19 2649.8,-675.34 2650.92,-682.25"/>
<text text-anchor="middle" x="1810.15" y="-716.49" font-family="Times,serif" font-size="14.00">methyl calls</text>
</g>
<!-- awkfilter1 -->
<g id="node17" class="node">
<title>awkfilter1</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="2314.65,-674.69 2110.54,-656.69 2314.65,-638.69 2518.75,-656.69 2314.65,-674.69"/>
<text text-anchor="middle" x="2314.65" y="-652.99" font-family="Times,serif" font-size="14.00">awk &#39;($4 == 100 || $4 == 0) { print }&#39;</text>
</g>
<!-- sorted methyl calls bedGraph\n(no headers)&#45;&gt;awkfilter1 -->
<g id="edge18" class="edge">
<title>sorted methyl calls bedGraph\n(no headers)&#45;&gt;awkfilter1</title>
<path fill="none" stroke="black" d="M1582.67,-751.01C1602.86,-744.5 1624.41,-736.64 1643.65,-727.69 1655.05,-722.39 1655.63,-716.41 1667.65,-712.69 1759.86,-684.17 2005.67,-704.98 2101.65,-694.69 2150.78,-689.43 2205.7,-679.8 2247.47,-671.68"/>
<polygon fill="black" stroke="black" points="2248.39,-675.06 2257.53,-669.7 2247.04,-668.2 2248.39,-675.06"/>
<text text-anchor="middle" x="1701.15" y="-716.49" font-family="Times,serif" font-size="14.00">methyl calls</text>
</g>
<!-- awkfilter2 -->
<g id="node18" class="node">
<title>awkfilter2</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="259.65,-674.69 47.5,-656.69 259.65,-638.69 471.8,-656.69 259.65,-674.69"/>
<text text-anchor="middle" x="259.65" y="-652.99" font-family="Times,serif" font-size="14.00">awk &#39;($4 != 100 &amp;&amp; $4 != 0) { print }&#39;</text>
</g>
<!-- sorted methyl calls bedGraph\n(no headers)&#45;&gt;awkfilter2 -->
<g id="edge19" class="edge">
<title>sorted methyl calls bedGraph\n(no headers)&#45;&gt;awkfilter2</title>
<path fill="none" stroke="black" d="M1396.73,-759.92C1323.59,-751.85 1225.75,-740.32 1139.65,-727.69 1101.26,-722.06 1092.24,-716.65 1053.65,-712.69 801.51,-686.85 735.9,-719.4 483.65,-694.69 431.6,-689.6 373.34,-679.85 329.29,-671.63"/>
<polygon fill="black" stroke="black" points="329.73,-668.16 319.26,-669.74 328.44,-675.03 329.73,-668.16"/>
<text text-anchor="middle" x="1173.15" y="-716.49" font-family="Times,serif" font-size="14.00">methyl calls</text>
</g>
<!-- sorted methyl calls bedGraph\n(no headers) with params -->
<g id="node49" class="node">
<title>sorted methyl calls bedGraph\n(no headers) with params</title>
<polygon fill="none" stroke="black" points="1584.16,-694.69 1292.23,-694.69 1217.14,-618.69 1509.06,-618.69 1584.16,-694.69"/>
<text text-anchor="middle" x="1400.65" y="-660.49" font-family="Times,serif" font-size="14.00">sorted methyl calls bedGraph</text>
<text text-anchor="middle" x="1400.65" y="-645.49" font-family="Times,serif" font-size="14.00">(no headers) with params</text>
</g>
<!-- sorted methyl calls bedGraph\n(no headers)&#45;&gt;sorted methyl calls bedGraph\n(no headers) with params -->
<g id="edge50" class="edge">
<title>sorted methyl calls bedGraph\n(no headers)&#45;&gt;sorted methyl calls bedGraph\n(no headers) with params</title>
<path fill="none" stroke="black" d="M1483.63,-746.01C1471.57,-733.17 1456.59,-717.23 1442.69,-702.44"/>
<polygon fill="black" stroke="black" points="1444.94,-699.72 1435.54,-694.83 1439.84,-704.51 1444.94,-699.72"/>
<text text-anchor="middle" x="1506.65" y="-716.49" font-family="Times,serif" font-size="14.00">add_params.py</text>
</g>
<!-- CpG&#45;&gt;sort &#45;k1,1 &#45;k2,2n_2 -->
<g id="edge8" class="edge">
<title>CpG&#45;&gt;sort &#45;k1,1 &#45;k2,2n_2</title>
<path fill="none" stroke="black" d="M1447.82,-999.07C1458.5,-988.96 1472.97,-975.26 1485,-963.87"/>
<polygon fill="black" stroke="black" points="1487.64,-966.2 1492.49,-956.78 1482.82,-961.11 1487.64,-966.2"/>
</g>
<!-- CHG&#45;&gt;sort &#45;k1,1 &#45;k2,2n_2 -->
<g id="edge9" class="edge">
<title>CHG&#45;&gt;sort &#45;k1,1 &#45;k2,2n_2</title>
<path fill="none" stroke="black" d="M1507.65,-996.25C1507.65,-988.22 1507.65,-978.48 1507.65,-969.5"/>
<polygon fill="black" stroke="black" points="1511.15,-969.46 1507.65,-959.46 1504.15,-969.46 1511.15,-969.46"/>
</g>
<!-- CHH&#45;&gt;sort &#45;k1,1 &#45;k2,2n_2 -->
<g id="edge10" class="edge">
<title>CHH&#45;&gt;sort &#45;k1,1 &#45;k2,2n_2</title>
<path fill="none" stroke="black" d="M1568.27,-999.07C1557.45,-988.96 1542.79,-975.26 1530.6,-963.87"/>
<polygon fill="black" stroke="black" points="1532.7,-961.05 1523,-956.78 1527.92,-966.16 1532.7,-961.05"/>
</g>
<!-- methylation calls in low\nmappability regions in ClinVar regions -->
<g id="node28" class="node">
<title>methylation calls in low\nmappability regions in ClinVar regions</title>
<ellipse fill="none" stroke="black" cx="912.65" cy="-540.82" rx="161.85" ry="26.74"/>
<text text-anchor="middle" x="912.65" y="-544.62" font-family="Times,serif" font-size="14.00">methylation calls in low</text>
<text text-anchor="middle" x="912.65" y="-529.62" font-family="Times,serif" font-size="14.00">mappability regions in ClinVar regions</text>
</g>
<!-- bedtools intersect &#45;u_2&#45;&gt;methylation calls in low\nmappability regions in ClinVar regions -->
<g id="edge28" class="edge">
<title>bedtools intersect &#45;u_2&#45;&gt;methylation calls in low\nmappability regions in ClinVar regions</title>
<path fill="none" stroke="black" d="M912.65,-618.41C912.65,-605.54 912.65,-591.14 912.65,-578.33"/>
<polygon fill="black" stroke="black" points="916.15,-578.01 912.65,-568.01 909.15,-578.01 916.15,-578.01"/>
</g>
<!-- methylation calls in low\nmappability regions -->
<g id="node29" class="node">
<title>methylation calls in low\nmappability regions</title>
<ellipse fill="none" stroke="black" cx="1847.65" cy="-540.82" rx="105.22" ry="26.74"/>
<text text-anchor="middle" x="1847.65" y="-544.62" font-family="Times,serif" font-size="14.00">methylation calls in low</text>
<text text-anchor="middle" x="1847.65" y="-529.62" font-family="Times,serif" font-size="14.00">mappability regions</text>
</g>
<!-- bedtools intersect &#45;u_2b&#45;&gt;methylation calls in low\nmappability regions -->
<g id="edge29" class="edge">
<title>bedtools intersect &#45;u_2b&#45;&gt;methylation calls in low\nmappability regions</title>
<path fill="none" stroke="black" d="M1847.65,-618.41C1847.65,-605.54 1847.65,-591.14 1847.65,-578.33"/>
<polygon fill="black" stroke="black" points="1851.15,-578.01 1847.65,-568.01 1844.15,-578.01 1851.15,-578.01"/>
</g>
<!-- methylation calls in ClinVar regions -->
<g id="node30" class="node">
<title>methylation calls in ClinVar regions</title>
<ellipse fill="none" stroke="black" cx="2766.65" cy="-540.82" rx="139.18" ry="18"/>
<text text-anchor="middle" x="2766.65" y="-537.12" font-family="Times,serif" font-size="14.00">methylation calls in ClinVar regions</text>
</g>
<!-- bedtools intersect &#45;u_2c&#45;&gt;methylation calls in ClinVar regions -->
<g id="edge30" class="edge">
<title>bedtools intersect &#45;u_2c&#45;&gt;methylation calls in ClinVar regions</title>
<path fill="none" stroke="black" d="M2766.65,-618.41C2766.65,-602.33 2766.65,-583.87 2766.65,-569.11"/>
<polygon fill="black" stroke="black" points="2770.15,-569.09 2766.65,-559.09 2763.15,-569.09 2770.15,-569.09"/>
</g>
<!-- pure methylation calls -->
<g id="node39" class="node">
<title>pure methylation calls</title>
<ellipse fill="none" stroke="black" cx="2314.65" cy="-540.82" rx="89.88" ry="18"/>
<text text-anchor="middle" x="2314.65" y="-537.12" font-family="Times,serif" font-size="14.00">pure methylation calls</text>
</g>
<!-- awkfilter1&#45;&gt;pure methylation calls -->
<g id="edge36" class="edge">
<title>awkfilter1&#45;&gt;pure methylation calls</title>
<path fill="none" stroke="black" d="M2314.65,-638.4C2314.65,-619.98 2314.65,-590.58 2314.65,-569.05"/>
<polygon fill="black" stroke="black" points="2318.15,-568.91 2314.65,-558.91 2311.15,-568.91 2318.15,-568.91"/>
</g>
<!-- mixed methylation calls -->
<g id="node40" class="node">
<title>mixed methylation calls</title>
<ellipse fill="none" stroke="black" cx="249.65" cy="-540.82" rx="96.68" ry="18"/>
<text text-anchor="middle" x="249.65" y="-537.12" font-family="Times,serif" font-size="14.00">mixed methylation calls</text>
</g>
<!-- awkfilter2&#45;&gt;mixed methylation calls -->
<g id="edge37" class="edge">
<title>awkfilter2&#45;&gt;mixed methylation calls</title>
<path fill="none" stroke="black" d="M258.13,-638.4C256.51,-619.98 253.93,-590.58 252.04,-569.05"/>
<polygon fill="black" stroke="black" points="255.51,-568.56 251.15,-558.91 248.54,-569.18 255.51,-568.56"/>
</g>
<!-- sorted low mappability regions\nin ClinVar annotations 2 -->
<g id="node19" class="node">
<title>sorted low mappability regions\nin ClinVar annotations 2</title>
<ellipse fill="none" stroke="black" cx="912.65" cy="-772.56" rx="131.54" ry="26.74"/>
<text text-anchor="middle" x="912.65" y="-776.36" font-family="Times,serif" font-size="14.00">sorted low mappability regions</text>
<text text-anchor="middle" x="912.65" y="-761.36" font-family="Times,serif" font-size="14.00">in ClinVar annotations</text>
</g>
<!-- sorted low mappability regions\nin ClinVar annotations 2&#45;&gt;bedtools intersect &#45;u_2 -->
<g id="edge23" class="edge">
<title>sorted low mappability regions\nin ClinVar annotations 2&#45;&gt;bedtools intersect &#45;u_2</title>
<path fill="none" stroke="black" d="M912.65,-745.44C912.65,-733.51 912.65,-719 912.65,-705.24"/>
<polygon fill="black" stroke="black" points="916.15,-704.92 912.65,-694.92 909.15,-704.92 916.15,-704.92"/>
<text text-anchor="middle" x="981.15" y="-716.49" font-family="Times,serif" font-size="14.00">low map. clinvar regions</text>
</g>
<!-- sorted low mappability regions 2 -->
<g id="node20" class="node">
<title>sorted low mappability regions 2</title>
<ellipse fill="none" stroke="black" cx="1847.65" cy="-772.56" rx="131.54" ry="26.74"/>
<text text-anchor="middle" x="1847.65" y="-776.36" font-family="Times,serif" font-size="14.00">sorted low mappability regions</text>
<text text-anchor="middle" x="1847.65" y="-761.36" font-family="Times,serif" font-size="14.00">in genome</text>
</g>
<!-- sorted low mappability regions 2&#45;&gt;bedtools intersect &#45;u_2b -->
<g id="edge24" class="edge">
<title>sorted low mappability regions 2&#45;&gt;bedtools intersect &#45;u_2b</title>
<path fill="none" stroke="black" d="M1847.65,-745.44C1847.65,-733.51 1847.65,-719 1847.65,-705.24"/>
<polygon fill="black" stroke="black" points="1851.15,-704.92 1847.65,-694.92 1844.15,-704.92 1851.15,-704.92"/>
<text text-anchor="middle" x="1895.65" y="-716.49" font-family="Times,serif" font-size="14.00">low map. regions</text>
</g>
<!-- sorted clinvar regions 3 -->
<g id="node21" class="node">
<title>sorted clinvar regions 3</title>
<ellipse fill="none" stroke="black" cx="2766.65" cy="-772.56" rx="98.08" ry="26.74"/>
<text text-anchor="middle" x="2766.65" y="-776.36" font-family="Times,serif" font-size="14.00">sorted ClinVar regions</text>
<text text-anchor="middle" x="2766.65" y="-761.36" font-family="Times,serif" font-size="14.00">in genome</text>
</g>
<!-- sorted clinvar regions 3&#45;&gt;bedtools intersect &#45;u_2c -->
<g id="edge25" class="edge">
<title>sorted clinvar regions 3&#45;&gt;bedtools intersect &#45;u_2c</title>
<path fill="none" stroke="black" d="M2766.65,-745.44C2766.65,-733.51 2766.65,-719 2766.65,-705.24"/>
<polygon fill="black" stroke="black" points="2770.15,-704.92 2766.65,-694.92 2763.15,-704.92 2770.15,-704.92"/>
<text text-anchor="middle" x="2807.65" y="-716.49" font-family="Times,serif" font-size="14.00">clinvar regions</text>
</g>
<!-- PREPROC_2 -->
<!-- PREPROC_2&#45;&gt;sorted low mappability regions\nin ClinVar annotations 2 -->
<g id="edge20" class="edge">
<title>PREPROC_2&#45;&gt;sorted low mappability regions\nin ClinVar annotations 2</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M912.65,-850.22C912.65,-838.91 912.65,-823.7 912.65,-809.81"/>
<polygon fill="black" stroke="black" points="916.15,-809.53 912.65,-799.53 909.15,-809.53 916.15,-809.53"/>
<text text-anchor="middle" x="965.65" y="-821.23" font-family="Times,serif" font-size="14.00">from preprocessing</text>
</g>
<!-- PREPROC_4 -->
<!-- PREPROC_4&#45;&gt;sorted low mappability regions 2 -->
<g id="edge21" class="edge">
<title>PREPROC_4&#45;&gt;sorted low mappability regions 2</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M1847.65,-850.22C1847.65,-838.91 1847.65,-823.7 1847.65,-809.81"/>
<polygon fill="black" stroke="black" points="1851.15,-809.53 1847.65,-799.53 1844.15,-809.53 1851.15,-809.53"/>
<text text-anchor="middle" x="1900.65" y="-821.23" font-family="Times,serif" font-size="14.00">from preprocessing</text>
</g>
<!-- PREPROC_5 -->
<!-- PREPROC_5&#45;&gt;sorted clinvar regions 3 -->
<g id="edge22" class="edge">
<title>PREPROC_5&#45;&gt;sorted clinvar regions 3</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M2766.65,-850.22C2766.65,-838.91 2766.65,-823.7 2766.65,-809.81"/>
<polygon fill="black" stroke="black" points="2770.15,-809.53 2766.65,-799.53 2763.15,-809.53 2770.15,-809.53"/>
<text text-anchor="middle" x="2819.65" y="-821.23" font-family="Times,serif" font-size="14.00">from preprocessing</text>
</g>
<!-- sorted low mappability ClinVar regions 2 -->
<g id="node25" class="node">
<title>sorted low mappability ClinVar regions 2</title>
<ellipse fill="none" stroke="black" cx="566.65" cy="-540.82" rx="150.57" ry="18"/>
<text text-anchor="middle" x="566.65" y="-537.12" font-family="Times,serif" font-size="14.00">sorted low mappability ClinVar regions</text>
</g>
<!-- bedtools intersect &#45;u_3b -->
<g id="node27" class="node">
<title>bedtools intersect &#45;u_3b</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="697.65,-462.95 437.65,-424.95 697.65,-386.95 957.65,-424.95 697.65,-462.95"/>
<text text-anchor="middle" x="697.65" y="-428.75" font-family="Times,serif" font-size="14.00">bedtools intersect &#45;wa &#45;sorted </text>
<text text-anchor="middle" x="697.65" y="-413.75" font-family="Times,serif" font-size="14.00">&#45;a &lt;ClinVar low map.&gt; &#45;b &lt;low map. calls&gt;</text>
</g>
<!-- sorted low mappability ClinVar regions 2&#45;&gt;bedtools intersect &#45;u_3b -->
<g id="edge27" class="edge">
<title>sorted low mappability ClinVar regions 2&#45;&gt;bedtools intersect &#45;u_3b</title>
<path fill="none" stroke="black" d="M586.26,-522.78C603.96,-507.39 630.64,-484.2 653.45,-464.37"/>
<polygon fill="black" stroke="black" points="655.77,-466.99 661.02,-457.79 651.18,-461.71 655.77,-466.99"/>
<text text-anchor="middle" x="681.65" y="-484.75" font-family="Times,serif" font-size="14.00">ClinVar low map.</text>
</g>
<!-- PREPROC_3 -->
<!-- PREPROC_3&#45;&gt;sorted low mappability ClinVar regions 2 -->
<g id="edge26" class="edge">
<title>PREPROC_3&#45;&gt;sorted low mappability ClinVar regions 2</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M553.08,-638.4C555.66,-619.98 559.8,-590.58 562.82,-569.05"/>
<polygon fill="black" stroke="black" points="566.32,-569.3 564.25,-558.91 559.39,-568.32 566.32,-569.3"/>
<text text-anchor="middle" x="613.65" y="-589.49" font-family="Times,serif" font-size="14.00">from preprocessing</text>
</g>
<!-- ClinVar regions for each low map. methylation call -->
<g id="node32" class="node">
<title>ClinVar regions for each low map. methylation call</title>
<ellipse fill="none" stroke="black" cx="547.65" cy="-298.48" rx="193.87" ry="18"/>
<text text-anchor="middle" x="547.65" y="-294.78" font-family="Times,serif" font-size="14.00">ClinVar regions for each low map. methylation call</text>
</g>
<!-- bedtools intersect &#45;u_3b&#45;&gt;ClinVar regions for each low map. methylation call -->
<g id="edge32" class="edge">
<title>bedtools intersect &#45;u_3b&#45;&gt;ClinVar regions for each low map. methylation call</title>
<path fill="none" stroke="black" d="M659.41,-392.22C633.47,-370.7 599.75,-342.71 576.16,-323.14"/>
<polygon fill="black" stroke="black" points="578.2,-320.28 568.27,-316.59 573.73,-325.67 578.2,-320.28"/>
</g>
<!-- methylation calls in low\nmappability regions in ClinVar regions&#45;&gt;bedtools intersect &#45;u_3b -->
<g id="edge31" class="edge">
<title>methylation calls in low\nmappability regions in ClinVar regions&#45;&gt;bedtools intersect &#45;u_3b</title>
<path fill="none" stroke="black" d="M865.73,-514.97C835.24,-498.82 795.09,-477.56 761.67,-459.86"/>
<polygon fill="black" stroke="black" points="762.92,-456.56 752.44,-454.97 759.64,-462.75 762.92,-456.56"/>
<text text-anchor="middle" x="865.15" y="-484.75" font-family="Times,serif" font-size="14.00">low map. calls</text>
</g>
<!-- awkaligner1 -->
<g id="node36" class="node">
<title>awkaligner1</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="1146.65,-442.95 976.56,-424.95 1146.65,-406.95 1316.73,-424.95 1146.65,-442.95"/>
<text text-anchor="middle" x="1146.65" y="-421.25" font-family="Times,serif" font-size="14.00">awk &#39;{ print $0&quot;\t&lt;aligner&gt;&quot; }&#39;</text>
</g>
<!-- methylation calls in low\nmappability regions in ClinVar regions&#45;&gt;awkaligner1 -->
<g id="edge38" class="edge">
<title>methylation calls in low\nmappability regions in ClinVar regions&#45;&gt;awkaligner1</title>
<path fill="none" stroke="black" d="M963.13,-515.26C1006.92,-493.95 1069.31,-463.59 1108.88,-444.33"/>
<polygon fill="black" stroke="black" points="1110.48,-447.45 1117.94,-439.92 1107.42,-441.15 1110.48,-447.45"/>
</g>
<!-- awkaligner2 -->
<g id="node37" class="node">
<title>awkaligner2</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="1847.65,-442.95 1677.56,-424.95 1847.65,-406.95 2017.73,-424.95 1847.65,-442.95"/>
<text text-anchor="middle" x="1847.65" y="-421.25" font-family="Times,serif" font-size="14.00">awk &#39;{ print $0&quot;\t&lt;aligner&gt;&quot; }&#39;</text>
</g>
<!-- methylation calls in low\nmappability regions&#45;&gt;awkaligner2 -->
<g id="edge41" class="edge">
<title>methylation calls in low\nmappability regions&#45;&gt;awkaligner2</title>
<path fill="none" stroke="black" d="M1847.65,-513.7C1847.65,-495.59 1847.65,-471.52 1847.65,-453.17"/>
<polygon fill="black" stroke="black" points="1851.15,-453.12 1847.65,-443.12 1844.15,-453.12 1851.15,-453.12"/>
</g>
<!-- bedtools intersect &#45;u_3c -->
<g id="node31" class="node">
<title>bedtools intersect &#45;u_3c</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="2739.65,-462.95 2558.65,-424.95 2739.65,-386.95 2920.65,-424.95 2739.65,-462.95"/>
<text text-anchor="middle" x="2739.65" y="-428.75" font-family="Times,serif" font-size="14.00">bedtools intersect &#45;wa &#45;sorted </text>
<text text-anchor="middle" x="2739.65" y="-413.75" font-family="Times,serif" font-size="14.00">&#45;a &lt;ClinVar&gt; &#45;b &lt;calls&gt;</text>
</g>
<!-- methylation calls in ClinVar regions&#45;&gt;bedtools intersect &#45;u_3c -->
<g id="edge40" class="edge">
<title>methylation calls in ClinVar regions&#45;&gt;bedtools intersect &#45;u_3c</title>
<path fill="none" stroke="black" d="M2762.55,-522.53C2759.31,-508.88 2754.64,-489.2 2750.37,-471.16"/>
<polygon fill="black" stroke="black" points="2753.75,-470.28 2748.04,-461.35 2746.94,-471.89 2753.75,-470.28"/>
<text text-anchor="middle" x="2788.15" y="-484.75" font-family="Times,serif" font-size="14.00">clinvar calls</text>
</g>
<!-- awkaligner3 -->
<g id="node38" class="node">
<title>awkaligner3</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="3214.65,-442.95 3044.56,-424.95 3214.65,-406.95 3384.73,-424.95 3214.65,-442.95"/>
<text text-anchor="middle" x="3214.65" y="-421.25" font-family="Times,serif" font-size="14.00">awk &#39;{ print $0&quot;\t&lt;aligner&gt;&quot; }&#39;</text>
</g>
<!-- methylation calls in ClinVar regions&#45;&gt;awkaligner3 -->
<g id="edge39" class="edge">
<title>methylation calls in ClinVar regions&#45;&gt;awkaligner3</title>
<path fill="none" stroke="black" d="M2827.07,-524.47C2913.79,-502.42 3072.45,-462.09 3157.39,-440.51"/>
<polygon fill="black" stroke="black" points="3158.3,-443.89 3167.13,-438.03 3156.57,-437.1 3158.3,-443.89"/>
</g>
<!-- ClinVar regions for each methylation call -->
<g id="node33" class="node">
<title>ClinVar regions for each methylation call</title>
<ellipse fill="none" stroke="black" cx="2711.65" cy="-298.48" rx="157.07" ry="18"/>
<text text-anchor="middle" x="2711.65" y="-294.78" font-family="Times,serif" font-size="14.00">ClinVar regions for each methylation call</text>
</g>
<!-- bedtools intersect &#45;u_3c&#45;&gt;ClinVar regions for each methylation call -->
<g id="edge33" class="edge">
<title>bedtools intersect &#45;u_3c&#45;&gt;ClinVar regions for each methylation call</title>
<path fill="none" stroke="black" d="M2731.7,-388.62C2727.26,-368.91 2721.83,-344.74 2717.74,-326.55"/>
<polygon fill="black" stroke="black" points="2721.11,-325.59 2715.5,-316.61 2714.28,-327.13 2721.11,-325.59"/>
</g>
<!-- ClinVar regions with count of low map. methylation calls -->
<g id="node34" class="node">
<title>ClinVar regions with count of low map. methylation calls</title>
<ellipse fill="none" stroke="black" cx="543.65" cy="-192" rx="215.56" ry="18"/>
<text text-anchor="middle" x="543.65" y="-188.3" font-family="Times,serif" font-size="14.00">ClinVar regions with count of low map. methylation calls</text>
</g>
<!-- ClinVar regions for each low map. methylation call&#45;&gt;ClinVar regions with count of low map. methylation calls -->
<g id="edge34" class="edge">
<title>ClinVar regions for each low map. methylation call&#45;&gt;ClinVar regions with count of low map. methylation calls</title>
<path fill="none" stroke="black" d="M546.99,-280.22C546.36,-263.97 545.42,-239.34 544.7,-220.4"/>
<polygon fill="black" stroke="black" points="548.19,-220.14 544.31,-210.28 541.19,-220.4 548.19,-220.14"/>
<text text-anchor="middle" x="585.15" y="-231.8" font-family="Times,serif" font-size="14.00">count_calls.py</text>
</g>
<!-- ClinVar regions with count of methylation calls -->
<g id="node35" class="node">
<title>ClinVar regions with count of methylation calls</title>
<ellipse fill="none" stroke="black" cx="2703.65" cy="-192" rx="179.27" ry="18"/>
<text text-anchor="middle" x="2703.65" y="-188.3" font-family="Times,serif" font-size="14.00">ClinVar regions with count of methylation calls</text>
</g>
<!-- ClinVar regions for each methylation call&#45;&gt;ClinVar regions with count of methylation calls -->
<g id="edge35" class="edge">
<title>ClinVar regions for each methylation call&#45;&gt;ClinVar regions with count of methylation calls</title>
<path fill="none" stroke="black" d="M2710.33,-280.22C2709.08,-263.97 2707.2,-239.34 2705.74,-220.4"/>
<polygon fill="black" stroke="black" points="2709.22,-219.98 2704.97,-210.28 2702.24,-220.52 2709.22,-219.98"/>
<text text-anchor="middle" x="2746.15" y="-231.8" font-family="Times,serif" font-size="14.00">count_calls.py</text>
</g>
<!-- ClinVar regions with count of low map. methylation calls with params -->
<g id="node51" class="node">
<title>ClinVar regions with count of low map. methylation calls with params</title>
<polygon fill="none" stroke="black" points="901.72,-123 290.75,-123 133.58,-87 744.55,-87 901.72,-123"/>
<text text-anchor="middle" x="517.65" y="-101.3" font-family="Times,serif" font-size="14.00">ClinVar regions with count of low map. methylation calls with params</text>
</g>
<!-- ClinVar regions with count of low map. methylation calls&#45;&gt;ClinVar regions with count of low map. methylation calls with params -->
<g id="edge52" class="edge">
<title>ClinVar regions with count of low map. methylation calls&#45;&gt;ClinVar regions with count of low map. methylation calls with params</title>
<path fill="none" stroke="black" d="M538.39,-173.8C534.79,-162.05 529.95,-146.24 525.86,-132.84"/>
<polygon fill="black" stroke="black" points="529.17,-131.71 522.9,-123.18 522.48,-133.76 529.17,-131.71"/>
<text text-anchor="middle" x="574.65" y="-144.8" font-family="Times,serif" font-size="14.00">add_params.py</text>
</g>
<!-- ClinVar regions with count of methylation calls with params -->
<g id="node52" class="node">
<title>ClinVar regions with count of methylation calls with params</title>
<polygon fill="none" stroke="black" points="3031.49,-123 2505.19,-123 2369.81,-87 2896.1,-87 3031.49,-123"/>
<text text-anchor="middle" x="2700.65" y="-101.3" font-family="Times,serif" font-size="14.00">ClinVar regions with count of methylation calls with params</text>
</g>
<!-- ClinVar regions with count of methylation calls&#45;&gt;ClinVar regions with count of methylation calls with params -->
<g id="edge53" class="edge">
<title>ClinVar regions with count of methylation calls&#45;&gt;ClinVar regions with count of methylation calls with params</title>
<path fill="none" stroke="black" d="M2703.04,-173.8C2702.63,-162.16 2702.08,-146.55 2701.61,-133.24"/>
<polygon fill="black" stroke="black" points="2705.1,-133.05 2701.25,-123.18 2698.11,-133.29 2705.1,-133.05"/>
<text text-anchor="middle" x="2743.65" y="-144.8" font-family="Times,serif" font-size="14.00">add_params.py</text>
</g>
<!-- methylation calls in low\nmappability regions in ClinVar regions\nwith identifier column -->
<g id="node42" class="node">
<title>methylation calls in low\nmappability regions in ClinVar regions\nwith identifier column</title>
<ellipse fill="none" stroke="black" cx="1154.65" cy="-298.48" rx="161.85" ry="37.45"/>
<text text-anchor="middle" x="1154.65" y="-309.78" font-family="Times,serif" font-size="14.00">methylation calls in low</text>
<text text-anchor="middle" x="1154.65" y="-294.78" font-family="Times,serif" font-size="14.00">mappability regions in ClinVar regions</text>
<text text-anchor="middle" x="1154.65" y="-279.78" font-family="Times,serif" font-size="14.00">with identifier column</text>
</g>
<!-- awkaligner1&#45;&gt;methylation calls in low\nmappability regions in ClinVar regions\nwith identifier column -->
<g id="edge43" class="edge">
<title>awkaligner1&#45;&gt;methylation calls in low\nmappability regions in ClinVar regions\nwith identifier column</title>
<path fill="none" stroke="black" d="M1147.74,-406.91C1148.75,-391.29 1150.28,-367.42 1151.65,-346.15"/>
<polygon fill="black" stroke="black" points="1155.15,-346.17 1152.3,-335.97 1148.17,-345.72 1155.15,-346.17"/>
</g>
<!-- methylation calls in low\nmappability regions\nwith identifier column -->
<g id="node41" class="node">
<title>methylation calls in low\nmappability regions\nwith identifier column</title>
<ellipse fill="none" stroke="black" cx="1847.65" cy="-298.48" rx="105.22" ry="37.45"/>
<text text-anchor="middle" x="1847.65" y="-309.78" font-family="Times,serif" font-size="14.00">methylation calls in low</text>
<text text-anchor="middle" x="1847.65" y="-294.78" font-family="Times,serif" font-size="14.00">mappability regions</text>
<text text-anchor="middle" x="1847.65" y="-279.78" font-family="Times,serif" font-size="14.00">with identifier column</text>
</g>
<!-- awkaligner2&#45;&gt;methylation calls in low\nmappability regions\nwith identifier column -->
<g id="edge42" class="edge">
<title>awkaligner2&#45;&gt;methylation calls in low\nmappability regions\nwith identifier column</title>
<path fill="none" stroke="black" d="M1847.65,-406.91C1847.65,-391.29 1847.65,-367.42 1847.65,-346.15"/>
<polygon fill="black" stroke="black" points="1851.15,-345.97 1847.65,-335.97 1844.15,-345.97 1851.15,-345.97"/>
</g>
<!-- methylation calls in ClinVar regions\nwith identifier column -->
<g id="node43" class="node">
<title>methylation calls in ClinVar regions\nwith identifier column</title>
<ellipse fill="none" stroke="black" cx="3267.65" cy="-298.48" rx="151.14" ry="26.74"/>
<text text-anchor="middle" x="3267.65" y="-302.28" font-family="Times,serif" font-size="14.00">methylation calls in ClinVar regions</text>
<text text-anchor="middle" x="3267.65" y="-287.28" font-family="Times,serif" font-size="14.00">with identifier column</text>
</g>
<!-- awkaligner3&#45;&gt;methylation calls in ClinVar regions\nwith identifier column -->
<g id="edge44" class="edge">
<title>awkaligner3&#45;&gt;methylation calls in ClinVar regions\nwith identifier column</title>
<path fill="none" stroke="black" d="M3221.58,-407.68C3229.42,-389.26 3242.43,-358.71 3252.62,-334.76"/>
<polygon fill="black" stroke="black" points="3255.96,-335.85 3256.66,-325.28 3249.52,-333.11 3255.96,-335.85"/>
</g>
<!-- pure methylation calls with params -->
<g id="node53" class="node">
<title>pure methylation calls with params</title>
<polygon fill="none" stroke="black" points="2513.98,-442.95 2196.89,-442.95 2115.32,-406.95 2432.41,-406.95 2513.98,-442.95"/>
<text text-anchor="middle" x="2314.65" y="-421.25" font-family="Times,serif" font-size="14.00">pure methylation calls with params</text>
</g>
<!-- pure methylation calls&#45;&gt;pure methylation calls with params -->
<g id="edge54" class="edge">
<title>pure methylation calls&#45;&gt;pure methylation calls with params</title>
<path fill="none" stroke="black" d="M2314.65,-522.53C2314.65,-504.11 2314.65,-474.71 2314.65,-453.18"/>
<polygon fill="black" stroke="black" points="2318.15,-453.04 2314.65,-443.04 2311.15,-453.04 2318.15,-453.04"/>
<text text-anchor="middle" x="2356.65" y="-484.75" font-family="Times,serif" font-size="14.00">add_params.py</text>
</g>
<!-- mixed methylation calls with params -->
<g id="node54" class="node">
<title>mixed methylation calls with params</title>
<polygon fill="none" stroke="black" points="419.44,-442.95 85.71,-442.95 -0.15,-406.95 333.59,-406.95 419.44,-442.95"/>
<text text-anchor="middle" x="209.65" y="-421.25" font-family="Times,serif" font-size="14.00">mixed methylation calls with params</text>
</g>
<!-- mixed methylation calls&#45;&gt;mixed methylation calls with params -->
<g id="edge55" class="edge">
<title>mixed methylation calls&#45;&gt;mixed methylation calls with params</title>
<path fill="none" stroke="black" d="M243.57,-522.53C237.04,-503.95 226.59,-474.18 219.01,-452.6"/>
<polygon fill="black" stroke="black" points="222.26,-451.31 215.65,-443.04 215.66,-453.63 222.26,-451.31"/>
<text text-anchor="middle" x="275.65" y="-484.75" font-family="Times,serif" font-size="14.00">add_params.py</text>
</g>
<!-- low map. methylation call data -->
<g id="node46" class="node">
<title>low map. methylation call data</title>
<ellipse fill="none" stroke="black" cx="1847.65" cy="-192" rx="121.58" ry="18"/>
<text text-anchor="middle" x="1847.65" y="-188.3" font-family="Times,serif" font-size="14.00">low map. methylation call data</text>
</g>
<!-- methylation calls in low\nmappability regions\nwith identifier column&#45;&gt;low map. methylation call data -->
<g id="edge47" class="edge">
<title>methylation calls in low\nmappability regions\nwith identifier column&#45;&gt;low map. methylation call data</title>
<path fill="none" stroke="black" stroke-dasharray="5,2" d="M1847.65,-260.94C1847.65,-247.56 1847.65,-232.66 1847.65,-220.23"/>
<polygon fill="black" stroke="black" points="1851.15,-220.05 1847.65,-210.05 1844.15,-220.05 1851.15,-220.05"/>
<text text-anchor="middle" x="1899.65" y="-231.8" font-family="Times,serif" font-size="14.00">Display in Tableau</text>
</g>
<!-- low map. clinvar methylation call data -->
<g id="node45" class="node">
<title>low map. clinvar methylation call data</title>
<ellipse fill="none" stroke="black" cx="1202.65" cy="-192" rx="148.37" ry="18"/>
<text text-anchor="middle" x="1202.65" y="-188.3" font-family="Times,serif" font-size="14.00">low map. clinvar methylation call data</text>
</g>
<!-- methylation calls in low\nmappability regions in ClinVar regions\nwith identifier column&#45;&gt;low map. clinvar methylation call data -->
<g id="edge46" class="edge">
<title>methylation calls in low\nmappability regions in ClinVar regions\nwith identifier column&#45;&gt;low map. clinvar methylation call data</title>
<path fill="none" stroke="black" stroke-dasharray="5,2" d="M1171.43,-260.94C1177.71,-247.28 1184.7,-232.06 1190.48,-219.48"/>
<polygon fill="black" stroke="black" points="1193.82,-220.6 1194.81,-210.05 1187.46,-217.68 1193.82,-220.6"/>
<text text-anchor="middle" x="1238.65" y="-231.8" font-family="Times,serif" font-size="14.00">Display in Tableau</text>
</g>
<!-- clinvar methylation call data -->
<g id="node44" class="node">
<title>clinvar methylation call data</title>
<ellipse fill="none" stroke="black" cx="3280.65" cy="-192" rx="112.38" ry="18"/>
<text text-anchor="middle" x="3280.65" y="-188.3" font-family="Times,serif" font-size="14.00">clinvar methylation call data</text>
</g>
<!-- methylation calls in ClinVar regions\nwith identifier column&#45;&gt;clinvar methylation call data -->
<g id="edge45" class="edge">
<title>methylation calls in ClinVar regions\nwith identifier column&#45;&gt;clinvar methylation call data</title>
<path fill="none" stroke="black" stroke-dasharray="5,2" d="M3270.89,-271.38C3272.82,-255.86 3275.27,-236.2 3277.23,-220.45"/>
<polygon fill="black" stroke="black" points="3280.74,-220.62 3278.5,-210.27 3273.79,-219.76 3280.74,-220.62"/>
<text text-anchor="middle" x="3327.65" y="-231.8" font-family="Times,serif" font-size="14.00">Display in Tableau</text>
</g>
<!-- clinvar methylation call data with params -->
<g id="node50" class="node">
<title>clinvar methylation call data with params</title>
<polygon fill="gray" stroke="gray" points="3514.18,-123 3144.27,-123 3049.11,-87 3419.02,-87 3514.18,-123"/>
<text text-anchor="middle" x="3281.65" y="-101.3" font-family="Times,serif" font-size="14.00">clinvar methylation call data with params</text>
</g>
<!-- clinvar methylation call data&#45;&gt;clinvar methylation call data with params -->
<g id="edge51" class="edge">
<title>clinvar methylation call data&#45;&gt;clinvar methylation call data with params</title>
<path fill="none" stroke="black" d="M3280.85,-173.8C3280.99,-162.16 3281.17,-146.55 3281.33,-133.24"/>
<polygon fill="black" stroke="black" points="3284.83,-133.22 3281.44,-123.18 3277.83,-133.13 3284.83,-133.22"/>
<text text-anchor="middle" x="3322.65" y="-144.8" font-family="Times,serif" font-size="14.00">add_params.py</text>
</g>
<!-- low map. clinvar methylation call data with params -->
<g id="node48" class="node">
<title>low map. clinvar methylation call data with params</title>
<polygon fill="gray" stroke="gray" points="1489.5,-123 1036.36,-123 919.79,-87 1372.93,-87 1489.5,-123"/>
<text text-anchor="middle" x="1204.65" y="-101.3" font-family="Times,serif" font-size="14.00">low map. clinvar methylation call data with params</text>
</g>
<!-- low map. clinvar methylation call data&#45;&gt;low map. clinvar methylation call data with params -->
<g id="edge49" class="edge">
<title>low map. clinvar methylation call data&#45;&gt;low map. clinvar methylation call data with params</title>
<path fill="none" stroke="black" d="M1203.05,-173.8C1203.33,-162.16 1203.69,-146.55 1204.01,-133.24"/>
<polygon fill="black" stroke="black" points="1207.51,-133.26 1204.24,-123.18 1200.51,-133.09 1207.51,-133.26"/>
<text text-anchor="middle" x="1246.65" y="-144.8" font-family="Times,serif" font-size="14.00">add_params.py</text>
</g>
<!-- low map. methylation call data with params -->
<g id="node47" class="node">
<title>low map. methylation call data with params</title>
<polygon fill="gray" stroke="gray" points="2093.36,-123 1702.48,-123 1601.93,-87 1992.81,-87 2093.36,-123"/>
<text text-anchor="middle" x="1847.65" y="-101.3" font-family="Times,serif" font-size="14.00">low map. methylation call data with params</text>
</g>
<!-- low map. methylation call data&#45;&gt;low map. methylation call data with params -->
<g id="edge48" class="edge">
<title>low map. methylation call data&#45;&gt;low map. methylation call data with params</title>
<path fill="none" stroke="black" d="M1847.65,-173.8C1847.65,-162.16 1847.65,-146.55 1847.65,-133.24"/>
<polygon fill="black" stroke="black" points="1851.15,-133.18 1847.65,-123.18 1844.15,-133.18 1851.15,-133.18"/>
<text text-anchor="middle" x="1889.65" y="-144.8" font-family="Times,serif" font-size="14.00">add_params.py</text>
</g>
<!-- POSTPROC_7 -->
<!-- low map. methylation call data with params&#45;&gt;POSTPROC_7 -->
<g id="edge62" class="edge">
<title>low map. methylation call data with params&#45;&gt;POSTPROC_7</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M1847.65,-86.8C1847.65,-75.16 1847.65,-59.55 1847.65,-46.24"/>
<polygon fill="black" stroke="black" points="1851.15,-46.18 1847.65,-36.18 1844.15,-46.18 1851.15,-46.18"/>
<text text-anchor="middle" x="1895.15" y="-57.8" font-family="Times,serif" font-size="14.00">to postprocessing</text>
</g>
<!-- POSTPROC_8 -->
<!-- low map. clinvar methylation call data with params&#45;&gt;POSTPROC_8 -->
<g id="edge63" class="edge">
<title>low map. clinvar methylation call data with params&#45;&gt;POSTPROC_8</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M1204.65,-86.8C1204.65,-75.16 1204.65,-59.55 1204.65,-46.24"/>
<polygon fill="black" stroke="black" points="1208.15,-46.18 1204.65,-36.18 1201.15,-46.18 1208.15,-46.18"/>
<text text-anchor="middle" x="1252.15" y="-57.8" font-family="Times,serif" font-size="14.00">to postprocessing</text>
</g>
<!-- POSTPROC_1 -->
<!-- sorted methyl calls bedGraph\n(no headers) with params&#45;&gt;POSTPROC_1 -->
<g id="edge56" class="edge">
<title>sorted methyl calls bedGraph\n(no headers) with params&#45;&gt;POSTPROC_1</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M1400.65,-618.41C1400.65,-602.33 1400.65,-583.87 1400.65,-569.11"/>
<polygon fill="black" stroke="black" points="1404.15,-569.09 1400.65,-559.09 1397.15,-569.09 1404.15,-569.09"/>
<text text-anchor="middle" x="1448.15" y="-589.49" font-family="Times,serif" font-size="14.00">to postprocessing</text>
</g>
<!-- POSTPROC_2 -->
<!-- clinvar methylation call data with params&#45;&gt;POSTPROC_2 -->
<g id="edge57" class="edge">
<title>clinvar methylation call data with params&#45;&gt;POSTPROC_2</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M3281.65,-86.8C3281.65,-75.16 3281.65,-59.55 3281.65,-46.24"/>
<polygon fill="black" stroke="black" points="3285.15,-46.18 3281.65,-36.18 3278.15,-46.18 3285.15,-46.18"/>
<text text-anchor="middle" x="3329.15" y="-57.8" font-family="Times,serif" font-size="14.00">to postprocessing</text>
</g>
<!-- POSTPROC_6 -->
<!-- ClinVar regions with count of low map. methylation calls with params&#45;&gt;POSTPROC_6 -->
<g id="edge61" class="edge">
<title>ClinVar regions with count of low map. methylation calls with params&#45;&gt;POSTPROC_6</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M517.65,-86.8C517.65,-75.16 517.65,-59.55 517.65,-46.24"/>
<polygon fill="black" stroke="black" points="521.15,-46.18 517.65,-36.18 514.15,-46.18 521.15,-46.18"/>
<text text-anchor="middle" x="565.15" y="-57.8" font-family="Times,serif" font-size="14.00">to postprocessing</text>
</g>
<!-- POSTPROC_5 -->
<!-- ClinVar regions with count of methylation calls with params&#45;&gt;POSTPROC_5 -->
<g id="edge60" class="edge">
<title>ClinVar regions with count of methylation calls with params&#45;&gt;POSTPROC_5</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M2700.65,-86.8C2700.65,-75.16 2700.65,-59.55 2700.65,-46.24"/>
<polygon fill="black" stroke="black" points="2704.15,-46.18 2700.65,-36.18 2697.15,-46.18 2704.15,-46.18"/>
<text text-anchor="middle" x="2748.15" y="-57.8" font-family="Times,serif" font-size="14.00">to postprocessing</text>
</g>
<!-- POSTPROC_3 -->
<!-- pure methylation calls with params&#45;&gt;POSTPROC_3 -->
<g id="edge58" class="edge">
<title>pure methylation calls with params&#45;&gt;POSTPROC_3</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M2314.65,-406.91C2314.65,-386.26 2314.65,-351.21 2314.65,-326.78"/>
<polygon fill="black" stroke="black" points="2318.15,-326.72 2314.65,-316.72 2311.15,-326.72 2318.15,-326.72"/>
<text text-anchor="middle" x="2362.15" y="-357.75" font-family="Times,serif" font-size="14.00">to postprocessing</text>
</g>
<!-- POSTPROC_4 -->
<!-- mixed methylation calls with params&#45;&gt;POSTPROC_4 -->
<g id="edge59" class="edge">
<title>mixed methylation calls with params&#45;&gt;POSTPROC_4</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M209.65,-406.91C209.65,-386.26 209.65,-351.21 209.65,-326.78"/>
<polygon fill="black" stroke="black" points="213.15,-326.72 209.65,-316.72 206.15,-326.72 213.15,-326.72"/>
<text text-anchor="middle" x="257.15" y="-357.75" font-family="Times,serif" font-size="14.00">to postprocessing</text>
</g>
</g>
</svg>