This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Centromere hg38

Good morning, I need centromere file of hg38. The one I download from ucsc is chunked with 110 lines instead of 24 lines. I wrote a script to merge them but if I take the positions and put them on ucsc genome browser it only falls on part of the centromere. I think then that I must have downloaded the wrong file.

Do you have an idea on where I can find the good file.

Thank you.

hg38 centromere
• 2,629 views • link
written by heureuse • heureuse

This one. enter image description here

For the script, for chr1 for example, if I have 3 rows, I will just take the smallest start and the greatest stop. That way I will have just one entry for chr1

• 0 views • link
written by heureuse • heureuse

You can use bedtools merge.

wget -q -O- https://genome.ucsc.edu/cgi-bin/hgTables\?hgsid\=1654900614_2cotpic8SktMadZRcx7tYZBbq6md\&boolshad.hgta_printCustomTrackHeaders\=0\&hgta_ctName\=tb_centromeres\&hgta_ctDesc\=table+browser+query+on+centromeres\&hgta_ctVis\=pack\&hgta_ctUrl\=\&fbQual\=whole\&fbUpBases\=200\&fbDownBases\=200\&hgta_doGetBed\=get+BED | 
sort -k1,1V -k2,2n | 
bedtools merge -d 1000000
• 0 views • link
written by barslmn • barslmn

0 answers

No answers yet.

Log in to answer this question.

More posts like this