View Single Post
Old 02-04-2011, 08:54 AM   #1
libiyagirl
Member
 
Location: canada

Join Date: Jan 2011
Posts: 27
Default Using Dindel: callback and output on realigned BAM

Hello, when I'm using dindel, I want to output the realigned BAM file for the genome viewer. So I used the option --outputRealignedBAM, and then I had many many files. Is it possible that I merge all the realigned BAM files into one? How can I do that?
Also, I used the callback script provided by the author:
#!/bin/bash
samtools view $1 >> $2.sam
echo $3 $4 $5 >> $2.windows
rm $1
But it's not working. Any idea to improve this?

Many thanks!
libiyagirl is offline   Reply With Quote