diff --git a/Snakemake/1000RNASeq_chicken/ASE/rules/STAR.smk b/Snakemake/1000RNASeq_chicken/ASE/rules/STAR.smk index 8df1b00bceefe000f49f35a935f2e91d269a062f..685105b93d03867f8e6f65269c9d87744c473568 100644 --- a/Snakemake/1000RNASeq_chicken/ASE/rules/STAR.smk +++ b/Snakemake/1000RNASeq_chicken/ASE/rules/STAR.smk @@ -12,7 +12,7 @@ wildcard_constraints: STAR_INDEX_FILE=['genomeParameters.txt', 'chrName.txt', 'chrLength.txt', 'chrStart.txt', 'chrNameLength.txt', 'exonGeTrInfo.tab', 'geneInfo.tab', 'transcriptInfo.tab', 'exonInfo.tab', 'sjdbList.fromGTF.out.tab', 'sjdbInfo.txt', 'sjdbList.out.tab', 'Genome', 'SA', 'SAindex'] -rule STAR_Index: +rule STAR_Index_1: input: fasta = 'Results/genomeMasked/' + os.path.splitext(os.path.basename(config['fasta_ref']))[0] + '_masked.fa', gtf = os.path.abspath(config["gtf_ref"]) @@ -157,7 +157,7 @@ def get_properly_paired_bam(wildcards): return "Results/STAR_Aln_2/" + wildcards.sample + "_rg_genomic.bam" -rule rmdup: +rule RmDuplicates: input: bam = get_properly_paired_bam output: