We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 3881415 commit 8780981Copy full SHA for 8780981
src/Snakefiles/8-hicPipeline.sm
@@ -103,7 +103,7 @@ cat > ./run_mashmap.sh <<EOF
103
set -e
104
echo "---Running MashMap"
105
{params.MASHMAP} -r ../{input.unitigs_hpc} -q ../{input.unitigs_hpc} -t {threads} -f none --pi {params.MASHMAP_DIV} -s 10000 -o mashmap.out
106
-cat mashmap.out |awk '{{if (\$NF > 99 && \$4-\$3 > 100000 && \$1 != \$6) print \$1"\\t"\$6"\\t"\$4-\$3}}'|sort |uniq > ../{output.unitigs_matches}
+cat mashmap.out |awk '{{if (\$4-\$3 > 100000 && \$1 != \$6) print \$1"\\t"\$6"\\t"\$4-\$3}}'|sort |uniq > ../{output.unitigs_matches}
107
EOF
108
109
chmod +x ./run_mashmap.sh
0 commit comments