Hi Thomas,
I've added support for direct access to an arbitrary subset of sequences
in FASTA files. It works in 3 steps:
Step 1: Call fasta.index() to get an index of what's in your files:
library(Biostrings)
filepaths <- c(path1, path2, path3, etc...)
fai <- fasta.index(filepaths,
Thanks Martin
This was thought as a feauture request/discussion of biostrings, which is why I
posted it here. Thought biostrings io capabilities was behind most other fasts
readers on bioconductor...
/Thomas
> Den 29/01/2015 kl. 15.45 skrev Martin Morgan :
>
>> On 01/29/2015 06:41 AM, Thomas
On 01/29/2015 06:41 AM, Thomas Lin Pedersen wrote:
Hi
I’m querying on whether there are any plans on supporting random access reading
of fasta files in the sense that it is possible to upfront specify the indexes
of sequences that should be read in.
I’m working on a package for comparative mi
