Schema for Pairwise Alignments - Human Genomes, Chain/Net pairwise alignments, as mapped by the HPRC project
  Database: hg38    Primary Table: chainHprcGCA_018503245v1    Row Count: 446   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 78276547double range score of chain
tName chr1varchar(255) values Target sequence name
tSize 248956422int(10) unsigned range Target sequence size
tStart 257673int(10) unsigned range Alignment start position in target
tEnd 1421757int(10) unsigned range Alignment end position in target
qName JAHEOQ010000255.1varchar(255) values Query sequence name
qSize 846340int(10) unsigned range Query sequence size
qStrand -char(1) values Query strand
qStart 8183int(10) unsigned range Alignment start position in query
qEnd 846340int(10) unsigned range Alignment end position in query
id 853395int(10) unsigned range chain id

Sample Rows
 
binscoretNametSizetStarttEndqNameqSizeqStrandqStartqEndid
978276547chr12489564222576731421757JAHEOQ010000255.1846340-8183846340853395
984682239chr124895642214309912325147JAHEOQ010000298.1895280+0895073852362
9158313807chr124895642223266844099965JAHEOQ010000212.11988138+01988137841664
1469648494chr124895642241004519053453JAHEOQ010000139.14978532+64978531784823
10107699447chr1248956422905385610187403JAHEOQ010000245.11142629+11142629850540
13926098745chr12489564221018812051906693JAHEOQ010000047.142454376+36342454376157666
851728150chr12489564221310425213122521JAHEOQ010000047.142454376-3954212639560396871034
127076419chr12489564221654832616954364JAHEOQ010000047.142454376-3527201935905184861504
711427504chr12489564221654832816552863JAHEOQ010000047.142454376+64487566453280871424
788512336chr12489564222664133126646733JAHEOQ010000047.142454376-2531232825317738871337

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