24小时热门版块排行榜    

Znn3bq.jpeg
查看: 988  |  回复: 7

denis2008

金虫 (正式写手)

[求助] 求助perl程序

各位编程高手,你们好,我目前刚学perl编程,现要编一个小程序,感觉有点困难,请求各位帮忙。
数据如下:
code1
CGTCTCGATc
TGCGCCGAGc
AACACTGGGc
CCCTGGGTCc
TACGTTTCAc
TGACGTGCAc
code2
AAAACAGGGc
CCGGGCTTGc
TAAATAGTTc
TACATCTTCc
code 3
GATC
ATCG
CGTA
输出格式:
code3+GGCACATTGATTTGGGAGTc+code1+code2+AGATCTGATGGCGCGAGGG
以及此序列的互补序列,分别输出到两个不同的文件中。
谢谢各位大侠帮忙!
回复此楼

» 收录本帖的淘帖专辑推荐

perl语言专栏

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

洗烧杯只是我的表面工作,其实我是研究僧。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mycaas

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
没看懂要做什么 code是代表序列还是只是简单的字符?
互补序列 tr reverse下就可以
2楼2012-09-28 16:07:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

denis2008

金虫 (正式写手)

引用回帖:
2楼: Originally posted by mycaas at 2012-09-28 16:07:04
没看懂要做什么 code是代表序列还是只是简单的字符?
互补序列 tr reverse下就可以

code代表的是它下面列出的各条序列。
洗烧杯只是我的表面工作,其实我是研究僧。
3楼2012-09-29 09:16:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scbible

木虫 (正式写手)

EPI the king of water

几年没用,忘记了,爱莫能助哈
个人BLOG:http://www.eoe.so
4楼2012-10-17 13:00:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木有鱼丸

铁杆木虫 (知名作家)

优秀版主优秀版主

【答案】应助回帖

#!usr/bin/perl
open(IN1,"E:\\code1.txt"||die"can't open the file:$!";
open(IN2,"E:\\code2.txt"||die"can't open the file:$!";
open(IN3,"E:\\code3.txt"||die"can't open the file:$!";
open(OUT1,">E:\\positive.txt"||die"can't open the file:$!";
open(OUT2,">E:\\negtive.txt"||die"can't open the file:$!";
$file1="GGCACATTGATTTGGGAGTC";
$file2="AGATCTGATGGCGCGAGGG";
@file1=split(//,$file1);
@file2=split(//,$file2);
foreach(@file1)
{
        if($_ eq "A"
        {
                $file3[$m]="T";
                $m++;
                }
                if($_ eq "T"
        {
                $file3[$m]="A";
                $m++;
                }
                if($_ eq "G"
        {
                $file3[$m]="C";
                $m++;
                }
                if($_ eq "C"
        {
                $file3[$m]="G";
                $m++;
                }
        }
        foreach(@file3)
        {
                $file3=$file3.$_;
                }
        foreach(@file2)
{
        if($_ eq "A"
        {
                $file4[$n]="T";
                $n++;
                }
                if($_ eq "T"
        {
                $file4[$n]="A";
                $n++;
                }
                if($_ eq "G"
        {
                $file4[$n]="C";
                $n++;
                }
                if($_ eq "C"
        {
                $file4[$n]="G";
                $n++;
                }
        }
                foreach(@file4)
        {
                $file4=$file4.$_;
                }
while(
{
        chomp($_);
        $code=$_;
        $code1=$code1.$_;
        }
        while(
{
        chomp($_);
        $code2=$code2.$_;
        }
                while(
{
        chomp($_);
        $code3=$code3.$_;
        }
print OUT1 "$code3$file1$code1$code2$file2\n";
print "$code3$file1$code1$code2$file2\n";
print OUT2 "$code3$file3$code1$code2$file4\n";
print "$code3$file3$code1$code4$file2\n";
你问我姓什么?姓章,章鱼丸。骗人?别逗了,我真姓章。。
5楼2012-10-22 21:46:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木有鱼丸

铁杆木虫 (知名作家)

)
6楼2012-10-22 21:48:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木有鱼丸

铁杆木虫 (知名作家)

优秀版主优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
denis2008: 金币+10, 有帮助 2012-10-23 11:20:58
我被小木虫打败了。那个表情是一个右括号。。。
你问我姓什么?姓章,章鱼丸。骗人?别逗了,我真姓章。。
7楼2012-10-22 21:48:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

denis2008

金虫 (正式写手)

引用回帖:
7楼: Originally posted by 木有鱼丸 at 2012-10-22 21:48:41
我被小木虫打败了。那个表情是一个右括号。。。

谢谢你哈
洗烧杯只是我的表面工作,其实我是研究僧。
8楼2012-10-23 11:20:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 denis2008 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 297,工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考博] 申博自荐 +4 Linxia林夏 2026-04-13 4/200 2026-04-19 19:55 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
信息提示
请填处理意见