analysis_preproc.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="1152pt" height="1736pt"
 viewBox="0.00 0.00 1152.00 1735.96" 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 1731.96)">
<title>workflow</title>
<polygon fill="white" stroke="transparent" points="-4,4 -4,-1731.96 1148,-1731.96 1148,4 -4,4"/>
<!-- ClinVar regions BED -->
<g id="node1" class="node">
<title>ClinVar regions BED</title>
<polygon fill="none" stroke="black" points="474.63,-1302.96 424.82,-1320.96 325.18,-1320.96 275.37,-1302.96 325.18,-1284.96 424.82,-1284.96 474.63,-1302.96"/>
<text text-anchor="middle" x="375" y="-1299.26" font-family="Times,serif" font-size="14.00">ClinVar regions BED</text>
</g>
<!-- expanded ClinVar regions BED -->
<g id="node32" class="node">
<title>expanded ClinVar regions BED</title>
<ellipse fill="none" stroke="black" cx="375" cy="-1215.96" rx="122.38" ry="18"/>
<text text-anchor="middle" x="375" y="-1212.26" font-family="Times,serif" font-size="14.00">expanded ClinVar regions BED</text>
</g>
<!-- ClinVar regions BED&#45;&gt;expanded ClinVar regions BED -->
<g id="edge34" class="edge">
<title>ClinVar regions BED&#45;&gt;expanded ClinVar regions BED</title>
<path fill="none" stroke="black" d="M375,-1284.76C375,-1273.12 375,-1257.51 375,-1244.2"/>
<polygon fill="black" stroke="black" points="378.5,-1244.14 375,-1234.14 371.5,-1244.14 378.5,-1244.14"/>
<text text-anchor="middle" x="451" y="-1255.76" font-family="Times,serif" font-size="14.00">bedtools slop &#45;l 2000 &#45;r 0 &#45;s</text>
</g>
<!-- grch38 fasta 2 -->
<g id="node2" class="node">
<title>grch38 fasta 2</title>
<polygon fill="none" stroke="black" points="884.09,-1709.96 848.05,-1727.96 775.95,-1727.96 739.91,-1709.96 775.95,-1691.96 848.05,-1691.96 884.09,-1709.96"/>
<text text-anchor="middle" x="812" y="-1706.26" font-family="Times,serif" font-size="14.00">reference fasta</text>
</g>
<!-- grch38 .fai file -->
<g id="node25" class="node">
<title>grch38 .fai file</title>
<ellipse fill="none" stroke="black" cx="812" cy="-1622.96" rx="63.09" ry="18"/>
<text text-anchor="middle" x="812" y="-1619.26" font-family="Times,serif" font-size="14.00">grch38 .fai file</text>
</g>
<!-- grch38 fasta 2&#45;&gt;grch38 .fai file -->
<g id="edge24" class="edge">
<title>grch38 fasta 2&#45;&gt;grch38 .fai file</title>
<path fill="none" stroke="black" d="M812,-1691.76C812,-1680.12 812,-1664.51 812,-1651.2"/>
<polygon fill="black" stroke="black" points="815.5,-1651.14 812,-1641.14 808.5,-1651.14 815.5,-1651.14"/>
<text text-anchor="middle" x="852.5" y="-1662.76" font-family="Times,serif" font-size="14.00">samtools faidx</text>
</g>
<!-- Bismap bedGraph -->
<g id="node3" class="node">
<title>Bismap bedGraph</title>
<ellipse fill="none" stroke="black" cx="637" cy="-1535.96" rx="75.29" ry="18"/>
<text text-anchor="middle" x="637" y="-1532.26" font-family="Times,serif" font-size="14.00">Bismap bedGraph</text>
</g>
<!-- sort &#45;k1,1 &#45;k2,2n_again -->
<g id="node4" class="node">
<title>sort &#45;k1,1 &#45;k2,2n_again</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="709,-1233.96 609.96,-1215.96 709,-1197.96 808.04,-1215.96 709,-1233.96"/>
<text text-anchor="middle" x="709" y="-1212.26" font-family="Times,serif" font-size="14.00">sort &#45;k1,1 &#45;k2,2n</text>
</g>
<!-- Bismap bedGraph&#45;&gt;sort &#45;k1,1 &#45;k2,2n_again -->
<g id="edge1" class="edge">
<title>Bismap bedGraph&#45;&gt;sort &#45;k1,1 &#45;k2,2n_again</title>
<path fill="none" stroke="black" d="M685.29,-1521.97C729.14,-1508.71 787,-1486.87 787,-1463.96 787,-1463.96 787,-1463.96 787,-1301.96 787,-1273.44 762.87,-1250.35 741.45,-1235.48"/>
<polygon fill="black" stroke="black" points="743.1,-1232.37 732.83,-1229.8 739.26,-1238.22 743.1,-1232.37"/>
</g>
<!-- bedtools complement -->
<g id="node23" class="node">
<title>bedtools complement</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="636,-1480.96 512.67,-1462.96 636,-1444.96 759.33,-1462.96 636,-1480.96"/>
<text text-anchor="middle" x="636" y="-1459.26" font-family="Times,serif" font-size="14.00">bedtools complement</text>
</g>
<!-- Bismap bedGraph&#45;&gt;bedtools complement -->
<g id="edge22" class="edge">
<title>Bismap bedGraph&#45;&gt;bedtools complement</title>
<path fill="none" stroke="black" d="M636.76,-1517.77C636.64,-1509.75 636.51,-1500.01 636.38,-1491.03"/>
<polygon fill="black" stroke="black" points="639.88,-1490.94 636.24,-1480.99 632.88,-1491.04 639.88,-1490.94"/>
</g>
<!-- Bismap complete bedGraph -->
<g id="node5" class="node">
<title>Bismap complete bedGraph</title>
<ellipse fill="none" stroke="black" cx="709" cy="-1128.96" rx="110.48" ry="18"/>
<text text-anchor="middle" x="709" y="-1125.26" font-family="Times,serif" font-size="14.00">Bismap complete bedGraph</text>
</g>
<!-- sort &#45;k1,1 &#45;k2,2n_again&#45;&gt;Bismap complete bedGraph -->
<g id="edge2" class="edge">
<title>sort &#45;k1,1 &#45;k2,2n_again&#45;&gt;Bismap complete bedGraph</title>
<path fill="none" stroke="black" d="M709,-1197.76C709,-1186.12 709,-1170.51 709,-1157.2"/>
<polygon fill="black" stroke="black" points="712.5,-1157.14 709,-1147.14 705.5,-1157.14 712.5,-1157.14"/>
</g>
<!-- BAM_PROC_1 -->
<!-- Bismap complete bedGraph&#45;&gt;BAM_PROC_1 -->
<g id="edge3" class="edge">
<title>Bismap complete bedGraph&#45;&gt;BAM_PROC_1</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M719.4,-1110.62C729.54,-1093.73 745.07,-1067.84 756.63,-1048.57"/>
<polygon fill="black" stroke="black" points="759.65,-1050.35 761.79,-1039.98 753.64,-1046.75 759.65,-1050.35"/>
<text text-anchor="middle" x="793" y="-1081.76" font-family="Times,serif" font-size="14.00">to BAM processing</text>
</g>
<!-- bedtools map &#45;o min -->
<g id="node9" class="node">
<title>bedtools map &#45;o min</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="519,-1059.96 352,-1021.96 519,-983.96 686,-1021.96 519,-1059.96"/>
<text text-anchor="middle" x="519" y="-1025.76" font-family="Times,serif" font-size="14.00">bedtools map &#45;o min </text>
<text text-anchor="middle" x="519" y="-1010.76" font-family="Times,serif" font-size="14.00">&#45;a &lt;ClinVar&gt; &#45;b &lt;Bismap&gt;</text>
</g>
<!-- Bismap complete bedGraph&#45;&gt;bedtools map &#45;o min -->
<g id="edge7" class="edge">
<title>Bismap complete bedGraph&#45;&gt;bedtools map &#45;o min</title>
<path fill="none" stroke="black" d="M688.3,-1111.28C674.97,-1100.96 657.01,-1087.85 640,-1077.96 623.18,-1068.18 604.32,-1058.88 586.58,-1050.8"/>
<polygon fill="black" stroke="black" points="587.98,-1047.59 577.43,-1046.69 585.12,-1053.97 587.98,-1047.59"/>
<text text-anchor="middle" x="683.5" y="-1081.76" font-family="Times,serif" font-size="14.00">Bismap</text>
</g>
<!-- bedtools map &#45;o mean -->
<g id="node11" class="node">
<title>bedtools map &#45;o mean</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="167,-1059.96 0,-1021.96 167,-983.96 334,-1021.96 167,-1059.96"/>
<text text-anchor="middle" x="167" y="-1025.76" font-family="Times,serif" font-size="14.00">bedtools map &#45;o mean </text>
<text text-anchor="middle" x="167" y="-1010.76" font-family="Times,serif" font-size="14.00">&#45;a &lt;ClinVar&gt; &#45;b &lt;Bismap&gt;</text>
</g>
<!-- Bismap complete bedGraph&#45;&gt;bedtools map &#45;o mean -->
<g id="edge9" class="edge">
<title>Bismap complete bedGraph&#45;&gt;bedtools map &#45;o mean</title>
<path fill="none" stroke="black" d="M652.6,-1113.41C608.56,-1102.45 545.7,-1087.7 490,-1077.96 425.16,-1066.62 407.93,-1070.74 343,-1059.96 315.29,-1055.36 285.33,-1049.42 258.11,-1043.66"/>
<polygon fill="black" stroke="black" points="258.73,-1040.21 248.22,-1041.55 257.27,-1047.06 258.73,-1040.21"/>
<text text-anchor="middle" x="584.5" y="-1081.76" font-family="Times,serif" font-size="14.00">Bismap</text>
</g>
<!-- awk_again -->
<g id="node27" class="node">
<title>awk_again</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="869,-500.74 715.5,-482.74 869,-464.74 1022.5,-482.74 869,-500.74"/>
<text text-anchor="middle" x="869" y="-479.04" font-family="Times,serif" font-size="14.00">awk &#39; if $4&lt;=0.01 { print }&#39;</text>
</g>
<!-- Bismap complete bedGraph&#45;&gt;awk_again -->
<g id="edge27" class="edge">
<title>Bismap complete bedGraph&#45;&gt;awk_again</title>
<path fill="none" stroke="black" d="M786.74,-1116.06C814.46,-1110.34 841.26,-1102.61 851,-1092.96 873.83,-1070.36 869,-1055.08 869,-1022.96 869,-1022.96 869,-1022.96 869,-563.61 869,-546.2 869,-526.66 869,-511.27"/>
<polygon fill="black" stroke="black" points="872.5,-510.88 869,-500.88 865.5,-510.88 872.5,-510.88"/>
</g>
<!-- 0 mappability regions bedGraph -->
<g id="node7" class="node">
<title>0 mappability regions bedGraph</title>
<ellipse fill="none" stroke="black" cx="634" cy="-1302.96" rx="125.38" ry="18"/>
<text text-anchor="middle" x="634" y="-1299.26" font-family="Times,serif" font-size="14.00">0 mappability regions bedGraph</text>
</g>
<!-- 0 mappability regions bedGraph&#45;&gt;sort &#45;k1,1 &#45;k2,2n_again -->
<g id="edge4" class="edge">
<title>0 mappability regions bedGraph&#45;&gt;sort &#45;k1,1 &#45;k2,2n_again</title>
<path fill="none" stroke="black" d="M649.18,-1284.76C660.63,-1271.78 676.44,-1253.86 688.94,-1239.7"/>
<polygon fill="black" stroke="black" points="691.86,-1241.68 695.85,-1231.86 686.61,-1237.05 691.86,-1241.68"/>
</g>
<!-- Bismap bigWig -->
<g id="node8" class="node">
<title>Bismap bigWig</title>
<polygon fill="none" stroke="black" points="714.1,-1622.96 675.55,-1640.96 598.45,-1640.96 559.9,-1622.96 598.45,-1604.96 675.55,-1604.96 714.1,-1622.96"/>
<text text-anchor="middle" x="637" y="-1619.26" font-family="Times,serif" font-size="14.00">Bismap bigWig</text>
</g>
<!-- Bismap bigWig&#45;&gt;Bismap bedGraph -->
<g id="edge5" class="edge">
<title>Bismap bigWig&#45;&gt;Bismap bedGraph</title>
<path fill="none" stroke="black" d="M637,-1604.76C637,-1593.12 637,-1577.51 637,-1564.2"/>
<polygon fill="black" stroke="black" points="640.5,-1564.14 637,-1554.14 633.5,-1564.14 640.5,-1564.14"/>
<text text-anchor="middle" x="693" y="-1575.76" font-family="Times,serif" font-size="14.00">bigWigToBedGraph</text>
</g>
<!-- ClinVar regions with\n min mappability bedGraph -->
<g id="node12" class="node">
<title>ClinVar regions with\n min mappability bedGraph</title>
<ellipse fill="none" stroke="black" cx="517" cy="-920.09" rx="118.59" ry="26.74"/>
<text text-anchor="middle" x="517" y="-923.89" font-family="Times,serif" font-size="14.00">ClinVar regions with</text>
<text text-anchor="middle" x="517" y="-908.89" font-family="Times,serif" font-size="14.00"> min mappability bedGraph</text>
</g>
<!-- bedtools map &#45;o min&#45;&gt;ClinVar regions with\n min mappability bedGraph -->
<g id="edge10" class="edge">
<title>bedtools map &#45;o min&#45;&gt;ClinVar regions with\n min mappability bedGraph</title>
<path fill="none" stroke="black" d="M518.26,-984.06C518.09,-975.37 517.9,-966.11 517.73,-957.44"/>
<polygon fill="black" stroke="black" points="521.22,-957.17 517.52,-947.24 514.22,-957.31 521.22,-957.17"/>
</g>
<!-- expanded sorted ClinVar regions BED -->
<g id="node10" class="node">
<title>expanded sorted ClinVar regions BED</title>
<ellipse fill="none" stroke="black" cx="375" cy="-1128.96" rx="146.47" ry="18"/>
<text text-anchor="middle" x="375" y="-1125.26" font-family="Times,serif" font-size="14.00">expanded sorted ClinVar regions BED</text>
</g>
<!-- expanded sorted ClinVar regions BED&#45;&gt;bedtools map &#45;o min -->
<g id="edge6" class="edge">
<title>expanded sorted ClinVar regions BED&#45;&gt;bedtools map &#45;o min</title>
<path fill="none" stroke="black" d="M398.14,-1111.09C417.99,-1096.61 447.24,-1075.29 471.99,-1057.24"/>
<polygon fill="black" stroke="black" points="474.17,-1059.98 480.19,-1051.26 470.05,-1054.32 474.17,-1059.98"/>
<text text-anchor="middle" x="464.5" y="-1081.76" font-family="Times,serif" font-size="14.00">ClinVar</text>
</g>
<!-- expanded sorted ClinVar regions BED&#45;&gt;bedtools map &#45;o mean -->
<g id="edge8" class="edge">
<title>expanded sorted ClinVar regions BED&#45;&gt;bedtools map &#45;o mean</title>
<path fill="none" stroke="black" d="M342.04,-1111.32C311.21,-1095.76 264.41,-1072.13 226.96,-1053.23"/>
<polygon fill="black" stroke="black" points="228.16,-1049.91 217.65,-1048.53 225,-1056.16 228.16,-1049.91"/>
<text text-anchor="middle" x="323.5" y="-1081.76" font-family="Times,serif" font-size="14.00">ClinVar</text>
</g>
<!-- ClinVar regions with\n mean mappability bedGraph -->
<g id="node13" class="node">
<title>ClinVar regions with\n mean mappability bedGraph</title>
<ellipse fill="none" stroke="black" cx="233" cy="-920.09" rx="124.4" ry="26.74"/>
<text text-anchor="middle" x="233" y="-923.89" font-family="Times,serif" font-size="14.00">ClinVar regions with</text>
<text text-anchor="middle" x="233" y="-908.89" font-family="Times,serif" font-size="14.00"> mean mappability bedGraph</text>
</g>
<!-- bedtools map &#45;o mean&#45;&gt;ClinVar regions with\n mean mappability bedGraph -->
<g id="edge11" class="edge">
<title>bedtools map &#45;o mean&#45;&gt;ClinVar regions with\n mean mappability bedGraph</title>
<path fill="none" stroke="black" d="M188.43,-988.54C195.5,-977.84 203.38,-965.92 210.52,-955.1"/>
<polygon fill="black" stroke="black" points="213.47,-956.99 216.06,-946.72 207.63,-953.13 213.47,-956.99"/>
</g>
<!-- paste -->
<g id="node14" class="node">
<title>paste</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="515,-856.22 474.61,-838.22 515,-820.22 555.39,-838.22 515,-856.22"/>
<text text-anchor="middle" x="515" y="-834.52" font-family="Times,serif" font-size="14.00">paste</text>
</g>
<!-- ClinVar regions with\n min mappability bedGraph&#45;&gt;paste -->
<g id="edge12" class="edge">
<title>ClinVar regions with\n min mappability bedGraph&#45;&gt;paste</title>
<path fill="none" stroke="black" d="M516.35,-892.94C516.13,-884.49 515.9,-875.11 515.69,-866.58"/>
<polygon fill="black" stroke="black" points="519.18,-866.28 515.43,-856.37 512.18,-866.45 519.18,-866.28"/>
</g>
<!-- ClinVar regions with\n mean mappability bedGraph&#45;&gt;paste -->
<g id="edge14" class="edge">
<title>ClinVar regions with\n mean mappability bedGraph&#45;&gt;paste</title>
<path fill="none" stroke="black" d="M306.36,-898.31C363.16,-882.22 439.17,-860.7 482.01,-848.56"/>
<polygon fill="black" stroke="black" points="482.98,-851.93 491.65,-845.84 481.07,-845.19 482.98,-851.93"/>
</g>
<!-- ClinVar regions with min and\n mean mappability bedGraph -->
<g id="node15" class="node">
<title>ClinVar regions with min and\n mean mappability bedGraph</title>
<ellipse fill="none" stroke="black" cx="515" cy="-756.35" rx="126.64" ry="26.74"/>
<text text-anchor="middle" x="515" y="-760.15" font-family="Times,serif" font-size="14.00">ClinVar regions with min and</text>
<text text-anchor="middle" x="515" y="-745.15" font-family="Times,serif" font-size="14.00"> mean mappability bedGraph</text>
</g>
<!-- paste&#45;&gt;ClinVar regions with min and\n mean mappability bedGraph -->
<g id="edge13" class="edge">
<title>paste&#45;&gt;ClinVar regions with min and\n mean mappability bedGraph</title>
<path fill="none" stroke="black" d="M515,-819.89C515,-812.16 515,-802.74 515,-793.58"/>
<polygon fill="black" stroke="black" points="518.5,-793.45 515,-783.45 511.5,-793.45 518.5,-793.45"/>
</g>
<!-- Mappability Distribution in ClinVar regions -->
<g id="node16" class="node">
<title>Mappability Distribution in ClinVar regions</title>
<polygon fill="gray" stroke="gray" points="795,-678.48 539,-678.48 539,-642.48 795,-642.48 795,-678.48"/>
<text text-anchor="middle" x="667" y="-656.78" font-family="Times,serif" font-size="14.00">Mappability Distribution in ClinVar regions</text>
</g>
<!-- ClinVar regions with min and\n mean mappability bedGraph&#45;&gt;Mappability Distribution in ClinVar regions -->
<g id="edge15" class="edge">
<title>ClinVar regions with min and\n mean mappability bedGraph&#45;&gt;Mappability Distribution in ClinVar regions</title>
<path fill="none" stroke="black" stroke-dasharray="5,2" d="M554.94,-730.69C578.39,-716.2 607.8,-698.04 630.61,-683.95"/>
<polygon fill="black" stroke="black" points="632.64,-686.81 639.31,-678.58 628.97,-680.85 632.64,-686.81"/>
<text text-anchor="middle" x="662" y="-700.28" font-family="Times,serif" font-size="14.00">Display in Tableau</text>
</g>
<!-- awk &#39;$11 &lt;= min_cutoff { print }&#39; -->
<g id="node19" class="node">
<title>awk &#39;$11 &lt;= min_cutoff { print }&#39;</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="511,-500.74 324.32,-482.74 511,-464.74 697.68,-482.74 511,-500.74"/>
<text text-anchor="middle" x="511" y="-479.04" font-family="Times,serif" font-size="14.00">awk &#39;$11 &lt;= min_cutoff { print }&#39;</text>
</g>
<!-- ClinVar regions with min and\n mean mappability bedGraph&#45;&gt;awk &#39;$11 &lt;= min_cutoff { print }&#39; -->
<g id="edge19" class="edge">
<title>ClinVar regions with min and\n mean mappability bedGraph&#45;&gt;awk &#39;$11 &lt;= min_cutoff { print }&#39;</title>
<path fill="none" stroke="black" d="M514.62,-729.31C513.86,-677.86 512.18,-563.66 511.4,-511.14"/>
<polygon fill="black" stroke="black" points="514.9,-510.86 511.25,-500.91 507.9,-510.96 514.9,-510.86"/>
</g>
<!-- low mappability cutoff value\n(min mappability &lt;= 0.01) -->
<g id="node18" class="node">
<title>low mappability cutoff value\n(min mappability &lt;= 0.01)</title>
<ellipse fill="none" stroke="black" cx="667" cy="-564.61" rx="123.99" ry="26.74"/>
<text text-anchor="middle" x="667" y="-568.41" font-family="Times,serif" font-size="14.00">low mappability cutoff value</text>
<text text-anchor="middle" x="667" y="-553.41" font-family="Times,serif" font-size="14.00">(min mappability &lt;= 0.01)</text>
</g>
<!-- Mappability Distribution in ClinVar regions&#45;&gt;low mappability cutoff value\n(min mappability &lt;= 0.01) -->
<g id="edge16" class="edge">
<title>Mappability Distribution in ClinVar regions&#45;&gt;low mappability cutoff value\n(min mappability &lt;= 0.01)</title>
<path fill="none" stroke="black" stroke-dasharray="5,2" d="M667,-642.26C667,-630.96 667,-615.75 667,-601.86"/>
<polygon fill="black" stroke="black" points="670.5,-601.58 667,-591.58 663.5,-601.58 670.5,-601.58"/>
<text text-anchor="middle" x="724.5" y="-613.28" font-family="Times,serif" font-size="14.00">Examine distribution</text>
</g>
<!-- Examine Data -->
<!-- low mappability cutoff value\n(min mappability &lt;= 0.01)&#45;&gt;awk &#39;$11 &lt;= min_cutoff { print }&#39; -->
<g id="edge17" class="edge">
<title>low mappability cutoff value\n(min mappability &lt;= 0.01)&#45;&gt;awk &#39;$11 &lt;= min_cutoff { print }&#39;</title>
<path fill="none" stroke="black" d="M620.19,-539.65C597.21,-527.88 569.9,-513.9 548.34,-502.86"/>
<polygon fill="black" stroke="black" points="549.68,-499.61 539.18,-498.17 546.49,-505.84 549.68,-499.61"/>
</g>
<!-- low mappability cutoff value\n(min mappability &lt;= 0.01)&#45;&gt;awk_again -->
<g id="edge28" class="edge">
<title>low mappability cutoff value\n(min mappability &lt;= 0.01)&#45;&gt;awk_again</title>
<path fill="none" stroke="black" d="M724.61,-540.83C756.88,-528.07 796.45,-512.43 826,-500.74"/>
<polygon fill="black" stroke="black" points="827.59,-503.88 835.61,-496.94 825.02,-497.37 827.59,-503.88"/>
</g>
<!-- low mappability ClinVar regions -->
<g id="node20" class="node">
<title>low mappability ClinVar regions</title>
<ellipse fill="none" stroke="black" cx="511" cy="-409.74" rx="126.98" ry="18"/>
<text text-anchor="middle" x="511" y="-406.04" font-family="Times,serif" font-size="14.00">low mappability ClinVar regions</text>
</g>
<!-- awk &#39;$11 &lt;= min_cutoff { print }&#39;&#45;&gt;low mappability ClinVar regions -->
<g id="edge18" class="edge">
<title>awk &#39;$11 &lt;= min_cutoff { print }&#39;&#45;&gt;low mappability ClinVar regions</title>
<path fill="none" stroke="black" d="M511,-464.55C511,-456.53 511,-446.79 511,-437.81"/>
<polygon fill="black" stroke="black" points="514.5,-437.77 511,-427.77 507.5,-437.77 514.5,-437.77"/>
</g>
<!-- sorted low mappability ClinVar regions -->
<g id="node21" class="node">
<title>sorted low mappability ClinVar regions</title>
<ellipse fill="none" stroke="black" cx="511" cy="-322.74" rx="150.57" ry="18"/>
<text text-anchor="middle" x="511" y="-319.04" font-family="Times,serif" font-size="14.00">sorted low mappability ClinVar regions</text>
</g>
<!-- low mappability ClinVar regions&#45;&gt;sorted low mappability ClinVar regions -->
<g id="edge20" class="edge">
<title>low mappability ClinVar regions&#45;&gt;sorted low mappability ClinVar regions</title>
<path fill="none" stroke="black" d="M511,-391.54C511,-379.9 511,-364.29 511,-350.98"/>
<polygon fill="black" stroke="black" points="514.5,-350.92 511,-340.92 507.5,-350.92 514.5,-350.92"/>
<text text-anchor="middle" x="557" y="-362.54" font-family="Times,serif" font-size="14.00">sort &#45;k1,1 &#45;k2,2n</text>
</g>
<!-- BAM_PROC_2 -->
<!-- sorted low mappability ClinVar regions&#45;&gt;BAM_PROC_2 -->
<g id="edge21" class="edge">
<title>sorted low mappability ClinVar regions&#45;&gt;BAM_PROC_2</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M510.5,-304.4C510.03,-287.96 509.32,-263.02 508.78,-243.93"/>
<polygon fill="black" stroke="black" points="512.27,-243.65 508.49,-233.76 505.27,-243.85 512.27,-243.65"/>
<text text-anchor="middle" x="617.5" y="-275.54" font-family="Times,serif" font-size="14.00">to BAM processing and postprocessing</text>
</g>
<!-- bedtools intersect &#45;u_0 -->
<g id="node30" class="node">
<title>bedtools intersect &#45;u_0</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="869,-253.74 594,-215.74 869,-177.74 1144,-215.74 869,-253.74"/>
<text text-anchor="middle" x="869" y="-219.54" font-family="Times,serif" font-size="14.00">bedtools intersect &#45;sorted </text>
<text text-anchor="middle" x="869" y="-204.54" font-family="Times,serif" font-size="14.00">&#45;a &lt;low map. regions&gt; &#45;b &lt;ClinVar low map.&gt;</text>
</g>
<!-- sorted low mappability ClinVar regions&#45;&gt;bedtools intersect &#45;u_0 -->
<g id="edge37" class="edge">
<title>sorted low mappability ClinVar regions&#45;&gt;bedtools intersect &#45;u_0</title>
<path fill="none" stroke="black" d="M619.66,-310.21C653.86,-304.93 691.45,-297.39 725,-286.74 753.2,-277.79 782.84,-264.18 807.85,-251.34"/>
<polygon fill="black" stroke="black" points="809.59,-254.38 816.84,-246.65 806.36,-248.17 809.59,-254.38"/>
<text text-anchor="middle" x="811" y="-275.54" font-family="Times,serif" font-size="14.00">ClinVar low map.</text>
</g>
<!-- 0 mappability regions BED -->
<g id="node24" class="node">
<title>0 mappability regions BED</title>
<ellipse fill="none" stroke="black" cx="634" cy="-1389.96" rx="108.58" ry="18"/>
<text text-anchor="middle" x="634" y="-1386.26" font-family="Times,serif" font-size="14.00">0 mappability regions BED</text>
</g>
<!-- bedtools complement&#45;&gt;0 mappability regions BED -->
<g id="edge23" class="edge">
<title>bedtools complement&#45;&gt;0 mappability regions BED</title>
<path fill="none" stroke="black" d="M635.52,-1444.77C635.29,-1436.75 635.02,-1427.01 634.76,-1418.03"/>
<polygon fill="black" stroke="black" points="638.26,-1417.89 634.48,-1407.99 631.26,-1418.08 638.26,-1417.89"/>
</g>
<!-- 0 mappability regions BED&#45;&gt;0 mappability regions bedGraph -->
<g id="edge30" class="edge">
<title>0 mappability regions BED&#45;&gt;0 mappability regions bedGraph</title>
<path fill="none" stroke="black" d="M634,-1371.76C634,-1360.12 634,-1344.51 634,-1331.2"/>
<polygon fill="black" stroke="black" points="637.5,-1331.14 634,-1321.14 630.5,-1331.14 637.5,-1331.14"/>
<text text-anchor="middle" x="694.5" y="-1342.76" font-family="Times,serif" font-size="14.00">awk &#39;{ print $0&quot;\t&quot;0 }&#39;</text>
</g>
<!-- grch38 .genome file -->
<g id="node26" class="node">
<title>grch38 .genome file</title>
<ellipse fill="none" stroke="black" cx="812" cy="-1535.96" rx="81.79" ry="18"/>
<text text-anchor="middle" x="812" y="-1532.26" font-family="Times,serif" font-size="14.00">grch38 .genome file</text>
</g>
<!-- grch38 .fai file&#45;&gt;grch38 .genome file -->
<g id="edge25" class="edge">
<title>grch38 .fai file&#45;&gt;grch38 .genome file</title>
<path fill="none" stroke="black" d="M812,-1604.76C812,-1593.12 812,-1577.51 812,-1564.2"/>
<polygon fill="black" stroke="black" points="815.5,-1564.14 812,-1554.14 808.5,-1564.14 815.5,-1564.14"/>
<text text-anchor="middle" x="837.5" y="-1575.76" font-family="Times,serif" font-size="14.00">cut &#45;f 1,2</text>
</g>
<!-- grch38 .genome file&#45;&gt;bedtools complement -->
<g id="edge26" class="edge">
<title>grch38 .genome file&#45;&gt;bedtools complement</title>
<path fill="none" stroke="black" d="M774.67,-1519.9C745.94,-1508.31 706.26,-1492.3 676.71,-1480.38"/>
<polygon fill="black" stroke="black" points="677.69,-1477 667.1,-1476.51 675.07,-1483.5 677.69,-1477"/>
</g>
<!-- &lt;=0.01 mappability regions bedGraph -->
<g id="node28" class="node">
<title>&lt;=0.01 mappability regions bedGraph</title>
<ellipse fill="none" stroke="black" cx="869" cy="-409.74" rx="147.57" ry="18"/>
<text text-anchor="middle" x="869" y="-406.04" font-family="Times,serif" font-size="14.00">&lt;=0.01 mappability regions bedGraph</text>
</g>
<!-- awk_again&#45;&gt;&lt;=0.01 mappability regions bedGraph -->
<g id="edge29" class="edge">
<title>awk_again&#45;&gt;&lt;=0.01 mappability regions bedGraph</title>
<path fill="none" stroke="black" d="M869,-464.55C869,-456.53 869,-446.79 869,-437.81"/>
<polygon fill="black" stroke="black" points="872.5,-437.77 869,-427.77 865.5,-437.77 872.5,-437.77"/>
</g>
<!-- sorted low mappability regions -->
<g id="node29" class="node">
<title>sorted low mappability regions</title>
<ellipse fill="none" stroke="black" cx="869" cy="-322.74" rx="120.78" ry="18"/>
<text text-anchor="middle" x="869" y="-319.04" font-family="Times,serif" font-size="14.00">sorted low mappability regions</text>
</g>
<!-- &lt;=0.01 mappability regions bedGraph&#45;&gt;sorted low mappability regions -->
<g id="edge31" class="edge">
<title>&lt;=0.01 mappability regions bedGraph&#45;&gt;sorted low mappability regions</title>
<path fill="none" stroke="black" d="M869,-391.54C869,-379.9 869,-364.29 869,-350.98"/>
<polygon fill="black" stroke="black" points="872.5,-350.92 869,-340.92 865.5,-350.92 872.5,-350.92"/>
<text text-anchor="middle" x="915" y="-362.54" font-family="Times,serif" font-size="14.00">sort &#45;k1,1 &#45;k2,2n</text>
</g>
<!-- sorted low mappability regions&#45;&gt;bedtools intersect &#45;u_0 -->
<g id="edge32" class="edge">
<title>sorted low mappability regions&#45;&gt;bedtools intersect &#45;u_0</title>
<path fill="none" stroke="black" d="M869,-304.4C869,-293.43 869,-278.68 869,-264.4"/>
<polygon fill="black" stroke="black" points="872.5,-264.15 869,-254.15 865.5,-264.15 872.5,-264.15"/>
<text text-anchor="middle" x="917" y="-275.54" font-family="Times,serif" font-size="14.00">low map. regions</text>
</g>
<!-- sorted low mappability regions\nin ClinVar annotations -->
<g id="node31" class="node">
<title>sorted low mappability regions\nin ClinVar annotations</title>
<ellipse fill="none" stroke="black" cx="869" cy="-113.87" rx="131.54" ry="26.74"/>
<text text-anchor="middle" x="869" y="-117.67" font-family="Times,serif" font-size="14.00">sorted low mappability regions</text>
<text text-anchor="middle" x="869" y="-102.67" font-family="Times,serif" font-size="14.00">in ClinVar annotations</text>
</g>
<!-- bedtools intersect &#45;u_0&#45;&gt;sorted low mappability regions\nin ClinVar annotations -->
<g id="edge33" class="edge">
<title>bedtools intersect &#45;u_0&#45;&gt;sorted low mappability regions\nin ClinVar annotations</title>
<path fill="none" stroke="black" d="M869,-177.55C869,-168.88 869,-159.65 869,-151.03"/>
<polygon fill="black" stroke="black" points="872.5,-150.87 869,-140.87 865.5,-150.87 872.5,-150.87"/>
</g>
<!-- BAM_PROC_3 -->
<!-- sorted low mappability regions\nin ClinVar annotations&#45;&gt;BAM_PROC_3 -->
<g id="edge36" class="edge">
<title>sorted low mappability regions\nin ClinVar annotations&#45;&gt;BAM_PROC_3</title>
<path fill="none" stroke="black" stroke-dasharray="1,5" d="M869,-86.97C869,-74.34 869,-59.14 869,-46.32"/>
<polygon fill="black" stroke="black" points="872.5,-46.22 869,-36.22 865.5,-46.22 872.5,-46.22"/>
<text text-anchor="middle" x="976.5" y="-57.8" font-family="Times,serif" font-size="14.00">to BAM processing and postprocessing</text>
</g>
<!-- expanded ClinVar regions BED&#45;&gt;expanded sorted ClinVar regions BED -->
<g id="edge35" class="edge">
<title>expanded ClinVar regions BED&#45;&gt;expanded sorted ClinVar regions BED</title>
<path fill="none" stroke="black" d="M375,-1197.76C375,-1186.12 375,-1170.51 375,-1157.2"/>
<polygon fill="black" stroke="black" points="378.5,-1157.14 375,-1147.14 371.5,-1157.14 378.5,-1157.14"/>
<text text-anchor="middle" x="421" y="-1168.76" font-family="Times,serif" font-size="14.00">sort &#45;k1,1 &#45;k2,2n</text>
</g>
</g>
</svg>