Hi,
I try to construct a BioMap structure from a BAM file. I use
info = baminfo(bamFilename,'ScanDictionary',true);bm = BioMap(bamFilename, 'SelectReference', info.ScannedDictionary{chromosomeNo});
Eventhough baminfo shows a valid reference (info.ScannedDictionary) and a number of mapped reads (info.ScannedDictionaryCount) when I create the bm object, this has no reads inside. Only 1 chromosome out of 16 gives me this problem, so I don't understand what's happening…
Any help is appreciated. Thanks!
Razvan
Best Answer