When running tophat_fusion_post with multiple threads, I found several of these errors in my log file:
Traceback (most recent call last):
File "/opt/wtsi-cgp/bin/tophat-fusion-post", line 2924, in <module>
sys.exit(main())
File "/opt/wtsi-cgp/bin/tophat-fusion-post", line 2901, in main
read_dist(params)
File "/opt/wtsi-cgp/bin/tophat-fusion-post", line 1491, in read_dist
parallel_work(pids, work)
File "/opt/wtsi-cgp/bin/tophat-fusion-post", line 1024, in parallel_work
work()
File "/opt/wtsi-cgp/bin/tophat-fusion-post", line 1486, in work
alignments_region(sample_name, bam_file_name, alignments_list)
File "/opt/wtsi-cgp/bin/tophat-fusion-post", line 1269, in alignments_region
chr1, chr2 = chr1.split('-')
ValueError: too many values to unpack
Strange enough, these errors do not stop the process, and I still have the output.
I suspect this has something to do with contig names in GRCh38, as we have many HLA-* contigs in the reference, and I do not see the error when using GRCh37.
I'm wondering if you can fix it, as I am not sure if these errors could also silently eliminate other fusion events from the output.
Hi,
When running tophat_fusion_post with multiple threads, I found several of these errors in my log file:
Strange enough, these errors do not stop the process, and I still have the output.
I suspect this has something to do with contig names in GRCh38, as we have many
HLA-*
contigs in the reference, and I do not see the error when using GRCh37.I'm wondering if you can fix it, as I am not sure if these errors could also silently eliminate other fusion events from the output.
The full log is below.
Thanks.