This track shows regions of this target genome (sheep - 2023-07-20 - University of Idaho) that has alignment to other query genomes ("chain" subtracks) or in synteny ("net" subtracks). The alignable parts are shown with thick blocks that look like exons. Non-alignable parts between these are shown like introns.
Other query genome assemblies aligning to this target genome assembly:
chains | syntenic | reciprocal best | lift over | common name | assembly |
---|---|---|---|---|---|
99.540 | 99.404 | 99.005 | 99.527 | sheep | GCF_016772045.1_ARS-UI_Ramb_v2.0 |
98.844 | 98.480 | 97.854 | 98.804 | Asiatic mouflon | GCA_014523465.1_CAU_Oori_1.0 |
97.846 | 96.627 | 96.015 | 97.672 | goat | GCF_001704415.2_ARS1.2 |
95.140 | 93.554 | 92.941 | n/a | Sheep | oviAri3 |
90.957 | 87.075 | 87.049 | 90.422 | cattle | GCF_002263795.3_ARS-UCD2.0 |
88.366 | 86.319 | 85.864 | 88.133 | Cow | bosTau9 |
26.058 | 24.397 | 24.824 | n/a | Mouse | mm10 |
The chain tracks shows alignments of the other genome assemblies to the sheep//hive/data/genomes/asmHubs/refseqBuild/GCF/016/772/045/GCF_016772045.2_ARS-UI_Ramb_v3.0/html/GCF_016772045.2_ARS-UI_Ramb_v3.0.names.tab/GCF_016772045.2/2023-07-20 genome using a gap scoring system that allows longer gaps than traditional affine gap scoring systems. It can also tolerate gaps in both query and target genomes 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. 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.
There could be four different types of chain tracks:
The alignment track shows the net derived from the chain data in the format of a pair-wise side by side alignment. The net file is converted to the MAF format for this display.
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.
At base level in full display mode, this track will show the sequence of query as it aligned to target. When the view is too large to show such detail, blocks of alignments will show corresponding alignments to other chromosomes with colors indicating other chromosomes.
The query genome was aligned to target genome with lastz. The resulting alignments were converted into axt format using the lavToAxt program. The axt alignments were fed into axtChain, which organizes all alignments between a single query chromosome and a single target chromosome into a group and creates a kd-tree out of the gapless subsections (blocks) of the alignments. A dynamic program was then run over the kd-trees to find the maximally scoring chains of these blocks.
Chains were derived from lastz alignments, using the methods described on the chain tracks description pages, and sorted with the highest-scoring chains in the genome ranked first. The program chainNet was then used to place the chains one at a time, trimming them as necessary to fit into sections not already covered by a higher-scoring chain. During this process, a natural hierarchy emerged in which a chain that filled a gap in a higher-scoring chain was placed underneath that chain. The program netSyntenic was used to fill in information about the relationship between higher- and lower-level chains, such as whether a lower-level chain was syntenic or inverted relative to the higher-level chain. The program netClass was then used to fill in how much of the gaps and chains contained Ns (sequencing gaps) in one or both species and how much was filled with transposons inserted before and after the two organisms diverged.
The resulting net file was converted to axt format via netToAxt, then converted to maf format via axtToMaf, then converted to the bigMaf format with mafToBigMaf and bedToBigBed
lastz was developed by Robert Harris, Pennsylvania State University.
The axtChain program was developed at the University of California at Santa Cruz by Jim Kent with advice from Webb Miller and David Haussler.
The browser display and database storage of the chains and nets were created by Robert Baertsch and Jim Kent.
The chainNet, netSyntenic, and netClass programs were developed at the University of California Santa Cruz by Jim Kent.
Harris, R.S. (2007) Improved pairwise alignment of genomic DNA Ph.D. Thesis, The Pennsylvania State University
Chiaromonte F, Yap VB, Miller W. Scoring pairwise genomic sequence alignments. Pac Symp Biocomput. 2002:115-26. PMID: 11928468
Kent WJ, Baertsch R, Hinrichs A, Miller W, Haussler D. Evolution's cauldron: duplication, deletion, and rearrangement in the mouse and human genomes. Proc Natl Acad Sci U S A. 2003 Sep 30;100(20):11484-9. PMID: 14500911; PMC: PMC208784
Schwartz S, Kent WJ, Smit A, Zhang Z, Baertsch R, Hardison RC, Haussler D, Miller W. Human-mouse alignments with BLASTZ. Genome Res. 2003 Jan;13(1):103-7. PMID: 12529312; PMC: PMC430961