Skip to content
Snippets Groups Projects
Verified Commit 0d79273c authored by Laurent Modolo's avatar Laurent Modolo
Browse files

emase: add diploid_genome process

parent 1db31f43
No related branches found
No related tags found
No related merge requests found
version = "0.10.16"
container_url = "lbmc/emase:${version}"
params.personalised_transcriptome = ""
process personalised_transcriptome {
params.diploid_genome = "-x"
params.diploid_genome_out = "-x"
process diploid_genome {
container = "${container_url}"
label "big_mem_mono_cpus"
tag "$file_id"
if (params.diploid_genome_out != "") {
publishDir "results/${params.diploid_genome_out}", mode: 'copy'
}
input:
tuple val(file_id), path(fasta)
tuple val(gtf_id), path(gtf)
tuple val(genome_a), path(fasta_a)
tuple val(genome_b), path(fasta_b)
output:
tuple val(file_id), path("${fasta.simpleName}.*"), emit: index
tuple val(file_id), path("*_bwa_report.txt"), emit: report
tuple val(file_id), path(".fa"), emit: index
script:
"""
prepare-emase ${personalised_transcriptome} -G ${REF_FASTA} -g ${REF_GTF} -o ${REF_DIR} -m --no-bowtie-index
// ${REF_DIR}/emase.transcriptome.fa
// ${REF_DIR}/emase.transcriptome.info
// ${REF_DIR}/emase.gene2transcripts.tsv
prepare-emase -G ${SAMPLE_DIR}/L.fa,${SAMPLE_DIR}/R.fa -s L,R -o ${SAMPLE_DIR}
prepare-emase -G ${fasta_a},${fasta_b} -s ${genome_a},${genome_b} ${params.diploid_genome}
"""
}
\ No newline at end of file
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