This is a test version of Biostars. For the public version, visit https://www.biostars.org.
known unmappable or poorly sequenced regions in hg38

Dear all,

is there any similar file to this, which known unmappable or poorly sequenced regions for hg38? Or any way to download it from ucsc table browser?

Thank you for any link or comments.

dna-seq bed ngs
• 4,945 views • link
updated by igor igor • written by Paul Paul

2 answers

The blacklist files are available from ENCODE directly (there is a drop-down menu to select the proper genome build): https://www.encodeproject.org/annotations/ENCSR636HFF/

I don't know if using LiftOver for these regions is a good idea, since their coordinates are uncertain in the first place.

• 0 views • link
written by igor igor

Thank you Igor. Did you validate this regions? I can see, that for build GRCh38 is there only 38 regions. For Hg19 it is 411 regions. Also when I check this version of BED, there are 1649 regions.

• 0 views • link
written by Paul Paul

You can use a combination of blacklisted regions from hg37 with the UCSC liftOver utility (or maybe someone else can recommend a blacklist for hg38).

• 0 views • link
written by harold.smith.tarheel harold.smith.tarheel

Thank you so much. I was thinking to use "wgEncodeDukeMapabilityRegionsExcludable.bed.gz" from here and then use UCSC liftOver, but is this good procedure? Isnt those regions "out of date" for HG38?

• 0 views • link
written by Paul Paul

Log in to answer this question.

More posts like this