Description

This track indicates any pair of exactly identical sequence for the 24 Jan 2018 Brachypodium distachyon/GCF_000005505.3_Brachypodium_distachyon_v3.0 genome assembly.

There may be two tracks in this composite collection:

The Gap Overlaps is thus a subset of the full Tandem Dups track.

This investigation began when an unusual number of paired sequences around gaps was noticed during the mouse strain sequencing project. This naturally raised the question, how common is this feature, and what type of assemblies can it be found in.

The Gap Overlaps track indicates any pair of exactly identical sequence on each side of gaps. Where a gap is any run of N's, including a single N. The end of an upstream sequence before the gap is duplicated exactly at the beginning of the downstream sequence following the gap in the assembly.
Data in track: .

The Tandem Dups track is a similar survey over the entire genome assembly. The separation gap between these paired sequences can range from 1 base up to 20,000 bases.
Data in track: Item count: 71,157; Bases covered: 67,682,473.

Methods

The Gap Overlap duplicate sequences were found by extracting 1,000 bases before and after each gap and aligned to each other with the blat command:

  blat -q=dna -minIdentity=95 -repMatch=10 upstreamContig.fa downstreamContig.fa
Filtering the PSL output for a perfect match, no mis-matches, and therefore of equal size matching sequence, where the alignment ends exactly at the end of the upstream sequence, and begins exactly at the start of the downstream sequence.

The Tandem Dups paired sequences were found with the following procedure:

The reason for starting with 29 base sized pairs and then selecting results of at least 30 base sized pairs results in a reasonable number of 30 base pairs. If the procedure starts with 30 base sized pairs, it produces way too many 30 base kmer pairs for a reasonable count.

See Also

Interactive tables of all results:

Credits

Thank you to Joel Armstrong and Benedict Paten of the Computational Genomics Lab at the U.C. Santa Cruz Genomics Institute for identifying this characteristic of genome assemblies.

The data and presentation of this track were prepared by Hiram Clawson, U.C. Santa Cruz Genomics Institute