Schema for Pairwise Alignments - Human Genomes, Chain/Net pairwise alignments, as mapped by the HPRC project
  Database: hg38    Primary Table: chainHprcGCA_018471105v1    Row Count: 348   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 9smallint(5) unsigned range Indexing field to speed chromosome range queries.
score 79325306double range score of chain
tName chr1varchar(255) values Target sequence name
tSize 248956422int(10) unsigned range Target sequence size
tStart 260780int(10) unsigned range Alignment start position in target
tEnd 1425449int(10) unsigned range Alignment end position in target
qName JAHALY010000158.1varchar(255) values Query sequence name
qSize 857007int(10) unsigned range Query sequence size
qStrand -char(1) values Query strand
qStart 847int(10) unsigned range Alignment start position in query
qEnd 857007int(10) unsigned range Alignment end position in query
id 747917int(10) unsigned range chain id

Sample Rows
 
binscoretNametSizetStarttEndqNameqSizeqStrandqStartqEndid
979325306chr12489564222607801425449JAHALY010000158.1857007-847857007747917
984260433chr124895642214302032317725JAHALY010000159.1888427-1888426746698
14671774989chr1248956422232514751906818JAHALY010000034.150401841-8050401840387767
851732569chr12489564221310425213122521JAHALY010000034.150401841+3938284339401112765577
7113334721chr12489564221654832616590456JAHALY010000034.150401841-1456857914634217765066
711164805chr12489564221656577516567539JAHALY010000034.150401841-1459697314598737765887
8813071919chr12489564221656577516746893JAHALY010000034.150401841-1444418615164688761024
7112707385chr12489564221661591016644329JAHALY010000287.128433-128432765226
11811426chr12489564221674689416862968JAHALY010000034.150401841-1497253415104427765501
788513918chr12489564222664133126646733JAHALY010000034.150401841+2530310125308502765781

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