Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error reading BUSCO files #9

Closed
ohudson1 opened this issue Dec 23, 2022 · 2 comments
Closed

Error reading BUSCO files #9

ohudson1 opened this issue Dec 23, 2022 · 2 comments

Comments

@ohudson1
Copy link

My problem when running the BUSCO_phylogenomics.py script is that it seems to look in all subdirectories within the sample directory (run_*) to have a folder called busco_sequences.
My input is:
python BUSCO_phylogenomics.py -i busco_runs/ -o BUSCO_tree_results --supermatrix --threads 8
Where inside busco_runs/ are 50 directories called run_samplename, for example run_FV-NE-02, and inside that are normal busco outputs: hmmer_output, busco_sequences, metaeuk_output, full_table.tsv

And the log reports that each run is 3 different runs because of their subdirectories:

run_FV-NE-02 busco_runs/run_FV-NE-02/hmmer_output
run_FV-NE-02 busco_runs/run_FV-NE-02/busco_sequences
run_FV-NE-02 busco_runs/run_FV-NE-02/metaeuk_output
No such file or directory: 'busco_runs/run_FV-NE-02/hmmer_output/busco_sequences/single_copy_busco_sequences'

I think that the script should be looking only in "busco_runs/run_FV-NE-02/busco_sequences" for single_copy_busco_sequences, and not the other directories. So instead of recognizing my 50 genome directories, it says there are 150.
How can I get it to recognize the busco_sequences folder as the appropriate one?

@ohudson1
Copy link
Author

Similar to the most recent issue raised, I was able to partially correct this error by making an extra directory layer.
I moved the "run_FV-NE-02" within a new directory that I labeled FV-NE-02_busco, within a parent directory test_busco/.

This worked well to initiate the .py script but failed to create the alignment files, so the error output was:
ERROR: Alignment not loaded: "/BUSCO_tree_results/supermatrix/alignments/50715at5125.aln" Check the file's content.

Please respond or comment if you have a corrected a similar problem.
Thanks!

@ohudson1
Copy link
Author

ohudson1 commented Jan 5, 2023

Resolved - needed an updated MUSCLE version - previous .py script versions include an option for old MUSCLE versions too.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant