Pindel: improved version for indels and structural variants

[KaiYe]

Hi KaiYe

I'm having problem running Pindel. Here's what I've done:
1) Download all files from http://www.ebi.ac.uk/~kye/pindel/v_0.2.0/
2) ran bam2pindel.pl on one paired-end samples (aligned using BWA). My bam file is sorted but it does not have the header expected by your program, so i used the -om to force the script to run.
a number of files is generated: e.g. myprefix.1.txt (chr1)
3) then I tried running pindel_x86_64, but i then got this error message: ./pindel_x86_64: /usr/lib64/libstdc++.so.6: version `GLIBCXX_3.4.9' not found (required by ./pindel_x86_64)
4) i tried upgrading some packages in my redhat linux, but still the same.
5) i then downloaded your source code from sourceforge (with svn) and compiled your pindel from scratch. It seems to work.
6) I find the "-i" parameter confusing as it says "-i, --config-file: the bam file later to be a config file;" in the script but "Input: the unmapped reads in a modified fastq format" in your powerpoint manual.
7) I assumed -i refers to the files generated by bam2pindel.pl, so i tested the command on some chromosomes. E.g.
pindel_64 -f hg19.fasta -i myprefix.1.txt -o otherprefix -c 1 -b empty.txt
8) but whichever chromosome i try, i always get "There are no reads for this chromosome":

Processing chromosome 1
Processing chromosome 2
Skip the rest of chromosomes.
1 249250621 269250621
26926 10000
BreakDancer events: 0
There are no reads for this chromosome.


What have i done wrong?

my email is jason.li @ petermac.org

Thanks
Jason

I will send you my source code via email.

[KaiYe]

Hi Kaiye,

interesting tool you got there. anyway, have you publish the method? I am curious about one thing, say you have 1 read, you will grow the pattern until you cannot get a match, then you find the rest of the read within the next 1-1M bps. What if there are several matches in the 1-1M bps region, which one do you use and what kind of consideration do you use to choose it?

Yes, Pindel has been published (http://www.ncbi.nlm.nih.gov/pubmed/19561018) and it was awarded best paper at ISMB 2009 Special Interest Group on Short Read Sequencing.


Only unique hit will be considered here.

[rwenang]

Hi Kaiye,

interesting tool you got there. anyway, have you publish the method? I am curious about one thing, say you have 1 read, you will grow the pattern until you cannot get a match, then you find the rest of the read within the next 1-1M bps. What if there are several matches in the 1-1M bps region, which one do you use and what kind of consideration do you use to choose it?

[jtjli]

help

Hi KaiYe

I'm having problem running Pindel. Here's what I've done:
1) Download all files from http://www.ebi.ac.uk/~kye/pindel/v_0.2.0/
2) ran bam2pindel.pl on one paired-end samples (aligned using BWA). My bam file is sorted but it does not have the header expected by your program, so i used the -om to force the script to run.
a number of files is generated: e.g. myprefix.1.txt (chr1)
3) then I tried running pindel_x86_64, but i then got this error message: ./pindel_x86_64: /usr/lib64/libstdc++.so.6: version `GLIBCXX_3.4.9' not found (required by ./pindel_x86_64)
4) i tried upgrading some packages in my redhat linux, but still the same.
5) i then downloaded your source code from sourceforge (with svn) and compiled your pindel from scratch. It seems to work.
6) I find the "-i" parameter confusing as it says "-i, --config-file: the bam file later to be a config file;" in the script but "Input: the unmapped reads in a modified fastq format" in your powerpoint manual.
7) I assumed -i refers to the files generated by bam2pindel.pl, so i tested the command on some chromosomes. E.g.
pindel_64 -f hg19.fasta -i myprefix.1.txt -o otherprefix -c 1 -b empty.txt
8) but whichever chromosome i try, i always get "There are no reads for this chromosome":

Processing chromosome 1
Processing chromosome 2
Skip the rest of chromosomes.
1 249250621 269250621
26926 10000
BreakDancer events: 0
There are no reads for this chromosome.


What have i done wrong?

my email is jason.li @ petermac.org

Thanks
Jason

[KaiYe]

I must be missing something.. it isn't producing any output for me, but it also isn't giving an error. I'm trying to convert my BAM file to the pindel format:

[heather@frankie (Mon Jul 26 13:33:04)]% bam2pindel.pl -i aln.sort.pindel.bam -o out.pindel -s retina -om -pi 150
[heather@frankie (Mon Jul 26 13:41:02)]% ls
. aln.sort.fix.bam.bai horse_genome_v2_all.fa.ann horse_genome_v2_all.fa.rsa
.. aln.sort.pindel.bam horse_genome_v2_all.fa.bwt horse_genome_v2_all.fa.sa
align.sort.bam aln_read1.sai horse_genome_v2_all.fa.fai tag_trim_6_1.fq
aln.bam aln_read2.sai horse_genome_v2_all.fa.pac tag_trim_6_2.fq
aln.sam horse_genome_v2_all.fa horse_genome_v2_all.fa.rbwt
aln.sort.fix.bam horse_genome_v2_all.fa.amb horse_genome_v2_all.fa.rpac

Would you please inform me your email address? I have cpp code to extract reads from sam files for Pindel.

Thanks.

[KaiYe]

Hi Kai,

Can Pintel call large structural variants (>1M) now?

Thanks.

Wu

Pindel can detect variants of any sizes as long as they are not inter-chromosome events. The only thing I worry about is speed. The runtime is linear to the maximum size of SVs.

I am currently testing a new version with the following additional functions:
1. Allow sequence errors/SNPs in the same reads containing INDELs/SVs
2. non-template sequence in deletions
3. inversions
4. tandem duplications
5. breakpoints of large insertions

Please send me an email for ask for it in case you want to test it.

Cheers,

Kai

[KaiYe]

Dear Kai Ye,
I've used the pindel software recently. I have heared that you have published a new version software, will you please give me the linkage please.




Best,
Cong chen
Wenzhou Medical College

Hi Cong Chen,

It seems to me that I have sent you my latest Pindel for test. Have you experienced any problems in using it?

Kai

[tinacai]

Dear Kai Ye,
I've used the pindel software recently. I have heared that you have published a new version software, will you please give me the linkage please.




Best,
Cong chen
Wenzhou Medical College

[wuhoucdc]

Hi Kai,

Can Pintel call large structural variants (>1M) now?

Thanks.

Wu

[raela]

I must be missing something.. it isn't producing any output for me, but it also isn't giving an error. I'm trying to convert my BAM file to the pindel format:

[heather@frankie (Mon Jul 26 13:33:04)]% bam2pindel.pl -i aln.sort.pindel.bam -o out.pindel -s retina -om -pi 150
[heather@frankie (Mon Jul 26 13:41:02)]% ls
. aln.sort.fix.bam.bai horse_genome_v2_all.fa.ann horse_genome_v2_all.fa.rsa
.. aln.sort.pindel.bam horse_genome_v2_all.fa.bwt horse_genome_v2_all.fa.sa
align.sort.bam aln_read1.sai horse_genome_v2_all.fa.fai tag_trim_6_1.fq
aln.bam aln_read2.sai horse_genome_v2_all.fa.pac tag_trim_6_2.fq
aln.sam horse_genome_v2_all.fa horse_genome_v2_all.fa.rbwt
aln.sort.fix.bam horse_genome_v2_all.fa.amb horse_genome_v2_all.fa.rpac

[sdvie]

Hi Sophia,

If you generate bam from sam directly after mapping with BWA, you don't have to do anything else.

Kai

cool, thanks.

[KaiYe]

Hi Sophia,

If you generate bam from sam directly after mapping with BWA, you don't have to do anything else.

Kai

[sdvie]

Hi Kai,

Just trying out pindel for the first time...one question on the bam2pindel step: in the user manual, the input for this script is described as "aln.NameSorted.MateFixed.bam". Does this mean I have to do something additional to the bam generated by samtools? If yes, what?

Thanks,
Sophia

[KaiYe]

hi all

Just put an improved Pindel on my website https://trac.nbic.nl/pindel/ with wiki, mail list, user manual.

An instruction using it from BWA mapping is provided. You can use it to detect indels and SVs at single-base resolution from SLX paired-end short reads.

Currently 1bp-1M bp deletions and 1bp-(read length -20)bp insertions can be detected. You can also find events of non-template insertion in deletions.

I am working on inversions and large insertions as well as using pindel for RNA-Seq data.
Please comment on Pindel and suggest additional functions.

Kai
[email protected]