Schema for Pairwise Alignments - Human Genomes, Chain/Net pairwise alignments, as mapped by the HPRC project
  Database: hg38    Primary Table: chainHprcGCA_018472595v1    Row Count: 377   Data last updated: 2023-10-14
Format description: Summary info about a chain of alignments
On download server: MariaDB table dump directory
fieldexampleSQL type info description
bin 590smallint(5) unsigned range Indexing field to speed chromosome range queries.
score 1421double range score of chain
tName chr1varchar(255) values Target sequence name
tSize 248956422int(10) unsigned range Target sequence size
tStart 683250int(10) unsigned range Alignment start position in target
tEnd 684060int(10) unsigned range Alignment end position in target
qName JAHBCB010000139.1varchar(255) values Query sequence name
qSize 765627int(10) unsigned range Query sequence size
qStrand +char(1) values Query strand
qStart 8915int(10) unsigned range Alignment start position in query
qEnd 8981int(10) unsigned range Alignment end position in query
id 783008int(10) unsigned range chain id

Sample Rows
 
binscoretNametSizetStarttEndqNameqSizeqStrandqStartqEndid
5901421chr1248956422683250684060JAHBCB010000139.1765627+89158981783008
969733225chr12489564226840651427495JAHBCB010000139.1765627+14072765344765808
985024120chr124895642214295942325147JAHBCB010000153.1900590+0899176764294
1736074381chr1248956422232724310186147JAHBCB010000093.18145367-08145366710591
12772238411chr12489564221018812039665701JAHBCB010000031.129788546+49829788540478171
7115688192chr12489564221655286816616314JAHBCB010000031.129788546+65239066598508780501
883282042chr12489564221661631416680665JAHBCB010000031.129788546+63800896477006781195
788513787chr12489564222664133126646733JAHBCB010000031.129788546-1304688213052283782467
11157187171chr12489564223966635151905625JAHBCB010000026.112235465-112235465690321
01697224332chr12489564225190715369907788JAHBCB010000063.118004768-018004768642663

Note: all start coordinates in our database are 0-based, not 1-based. See explanation here.

Pairwise Alignments (hprcChainNet) Track Description
 

Description

This track shows regions of the human genome that are alignable to other Homo sapiens genomes. The alignable parts are shown with thick blocks that look like exons. Non-alignable parts between these are shown with thin lines like introns. More description on this display can be found below.

Other assemblies included in this track are from the HPRC project.

Display Conventions and Configuration

Chain Track

The chain track shows alignments of the human genome to other Homo sapiens genomes using a gap scoring system that allows longer gaps than traditional affine gap scoring systems. It can also tolerate gaps in both source and target assemblies simultaneously. These "double-sided" gaps can be caused by local inversions and overlapping deletions in both species.

The chain track displays boxes joined together by either single or double lines. The boxes represent aligning regions. Single lines indicate gaps that are largely due to a deletion in the query assembly or an insertion in the target assembly. assembly. Double lines represent more complex gaps that involve substantial sequence in both species. This may result from inversions, overlapping deletions, an abundance of local mutation, or an unsequenced gap in one species. In cases where multiple chains align over a particular region of the target genome, the chains with single-lined gaps are often due to processed pseudogenes, while chains with double-lined gaps are more often due to paralogs and unprocessed pseudogenes.

In the "pack" and "full" display modes, the individual feature names indicate the chromosome, strand, and location (in thousands) of the match for each matching alignment.

By default, the chains to chromosome-based assemblies are colored based on which chromosome they map to in the aligning organism. To turn off the coloring, check the "off" button next to: Color track based on chromosome.

To display only the chains of one chromosome in the aligning organism, enter the name of that chromosome (e.g. chr4) in box next to: Filter by chromosome.

Methods

The bigChain files were obtained from the HPRC S3 bucket (Amazon Web Services). For more information about how the bigChain files were generated, please refer to the HPRC publication below.

Credits

Thank you to Glenn Hickey for providing the HAL file from the HPRC project.

References

Liao WW, Asri M, Ebler J, Doerr D, Haukness M, Hickey G, Lu S, Lucas JK, Monlong J, Abel HJ et al. A draft human pangenome reference. Nature. 2023 May;617(7960):312-324. DOI: 10.1038/s41586-023-05896-x; PMID: 37165242; PMC: PMC10172123

Hickey G, Monlong J, Ebler J, Novak AM, Eizenga JM, Gao Y, Human Pangenome Reference Consortium, Marschall T, Li H, Paten B. Pangenome graph construction from genome alignments with Minigraph-Cactus. Nat Biotechnol. 2023 May 10;. DOI: 10.1038/s41587-023-01793-w; PMID: 37165083; PMC: PMC10638906

Armstrong J, Hickey G, Diekhans M, Fiddes IT, Novak AM, Deran A, Fang Q, Xie D, Feng S, Stiller J et al. Progressive Cactus is a multiple-genome aligner for the thousand-genome era. Nature. 2020 Nov;587(7833):246-251. DOI: 10.1038/s41586-020-2871-y; PMID: 33177663; PMC: PMC7673649

Paten B, Earl D, Nguyen N, Diekhans M, Zerbino D, Haussler D. Cactus: Algorithms for genome multiple sequence alignment. Genome Res. 2011 Sep;21(9):1512-28. DOI: 10.1101/gr.123356.111; PMID: 21665927; PMC: PMC3166836