millironx/taxprofiler
1
0
Fork 0
mirror of https://github.com/MillironX/taxprofiler.git synced 2024-11-22 10:29:54 +00:00
Code Issues Releases Packages Activity

Fix MALT save_alignemnts reads and some clean up

Browse source
This commit is contained in:
James Fellows Yates 2022-06-03 22:29:04 +02:00
parent d48b3be5a7
commit 9462032d00
4 changed files with 9 additions and 8 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
conf/test.config
View file

@ -47,4 +47,7 @@ process {
withName: MALT_RUN {
maxForks = 1
}
withName: MEGAN_RMA2INFO {
maxForks = 1
}
}

2
nextflow.config
View file

@ -196,7 +196,7 @@ profiles {
test { includeConfig 'conf/test.config' }
test_full { includeConfig 'conf/test_full.config' }
test_noprofiling { includeConfig 'conf/test_noprofiling.config' }
test_nopreprocessing { includeConfig 'conf/test_preprocessing.config' }
test_nopreprocessing { includeConfig 'conf/test_nopreprocessing.config' }
}
// Load igenomes.config if required

2
subworkflows/local/longread_preprocessing.nf
View file

@ -48,7 +48,7 @@ workflow LONGREAD_PREPROCESSING {
}
FASTQC_PROCESSED ( ch_processed_reads.dump(tag: "filtlong") )
FASTQC_PROCESSED ( ch_processed_reads )
ch_multiqc_files = ch_multiqc_files.mix( FASTQC_PROCESSED.out.zip )
emit:

10
subworkflows/local/profiling.nf
View file

@ -66,19 +66,17 @@ workflow PROFILING {
.filter { it[0]['instrument_platform'] == 'ILLUMINA' }
.map {
meta, reads, db_meta, db ->
def sam_format = params.malt_save_reads ? ' --alignments' : ""
// TODO No MALT SAM?
// TODO check all aligned reads published
// TODO try turning on/off aligned reads
def sam_format = params.malt_save_reads ? ' --alignments ./ -za false' : ""
// TODO wut? [9a/a441d6] Submitted process > NFCORE_TAXPROFILER:TAXPROFILER:PROFILING:MALT_RUN (null)
def temp_meta = [ id: meta['db_name'] ]
def new_db_meta = db_meta.clone()
new_db_meta['db_params'] = db_meta['db_params'] + sam_format
def new_meta = temp_meta + new_db_meta
new_meta['id'] = new_meta['db_name']
[ new_meta, reads, db ]
}
.groupTuple(by: [0,2])
.dump(tag: "into_malt")
.multiMap {
it ->
reads: [ it[0], it[1].flatten() ]
@ -192,7 +190,7 @@ workflow PROFILING {
// this will replace output file!
ch_diamond_reads_format = params.diamond_save_reads ? 'sam' : params.diamond_output_format
DIAMOND_BLASTX ( ch_input_for_diamond.reads, ch_input_for_diamond.db, params.diamond_output_format, [] )
DIAMOND_BLASTX ( ch_input_for_diamond.reads, ch_input_for_diamond.db, ch_diamond_reads_format , [] )
ch_versions = ch_versions.mix( DIAMOND_BLASTX.out.versions.first() )
ch_raw_profiles = ch_raw_profiles.mix( DIAMOND_BLASTX.out.tsv )