Re: [Bioc-devel] Random access to sequences in fasta files

2015-02-20 Thread Hervé Pagès
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,

Re: [Bioc-devel] Random access to sequences in fasta files

2015-01-29 Thread Thomas Dybdal Pedersen
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

Re: [Bioc-devel] Random access to sequences in fasta files

2015-01-29 Thread Martin Morgan
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

[Bioc-devel] Random access to sequences in fasta files

2015-01-29 Thread Thomas Lin Pedersen
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 microbial genomics and it would be a huge speed improv