funannotate: CMD ERROR during update using fun annotate-docker

I am trying to run update following the tutorial, but I ran into a CMD ERROR as seen below, which didn’t explicitly state the issue occurring:

$ funannotate-docker update -i ${WDIR}/genomes/04_ARF-L_fun --cpus 19
-------------------------------------------------------
[Jan 09 10:11 PM]: OS: Debian GNU/Linux 10, 20 cores, ~ 132 GB RAM. Python: 3.8.12
[Jan 09 10:11 PM]: Running 1.8.16
[Jan 09 10:11 PM]: No NCBI SBT file given, will use default, for NCBI submissions pass one here '--sbt'
[Jan 09 10:11 PM]: Found relevant files in /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training, will re-use them:
        GFF3: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/predict_results/Hyalorbilia_cf._ulicicola_ARF-L.gff3
        Genome: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/predict_results/Hyalorbilia_cf._ulicicola_ARF-L.scaffolds.fa
        Single reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/single.fq.gz
        Forward reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/left.fq.gz
        Reverse reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/right.fq.gz
        Forward Q-trimmed reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/trimmomatic/trimmed_left.fastq.gz
        Reverse Q-trimmed reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/trimmomatic/trimmed_right.fastq.gz
        Single Q-trimmed reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/trimmomatic/trimmed_single.fastq.gz
        Forward normalized reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/normalize/left.norm.fq
        Reverse normalized reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/normalize/right.norm.fq
        Single normalized reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/normalize/single.norm.fq
        Trinity results: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/funannotate_train.trinity-GG.fasta
        PASA config file: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/pasa/alignAssembly.txt
        BAM alignments: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/funannotate_train.coordSorted.bam
        StringTie GTF: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/funannotate_train.stringtie.gtf
[Jan 09 10:11 PM]: Reannotating Hyalorbilia cf. ulicicola, NCBI accession: None
[Jan 09 10:11 PM]: Previous annotation consists of: 16,589 protein coding gene models and 85 non-coding gene models
[Jan 09 10:11 PM]: Existing annotation: locustag=FUN_ genenumber=16674
[Jan 09 10:11 PM]: Converting transcript alignments to GFF3 format
[Jan 09 10:11 PM]: Converting Trinity transcript alignments to GFF3 format
[Jan 09 10:11 PM]: PASA database is SQLite: Hyalorbilia_cf__ulicicola_ARF_L_pasa
[Jan 09 10:11 PM]: Running PASA annotation comparison step 1
[Jan 09 10:11 PM]: CMD ERROR: /venv/opt/pasa-2.4.1/Launch_PASA_pipeline.pl -c /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/update_misc/pasa/annotCompare.txt -g /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/update_misc/genome.fa -t /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/update_misc/trinity.fasta.clean -A -L --CPU 19 --annots /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/update_misc/genome.gff3
--------------

Previously I’d run into issues and ended up using the mysql database when I ran the train command, and it completed successfully. So I added the --pasa_db option on here as well, hoping it would handle the error, but I got a different one:

$ funannotate-docker update -i ${WDIR}/genomes/04_ARF-L_fun \
    --pasa_db mysql --cpus 19
-------------------------------------------------------
[Jan 09 10:17 PM]: OS: Debian GNU/Linux 10, 20 cores, ~ 132 GB RAM. Python: 3.8.12
[Jan 09 10:17 PM]: Running 1.8.16
[Jan 09 10:17 PM]: No NCBI SBT file given, will use default, for NCBI submissions pass one here '--sbt'
[Jan 09 10:17 PM]: Found relevant files in /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training, will re-use them:
        GFF3: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/predict_results/Hyalorbilia_cf._ulicicola_ARF-L.gff3
        Genome: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/predict_results/Hyalorbilia_cf._ulicicola_ARF-L.scaffolds.fa
        Single reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/single.fq.gz
        Forward reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/left.fq.gz
        Reverse reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/right.fq.gz
        Forward Q-trimmed reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/trimmomatic/trimmed_left.fastq.gz
        Reverse Q-trimmed reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/trimmomatic/trimmed_right.fastq.gz
        Single Q-trimmed reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/trimmomatic/trimmed_single.fastq.gz
        Forward normalized reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/normalize/left.norm.fq
        Reverse normalized reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/normalize/right.norm.fq
        Single normalized reads: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/normalize/single.norm.fq
        Trinity results: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/funannotate_train.trinity-GG.fasta
        PASA config file: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/pasa/alignAssembly.txt
        BAM alignments: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/funannotate_train.coordSorted.bam
        StringTie GTF: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/training/funannotate_train.stringtie.gtf
[Jan 09 10:17 PM]: Reannotating Hyalorbilia cf. ulicicola, NCBI accession: None
[Jan 09 10:17 PM]: Previous annotation consists of: 16,589 protein coding gene models and 85 non-coding gene models
[Jan 09 10:17 PM]: Existing annotation: locustag=FUN_ genenumber=16674
[Jan 09 10:17 PM]: Existing BAM alignments found: /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/update_misc/trinity.alignments.bam, /mnt/data1/PN97_fungal_genomes/genomes/04_ARF-L_fun/update_misc/transcript.alignments.bam
Traceback (most recent call last):
  File "/venv/bin/funannotate", line 8, in <module>
    sys.exit(main())
  File "/venv/lib/python3.8/site-packages/funannotate/funannotate.py", line 717, in main
    mod.main(arguments)
  File "/venv/lib/python3.8/site-packages/funannotate/update.py", line 3278, in main
    runPASA(
  File "/venv/lib/python3.8/site-packages/funannotate/update.py", line 846, in runPASA
    if not getPASAinformation(configFile, DataBaseName, folder, genome):
  File "/venv/lib/python3.8/site-packages/funannotate/update.py", line 741, in getPASAinformation
    with open(pasaconf_file, "r") as pasaconf:
FileNotFoundError: [Errno 2] No such file or directory: '/venv/opt/pasa-2.4.1/pasa_conf/conf.txt'

About this issue

  • Original URL
  • State: closed
  • Created 6 months ago
  • Comments: 15 (7 by maintainers)

Most upvoted comments

I usually just use HPC if I need more resources than my current computer has. I selected mysql primarily because the person working on this project before me did and I was just referencing their commands while trying to re-run the process. Assuming sqlite takes much more time I will probably go a different route! I’ll give it a try and see what happens. Thanks again!