realign.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="795pt" height="564pt"
 viewBox="0.00 0.00 794.59 564.00" 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 560)">
<title>workflow</title>
<polygon fill="white" stroke="transparent" points="-4,4 -4,-560 790.59,-560 790.59,4 -4,4"/>
<!-- Reference FASTA -->
<g id="node1" class="node">
<title>Reference FASTA</title>
<polygon fill="none" stroke="black" points="577.62,-358 533.31,-376 444.69,-376 400.38,-358 444.69,-340 533.31,-340 577.62,-358"/>
<text text-anchor="middle" x="489" y="-354.3" font-family="Times,serif" font-size="14.00">Reference FASTA</text>
</g>
<!-- Bismark -->
<g id="node5" class="node">
<title>Bismark</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="203,-196 10.79,-178 203,-160 395.21,-178 203,-196"/>
<text text-anchor="middle" x="203" y="-174.3" font-family="Times,serif" font-size="14.00">bismark &#45;1 &lt;read 1s&gt; &#45;2 &lt;read 2s&gt;</text>
</g>
<!-- Reference FASTA&#45;&gt;Bismark -->
<g id="edge10" class="edge">
<title>Reference FASTA&#45;&gt;Bismark</title>
<path fill="none" stroke="black" d="M484,-339.89C478.88,-324.67 469.73,-302.48 456,-287 409.45,-234.49 332.99,-206.68 276.29,-192.51"/>
<polygon fill="black" stroke="black" points="276.85,-189.05 266.31,-190.11 275.21,-195.85 276.85,-189.05"/>
</g>
<!-- bwa&#45;meth -->
<g id="node7" class="node">
<title>bwa&#45;meth</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="600,-196 413.32,-178 600,-160 786.68,-178 600,-196"/>
<text text-anchor="middle" x="600" y="-174.3" font-family="Times,serif" font-size="14.00">bwameth.py &lt;read 1s&gt; &lt;read 2s&gt;</text>
</g>
<!-- Reference FASTA&#45;&gt;bwa&#45;meth -->
<g id="edge11" class="edge">
<title>Reference FASTA&#45;&gt;bwa&#45;meth</title>
<path fill="none" stroke="black" d="M490.81,-339.87C494.21,-314.48 503.47,-266.27 527,-233 536.63,-219.38 550.78,-207.79 564.05,-198.89"/>
<polygon fill="black" stroke="black" points="566.06,-201.76 572.6,-193.43 562.29,-195.86 566.06,-201.76"/>
</g>
<!-- Bismark Index -->
<g id="node10" class="node">
<title>Bismark Index</title>
<ellipse fill="none" stroke="black" cx="297" cy="-251" rx="63.89" ry="18"/>
<text text-anchor="middle" x="297" y="-247.3" font-family="Times,serif" font-size="14.00">Bismark Index</text>
</g>
<!-- Reference FASTA&#45;&gt;Bismark Index -->
<g id="edge8" class="edge">
<title>Reference FASTA&#45;&gt;Bismark Index</title>
<path fill="none" stroke="black" d="M437.7,-342.8C379.81,-326.76 293.57,-302.75 293,-302 288.04,-295.46 287.28,-286.96 288.29,-278.83"/>
<polygon fill="black" stroke="black" points="291.74,-279.41 290.34,-268.91 284.89,-277.99 291.74,-279.41"/>
<text text-anchor="middle" x="374.5" y="-290.8" font-family="Times,serif" font-size="14.00">bismark_genome_preparation</text>
</g>
<!-- bwa&#45;meth Index -->
<g id="node11" class="node">
<title>bwa&#45;meth Index</title>
<ellipse fill="none" stroke="black" cx="605" cy="-251" rx="68.79" ry="18"/>
<text text-anchor="middle" x="605" y="-247.3" font-family="Times,serif" font-size="14.00">bwa&#45;meth Index</text>
</g>
<!-- Reference FASTA&#45;&gt;bwa&#45;meth Index -->
<g id="edge9" class="edge">
<title>Reference FASTA&#45;&gt;bwa&#45;meth Index</title>
<path fill="none" stroke="black" d="M507.9,-339.9C527.34,-322.3 557.75,-294.77 579.35,-275.22"/>
<polygon fill="black" stroke="black" points="581.81,-277.71 586.88,-268.4 577.11,-272.52 581.81,-277.71"/>
<text text-anchor="middle" x="615" y="-290.8" font-family="Times,serif" font-size="14.00">bwameth.py index</text>
</g>
<!-- Input BAMs -->
<g id="node2" class="node">
<title>Input BAMs</title>
<polygon fill="none" stroke="black" points="237.58,-538 205.29,-556 140.71,-556 108.42,-538 140.71,-520 205.29,-520 237.58,-538"/>
<text text-anchor="middle" x="173" y="-534.3" font-family="Times,serif" font-size="14.00">Input BAMs</text>
</g>
<!-- Collated BAMs -->
<g id="node8" class="node">
<title>Collated BAMs</title>
<ellipse fill="none" stroke="black" cx="173" cy="-451" rx="66.89" ry="18"/>
<text text-anchor="middle" x="173" y="-447.3" font-family="Times,serif" font-size="14.00">Collated BAMs</text>
</g>
<!-- Input BAMs&#45;&gt;Collated BAMs -->
<g id="edge1" class="edge">
<title>Input BAMs&#45;&gt;Collated BAMs</title>
<path fill="none" stroke="black" d="M173,-519.8C173,-508.16 173,-492.55 173,-479.24"/>
<polygon fill="black" stroke="black" points="176.5,-479.18 173,-469.18 169.5,-479.18 176.5,-479.18"/>
<text text-anchor="middle" x="217.5" y="-490.8" font-family="Times,serif" font-size="14.00">samtools collate</text>
</g>
<!-- Final Bismark BAMs -->
<g id="node3" class="node">
<title>Final Bismark BAMs</title>
<polygon fill="none" stroke="black" points="331.79,-36 126.91,-36 74.21,0 279.09,0 331.79,-36"/>
<text text-anchor="middle" x="203" y="-14.3" font-family="Times,serif" font-size="14.00">Final Bismark BAMs</text>
</g>
<!-- Final bwa&#45;meth BAMs -->
<g id="node4" class="node">
<title>Final bwa&#45;meth BAMs</title>
<polygon fill="none" stroke="black" points="736.04,-123 519.63,-123 463.96,-87 680.37,-87 736.04,-123"/>
<text text-anchor="middle" x="600" y="-101.3" font-family="Times,serif" font-size="14.00">Final bwa&#45;meth BAMs</text>
</g>
<!-- Realigned Bismark BAMs -->
<g id="node12" class="node">
<title>Realigned Bismark BAMs</title>
<ellipse fill="none" stroke="black" cx="203" cy="-105" rx="105.08" ry="18"/>
<text text-anchor="middle" x="203" y="-101.3" font-family="Times,serif" font-size="14.00">Realigned Bismark BAMs</text>
</g>
<!-- Bismark&#45;&gt;Realigned Bismark BAMs -->
<g id="edge12" class="edge">
<title>Bismark&#45;&gt;Realigned Bismark BAMs</title>
<path fill="none" stroke="black" d="M203,-159.81C203,-151.79 203,-142.05 203,-133.07"/>
<polygon fill="black" stroke="black" points="206.5,-133.03 203,-123.03 199.5,-133.03 206.5,-133.03"/>
</g>
<!-- samtools fastq -->
<g id="node6" class="node">
<title>samtools fastq</title>
<polygon fill="#d9d9d9" stroke="#d9d9d9" points="173,-396 0,-358 173,-320 346,-358 173,-396"/>
<text text-anchor="middle" x="173" y="-361.8" font-family="Times,serif" font-size="14.00">samtools fastq &#45;1 &lt;read 1s&gt; </text>
<text text-anchor="middle" x="173" y="-346.8" font-family="Times,serif" font-size="14.00">&#45;2 &lt;read 2s&gt; &lt;BAM file&gt;</text>
</g>
<!-- FASTQs -->
<g id="node9" class="node">
<title>FASTQs</title>
<ellipse fill="none" stroke="black" cx="173" cy="-251" rx="42.49" ry="18"/>
<text text-anchor="middle" x="173" y="-247.3" font-family="Times,serif" font-size="14.00">FASTQs</text>
</g>
<!-- samtools fastq&#45;&gt;FASTQs -->
<g id="edge3" class="edge">
<title>samtools fastq&#45;&gt;FASTQs</title>
<path fill="none" stroke="black" d="M173,-319.98C173,-306.56 173,-291.65 173,-279.22"/>
<polygon fill="black" stroke="black" points="176.5,-279.04 173,-269.04 169.5,-279.04 176.5,-279.04"/>
</g>
<!-- bwa&#45;meth&#45;&gt;Final bwa&#45;meth BAMs -->
<g id="edge14" class="edge">
<title>bwa&#45;meth&#45;&gt;Final bwa&#45;meth BAMs</title>
<path fill="none" stroke="black" d="M600,-159.81C600,-151.79 600,-142.05 600,-133.07"/>
<polygon fill="black" stroke="black" points="603.5,-133.03 600,-123.03 596.5,-133.03 603.5,-133.03"/>
</g>
<!-- Collated BAMs&#45;&gt;samtools fastq -->
<g id="edge2" class="edge">
<title>Collated BAMs&#45;&gt;samtools fastq</title>
<path fill="none" stroke="black" d="M173,-432.88C173,-425.25 173,-415.83 173,-406.31"/>
<polygon fill="black" stroke="black" points="176.5,-406.02 173,-396.02 169.5,-406.02 176.5,-406.02"/>
</g>
<!-- FASTQs&#45;&gt;Bismark -->
<g id="edge4" class="edge">
<title>FASTQs&#45;&gt;Bismark</title>
<path fill="none" stroke="black" d="M180.11,-233.17C183.66,-224.77 188.03,-214.42 192,-205.03"/>
<polygon fill="black" stroke="black" points="195.27,-206.28 195.94,-195.71 188.83,-203.55 195.27,-206.28"/>
</g>
<!-- FASTQs&#45;&gt;bwa&#45;meth -->
<g id="edge5" class="edge">
<title>FASTQs&#45;&gt;bwa&#45;meth</title>
<path fill="none" stroke="black" d="M204.03,-238.67C210.57,-236.57 217.46,-234.55 224,-233 316.28,-211.18 423.53,-196.83 499.81,-188.45"/>
<polygon fill="black" stroke="black" points="500.32,-191.92 509.89,-187.36 499.57,-184.96 500.32,-191.92"/>
</g>
<!-- Bismark Index&#45;&gt;Bismark -->
<g id="edge6" class="edge">
<title>Bismark Index&#45;&gt;Bismark</title>
<path fill="none" stroke="black" d="M275.67,-233.89C262.54,-223.97 245.55,-211.14 231.3,-200.38"/>
<polygon fill="black" stroke="black" points="233.37,-197.55 223.28,-194.32 229.15,-203.14 233.37,-197.55"/>
</g>
<!-- bwa&#45;meth Index&#45;&gt;bwa&#45;meth -->
<g id="edge7" class="edge">
<title>bwa&#45;meth Index&#45;&gt;bwa&#45;meth</title>
<path fill="none" stroke="black" d="M603.79,-232.81C603.22,-224.79 602.54,-215.05 601.91,-206.07"/>
<polygon fill="black" stroke="black" points="605.39,-205.76 601.2,-196.03 598.41,-206.25 605.39,-205.76"/>
</g>
<!-- Realigned Bismark BAMs&#45;&gt;Final Bismark BAMs -->
<g id="edge13" class="edge">
<title>Realigned Bismark BAMs&#45;&gt;Final Bismark BAMs</title>
<path fill="none" stroke="black" d="M203,-86.8C203,-75.16 203,-59.55 203,-46.24"/>
<polygon fill="black" stroke="black" points="206.5,-46.18 203,-36.18 199.5,-46.18 206.5,-46.18"/>
<text text-anchor="middle" x="260" y="-57.8" font-family="Times,serif" font-size="14.00">deduplicate_bismark</text>
</g>
</g>
</svg>