Hey Guys
I am wondering if there is a slick way access all the possible alignments for a read present in sam or bam file given the read header. Since the existing codebase is in perl I would prefer something which can be done in/via perl.
By default BAM's are indexed by location so the inbuilt samtools indexing wont work I guess.
I should also say the input bam file will have in the order of 500 million total alignments and many reads are expected to be aligned to more than one place in the genome. Given the size of the data loading it all in one big hash is not turning out to be memory friendly.
Thanks!
-Abhi
Have a look at Bio::DB::Sam. This does not ship with the core BioPerl, and to compile this you will need to get the Kent source tree a recent version of SAMtools. Building this can be tricky depending on your OS, but I'm sure there are plenty of people that can offer advice.
Actually, you don't need the Kent source tree (that would be Bio::DB::BigFile), but you do need a recent version of samtools.
Actually, you don't need the Kent source tree (that would be needed for Bio::DB::BigFile), but you do need a recent version of samtools
•
written
by
Chris Fields
Log in to answer this question.
