This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Pindel Error Message: Coordinate Of The Close End Is Bad

Hi,

I'm running the latest version of pindel (0.2.4t) on my BAM of just chromosome 2 reads and I'm getting this error:

$ pindel -i bam_config_file.txt -f /path/to/reference/ -o output_prefix
...
...
Looking at chromosome chr2 bases 10000000 to 20000000.
getReads chr2 263199373
Insertsize in bamreads: 250
Number of reads in current window:                      7490, + 3804 - 3686
Number of reads where the close end could be mapped:    4146, + 2400 - 1738
Percentage of reads which could be mapped: + 63.09% - 47.15%`
BAM file index  0       4146
There are 4146 reads for this chromosome region.
Coordinate of the close end 19994552 is bad.

As a sanity check, the reference I'm passing into pindel is the same same reference used to generate the BAM.

I see that a similar question was posted here, but the response from the pindel author (Kai Ye) was to use the latest version, which I am already: https://trac.nbic.nl/pipermail/pindel-users/2013-April/000247.html

Thank you for any suggestions!

pindel

2 answers

The latest version of pindel is actually 0.2.4w: https://github.com/genome/pindel.

See also: https://trac.nbic.nl/pipermail/pindel-users/2013-May/000271.html

it just turned 0.2.4x. check https://github.com/genome/pindel.

Kai

Log in to answer this question.