Problem 969. Genome Sequence 003: DNA Sequence with random positioned segments
This Challenge series will evolve the complexity of Genome DNA Sequencing. DNA Sequencing and the Shot Gun Method will be naively simplified into Cody Challenges. Genome sizes is another interesting wiki page.
DNA is represented by symbols ACGT, which for Matlab will be encoded as 0123. The basic goal is to reconstruct the original serial string of ACGT given multiple short segments. Segments are gleaned from multiple copies of the Virus/Bacteria/Chromosome thus there are overlapping, duplicative, and flipped segments. There are potential errors and duplicative stretches in the created segments. Chromosome 20 in its 59,187,298 base pairs has a segment of 820 that is repeated in at least two locations. The data being non-random largely increase lengths of duplicative stretches.
Example: Genome = ACGTCGGCCATGGACATTACG
Given three overlapping pieces, ACGTCGGCCA,CTACAGGTACCG, and GACATTACG these can be readily seen to overlap and create the original if the middle is recognized as having asymetric overlaps with its adjacent segments.
ACGTCGGCCATGGACATTACG
ACGTCGGCCAsssssssssss
ssssCTACAGGTACCGsssss Middle
ssssssssssssGACATTACG
Genome_003 Challenge is to reconstruct a genome under near ideal error free segment creation conditions. No segments are reversed. The segments start at random locations.
- Segments start at random positions (Genome_003 change from 001)
- Genome length is unconstrained
- Length of each segment - 48
- All segments may overlap by 16 to 47 characters
- No errors in the segments
- Genome is random (No duplicate starts or ends for 16 symbols of segments)
- Segment order is scrambled
Input: segs, Array of M rows of 48 value segments. Values are [0, 1, 2, 3].
Output: Gout, Genome vector of values [0,1,2,3]
Example: [0 1 2 3 2; 1 2 3 2 2; 2 2 1 1 2] creates [0 1 2 3 2 2 1 1 2] W=5, Overlap=varies
Future: Flipped segments(002), Random Position of Segment start locations(003), Extra Segments, Phage Phi X174, Parallel Processing Simulation(Shot Gun Approach), Haemophilus Influenza, Sequence with Segment Errors, and Chromosome 20 with its 59M length using 100K 4K-segments
Solution Stats
Problem Comments
-
1 Comment
Are Mjaavatten
on 19 Jul 2021
This is the hardest of the four sequencing problems by @Richard Zapor
Solution Comments
Show commentsProblem Recent Solvers3
Suggested Problems
-
Remove the small words from a list of words.
1489 Solvers
-
13358 Solvers
-
Convert a numerical matrix into a cell array of strings
1764 Solvers
-
middleAsColumn: Return all but first and last element as a column vector
601 Solvers
-
279 Solvers
More from this Author308
Problem Tags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!