Skip to content

trees by Fred Lemoine : replace the code inside main.nf and debugged

This

    fasta=Reformat(aa_tsv_ch2)
    fastagroups=DefineGroups(fasta,parseDb_filtering.out.select_ch).flatten()
    align=Align(fastagroups,phylo_tree_heavy)
    // Keeps only alignments with at least 3 sequences
    filtered=NbSequences(align).filter{it[0].toInteger()>=3}.map{it->it[1]}
    tree=Tree(filtered,phylo_tree_model_file)
    itolmeta=ProcessMeta(meta_file)
    ITOL(tree[0],itolmeta,phylo_tree_itolkey)

has been removed from main.nf and must be replaced just before

    backup(
        config_file, 
        log_file
    )