作者alpe (薛丁格的猫)
看板BioMedInfo
标题Re: [问题] 请问有无方便的程式或方法可帮我解决难 …
时间Sat Aug 2 01:09:03 2008
※ 引述《enisx (东方有比目鱼,不比不行)》之铭言:
: 请教各位高手,
: 我手上有外送回来定序的DNA序列文字档(.txt)
: 里头大约有约700-1000条,档案已经转成FASTA格式了,
: 但是,我想从每条DNA序列中"撷取"出我想要的部份,
: 举例如下:
: -------------------------------------------------------------
: >123456
: TACCTCTGTGCTGTGGGGCGGATCTGAAAAGCTGGTCTTT
: → hgsfhevil:c++也可以用这样的东西去做吗? 08/01 14:03
: → alpe:晚上来 release win 下的执行档好了 08/01 18:13
http://www.uploadpedia.com/19Q94XWLMLGH/seq-rar.html
seq.rar => seq.exe
test
>123456
TCAXXXXXXXXXXXXXXXTTTTCAXXXXXXXXXXXXXXXTTTTCAXXXXXXXXXXXXXXXTTT
TCAXXXXXXXXXXXXXXXTTT
test2
>123456
TCAzzzzzzzzzzzzzzzTTTTCAzzzzzzzzzzzzzzzTTTTCAzzzzzzzzzzzzzzzTTT
TCAzzzzzzzzzzzzzzzTTT
D:\>seq '(TCA.*?TTT)' test test2
Create By Mervyn Wang 2008.07.31.
(TCA.*?TTT)
File: test
-------------------------
1 : TCAXXXXXXXXXXXXXXXTTT
2 : TCAXXXXXXXXXXXXXXXTTT
3 : TCAXXXXXXXXXXXXXXXTTT
4 : TCAXXXXXXXXXXXXXXXTTT
File: test2
-------------------------
1 : TCAzzzzzzzzzzzzzzzTTT
2 : TCAzzzzzzzzzzzzzzzTTT
3 : TCAzzzzzzzzzzzzzzzTTT
4 : TCAzzzzzzzzzzzzzzzTTT
好像应该要写一个 output 到file 的 参数的... ...
--
Exactly. For that one fraction of a second, you were open to options
you had never considered. THAT is the exploration that awaits you:
not mapping stars and studying nebulae,but
charting the unknown possibilities of existence.
Star Trek S7E26 "All Good Thing"
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 59.121.121.35
※ 编辑: alpe 来自: 59.121.121.35 (08/02 01:19)
1F:→ enisx:感谢您的方法! 08/03 18:22