4 ms·

Who knows the real story. This could be a case of the author taking journalistic liberties, or simply misunderstanding the process. My understanding is that a

Who knows the real story. This could be a case of the author taking journalistic liberties, or simply misunderstanding the process. My understanding is that a simple diff would take care of 90% of the work by eliminating all common base pairs. Then, from the results, his job would be to cross-reference each of the mutations against the gene database. Any unknown mutation could be the one he's looking for. Perhaps this is actually what he did, and the "no feasible way to automate the hunt" refers to the cross-referencing step, not the differencing step. Still though, reading that part almost made me email him saying, "I'll find you a way to automate it!".

a-priori 18y ago
So I just emailed Dr. Rienhoff, linking to this page, and encouraging him to post any problems he encounters to HN.
etal 18y ago
Here's his own telling of the story: http://mydaughtersdna.org/Members/hyrjr/the-reasoning-behind-the-dna-analysis http://mydaughtersdna.org/Members/hyrjr/the-reasoning-behind... An example of the zigzag pattern he probably received from the sequencing lab (his would have been much larger): http://www.ensembl.org/Homo_sapiens/Transcript/Summary?db=core;g=ENSG00000135503;r=12:50631753-50677124;t=ENST00000257963 http://www.ensembl.org/Homo_sapiens/Transcript/Summary?db=co... The reference sequence for the gene (this is what he was looking at): http://www.ensembl.org/Homo_sapiens/Transcript/Exons?db=core;g=ENSG00000135503;r=12:50631753-50677124;t=ENST00000257963 http://www.ensembl.org/Homo_sapiens/Transcript/Exons?db=core... and his manual diff (at his site): http://mydaughtersdna.org/Members/hyrjr/acvr1b http://mydaughtersdna.org/Members/hyrjr/acvr1b The usual way to do a gene diff is by aligning two sequences with a program like ClustalW (free, and packaged with Debian now). That shows the differences and positions in a format other programs can read. But I don't know how much annotation the lab gave him, and annotating a pile of introns and exons does take a long time, even with the right tools -- the diffs aren't necessarily helpful until the rest of that information is in place.

a-priori 18y ago

So I just emailed Dr. Rienhoff, linking to this page, and encouraging him to post any problems he encounters to HN.

etal 18y ago

Here's his own telling of the story: http://mydaughtersdna.org/Members/hyrjr/the-reasoning-behind-the-dna-analysis http://mydaughtersdna.org/Members/hyrjr/the-reasoning-behind... An example of the zigzag pattern he probably received from the sequencing lab (his would have been much larger): http://www.ensembl.org/Homo_sapiens/Transcript/Summary?db=core;g=ENSG00000135503;r=12:50631753-50677124;t=ENST00000257963 http://www.ensembl.org/Homo_sapiens/Transcript/Summary?db=co... The reference sequence for the gene (this is what he was looking at): http://www.ensembl.org/Homo_sapiens/Transcript/Exons?db=core;g=ENSG00000135503;r=12:50631753-50677124;t=ENST00000257963 http://www.ensembl.org/Homo_sapiens/Transcript/Exons?db=core... and his manual diff (at his site): http://mydaughtersdna.org/Members/hyrjr/acvr1b http://mydaughtersdna.org/Members/hyrjr/acvr1b The usual way to do a gene diff is by aligning two sequences with a program like ClustalW (free, and packaged with Debian now). That shows the differences and positions in a format other programs can read. But I don't know how much annotation the lab gave him, and annotating a pile of introns and exons does take a long time, even with the right tools -- the diffs aren't necessarily helpful until the rest of that information is in place.