chloroplast is OK but mitochondrial genome always failed #375
Unanswered
lucky-riverstar
asked this question in
Q&A
Replies: 0 comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Uh oh!
There was an error while loading. Please reload this page.
-
I am sure the data is OK. Because based on the same data , GetOrganelle assembly the chloroplast successfully, while mitochondrial genome always fail. Code uesd as follow ,respectivelly:
chloroplast :
get_organelle_from_reads.py -1 SRR14309715_1_clean.fq.gz -2 SRR14309715_2_clean.fq.gz -o SRR14309715_pt2 -R 15 -t 20 -k 21,45,65,85,105 -F embplant_pt
mitochondrial :
get_organelle_from_reads.py -1 SRR14309715_1_clean.fq.gz -2 SRR14309715_2_clean.fq.gz -o SRR14309715_mt -R 20 -t 20 -k 21,45,65,85,105 -F embplant_mt
mitochondrial log:
2025-06-03 22:59:06,964 - INFO: Pre-reading fastq ...
2025-06-03 22:59:06,965 - INFO: Estimating reads to use ... (to use all reads, set '--reduce-reads-for-coverage inf --max-reads inf')
2025-06-03 22:59:07,038 - INFO: Tasting 100000+100000 reads ...
2025-06-03 22:59:15,631 - INFO: Tasting 500000+500000 reads ...
2025-06-03 22:59:31,203 - INFO: Tasting 2500000+2500000 reads ...
2025-06-03 23:00:20,804 - INFO: Estimating reads to use finished.
2025-06-03 23:00:20,804 - INFO: Unzipping reads file: SRR14309715_1_clean.fq.gz (16978450779 bytes)
2025-06-03 23:07:51,497 - ERROR:
Traceback (most recent call last):
File "/ifs1/User/01_biosoft/08_population_soft/GetOrganelle/get_organelle_from_reads.py", line 4310, in main
unzip(read_file, target_fq, 4 * all_read_nums[file_id],
File "/ifs1/User/01_biosoft/08_population_soft/GetOrganelle/GetOrganelleLib/pipe_control_func.py", line 840, in unzip
raise NotImplementedError("unzipping failed!")
NotImplementedError: unzipping failed!
So I do not understand what is happened and how to deal with it?
Besides ,the result of chloroplast is circular, but is it right about the place where red arrows pointed in the picture? If not ,how should I deal with it?

Thank you for your reply.
Beta Was this translation helpful? Give feedback.
All reactions