Skip to content

Commit

Permalink
fix bug for contig stats final
Browse files Browse the repository at this point in the history
  • Loading branch information
SilasK committed Jul 20, 2023
1 parent be4d888 commit 6affeb7
Showing 1 changed file with 21 additions and 18 deletions.
39 changes: 21 additions & 18 deletions workflow/rules/assemble.smk
Original file line number Diff line number Diff line change
Expand Up @@ -492,24 +492,6 @@ rule rename_contigs:
" minscaf={params.minlength} &> {log} "


rule calculate_contigs_stats:
input:
"{sample}/assembly/{sample}_contigs.fasta",
output:
"{sample}/assembly/contig_stats/final_contig_stats.txt",
conda:
"../envs/required_packages.yaml"
log:
"{sample}/logs/assembly/post_process/contig_stats_final.log",
threads: 1
resources:
mem=1,
time=config["runtime"]["simplejob"],
shell:
"stats.sh in={input} format=3 out={output} &> {log}"




if config["filter_contigs"]:

Expand Down Expand Up @@ -623,6 +605,27 @@ rule finalize_contigs:
os.symlink(os.path.relpath(input[0], os.path.dirname(output[0])), output[0])




rule calculate_contigs_stats:
input:
"{sample}/{sample}_contigs.fasta",
output:
"{sample}/assembly/contig_stats/final_contig_stats.txt",
conda:
"../envs/required_packages.yaml"
log:
"{sample}/logs/assembly/post_process/contig_stats_final.log",
threads: 1
resources:
mem=1,
time=config["runtime"]["simplejob"],
shell:
"stats.sh in={input} format=3 out={output} &> {log}"




# generalized rule so that reads from any "sample" can be aligned to contigs from "sample_contigs"
rule align_reads_to_final_contigs:
input:
Expand Down

0 comments on commit 6affeb7

Please sign in to comment.