+ -
当前位置:首页 → 问答吧 → 求助:如何用perl来连接字符串

求助:如何用perl来连接字符串

时间:2011-03-30

来源:互联网

我有306个文件,每个文件里面的内容类似于:
  1. #文件1:
  2. >seqAn1
  3. actgagcatgacagtagac
  4. agcatag
  5. >seqBn1
  6. actagacagatcagagatc
  7. acagat
  8. >seqCn1
  9. acagtagacgatgacagtc
  10. acatcg
复制代码
  1. #文件2:
  2. >seqAn2
  3. acagatgacagatagacagat
  4. acagatca
  5. >seqBn2
  6. acagatgacagatacacagat
  7. acacgatagacacagatagga
  8. cacgt
  9. >seqCn2
  10. acagatgacacagatgacaca
  11. acagatgacaacacagatgaa
  12. caccagagag--agcaagta
复制代码
  1. #文件306:
  2. >seqAn306
  3. acagatgacacagatgacagatagaca
  4. acagatgacagatagacagtagacagt
  5. acaagatgac
  6. >seqBn306
  7. acagatgacacagatgacacagatgac
  8. acagatgacacagtta
  9. >seqCn306
  10. acagatgaca---cagatgaca-acaga
  11. acacag
复制代码
我想达到这样的目的,将seq相同的序列合在一起,如:将>seqAn1,>seqAn2…>seqAn306下面的那些字符序列按顺序合在一起;

作者: xp198766   发布时间: 2011-03-30

利用 HASH了

作者: liyangole   发布时间: 2011-03-30

my ($title,%seq,@array);
while (<FH>) {
   chomp;
   if (/^>/) {
     $title = $_;
     print
   }
   else {
     $seq{$title} .= $_;
     print
   }
}

作者: leigh111   发布时间: 2011-03-30

回复 liyangole


    我是新手,也想到用hash,感觉不知道从哪里下手

作者: xp198766   发布时间: 2011-03-30

回复 xp198766


    请问,我的文件非常多,怎么都输入perl里面呢?一个一个点吗?

作者: xp198766   发布时间: 2011-03-30

为什么这么多搞基因的连门都没入就开始做事……

作者: zhlong8   发布时间: 2011-03-30

回复 zhlong8


    生物出身,想学习来着,非常初级的还行,稍微复杂点的就完不成了……还请各位包涵……

作者: xp198766   发布时间: 2011-03-30

热门下载

更多