Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Blaise LI
bioinfo_utils
Commits
783a77ac
Commit
783a77ac
authored
Feb 10, 2020
by
Blaise Li
Browse files
Compute TPM and "efficiencies" on protein_coding.
parent
23c1b898
Changes
4
Hide whitespace changes
Inline
Side-by-side
Degradome-seq/Degradome-seq.snakefile
View file @
783a77ac
...
...
@@ -130,7 +130,7 @@ rule all:
# lib=LIBS, biotype=["alltypes"], orientation=ORIENTATIONS),
expand(
OPJ(counts_dir, "diff_%s" % eff_name, "{contrast}", "{orientation}_{biotype}", "{contrast}_diff_%s.txt" % eff_name),
contrast=DD_CONTRASTS, orientation=["fwd"], biotype=["alltypes"]),
contrast=DD_CONTRASTS, orientation=["fwd"], biotype=["alltypes"
, "protein_coding"
]),
include: ensure_relative(irules["link_raw_data"], workflow.basedir)
...
...
PRO-seq/PRO-seq.snakefile
View file @
783a77ac
...
...
@@ -357,7 +357,7 @@ rule all:
expand(OPJ(
"{trimmer}", aligner, "mapped_C_elegans", "{counter}",
"all_on_C_elegans", "{biotype}_{orientation}_TPM.txt"),
trimmer=TRIMMERS, counter=COUNTERS, biotype=["alltypes"], orientation=ORIENTATIONS),
trimmer=TRIMMERS, counter=COUNTERS, biotype=["alltypes"
, "protein_coding"
], orientation=ORIENTATIONS),
#expand(OPJ("{trimmer}", "figures", aligner, "{lib}_mean", "{orientation}_on_merged_isolated_%d_{biotype}_min_%d_meta_profile.pdf" % (MIN_DIST, META_MIN_LEN)), trimmer=TRIMMERS, lib=LIBS, orientation=["all"], biotype=["protein_coding"]),
#expand(OPJ("{trimmer}", "figures", aligner, "{lib}_mean", "{orientation}_on_merged_isolated_%d_{biotype}_min_%d_meta_profile.pdf" % (MIN_DIST, META_MIN_LEN)), trimmer=TRIMMERS, lib=LIBS, orientation=["all"], biotype=METAGENE_BIOTYPES),
# TODO: Add metagene profiles similar to small RNA-seq
...
...
RNA_Seq_Cecere/RNA-seq.snakefile
View file @
783a77ac
...
...
@@ -318,7 +318,7 @@ counts_files = [
OPJ(aligner, f"mapped_{genome}", "{counter}",
"all_{mapped_type}", "{biotype}_{orientation}_TPM.txt"),
counter=COUNTERS, mapped_type=[f"on_{genome}", f"unique_on_{genome}"],
biotype=["alltypes"], orientation=ORIENTATIONS),
biotype=["alltypes"
, "protein_coding"
], orientation=ORIENTATIONS),
expand(
OPJ(aligner, f"mapped_{genome}", "{counter}",
"all_{mapped_type}", "{biotype}_{orientation}_RPK.txt"),
...
...
Ribo-seq/Ribo-seq.snakefile
View file @
783a77ac
...
...
@@ -579,7 +579,7 @@ rule all:
# translation_efficiency
expand(
OPJ(counts_dir, "{lib}_mean_{read_type}_on_%s" % genome, "{lib}_{biotype}_{orientation}_%s.txt" % eff_name),
lib=EFF_LIBS, read_type=["RPF"], biotype=["alltypes"], orientation=["fwd"]),
lib=EFF_LIBS, read_type=["RPF"], biotype=["alltypes"
, "protein_coding"
], orientation=["fwd"]),
# DESeq2 results
expand(
OPJ(
...
...
@@ -599,7 +599,7 @@ rule all:
counts_dir,
"diff_%s_{read_type}" % eff_name,
"{contrast}", "{orientation}_{biotype}", "{contrast}_diff_%s.txt" % eff_name),
read_type=["RPF"], contrast=DT_CONTRASTS, orientation=["fwd"], biotype=["alltypes"]),
read_type=["RPF"], contrast=DT_CONTRASTS, orientation=["fwd"], biotype=["alltypes"
, "protein_coding"
]),
# rule future_all:
# meta_profiles,
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment