cobra
cobra copied to clipboard
Missing `COBRA_retrieved_for_joining` contig file
Hello Cobra authors!
Thanks for building this tool.
I keep getting this error of a contig file not being created.
Traceback (most recent call last):
File "/home/mall0133/miniconda3/envs/cobra/bin/cobra-meta", line 10, in <module>
sys.exit(main())
File "/home/mall0133/miniconda3/envs/cobra/lib/python3.8/site-packages/cobra.py", line 1747, in main
'\t'.join([contig, str(header2len[contig]), summarize(contig), query2current[contig]]) + '\n')
File "/home/mall0133/miniconda3/envs/cobra/lib/python3.8/site-packages/cobra.py", line 575, in summarize
b = count_seq('COBRA_retrieved_for_joining/{0}_retrieved.fa'.format(item)) # number of retrieved contigs
File "/home/mall0133/miniconda3/envs/cobra/lib/python3.8/site-packages/cobra.py", line 482, in count_seq
a = open(fasta_file, 'r')
FileNotFoundError: [Errno 2] No such file or directory: 'COBRA_retrieved_for_joining/NODE_1990_length_17225_cov_13.396421_retrieved.fa'
Every time I try re-running, it fails at this stage with a FileNotFoundError
for a different contig.
This is my Cobra command. I'm using the latest version (version 1.2.3).
cobra-meta -f contigs.fasta -q query_contigs.fasta -c coverage.tsv -m sorted_reads.bam -a metaspades -mink 21 -maxk 127
My query_contigs.fasta
file contains about 2300 contig sequences. I've also attached the log file of the run.
log.txt
Any advice on how to fix this error and get Cobra running will be appreciated.
Thanks!