diff --git a/Min_test_data/nextflow_results/MinP1831_erik_protein.fa b/Min_test_data/nextflow_results/MinP1831_erik_protein.fa new file mode 120000 index 0000000000000000000000000000000000000000..55f9c3b6e2a11129eada9ff014f51b917faa6b54 --- /dev/null +++ b/Min_test_data/nextflow_results/MinP1831_erik_protein.fa @@ -0,0 +1 @@ +/home/nibio/Documents/Git/DivGene-genome-annotation/work/b2/4d0430cc99d9679f13bbc44fdb2b14/MinP1831_erik_protein.fa \ No newline at end of file diff --git a/Min_test_data/nextflow_results/emapper/MinP1831_erik_em.emapper.annotations b/Min_test_data/nextflow_results/emapper/MinP1831_erik_em.emapper.annotations new file mode 120000 index 0000000000000000000000000000000000000000..a16793174a65face9e9e1e43408a72ea007cca01 --- /dev/null +++ b/Min_test_data/nextflow_results/emapper/MinP1831_erik_em.emapper.annotations @@ -0,0 +1 @@ +/home/nibio/Documents/Git/DivGene-genome-annotation/work/3c/48d999464a895324d779655902e4b4/MinP1831_erik_em.emapper.annotations \ No newline at end of file diff --git a/Min_test_data/nextflow_results/fixed.gff b/Min_test_data/nextflow_results/fixed.gff new file mode 120000 index 0000000000000000000000000000000000000000..7530602fd877bb9011ff574c5d2085c893e899ce --- /dev/null +++ b/Min_test_data/nextflow_results/fixed.gff @@ -0,0 +1 @@ +/home/nibio/Documents/Git/DivGene-genome-annotation/work/e3/b57eac02db62ca9e83c859e498918e/fixed.gff \ No newline at end of file diff --git a/Min_test_data/nextflow_results/interpro/MinP1831_erik_interpro.xml b/Min_test_data/nextflow_results/interpro/MinP1831_erik_interpro.xml new file mode 120000 index 0000000000000000000000000000000000000000..8def5190be5819292300f2ed3f2769c84a8aca60 --- /dev/null +++ b/Min_test_data/nextflow_results/interpro/MinP1831_erik_interpro.xml @@ -0,0 +1 @@ +/home/nibio/Documents/Git/DivGene-genome-annotation/work/d9/759832a0f6ff7f32cf96c6610b9246/MinP1831_erik_interpro.xml \ No newline at end of file diff --git a/effector_annot.r b/effector_annot.r new file mode 100644 index 0000000000000000000000000000000000000000..d57afb19612dba58393ab17aaf69aa960424a312 --- /dev/null +++ b/effector_annot.r @@ -0,0 +1,8 @@ +library(effectR) + +# fasta_prot <- seqinr::read.fasta(commandArgs(trailingOnly = TRUE)[1]) +test <- seqinr::read.fasta("Mid_test_data/nextflow_results/MidP1831_erik_protein.fa") +fasta_prot <- test + +rxlr_cand <- regex.search(fasta_prot) +crn_cand <- regex.search(fasta_prot, motif = "CRN") \ No newline at end of file diff --git a/nextflow.config b/nextflow.config index 50bbaf22bf27e4e11b3570395e6b46c82a59b0c9..2dec088b1bc9eafc7467bb79fcfc881f2a129fe3 100644 --- a/nextflow.config +++ b/nextflow.config @@ -1,5 +1,5 @@ -params.genome = "MidP1831_erik" -params.project = "/home/nibio/Documents/Git/DivGene-genome-annotation/Mid_test_data" +params.genome = "MinP1831_erik" +params.project = "/home/nibio/Documents/Git/DivGene-genome-annotation/Min_test_data" profiles { conda { diff --git a/pinf_genome_annotation.nf b/pinf_genome_annotation.nf index 1ae7ee7b19d7e85fd877756f832ffc55e5b3b1fd..4e25b99adf0ff8b2d7ba48e6714d656388dcbedd 100644 --- a/pinf_genome_annotation.nf +++ b/pinf_genome_annotation.nf @@ -80,7 +80,7 @@ publishDir "${params.project}/nextflow_results/interpro", mode: 'symlink' ~/interproscan-5.68-100.0/interproscan.sh \ -i $proteinFasta \ -f xml \ - -b "${params.genome}_interpro.xml" + -b "${params.genome}_interpro" """ } @@ -113,9 +113,10 @@ funannotate annotate \ workflow { def gff_ch = Channel.fromPath(params.genomeGFF) + def sp_ch = Channel.value(params.speciesStrain) fixGFF(gff_ch) extractProteins(params.genomeFasta, fixGFF.out) searchAndAnnotateEggnogg(extractProteins.out) searchAndAnnotateInterPro(extractProteins.out) - funannotateIntegration(fixGFF.out, params.genomeFasta, searchAndAnnotateEggnogg.out, searchAndAnnotateInterPro.out, params.sbt, params.speciesStrain) + funannotateIntegration(fixGFF.out, params.genomeFasta, searchAndAnnotateEggnogg.out, searchAndAnnotateInterPro.out, params.sbt, "Pinf") } \ No newline at end of file