| 查看: 2263 | 回复: 7 | |||
amo11000金虫 (小有名气)
|
[求助]
编写一个程序处理文本内容，删除长度<10，开头的，不以CTG结尾的那些序列
|
|
我有一个文本文件1.txt，内容如下：>序列名字1（序列名字可以任意长度）ATGgcggggccccATAGGGGGGTTCCUTG>序列名字2ATGgcggggccccAgTAGGGGGGTTCCUTG>序列名字3CTGgggggggccccccgggggutg。。。。。。后面还有很多序列 要求删除那些长度<10，不以ATG开头的，不以CTG结尾的那些序列。名字不包括在序列的长度内，最后得到的文件内容应该是： >序列名字1（序列名字可以任意长度）ATGgcggggccccATAGGGGGGTTCCUTG>序列名字2ATGgcggggccccAgTAGGGGGGTTCCUTG。。。。。。后面还有很多序列。 由于序列3的开头字符不是ATG，而是CTG，所以被删除掉了。 谁能帮忙编写一个vb程序处理上面的问题啊？谢谢了 |
回复此楼» 猜你喜欢
filecode=后面第一个是大写字母
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有81人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
感觉是下周放榜了
已经有7人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 删除vip讨论子版的建议
已经有5人回复
- 自科青基中（三）经费申请说明 可否删除？
已经有5人回复
- HPLC延长柱前管路长度会对谱图有何影响？
已经有13人回复
- 修稿中审稿专家要求considerly shorten论文，修改时可以删除一部分内容吗？
已经有4人回复
- 青年基金申请书正文行间距及自动生成的各章节括号内说明内容能否删除等问题
已经有8人回复
lt292
金虫 (正式写手)
Holy Shit
- 应助: 21 (小学生)
- 金币: 1256.1
- 散金: 50
- 红花: 3
- 帖子: 591
- 在线: 36.4小时
- 虫号: 668772
- 注册: 2008-12-05
- 性别: GG
2楼2012-02-28 23:14:03
lurencyj
木虫 (著名写手)
- 应助: 159 (高中生)
- 金币: 2869.2
- 散金: 520
- 红花: 8
- 沙发: 10
- 帖子: 1244
- 在线: 148.3小时
- 虫号: 888093
- 注册: 2009-10-29
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2012-02-28 23:15:32
4楼2012-02-29 14:37:40
5楼2012-02-29 15:28:50
jackie1179
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1929.2
- 红花: 1
- 帖子: 460
- 在线: 47.2小时
- 虫号: 1119284
- 注册: 2010-10-11
- 专业: 生物信息学
6楼2012-03-01 09:49:25
7楼2012-03-01 10:48:33
christina1213
木虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 4181.9
- 散金: 20
- 帖子: 658
- 在线: 118小时
- 虫号: 599542
- 注册: 2008-09-11
- 专业: 生物信息学
【答案】应助回帖
★
感谢参与，应助指数 +1
jjdg(金币+1): 感谢参与 2012-03-09 03:03:14
感谢参与，应助指数 +1
jjdg(金币+1): 感谢参与 2012-03-09 03:03:14
|
用perl脚本吧，很方便，我写了一个大概，要根据自己的序列格式更改一下（也可能有一些小错误），你的格式应该是fasta格式的测序序列，>文件名应该是单独一行。 #!/usr/bin/perl use strict; use warnings; open (INFILE,"<1.txt"  ;open (OUTFILE,">>out.txt" ;while (my $line=<> { my $num=length($line); if (($num>9)&&(~m/^ATG/)&&(CTG$)) { print OUTFILE $line; } if (~m/^>/) { print OUTFILE $line; #这个地方要加一个怎样判断序列名称和你输出序列 的对应的语句。不然每个文件名都会输出，时间有限，我就这样简写了，抱歉了…… } } close INFILE; close OUTFILE; |
8楼2012-03-08 20:08:59
25