Question: Problems With Sam2Pindel
0
gravatar for jason.r.gallant
4.6 years ago by
jason.r.gallant10 wrote:

Hi all,

I'm trying to operate the PINDEL pipeline on some resequencing data that was aligned with Bowtie2. I am starting with the sam2pindel program, and things seem to proceed nicely for several hundred megabytes of output... however the program mysteriously hangs suddenly. Has anyone observed this problem before, or have any advice how better to troubleshoot the issue?

Many thanks, Jason Gallant

pindel • 1.2k views
ADD COMMENTlink modified 16 months ago by Biostar ♦♦ 20 • written 4.6 years ago by jason.r.gallant10

Does it tell anything? Any output or any log file.

ADD REPLYlink written 4.6 years ago by Ashutosh Pandey11k

Nope, no output or log, just hangs.

ADD REPLYlink written 4.6 years ago by jason.r.gallant10

Have you tried bam2pindel.pl?

ADD REPLYlink written 4.6 years ago by PoGibas4.7k
0
gravatar for liangkaiye
4.6 years ago by
liangkaiye250
United States
liangkaiye250 wrote:

hi Jason,

can you tell me how you run sam2pindel? did you pipe samtools view (or other ways) or use sam file as input? can you cut the bam where it went wrong and send to me k.ye@genome.wustl.edu)? sam2pindel is rather simple way of selecting reads for pindel run. have you try -i and feed bowtie BAM files directly to Pindel?

I do not have bowtie bam here but would spend sometime to study bowtie sam output to optimize read selection, if this is requested by Pindel users.

Kai

ADD COMMENTlink written 4.6 years ago by liangkaiye250

Hi Kai,

Thanks for your note. I am indeed piping samtoolsview into sam2pindel. I wasn't brave enough to just try feeding the bowtie2 bam files directly into pindel :)

I will try this now, and see how it goes. I will peek into the bamfile and see if I can send it to you directly.

Cheers, Jason

ADD REPLYlink written 4.6 years ago by jason.r.gallant10

great, it might be ideal to have the same data aligned with bwa and bowtie, then pass the bam files in one run to pindel. then you will see the difference clearly.

use -i option as input in the config file (assume 500 insert size): /path/to/bwa.bam 500 bwa /path/to/bowtie.bam 500 bowtie

Kai

ADD REPLYlink written 4.6 years ago by liangkaiye250
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1273 users visited in the last hour