mossmatters / HybPiper

Recovering genes from targeted sequence capture data
GNU General Public License v3.0
111 stars 45 forks source link

get_seq_lengths.py and hybpiper_stats.py issues #65

Closed perla-far closed 1 year ago

perla-far commented 3 years ago

Hello, when I run the script get_seq_lengths.py, it works but it give a txt file for which at the end of the file, there are 3 rows with the value "0". These 3 rows are not assigned to any sample (screenshot below). I am not sure why I have this result. please if you can help in this issue it would be great. As well, I have an error when I run the hybpiper_stats.py (I have tried using the output text file from the script get_seq_lengths.py including and excluding the values "0" at the end of the file). This error indicate that the .Bam file could not be found (screenshot below). I have followed the same folder' hierarchy of the "test_dataset". please let me know how can I fix these issues

Thank you so much

image image

chrisjackson-pellicle commented 2 years ago

Hi @perla-far,

Are you using a namelist.txt file that has empty lines at the end? If so, that's the likely cause of what you're seeing. If you're running HybPiper version 1.3, you can simply delete the empty lines and run get_seq_lengths.py and hybpiper_stats.py again. Alternatively, this has been fixed in HybPiper2 version 2.0.1rc build 12 - empty lines in your namelist.txt file are now ignored.

Cheers,

Chris

chrisjackson-pellicle commented 1 year ago

Hi @perla-far - I'm closing this issue now, but feel free to re-open it if you still have a problem with this.