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

Sample Rows
 
binscoretNametSizetStarttEndqNameqSizeqStrandqStartqEndid
971992464chr12489564222607831356910JAHBBZ010000189.1803069-12731803069732620
985029734chr124895642214286672325086JAHBBZ010000201.1897271-8897259730927
7531627463chr124895642223251472656440JAHBBZ010000181.1465475+1313341483736862
1704510269chr1248956422275282610188049JAHBBZ010000082.17546380-582897546378663308
13928404109chr12489564221019132851906853JAHBBZ010000033.142227865+142227157221471
7111286715chr12489564221655286816567539JAHBBZ010000033.142227865+65381746559355745378
8810069830chr12489564221656577516680665JAHBBZ010000033.142227865+62857586494621742453
788513918chr12489564222664133126646733JAHBBZ010000033.142227865-2530163525307036745855
03278740027chr12489564225190835886797362JAHBBZ010000076.134843587+034843369325603
05183394290chr124895642286800940168472103JAHBBZ010000075.173274481+117673274446113928

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