Entering edit mode
Ruppert Valentino
▴
270
@ruppert-valentino-1376
Last seen 10.1 years ago
Hello,
I am trying to access next gen sequencing data remotely via
R/bioconductor but I can't seem to send queries to it like using
biomaRt. I tried Rsamtools but even with that there is no way to query
the sequence file directly.
What I am trying to do is to get sequence data for specific regions
e.g. chrom5 150100000 to 150101000 from http://www.1000genomes.org/
cases e.g. NA19240, however there doesn't seem to be any tool to this
easily.
In the Rsamtools they mention that initially they downloaded this
using samtools view bamfile
Does anyone know of a way to access next gen sequence data remotely
without having to download them locally, if so I would appreciate it
if they email me the R script for that.
Thanks
[[alternative HTML version deleted]]