Skip to content

Commit

Permalink
use script to prep SRA metadata for EAGER tsv
Browse files Browse the repository at this point in the history
  • Loading branch information
ktmeaton committed May 6, 2020
1 parent 58eaf5b commit f82390f
Showing 1 changed file with 8 additions and 8 deletions.
16 changes: 8 additions & 8 deletions pipeline.nf
Original file line number Diff line number Diff line change
Expand Up @@ -258,7 +258,7 @@ if( (params.sqlite || ( params.ncbimeta_update && params.ncbimeta_annot) ) && !p
file sqlite from ch_sqlite
output:
file params.file_assembly_for_download_ftp into ch_assembly_for_download_ftp
file params.file_sra_for_download_acc into ch_sra_for_download_acc
file params.eager_tsv into ch_sra_tsv_for_eager

// Shell script to execute
script:
Expand All @@ -272,13 +272,13 @@ if( (params.sqlite || ( params.ncbimeta_update && params.ncbimeta_annot) ) && !p
echo \$asm_ftp >> ${params.file_assembly_for_download_ftp}
fi;
done;
# Select the SRA Run Accessions
sqlite3 ${sqlite} ${params.sqlite_select_command_sra} | grep . | head -n ${params.max_datasets} | sed -E -e 's/ |;/\\n/g' | while read line;
do
if [[ ! -z \$line ]]; then
echo \$line >> ${params.file_sra_for_download_acc}
fi;
done
# Extract SRA Metadata for EAGER tsv
${params.scriptdir}/sqlite_EAGER_tsv.py \
--database ${sqlite} \
--query ${params.sqlite_select_command_sra} \
--organism ${params.eager_organism} \
--max-datasets ${params.max_datasets} \
--output ${params.eager_tsv}
"""
}

Expand Down

0 comments on commit f82390f

Please sign in to comment.