Question: How Can I Get The Human Chromosome Centromere Position And Chromosome Length In Grch37/Hg19
22
gravatar for Jdk
9.8 years ago by
Jdk220
Jdk220 wrote:

This site (https://lists.soe.ucsc.edu/pipermail/genome/2003-November/003423.html) used to give a good answer, but it does not work now. I can not find ChrN_gap under the positional tables drop box. Does anyone know how to get centromere position in Table Browser?

position ucsc • 29k views
ADD COMMENTlink modified 20 months ago by Shixiang60 • written 9.8 years ago by Jdk220
23
gravatar for Pierre Lindenbaum
9.8 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum129k wrote:

The positions of the cytobands are available in : http://hgdownload.cse.ucsc.edu/goldenPath/hg18/database/cytoBand.txt.gz

curl -s "http://hgdownload.cse.ucsc.edu/goldenPath/hg18/database/cytoBand.txt.gz" | gunzip -c | grep acen | head
chr1    121100000    124300000    p11.1    acen
chr1    124300000    128000000    q11    acen
chr10    38800000    40300000    p11.1    acen
chr10    40300000    42100000    q11.1    acen
chr11    51400000    52900000    p11.11    acen
chr11    52900000    56400000    q11    acen
chr12    33200000    35400000    p11.1    acen
chr12    35400000    36500000    q11    acen
chr13    13500000    16000000    p11.1    acen
chr13    16000000    18400000    q11    acen
ADD COMMENTlink modified 19 months ago by RamRS27k • written 9.8 years ago by Pierre Lindenbaum129k

I used h18/build36, but of course you should use h19/build37 for your probem

ADD REPLYlink written 9.8 years ago by Pierre Lindenbaum129k

Wow! Thank you very much!

ADD REPLYlink written 9.8 years ago by Jdk220
12
gravatar for Jeremy Leipzig
9.8 years ago by
Philadelphia, PA
Jeremy Leipzig19k wrote:

http://genome.ucsc.edu/cgi-bin/hgTables

Just choose All Tables -> table:gap

Centromeres are treated as a gap "type"

#bin    chrom    chromStart    chromEnd    ix    n    size    type    bridge
23    chr1    121535434    124535434    1270    N    3000000    centromere    no

Chromosome length can be accessed from chromInfo

#chrom    size    fileName
chr1    249250621    /gbdb/hg19/hg19.2bit
chr2    243199373    /gbdb/hg19/hg19.2bit
chr3    198022430    /gbdb/hg19/hg19.2bit
chr4    191154276    /gbdb/hg19/hg19.2bit
chr5    180915260    /gbdb/hg19/hg19.2bit
chr6    171115067    /gbdb/hg19/hg19.2bit
ADD COMMENTlink modified 19 months ago by RamRS27k • written 9.8 years ago by Jeremy Leipzig19k
1

However, hg38 version cannot find centromeres in this way

ADD REPLYlink written 20 months ago by Shixiang60
1

As for hg38 data, you can download it through FTP annotation database directly.

ADD REPLYlink written 13 months ago by Timze W60

It works very well! Thank you!

ADD REPLYlink written 9.8 years ago by Jdk220
5
gravatar for Hypotheses
7.2 years ago by
Hypotheses90
Bangkok, Thailand
Hypotheses90 wrote:

http://genome.ucsc.edu/FAQ/FAQtracks.html#tracks20

Question:
"How do I find the positions of the centromeres and telomeres in a particular assembly?"

Response:
This information can be found in the "gap" database table. Use the Table Browser to extract it. To do this, select your assembly and the gap table, then click the "filter Create" button. Set the "type" field to centromere telomere (separated by a space). For help using the Table Browser, visit the User's Guide.

ADD COMMENTlink modified 19 months ago by RamRS27k • written 7.2 years ago by Hypotheses90
3
gravatar for Shixiang
20 months ago by
Shixiang60
Shanghai
Shixiang60 wrote:

I am also looking for this, and I find all what we need are included in package GWASTools https://www.rdocumentation.org/packages/GWASTools/versions/1.18.0/topics/centromeres

ADD COMMENTlink modified 20 months ago • written 20 months ago by Shixiang60
1
gravatar for 516356412
5.3 years ago by
51635641210
51635641210 wrote:

I have tried

curl -s "http://hgdownload.cse.ucsc.edu/goldenPath/dm3/database/cytoBand.txt.gz" | gunzip -c | grep acen | head

, but it didn't show any results.I also tried

http://genome.ucsc.edu/cgi-bin/hgTables

just choose All Tables -> table:gap

Centromeres are treated as a gap "type". The result is different from human genome, the type of the reads are contig not centromere.

#bin    chrom    chromStart    chromEnd    ix    n    size    type    bridge
594    chr4    1221288    1221388    138    N    100    contig    yes
587    chrU    325170    325270    141    N    100    contig    yes
588    chrU    429354    442354    143    N    13000    contig    yes
589    chrU    543815    543915    145    N    100    contig    yes
589    chrU    639309    639409    147    N    100    contig    yes

Lastly I tried the third method, it showed the message as follows:

#filter: (gap.type = 'telomere' OR gap.type = 'centromere')
#bin    chrom    chromStart    chromEnd    ix    n    size    type    bridge
# No results passing filter.
ADD COMMENTlink modified 19 months ago by RamRS27k • written 5.3 years ago by 51635641210
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1208 users visited in the last hour