Breakdancer And Bam2Cfg.Pl
1
0
Entering edit mode
12.2 years ago
madkitty ▴ 690

We have whole genome sequencing file aligned into a bam format file. Each sample has at least 2 bam files.

I want to find CNV Gain and Deletion from those bam files with Breakdancer, but I need to create a configuration file with bam2cfg.pl

They said I have to type : bam2cfg.pl -q [BDmapq] -n [readsToSample] bamfile > bd.config

How do I know BDmapq? How do I know readsToSample?

What's "bamfile" ? Is it the path to my bam file ? What's "bd.config" ??

Let's I have the following bam file /home/myfile.bam /home/myfile.bai

What do I do from here ??

cnv breakdancer • 6.2k views
ADD COMMENT
2
Entering edit mode
12.2 years ago
Mike Smith ★ 2.1k

"BDmapq" is an integer specifying the minimum mapping quality for a read to be used by the algorithm.

"readsToSample" is also an integer, this time giving the number of read pairs that should be used to estimate the insert size of your library.

However, you can omit both and Breakdancer will use its default values. You can find a few more details of the arguments here: http://breakdancer.sourceforge.net/bam2cfg.html

The "bamfile" argument is a space separated list of the bam files you want included. You say you have at least two bam files per sample, so you probably want something like: /home/myfile1.bam /home/myfile2.bam

The final part, "bd.config", is the name of the config file to be created. if you don't include it the resulting config information will be printed to the screen. Using ">" will direct that to the specified file, rather than the screen.

ADD COMMENT
0
Entering edit mode

Thanks for your answer. I tried with all the values and now it says "Can't locate GD/Graph/histogram.pm in @INC" ..

ADD REPLY
1
Entering edit mode

Hi,

Then you have to install that module. You can find it (and plenty more of them) in the CPAN Repository. The one you're lookin for is this

http://search.cpan.org/~lds/GD-2.46/GD.pm

There are some alternative (and faster) ways to install a CPAN Module but if you haven't done it before I'd suggest the manual installation. Googling you can find lots of help for doing this. Like this link.

http://www.thegeekstuff.com/2008/09/how-to-install-perl-modules-manually-and-using-cpan-command/

A common problem in not having permissions to install modules on the system so you have to install it locally in your directory and then tell the system where to find that module. In oder words, include the folder where you installed the module in your $PATH. Some guidelines here

http://www.perlmonks.org/index.pl?node_id=128077

ADD REPLY

Login before adding your answer.

Traffic: 1154 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6