From 978feeaa394fc48b08fe63849b3e1a2bd3e6a28d Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?=E2=80=98ens-ftricomi=E2=80=99?= Date: Thu, 15 Aug 2024 15:26:17 +0100 Subject: [PATCH] import changes --- .../transcriptomic_alignment/index_genome.nf | 22 ++++++++++++------- 1 file changed, 14 insertions(+), 8 deletions(-) diff --git a/pipelines/nextflow/modules/transcriptomic_alignment/index_genome.nf b/pipelines/nextflow/modules/transcriptomic_alignment/index_genome.nf index af203c2..45ddebc 100644 --- a/pipelines/nextflow/modules/transcriptomic_alignment/index_genome.nf +++ b/pipelines/nextflow/modules/transcriptomic_alignment/index_genome.nf @@ -47,11 +47,15 @@ process INDEX_GENOME { // Initialize variables def numberOfReferences = 0 def genomeLength = 0 - // Read the FASTA file line by line - genomefilePath.eachLine { line -> - if (line.startsWith('>')) { - numberOfReferences++ - } else { + def retryCount = 0 + def maxRetries = 3 + def filesValid = false + while (!filesValid && retryCount < maxRetries) { + // Read the FASTA file line by line + genomefilePath.eachLine { line -> + if (line.startsWith('>')) { + numberOfReferences++ + } else { genomeLength += line.trim().length() } } @@ -66,6 +70,10 @@ process INDEX_GENOME { // Calculate the genome length (excluding header lines) //def genomeLength = fastaContent.split('\n').findAll { !it.startsWith('>') }.join('').length() + if (genomeLength != 0) { + filesValid = true + } + } // Define read length (you may need to adjust this based on your data) def readLength = 100 @@ -97,6 +105,4 @@ process INDEX_GENOME { """ } - } - - +} \ No newline at end of file