7 ms·
https://youtu.be/fCd6B5HRaZ8 https://youtu.be/fCd6B5HRaZ8 is the best visualization of how the most popular type of DNA sequencer works (that I've found). Imag
by shpx 9y ago
https://youtu.be/fCd6B5HRaZ8 https://youtu.be/fCd6B5HRaZ8 is the best visualization of how the most popular type of DNA sequencer works (that I've found).
Imagine you have a string of length 3 billion made by randomly choosing from 4 characters. Like this
dna = ''.join(random.choices('atgc', weights=[30.9, 29.4, 19.9, 19.8], k=3_234_830_000))
you get to randomly sample 1 billion[3, page 7] overlapping substrings of length 200[3, page 7] with .1% of the characters randomly changed[3, page 8]. Trying to find the original string from this is technically an undecidable problem. If there's a sequence 400 characters long that repeats multiple times, how could you know if it repeats 5 times or 50 times? (this would be unlikely to happen with random.choices() but DNA isn't random). This is called sequence alignment and it's one of the hard problems in bioinformatics[4].
[0] https://docs.python.org/3/library/random.html https://docs.python.org/3/library/random.html random.choices() was added in Python 3.6
[1] http://www.biology-pages.info/B/BasePairing.html http://www.biology-pages.info/B/BasePairing.html source for `weights`
[2] https://en.wikipedia.org/wiki/Human_genome https://en.wikipedia.org/wiki/Human_genome source for 3_234_830_000 (python ignores underscores in numbers)
[3] http://sci-hub.io/10.1111/j.1755-0998.2011.03024.x http://sci-hub.io/10.1111/j.1755-0998.2011.03024.x illumina is the most popular producer of genome sequencers
[4] https://news.ycombinator.com/item?id=5123377 https://news.ycombinator.com/item?id=5123377
- nether 9y agoSide note: your comment is one of the few times I've seen the code tag used for actual code on HN, rather than quotations or just indentation.
- Bromskloss 9y agoYeah, why does Hacker News not have a real way to indent things? It wouldn't make loading the page any less light-weight.
- labster 9y agoSeriously. I can't go a week in HN without reading someone's complaint about not being able to read a quote, probably on mobile. If users keep complaining regularly, it's a problem with the software interface, not the users.
- nether 9y agoI emailed dang about it a while (years) ago. He said he's added it to the fix list, after initially saying that having a quote tag would mess up the site's "character." You may want to email him too.
- lolreally1 9y agoI fail to see how a quote tag would drive away the throngs of Silicon Valley elites, people who begin sentences with “I’m not an expert, but...”, hardcore Randians, and Paul Graham sycophants.
- sspiff 9y agoAnd yet here you are, on the site whose users you are attacking, commenting on an off-topic comment thread. Perhaps you should rethink your life choices?
- technofiend 9y agoYou forgot the folks who turn every article on the latest Google phone into a discussion on Apple products. See the latest posting about the Pixel 2 as an example.
- cup-of-tea 9y agoNobody said anything about mobile. Am I missing something?
- deleted 9y ago[deleted]
- deleted 9y ago[deleted]
- baxtr 9y agoTrue, and unfortunately it does not render it properly, it‘s cut after 20 chars or so
- deleted 9y ago[deleted]
- neikos 9y agoYou should be able to scroll it.
- baxtr 9y agoThanks, you’re right. I wasn’t aware of that. The experience is rather crappy though. Cumbersome to scroll on a smartphone screen and I’d love to see to entire statement not parts of only at a time.
- KGIII 9y ago> dna = ''.join(random.choices('atgc', weights=[30.9, 29.4, 19.9, 19.8], k=3_234_830_000)) On a tablet, I highlight and pull down. It grabs it all and I can then read it. It is a very awful workaround that impacts the site's usability.
- bsmrea 9y agoThis problem is called sequence assembly or genome assembly, not sequence alignment (which is a related problem).
- deleted 9y ago[deleted]
- teekert 9y agoVery nice movie indeed. This Sequencing by Synthesis is indeed the most popular method, the words in this type of DNA sequencing are usually 150 base bairs long, but by using fragments that can be much longer (say 2500 base pairs) and sequencing from both ends (as the movie explains) one can derived information on the position of 2 150 base pair long "words" relative to each other. The alternative technologies, sometimes called third generation sequencing are nanopore based (PacBio and Oxford Nanopore technologies), they work by pulling long DNA molecules through tiny holes and read the sequence directly, more specifically, they read a couple of base pairs at a time and infer the sequence from the current between the sides of the pores. This is very error prone and the (around 1 in 10?) random errors are corrected by sequencing every part of the genomic area of interest multiple times. To be clear, most sequencing experiments (in research and in the hospital) are designed to look at variations we can understand and thus the "words" (we call them reads) are aligned against "the" reference genome (as opposed to assembling de novo), which is the genome for as far as we understand it (the last official version is from 2013 but is updated regularly with new parts). Typically, a significant portion of the reads a sequencer produces cannot be mapped to the reference, probably as a results of unmapped regions indeed or perhaps deviations between the patient and the reference genome, this is of course even worse in patients with tumors with unstable genomes. This field still has a lot of challenges to tackle.
- cup-of-tea 9y ago> with .1% of the characters randomly changed Nope, not randomly. Having a random error profile is the best case. In reality these machines suffer from systematic errors, especially the Illumina ones. That means that they always fail to read through certain pieces of DNA, or make the same error every time.