Question: How Can I Get The Human Chromosome Centromere Position And Chromosome Length In Grch37/Hg19
13
gravatar for Jdk
6.9 years ago by
Jdk130
Jdk130 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 • 19k views
ADD COMMENTlink modified 13 months ago by Zeinab 0 • written 6.9 years ago by Jdk130
15
gravatar for Pierre Lindenbaum
6.9 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum96k 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 written 6.9 years ago by Pierre Lindenbaum96k

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

ADD REPLYlink written 6.9 years ago by Pierre Lindenbaum96k

Wow! Thank you very much!

ADD REPLYlink written 6.9 years ago by Jdk130
9
gravatar for Jeremy Leipzig
6.9 years ago by
Philadelphia, PA
Jeremy Leipzig17k 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 6.9 years ago • written 6.9 years ago by Jeremy Leipzig17k

It works very well! Thank you!

ADD REPLYlink written 6.9 years ago by Jdk130
4
gravatar for Hypotheses
4.2 years ago by
Hypotheses60
Bangkok, Thailand
Hypotheses60 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 written 4.2 years ago by Hypotheses60
1
gravatar for 516356412
2.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 written 2.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: 805 users visited in the last hour