[ensembl-dev] Display BAM with respect to Ensembl features
th3 at sanger.ac.uk
Thu Aug 22 10:21:11 BST 2013
The coordinates of your reads in the bam file should be genomic so there is no easy way. You will need to a write a Perl script to modify your bam file.
As you used Ensembl gene models the transformation of your coordinates should be pretty easy as the reference name in your bam file should be the stable id of your feature. Then you can get the genomic coordinates of the feature with the Perl API and calculate the genomic coordinates of your read.
Hope this help
On 22 Aug 2013, at 01:30, Ben Warren <Ben.Warren at plantandfood.co.nz> wrote:
> I have a RNA-Seq alignment, in the form of a BAM file, where the reads are aligned to gene models, which correspond to the stable ids of gene features in my local Ensembl instance.
> Is there any way to display the BAM in Ensembl? i.e. is it possible to display this BAM which references Ensembl gene features, not the scaffolds or contigs in a Ensembl assembly/coordinate-system?
> Ben Warren
> Research Technologist
> F: +64 9 925 7001
> ben.warren at plantandfood.co.nz
> The New Zealand Institute for Plant & Food Research Limited
> Postal Address: Plant & Food Research Mt Albert
> Private Bag 92169, Auckland, 1142, New Zealand
> Physical Address: Plant & Food Research Mt Albert
> 120 Mt Albert Road, Sandringham, Auckland, 1025, New Zealand
> The contents of this e-mail are confidential and may be subject to legal privilege.
> If you are not the intended recipient you must not use, disseminate, distribute or
> reproduce all or any part of this e-mail or attachments. If you have received this
> e-mail in error, please notify the sender and delete all material pertaining to this
> e-mail. Any opinion or views expressed in this e-mail are those of the individual
> sender and may not represent those of The New Zealand Institute for Plant and
> Food Research Limited.
> Dev mailing list Dev at ensembl.org
> Posting guidelines and subscribe/unsubscribe info: http://lists.ensembl.org/mailman/listinfo/dev
> Ensembl Blog: http://www.ensembl.info/
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Dev