Schema for Problematic Regions - Problematic Regions for NGS or Sanger sequencing or very variable regions
  Database: hg19    Primary Table: getRmNgsDeadZone Data last updated: 2020-04-23
Big Bed File Download: /gbdb/hg19/bbi/problematic/deadZone.bb
Item Count: 4,262
The data is stored in the binary BigBed format.

Format description: Browser extensible data (<=12 fields)
fieldexampledescription
chromchr1Chromosome (or contig, scaffold, etc.)
chromStart206139240Start position in chromosome
chromEnd206139521End position in chromosome
start206139305 start
end206139456 end
geneFAM72A gene
plusminus65bpexonlength282 plus minus 65bp exon length
exonlength152 exon length
averagemappabilityl250m00.481678450355 average mappability l250m0
exonmappbelow1l250m01yes0noTrue exon mapp below 1 l250m0 (1=yes, 0=no)
positionsmappbelow1l250m0282 positions mapp below 1 l250m0
percpositionsmappbelow1l250m01 %positions mapp below 1 l250m0
maxcontigousbasesl250m0282 max contigous bases l250m0

Sample Rows
 
chromchromStartchromEndstartendgeneplusminus65bpexonlengthexonlengthaveragemappabilityl250m0exonmappbelow1l250m01yes0nopositionsmappbelow1l250m0percpositionsmappbelow1l250m0maxcontigousbasesl250m0
chr1206139240206139521206139305206139456FAM72A2821520.481678450355True2821282
chr1206139240206139521206139305206139456FAM72C2821520.481678450355True2821282
chr1206141381206141588206141446206141523FAM72A208780.354166375True2081208
chr1206141381206141588206141446206141523FAM72C208780.354166375True2081208
chr1206574715206575069206574780206575004SRGAP23552250.5True313
chr1206578543206578774206578608206578709SRGAP22321020.5True2321232
chr1206578543206578774206578608206578709SRGAP2B2321020.5True2321232
chr1206578543206578774206578608206578709SRGAP2C2321020.5True2321232
chr1206578543206578774206578608206578709SRGAP2D2321020.5True2321232
chr1207700301207700533207700366207700468CR12331030.5True2331233

Problematic Regions (problematic) Track Description
 

Description

This track helps call out sections of the genome that often cause problems for bioinformaticians. The 12 subtracks identify genomic regions known to cause analysis artifacts for common sequencing downstream computations, such as alignment, variant calling, or peak calling. The underlying data was imported from the NCBI GeT-RM, the Genome-in-a-Bottle, and Anshul Kundaje's ENCODE Blacklist projects.

The only exception is the UCSC Unusual Regions subtrack, which contains annotations of a few special gene clusters (IGH, IGL, PAR1/2, TCRA, TCRB, etc) and fixed sequences, alternate haplotypes, unplaced contigs, pseudo-autosomal regions, and mitochondria. These loci can yield alignments with low-quality mapping scores and discordant read pairs. This data set was manually curated, based on the Genome Browser's assembly description, the FAQs about assembly, and the NCBI RefSeq "other" annotations track data.

The ENCODE Blacklist subtrack contains a comprehensive set of regions which are troublesome for high-throughput Next-Generation Sequencing (NGS) aligners. These regions tend to have a very high ratio of multi-mapping to unique mapping reads and high variance in mappability due to repetitive elements such as satellite, centromeric and telomeric repeats.

The Genome-In-A-Bottle (GIAB) track set contains defined regions where it is difficult to make a confident call, due to low coverage, systematic sequencing errors, and local alignment problems. These regions were identified from sequencing data generated by multiple technologies.

The NCBI GeT-RM, Genetic Testing Reference Materials, track set contains highly homologous gene- and exon-level regions difficult or impossible to analyze with standard Sanger or short-read NGS approaches and are relevant to current clinical testing.

Display Conventions and Configuration

Each track contains a set of regions of varying length with no special configuration options. The UCSC Unusual Regions track has a mouse-over description, all other tracks have at most a name field, which can be shown in pack mode. The tracks are usually kept in dense mode.

The Hide empty subtracks control hides subtracks with no data in the browser window. Changing the browser window by zooming or scrolling may result in the display of a different selection of tracks.

Data access

The raw data can be explored interactively with the Table Browser or the Data Integrator.

For automated download and analysis, the genome annotation is stored in bigBed files that can be downloaded from our download server. Individual regions or the whole genome annotation can be obtained using our tool bigBedToBed which can be compiled from the source code or downloaded as a precompiled binary for your system. Instructions for downloading source code and binaries can be found here. The tool can also be used to obtain only features within a given range, e.g. bigBedToBed http://hgdownload.soe.ucsc.edu/gbdb/hg19/bbi/problematic/deadZone.bb -chrom=chr21 -start=0 -end=100000000 stdout

Methods

Files were downloaded from the respective databases and converted to bigBed format. The procedure is documented in our hg19 makeDoc file (search problematic).

Credits

Thanks to Anna Benet-Pages, Max Haeussler, Angie Hinrichs, and Daniel Schmelter at the UCSC Genome Browser for planning, building, and testing these tracks. The underlying data comes from the ENCODE Blacklist, the GeT-RM, and the Genome-in-a-Bottle projects.

References

Amemiya HM, Kundaje A, Boyle AP. The ENCODE Blacklist: Identification of Problematic Regions of the Genome. Sci Rep. 2019 Jun 27;9(1):9354. PMID: 31249361; PMC: PMC6597582

Zook JM, Chapman B, Wang J, Mittelman D, Hofmann O, Hide W, Salit M. Integrating human sequence data sets provides a resource of benchmark SNP and indel genotype calls. Nat Biotechnol. 2014 Mar;32(3):246-51. PMID: 24531798

Mandelker D, Schmidt RJ, Ankala A, McDonald Gibson K, Bowser M, Sharma H, Duffy E, Hegde M, Santani A, Lebo M et al. Navigating highly homologous genes in a molecular diagnostic setting: a resource for clinical next- generation sequencing. Genet Med. 2016 Dec;18(12):1282-1289. PMID: 27228465