gatb-core icon indicating copy to clipboard operation
gatb-core copied to clipboard

Reading BAM files

Open mikael-s opened this issue 2 years ago • 4 comments

Hi!

Is there any possibility to read a BAM file as a bank of sequences? Do you know of any software relying on GATB which would have done that?

Thanks!

mikael-s avatar Jul 19 '21 07:07 mikael-s

Hi,

I did that some time ago using htslib, but did not push it to GATB master branch because of the added dependency to htslib that I thought would bother some people. Let me see if I can put that in a branch or fork.

rizkg avatar Jul 19 '21 16:07 rizkg

another way would be a wrapper script with named pipes? e.g. https://github.com/COMBINE-lab/salmon/issues/29#issuecomment-168534703

rchikhi avatar Jul 20 '21 13:07 rchikhi

I did that some time ago using htslib, but did not push it to GATB master branch because of the added dependency to htslib that I thought would bother some people. Let me see if I can put that in a branch or fork.

Awesome, thanks!

another way would be a wrapper script with named pipes? e.g. COMBINE-lab/salmon#29 (comment)

Named pipes are nice but unfortunately I may need to read the input file more than once.

mikael-s avatar Jul 21 '21 08:07 mikael-s

yes exactly, I came back to post this, sorry this pipe solution can't work in gatb as it already reads the file more than once internally

rchikhi avatar Jul 21 '21 09:07 rchikhi