4 answers
In Ensembl, PAR regions are stored as Bio::EnsEMBL::AssemblyExceptionFeatures. You can fetch these via a Bio::EnsEMBL::DBSQL::AssemblyExceptionFeatureAdaptor.
# get DB adaptor $dba by usual means
my $aefa = $dba->get_AssemblyExceptionFeatureAdaptor();
# get $slice e.g. per chromosome
my @aefs = @{$aefa->fetch_all_by_Slice($slice)};
AssemblyExceptionFeatures also include haplotypes, so you may need to filter the features.
•
written
by
biobot 0.0.77.a.1099
There is a public Ensembl MySQL database at http://www.ensembl.org/info/data/mysql.html and file downloads at http://www.ensembl.org/info/data/ftp/index.html . If you use the MySQL instance, you will need to understand the database schema (which is documented). The Perl API hides the schema details.
There is a public Ensembl MySQL database described at ensembl.org/info/data/mysql.html and file downloads at ensembl.org/info/data/ftp/index.html . If you use the MySQL instance, you will need to understand the database schema (which is documented). The Perl API hides the schema details. – Keith James 0 secs ago
There is a public Ensembl MySQL database described at ensembl.org/info/data/mysql.html and file downloads at ensembl.org/info/data/ftp/index.html . If you use the MySQL instance, you will need to understand the database schema (which is documented). The Perl API hides the schema details.
There is a public Ensembl MySQL database described at http://www.ensembl.org/info/data/mysql.html and file downloads at http://www.ensembl.org/info/data/ftp/index.html . If you use the MySQL instance, you will need to understand the database schema (which is documented). The Perl API hides the schema details.
If you're just doing one-off queries, why not just go directly to Ensembl? For example, for my favorite test gene CDK2:
http://uswest.ensembl.org/Homo_sapiens/Gene/Summary?g=ENSG00000123374
Clearly shows the genomic position at the top:
Chromosome 12: 56,360,556-56,366,565 forward strand
Now, the hard part (for me) is knowing the Ensembl gene IDs for PAR1 and PAR2, since PAR2 is apparently not an official gene human gene symbol, and the PAR1 in Entrez Gene doesn't have a match in Ensembl. But presumably you can find your Ensembl Gene ID more easily?
(Incidentally, NCBI Entrez Gene also will give genomic coordinates, and they use the same assembly as Ensembl...)
Thanks Andrew, but when I said PAR1 and PAR2 I was referring to the Pseudo-Autosomal Regions rather than genes (http://en.wikipedia.org/wiki/Pseudoautosomal_region). Apologies for being unclear!
Log in to answer this question.
-
How come on T2T human genome assembly the PAR regions of chrY are bigger than the corresponding r…
written by Duarte Molha
-
How to remove XTR region from my vcf file.
written by MHSATO •
-
Pseudoautosomal regions in cow
written by misterie
-
BWA mem of a heavily masked reference index takes more time than the normal reference?
written by Duarte Molha
-
telomeric ends fasta
written by 2nelly
-
Blacklisted regions for hg38?
written by Endre Bakken Stovner
-
Ensembl treatment of pseudoautosomal regions (PAR)- give us your input
written by Giulietta - Ensembl Helpdesk
-
How Are Human Sex Chromosomes Handled When Mapping Reads?
written by robert •
-
Getting Pairwise Genome Alignments/Syntenic Regions From Ensembl Or Ucsc?
written by User 9334 •
-
Why Does Ensembl 59 Annotate A Gene In The Par-1 Region Of Chromosome Y?
written by Bio_X2Y
Thanks to the folks below for answers. The PAR regions (as well as haplotypes) are defined in the assembly_exception table of the Ensembl human core mysql database.