News for hackers. Stuff that matters.
Re:No human genome has ever been completely sequenced by shpx on Monday October 23, 2017 @01:06PM (#15534731)
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 random.choices() was added in Python 3.6

[1] http://www.biology-pages.info/B/BasePairing.html source for `weights`

[2] 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 illumina is the most popular producer of genome sequencers

[4] https://news.ycombinator.com/item?id=5123377

Re:No human genome has ever been completely sequenced by bsmrea on Monday October 23, 2017 @03:22PM (#15536158)
This problem is called sequence assembly or genome assembly, not sequence alignment (which is a related problem).
Re:No human genome has ever been completely sequenced by teekert on Tuesday October 24, 2017 @01:58AM (#15539545)
Very 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.

Re:No human genome has ever been completely sequenced by cup-of-tea on Tuesday October 24, 2017 @07:18AM (#15540755)
> 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.

Re:No human genome has ever been completely sequenced by nether on Monday October 23, 2017 @02:06PM (#15535414)
Side 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.
Re:No human genome has ever been completely sequenced by Bromskloss on Monday October 23, 2017 @03:11PM (#15536054)
Yeah, why does Hacker News not have a real way to indent things? It wouldn't make loading the page any less light-weight.
Re:No human genome has ever been completely sequenced by labster on Monday October 23, 2017 @03:43PM (#15536334)
Seriously. 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.
Re:No human genome has ever been completely sequenced by nether on Monday October 23, 2017 @07:52PM (#15537974)
I 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.
[dead]
Re:No human genome has ever been completely sequenced by technofiend on Tuesday October 24, 2017 @02:04AM (#15539573)
You 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.
Re:No human genome has ever been completely sequenced by sspiff on Tuesday October 24, 2017 @01:24AM (#15539422)
And 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?
Re:No human genome has ever been completely sequenced by cup-of-tea on Tuesday October 24, 2017 @07:21AM (#15540768)
Nobody said anything about mobile. Am I missing something?
Re:No human genome has ever been completely sequenced by baxtr on Monday October 23, 2017 @02:20PM (#15535542)
True, and unfortunately it does not render it properly, it‘s cut after 20 chars or so
Re:No human genome has ever been completely sequenced by neikos on Monday October 23, 2017 @02:40PM (#15535744)
You should be able to scroll it.
Re:No human genome has ever been completely sequenced by baxtr on Monday October 23, 2017 @03:54PM (#15536406)
Thanks, 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.
Re:No human genome has ever been completely sequenced by KGIII on Monday October 23, 2017 @04:25PM (#15536666)
> 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.