Skip to content

Commit

Permalink
only start purge_dups2 after succesfully completed purge_dups
Browse files Browse the repository at this point in the history
  • Loading branch information
LiaOb21 committed Jan 4, 2024
1 parent 099d61b commit ca943d5
Showing 1 changed file with 5 additions and 1 deletion.
6 changes: 5 additions & 1 deletion workflow/rules/purge_dups2.smk
Original file line number Diff line number Diff line change
@@ -1,6 +1,10 @@
# This rule uses purge_dups purge haplotigs and overlaps from the alternate assembly produced by hifiasm

rule run_purge_dups_alt:
depends:
"results/purge_dups_alt/hifiasm_p_purged.fa",
"results/purge_dups/hap.fa",
"results/purge_dups/hist.out.png"
input:
reads = "results/reads/hifi/hifi.fastq.gz",
fasta = "results/hifiasm/hifiasm.asm.a_ctg.fa",
Expand All @@ -26,7 +30,7 @@ rule run_purge_dups_alt:
"logs/purge_dups_alt.log"
shell:
"""
cat {input.fasta} {input.hap_fa_in} > merged.fa >> {log} 2>&1
cat {input.fasta} {input.hap_fa_in} > merged.fa
minimap2 -xasm20 merged.fa {input.reads} -t {config[minimap2][t]} | gzip -c - > hifi_vs_hifiasm_contigs.paf.gz >> {log} 2>&1
pbcstat hifi_vs_hifiasm_contigs.paf.gz >> {log} 2>&1
calcuts PB.stat > cutoffs 2>calcults.log
Expand Down

0 comments on commit ca943d5

Please sign in to comment.