This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problem for pindel2vcf to create VCF file from pindel created _D files

Hi

I am new pindel user, and use ../pindel -f hg19.fasta -i p1_config.txt -c ALL -o P2_Index_95_P319 to create P2_Index_95_P319 with _SI, _D, and so on. When I use pindel2vcf to create vcf file with _SI. I got it. and also found some interesting insertion. But when I create vcf filw with_D, using ../pindel2vcf -p P2_Index_95_P319_D -r hg19.fasta -R hg19-GRCh37 -d 200902 -v P2_Index_95_P319_D.vcf, some error is as following:

Exiting reference scanning.
Processing chromosome chr1
Reading region 0-300000000
Reading chromosome chr1 into memory.
terminate called after throwing an instance of 'std::out_of_range'
what(): basic_string::erase
Aborted (core dumped)

It will be appreciated to give suggestion from anybody. I try other samples, All _SI get its vcf, but _D gives same errors.

software-error pindel pindel2vcf
• 7,428 views • link
updated by Ram Ram • written by hongbing.li • hongbing.li

chr1 is ~250E6 bp but your tool is trying to read 300E6 bp ? did you write a config file with this '300000000' ?

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Can you communicate with me by email on this? kye@genome.wustl.edu

I need the pindel output lines causing the error.

• 0 views • link
updated by Ram Ram • written by liangkaiye liangkaiye

Hi there, I am having a similar problem with pindel2vcf. Here is my command:

/path/to/pindel2vcf -p /path/to/pindel_chrXsubset_D -r /path/to/canFam3.fa -R canFam31 -d 20141204 -v /path/to/pindel_chrXsubset_D.vcf -G

And the output I'm getting is:

terminate called after throwing an instance of 'std::out_of_range'
  what():  basic_string::erase

Thanks for any thoughts you might have.

Steve

• 0 views • link
updated by Ram Ram • written by steven_friedenberg • steven_friedenberg

Problem has been resolved by Lingkaiye.

Thanks for all of you!

• 0 views • link
written by hongbing.li • hongbing.li

Hi.

I experience the same sort of problem. Could you post about how to fix it please?

/winnie

My Error message:

Exiting reference scanning.
Processing chromosome P0116_N055-10_trimmed_(paired)_contig_1_consensus
Reading region 0-300000000
Reading chromosome P0116_N055-10_trimmed_(paired)_contig_1_consensus into memory.
terminate called after throwing an instance of 'std::out_of_range'
  what():  basic_string::erase: __pos (which is 60) > this->size() (which is 3)
Abort trap: 6
• 0 views • link
updated by Ram Ram • written by winnie • winnie

Hi winnie,

Following is resolved by Liangkeye for me

grep ChrID P2_Index_13_P375_D > P2_Index_13_P375_D.head.txt
../pindel2vcf -p P2_Index_13_P375_D.head.txt -r hg19.fasta -R hg19-GRCh37 -d 200902 -v P2_Index_13_P375_D.head.txt.vcf
• 0 views • link
updated by Ram Ram • written by hongbing.li • hongbing.li

I will give a go. Thank you.

best regards,

Winnie

• 0 views • link
updated by Ram Ram • written by winnie • winnie

1 answer

Thanks for all your response, and I have sent files to kye@genome.wustl.edu

• 0 views • link
updated by Ram Ram • written by hongbing.li • hongbing.li

Log in to answer this question.

More posts like this