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]]