How to open a BAM file?

16 次查看（过去 30 天）
Razvan
Razvan 2012-11-10
Hi, I have a problem opening a BAM file (the one from here: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM754388 the link is at the bottom of the page)
Usually I open BAM files like this:
bamFilename = 'GSM754388_P2_5_20110519_7.bam';
info = baminfo(bamFilename,'ScanDictionary',true);
bm = BioMap(bamFilename, 'SelectReference', info.ScannedDictionary{1});
but this time it doesn't work, and some error pops up...
I tried also to open it with bamread but with no success:
bamFilename = 'GSM754388_P2_5_20110519_7.bam';
info = baminfo(bamFilename);
RefSeq = info.SequenceDictionary;
data = bamread(bamFilename, RefSeq(1).SequenceName, [1, RefSeq(1).SequenceLength]);
Is there any other way to access this BAM file in MATLAB?
Thanks,
Razvan

请先登录，再回答此问题。

采纳的回答

Lucio Cetto
Lucio Cetto 2012-11-14
It may be possible that the BAM file is not ordered. The Bioinfomatics Toolbox (12b) automatically indexes the BAM file, but for this it must be ordered. You can order the BAM file using SAMTOOLS or by using the following undocummented option in MATLAB:
bioinfoprivate.bamaccessmex('bamsort', 'unorderedname.bam', 'orderedname')

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Bioinformatics Toolbox 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by