PBJelly gap file
I am trying to run PBJelly from PBsuite using Pacbio long reads, I wonder if anyone knows the format of the gapInfo.bed file it uses, as a standard bed file of chromosome and start and end position doesn't appear to be accepted.
The error is
Traceback (most recent call last):
File "Jelly.py", line 270, in <module>
prog.run()
File "Jelly.py", line 256, in run
self.options.extras)
File "Stages.py", line 149, in assembly
gapInfo = GapInfoFile(gapInfoFn)
File "FileHandlers.py", line 228, in __init__
curGap = Gap(*line.strip().split('\t'))
File "FileHandlers.py", line 267, in __init__
ref, lcontig, rcontig = name.split('_')
ValueError: need more than 1 value to unpack
and head of my bed file is (tab delimited)
chr1 148 326
chr1 1448 1474
chr1 2752 2850
I tried looking for an answer to this on source forge and the ticket system that PBsuite has but I couldn't find an answer.
And I made this gap file in python by searching for strings of Ns in the fasta file.
• 435 views
•
link
0 answers
No answers yet.
Log in to answer this question.