diff --git a/Snakemake/1000RNASeq_chicken/ASE/rules/STAR.smk b/Snakemake/1000RNASeq_chicken/ASE/rules/STAR.smk index 315f8d734214eee7caaccf74dc16820f0fe66dbd..053841c211ff08cb2679abcea287cc43f7a5e003 100644 --- a/Snakemake/1000RNASeq_chicken/ASE/rules/STAR.smk +++ b/Snakemake/1000RNASeq_chicken/ASE/rules/STAR.smk @@ -188,7 +188,8 @@ rule SplitNCigarReads: fasta = 'Results/genomeMasked/' + os.path.splitext(os.path.basename(config['fasta_ref']))[0] + '_masked.fa', idx = 'Results/genomeMasked/' + os.path.splitext(os.path.basename(config['fasta_ref']))[0] + '_masked.fa.fai', dict = 'Results/genomeMasked/' + os.path.splitext(os.path.basename(config['fasta_ref']))[0] + '_masked.dict', - bam = "Results/STAR_Aln_2/{sample}_rg_genomic{properlyPaired}rmdup_uniq.bam" + bam = "Results/STAR_Aln_2/{sample}_rg_genomic{properlyPaired}rmdup_uniq.bam", + bai = "Results/STAR_Aln_2/{sample}_rg_genomic{properlyPaired}rmdup_uniq.bai" output: bam = "Results/STAR_Aln_2/{sample}_rg_genomic{properlyPaired}rmdup_uniq_split.bam" params: diff --git a/Snakemake/1000RNASeq_chicken/ASE/rules/basics.smk b/Snakemake/1000RNASeq_chicken/ASE/rules/basics.smk index 3a7d39750726b5937600c5ab6901c879293eee7b..f37be0839f8295cb05ab04d52f67d7210fd702ce 100644 --- a/Snakemake/1000RNASeq_chicken/ASE/rules/basics.smk +++ b/Snakemake/1000RNASeq_chicken/ASE/rules/basics.smk @@ -52,5 +52,5 @@ rule indexFasta: shell: """ samtools faidx {input.ref} - java -Xmx{params.mem} -jar {params.jar} R={input.ref} O={output.dict} + java -Xmx{params.mem} -jar {params.jar} CreateSequenceDictionary R={input.ref} O={output.dict} """