This is a test version of Biostars. For the public version, visit https://www.biostars.org.
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.

pacbio pbjelly

0 answers

No answers yet.

Log in to answer this question.