作者plankton (我是没人要的小孩)
看板perl
标题Re: [请益] 打开档案作样式比对
时间Wed Feb 28 13:06:28 2007
※ 引述《tear2001 ( Play ball !!)》之铭言：
: 冒昧问一个问题 由於小弟是新手中的新手 太笨的话各位别生气 Q_Q
: 目前我写出一个比对用的script
: 才再测试阶段 但是不管怎麽比对都....orz
: $File = $ARGV[0];
: open( FHD, $File ) || die "$!\n"; #打开档案为唯读型式
: foreach $line (<FHD>) {
: if( $line =~ m/"Mycoplasma pneumoniae M129"\i/ ) {
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
我觉得你比对的地方怪怪的
: #对打开的档案内容和 Mycoplasma pneumoniae M129 进行比对
: print "bingo";
: }
: else {
: print "fuck";
: }
: }
这样因该可以
if ($line =~ /Mycoplasma pneumoniae M129/i)
1.你加""他也会比对""..你的档案里似乎没有这个东西
2. \i ?? 你是要大小写都可以吗? 那应该是写在 /xxpkfkdfk/ 之後
我也很逊 希望可以帮到你
: #比对成功则印出bingo 不然就fuck
: 我输入的档案长的是这样子:
: Query= gi|12045217|ref|NP_073028.1| acetate/propionate kinase
: [Mycoplasma genitalium G37]
: (393 letters)
: Database: refseq_gn
: 1,370,603 sequences; 432,705,335 total letters
: Searching..................................................done
: Score E
: Sequences producing significant alignments: (bits) Value
: acetate/propionate kinase [Mycoplasma pneumoniae M129] 600 e-170
: acetate/propionate kinase [Mycoplasma gallisepticum R] 425 e-118
: acetate/propionate kinase [Mycoplasma penetrans HF-2] 356 3e-097
: acetate kinase [Mycoplasma mycoides subsp. mycoides SC str. PG1] 316 3e-085
: acetate kinase [Mesoplasma florum L1] 315 8e-085
: .
: .
: .
: .
: 我原本预计会出现bingo的....但是却.......一直fuck....orz
: 实在找不差哪里错的 应该是我的脑子吧orz
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.129.160.66
1F：推 tear2001:谢谢 <(_ _)> 02/28 13:14
2F：推 tear2001:我再试试好了 02/28 13:17