Skip to content
Snippets Groups Projects
Commit 44e79e5a authored by Valentina Galata's avatar Valentina Galata
Browse files

updated targets in mapping step (issue #66)

parent d1e53f46
No related branches found
No related tags found
1 merge request!76Merge "cleanup" branch with "master" branch
......@@ -11,7 +11,7 @@ rule MAPPING:
),
metat=expand(
os.path.join(RESULTS_DIR, "mapping/metat/{rtype_tool}.bam"),
rtype_tool=["{r}/{t}/ASSEMBLY_sr".format(r=r, t=t) for r, t in READ_ASSEMBLERS]
rtype_tool=["{r}/{t}/ASSEMBLY.FILTERED.sr".format(r=r, t=t) for r, t in READ_ASSEMBLERS]
) if "metat" in META_TYPES else []
output:
touch("status/mapping.done")
......@@ -24,7 +24,7 @@ rule COVERAGE:
),
metat=expand(
os.path.join(RESULTS_DIR, "mapping/metat/{rtype_tool}.cov.ave.txt"),
rtype_tool=["{r}/{t}/ASSEMBLY_sr".format(r=r, t=t) for r, t in READ_ASSEMBLERS]
rtype_tool=["{r}/{t}/ASSEMBLY.FILTERED.sr".format(r=r, t=t) for r, t in READ_ASSEMBLERS]
) if "metat" in META_TYPES else []
output:
touch("status/coverage.done")
......@@ -33,12 +33,12 @@ rule MAPPABILITY:
input:
metag=expand(
os.path.join(RESULTS_DIR, "mapping/metag/{rtype_tool}.{ext}"),
rtype_tool=["{r}/{t}/ASSEMBLY_{r}".format(r=r, t=t) for r, t in READ_ASSEMBLERS],
rtype_tool=["{r}/{t}/ASSEMBLY.FILTERED.{r}".format(r=r, t=t) for r, t in READ_ASSEMBLERS],
ext=["flagstat.txt", "idxstats.txt", "unique.txt"]
),
metat=expand(
os.path.join(RESULTS_DIR, "mapping/metat/{rtype_tool}.{ext}"),
rtype_tool=["{r}/{t}/ASSEMBLY_sr".format(r=r, t=t) for r, t in READ_ASSEMBLERS],
rtype_tool=["{r}/{t}/ASSEMBLY.FILTERED.sr".format(r=r, t=t) for r, t in READ_ASSEMBLERS],
ext=["flagstat.txt", "idxstats.txt", "unique.txt"]
) if "metat" in META_TYPES else []
output:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment