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

Probleme with MitoHifi #145

Open
LucileSol opened this issue Dec 16, 2024 · 0 comments
Open

Probleme with MitoHifi #145

LucileSol opened this issue Dec 16, 2024 · 0 comments

Comments

@LucileSol
Copy link
Contributor

Describe the bug (12/12)
I am running the pipeline on Drosophila histrio from Darwin tree of life :
https://portal.darwintreeoflife.org/data/root/details/Drosophila%20histrio

I got the following "Note" : [05/75069e] NOTE: Process ASSEMBLE_ORGANELLES:MITOHIFI_MITOHIFI (hifiasm-raw-default) terminated with an error exit status (1) -- Error is ignored

It seems there is no mitochondrial contig found but drosophila can have NUMTS

.command.log

2024-12-12 09:32:26 [INFO] Welcome to MitoHifi v2. Starting pipeline...
2024-12-12 09:32:26 [INFO] Length of related mitogenome is: 17535 bp
2024-12-12 09:32:26 [INFO] Number of genes on related mitogenome: 37
2024-12-12 09:32:26 [INFO] Running MitoHifi pipeline in contigs mode...
2024-12-12 09:32:26 [INFO] 1. Fixing potentially conflicting FASTA headers
2024-12-12 09:32:30 [INFO] 2. Let's run the blast of the contigs versus the close-related mitogenome
2024-12-12 09:32:30 [INFO] 2.1. Creating BLAST database:
2024-12-12 09:32:30 [INFO] makeblastdb -in OP035243.1.fasta -dbtype nucl
2024-12-12 09:32:34 [INFO] Makeblastdb done.
2024-12-12 09:32:34 [INFO] 2.2. Running blast of contigs against close-related mitogenome:
2024-12-12 09:32:34 [INFO] blastn -query fixed_header_contigs.fasta -db OP035243.1.fasta -num_threads 12 -out contigs.blastn -outfmt 6 std qlen slen
2024-12-12 09:33:32 [INFO] Blast done.
2024-12-12 09:33:32 [INFO] 3. Filtering BLAST output to select target sequences
2024-12-12 09:33:32 [INFO] Filtering thresholds applied:
2024-12-12 09:33:32 [INFO] Minimum query percentage = 50
2024-12-12 09:33:32 [INFO] Minimum query length = 80% subject length
2024-12-12 09:33:32 [INFO] Maximum query length = 5 times subject length
Attention!
'parsed_blast.txt' and 'parsed_blast_all.txt' files are empty.
The pipeline has stopped !! You need to run further scripts to check if you have mito reads pulled to a large NUMT!

It seems to be a MitoHifi problem

update (16/12): the pipeline seems to be stuck so I killed it

Environment (please complete the following information):

  • OS: nac
@mahesh-panchal mahesh-panchal transferred this issue from NBISweden/assembly-project-template Jan 13, 2025
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