support for more genomes [igvR]
1
2
Entering edit mode
Peng ▴ 20
@peng-24004
Last seen 4.0 years ago

I'm a new user of igvR and finds it a perfect tools to visualize my interested regions in R! Just a follow-up on the plan to support custom genomes, it would be great if the following method can be added


setCustomGenome(igv, name="TAIR10",
                fastaURL="https://myHost.net/tair10/Arabidopsis_thaliana.TAIR10.dna.toplevel.fa",
                indexURL="https://myHost.net/tair10/Arabidopsis_thaliana.TAIR10.dna.toplevel.fa.fai",
                aliasURL: "https://myHost.net/tair10/chromosomeAliases.txt",
                annotationURL: "https://myHost.net/TAIR10_genes.sorted.chrLowered.gff3.gz")

I'd be definitely interested in giving it a try if possible!

Thanks!

igvR • 1.2k views
ADD COMMENT
1
Entering edit mode

Working on this now. thanks for the nudge, Peng.

ADD REPLY
4
Entering edit mode
Paul Shannon ▴ 470
@paul-shannon-5944
Last seen 2.5 years ago
United States

igvR in devel, version 1.11.1, has the new method, setCustomGenome. A vignette is coming soon. Man page for the function is in the package, and two examples are in inst/unitTests/test_setCustomGenome, of which this is the all-parameters-specified version, using hg38 files offered on amazon by the Broad:

 setCustomGenome(igv,
                 id="hg38",
                 genomeName="Human (GRCh38/hg38)",
                 fastaURL="https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg38/hg38.fa",
                 fastaIndexURL="https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg38/hg38.fa.fai",
                 cytobandURL="https://s3.amazonaws.com/igv.broadinstitute.org/annotations/hg38/cytoBandIdeo.txt",
                 chromosomeAliasURL=NA,
                 geneAnnotationName="Refseq Genes",
                 geneAnnotationURL="https://s3.amazonaws.com/igv.org.genomes/hg38/refGene.txt.gz",
                 geneAnnotationTrackHeight=500,
                 geneAnnotationTrackColor="red",
                 initialLocus="chr5:88,621,308-89,001,037",
                 visibilityWindow=5000000)
ADD COMMENT
2
Entering edit mode

igvR version 1.11.2 in devel has this feature, as demonstrated in a new vignette: 1

The vignette also explains how to start a simple local webserver, with byte ranges and CORS headers.

ADD REPLY
0
Entering edit mode

That is great! Thanks a bunch Paul!

ADD REPLY

Login before adding your answer.

Traffic: 919 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6