how to make BreakDancer's bam2cfg.pl work?

[cliff]

Dear All,

I hope someone here can help me figure out how to make BreakDancer's bam2cfg.pl work. I just downloaded the latest version of BreakDancer from http://breakdancer.sourceforge.net/index.html

From this page http://breakdancer.sourceforge.net/moreperl.html, it looks like we need download two perl modules. So, I downloaded them and put them into a library directory and added

use lib 'myDirectory/lib/';

at the beginning of bam2cfg.pl.

But after I ran bam2cfg.pl, I still got this error

Use of uninitialized value in printf at /myDirectory/bam2cfg.pl line 223.

Can anyone share your experience of doing this?

All your help will be greatly appreciated.

-c

[KevinLam]

actually on my CentOS system i found that i needed 3 modules

this are my own notes
#Perl version requires below install using cpan as root
cpan
install Statistics:escriptive;
install Math::CDF;

install GD::Graph::histogram

if you have root access and install using CPAN you do not need to fuss with the use lib option.


but your error seems like it has another reason. can you print the line 223 here?

for my version i only see

}

[natstreet]

I downloaded breakdancer a minute ago and tried running bam2cfg.pl and get the error message

Use of uninitialized value in printf at perl/bam2cfg.pl line 222

This is line 222

      printf "%d",$flagHgram{$rg}{all};

I installed the three libraries required using CPAN.

[KevinLam]

Hmmm I am using version
breakdancer-1.0_2010_06_24

but I see the same code.
will upgrade and try tomorrow.
$rg is the readgroup. perhaps the read group wasn't defined?

[cliff]

I downloaded breakdancer a minute ago and tried running bam2cfg.pl and get the error message

Use of uninitialized value in printf at perl/bam2cfg.pl line 222

This is line 222

      printf "%d",$flagHgram{$rg}{all};

I installed the three libraries required using CPAN.

my error is the same. I added

use lib 'myDirectory/lib';

at the beginning of bam2cfg.pl so my line 223 is your line 222.

Have you changed line 61 or 60 to

open(BAM,"yourDirectory/samtools view -h $fbam |") || die "unable to open $fbam\n";

I installed three modules and changed line 61 but still got the same error...

[Xian]

The error message is probably due to that the readgroup hasn't shown up as the key for flagHgram. Please change line 222 to

if(defined $flagHgram{$rg} && defined $flagHgram{$rg}{all}){
printf "%d", $flagHgram{$rg}{all};
}

Basically just adding a condition to print out the histogram number.

If the same error still appears after adding the condition, please try not using the option -g.