In [1]:
from ga4gh.client import client
c = client.HttpClient("http://1kgenomes.ga4gh.org")
In [3]:
len(list(c.searchReads(referenceId="WyJOQ0JJMzciLCIxIl0", readGroupIds=[ 'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSMDIzMzAwIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSMDMyODU0Il0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSMDMyODU1Il0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjIyIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjIzIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjMyIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjQzIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjUxIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjU0Il0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjYwIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjYzIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjcwIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjgxIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4Njg0Il0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4Njg2Il0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4Njg4Il0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NjkwIl0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4Njk3Il0',
'WyIxa2dlbm9tZXMiLCJyZ3MiLCJOQTE5MTAyIiwiU1JSNzg4NzEwIl0' ], start=10100, end=10111)))
Out[3]:
In [4]:
import pysam
In [5]:
url = "http://s3.amazonaws.com/1000genomes/phase3/data/NA19102/alignment/NA19102.mapped.ILLUMINA.bwa.YRI.low_coverage.20130415.bam"
In [7]:
amazonalignments = pysam.AlignmentFile(url)
In [9]:
len(list(amazonalignments.fetch("1", 10100, 10111)))
Out[9]:
In [ ]: