[glede]

I think you can refer to the manual at https://trac.nbic.nl/pindel/wiki/UserManual

[mig174]

I ran into some problems while trying to use pindel:
1. is there a 32bit-version?
2. while trying to compile version 0.2.4 there was an error with TIME_MINE_E, which is only set once to time(null) and never used in the rest of the code; as i did not see a possibility to fix that without writing around in your code, i decided to download version 0.2.0
3. after downloading version 0.2.0 and trying to run it with the parameters described on https://trac.nbic.nl/pindel/wiki/UserManual i only got the following usage information:

Welcome to Pindel, developed by Kai Ye, [email protected]

7 parameters are required here:
1. Input: the reference genome sequences in fasta format;
2. Input: the unmapped reads in a modified fastq format;
Better to use bam2pindel.pl to convert BAM files to Pindel input.
If the perl script fails for some reasons, please use the provided
sam2pindel.cpp to extract reads from sam files.
Compile cpp file first: g++ sam2pindel.cpp -o sam2pindel -O3
3. Output folder
4. BreakDancer result:
ChrA LocA stringA ChrB LocB stringB others
If you don't have BreakDancer result, please provide an empty file here.
5. Maximum event size index. 5 is recommended
2: 128
3: 512
4: 2,048
5: 8,092
6: 32,368
7: 129,472
8: 517,888
9: 2,071,552
10: 8,286,208
11: 33,144,832
12: 132,579,328
6. Number of threads
7. Which chr/fragment
Pindel will process reads for one chr each time
ChrName must be the same as in reference sequence and in read file

so it seems as if the information on the homepage is no longer accurate

finally i got it to run and it seems to work perfectly (thanks for that :-) ) but the way of getting there was a little painful.

[KaiYe]

I ran into some problems while trying to use pindel:
1. is there a 32bit-version?
2. while trying to compile version 0.2.4 there was an error with TIME_MINE_E, which is only set once to time(null) and never used in the rest of the code; as i did not see a possibility to fix that without writing around in your code, i decided to download version 0.2.0
3. after downloading version 0.2.0 and trying to run it with the parameters described on https://trac.nbic.nl/pindel/wiki/UserManual i only got the following usage information:

so it seems as if the information on the homepage is no longer accurate

finally i got it to run and it seems to work perfectly (thanks for that :-) ) but the way of getting there was a little painful.

1. Can you indicate which linux distribution you are using? I would install it on my PC to test it.
2. There are many nice features and improvement in new version. It is strongly advised to use the new version. And current stable version is 0.2.4h. I would like to you help you compiling the latest version. Please communicate with me by email, [email protected]
3. Current user manual is for the versions after 0.2.3, which has different but a more user friendly interface.
4. We support both BWA and mosaik BAMs now.

Please contact me by email ([email protected]) or phone (+31 71 526 9745). We will find a solution for you to run the latest version of Pindel.

Cheers,

Kai

[fyusufi]

Hi Kai,

I was wondering if there were any special parameters required while running BWA, or are the default values ok?

Also, if I have sequence libraries with different insert sizes should I map them into separate bam files?

Thanks!

[KaiYe]

Hi Kai,

I was wondering if there were any special parameters required while running BWA, or are the default values ok?

Also, if I have sequence libraries with different insert sizes should I map them into separate bam files?

Thanks!

Default parameters of BWA is fine.

If you have multiple insert sizes, it is better to map them separately. It is possible to put them in one BAM if the insert sizes don't differ very much, say 400 vs 500. Pindel will internally double the number so that slightly underestimation is fine.

Kai

[fyusufi]

Thanks for the quick reply!

[libiyagirl]

Hey Kai,

I am just wondering about the power of Pindel to detect long insertions. When I was running Pindel on a single chromosome, I didn't find any large insertions. But actually there are some duplications in that chromosome identified by breakdancer. Did I miss something? The "zero" LI looks very strange to me.

thanks a lot,

[KaiYe]

Hey Kai,

I am just wondering about the power of Pindel to detect long insertions. When I was running Pindel on a single chromosome, I didn't find any large insertions. But actually there are some duplications in that chromosome identified by breakdancer. Did I miss something? The "zero" LI looks very strange to me.

thanks a lot,

Please check whether you switch on long insertion option. In some versions, this function is turned off by default.

Kai

[ralonso]

Hello,

I am about to use pindel and I am wondering some questions.
I have a file mapped with bwa and illumnia reads.
As I see in the pindel webpage https://trac.nbic.nl/pindel/wiki/UserManual I should use "unmappable read". I have several bams, so I don't know which to use, those are:
1. sample.bam: mapped and unmmapped reads,it is not sorted
2. sample_mapped.bam: just mapped reads, it is not sorted, it has not unmapped
3. sample_mapped_singlehit_sorted.bam: just mapped reads, it is sorted, it has only single hit and it is not unmapped

should I use one of this files? or maybe another with the unmappable reads?

thank you very much!

[KaiYe]

Hello,

I am about to use pindel and I am wondering some questions.
I have a file mapped with bwa and illumnia reads.
As I see in the pindel webpage https://trac.nbic.nl/pindel/wiki/UserManual I should use "unmappable read". I have several bams, so I don't know which to use, those are:
1. sample.bam: mapped and unmmapped reads,it is not sorted
2. sample_mapped.bam: just mapped reads, it is not sorted, it has not unmapped
3. sample_mapped_singlehit_sorted.bam: just mapped reads, it is sorted, it has only single hit and it is not unmapped

should I use one of this files? or maybe another with the unmappable reads?

thank you very much!

sorted sample.bam will be fine. Pindel reads bam files directly and get useful reads by itself.

You need to use config file to tell Pindel where are your files, insert size and sample names.

Kai

[ralonso]

I don't have it sorted, could I use it even this fact? or is it totally necessary to sort it?

[KaiYe]

I don't have it sorted, could I use it even this fact? or is it totally necessary to sort it?

It is better sorted and then let Pindel directly read your sorted bam, -i config file.

Otherwise, you may use sam2pindel to extract reads (then use -p extracted_reads), but it is not recommended, as this takes additional space and runtime.

Kai

[sam.a]

hi, KaiYe
i'm wondering - does Pindel works only with human genome?

[KaiYe]

hi, KaiYe
i'm wondering - does Pindel works only with human genome?

Of course, as long as there is a reference genome or contigs.

[mikhmv]

Hi KaiYe,

I used MarkDuplicates from Picard and mark duplicates but did not removed them.
I have a question does pindel need file with removed duplicates or can work with Marked ones?

Thank you in an advance.