发布时间：2014-09-30 08:38:19

系谱分析（pedigree analysis）是了解遗传病的一个常用的方法。先对某家族各成员出现的某种遗传病的情况进行详细的调查，再以特定的符号和格式绘制成反映家族各成员相互关系和发生情况的图解，然后根据孟德尔定律对各成员的表现型和基因型进行分析。通过分析，可以判断某种性状或遗传病是属于哪一种遗传方式（单基因.........【阅读全文】

发布时间：2014-09-29 15:37:39

问题说明：假设有文件夹data，里面含有很多.TXT域名的文件，而文件里有特定字符需要匹配，并输出该行。#!/usr/bin/perl -wuse strict;opendir IN,"data" || die "can't open : $!\n";open OUT,">>outfile";my @file=grep (/\.TXT/,readdir IN);foreach my $list(@file){    open TWO,"data/$list" o.........【阅读全文】

发布时间：2014-09-28 16:51:47

#!/usr/bin/perl -wuse strict;system ("../samtools view -b -S 1a.sam >1a.bam");system ("../samtools view -b -S 1b.sam >1b.bam");system ("../samtools view -b -S 1c.sam >1c.bam");system ("../samtools view -b -S 1d.sam >1d.bam");system ("../samtools view -b -S 2a.sam >2a.bam");system ("../samtool.........【阅读全文】

发布时间：2014-09-28 16:49:52

#!/usr/bin/perl -wuse strict;print "input your file,then press Enter!\n";my $file=<STDIN>;my $key;my %seq;open IN,"<$file";open OUT,">resulta";   while (<IN>){     if (/>/){ $key=$_;}else{$seq{$key}.=$_;}   foreach my $target(values %seq){ &nb.........【阅读全文】

发布时间：2014-09-28 16:49:01

#!/usr/bin/perl -w#use strict;open OUT,">mac";print "please input your file name(match target),then press Enter.\n";while (){    chomp;          if (/>/){ $key=$_;       my       @m.........【阅读全文】

发布时间：2014-09-28 15:55:44

建库：makeblastdb -in my.fasta -dbtype nucl -title db -out db比对：blast* -query db -db db -out db -outfmt {1..10}......【阅读全文】

发布时间：2014-09-28 15:51:04

#!/usr/bin/perl -w$|=1;print "input your genome sequence,then press Enter!\n";my $seq=<STDIN>;open FASTA,"$seq" or die "can not open FASTA file,$!";print "input your gff file name,then press Enter!\n";my $gff=<STDIN>;open GFF,"$gff" or die "can not open GFF file,$!";open OUT,">gene_sequence".........【阅读全文】

发布时间：2014-09-28 15:45:58

文件说明：数据如下，我想提取第三列含gene的行的特定##gff-version 3Chr1    phytozome9_0    gene    3631    5899    .    +    .    ID=AT1G01010;Name=AT1G01010Chr1    .........【阅读全文】

发布时间：2014-09-28 14:43:46

#!/usr/bin/perl -w  #使用方法：perl transfer.pl file.gb ##文件的格式有三种$inputfile="<id.gb";$format1="genbank";$outputfile=">id.fasta";$format2="fasta";use Bio::SeqIO;         $in = Bio::SeqIO->new(-file => "$inputfile" ,    .........【阅读全文】

发布时间：2014-09-27 21:42:21

问题描述：比如我的数据（如下），我现在想提取表头里第二列数字小于99的序列，以下perl代码可以完成这个工作。>example 100 other tail  1   NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN >example 99 other tail  2   NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN.........【阅读全文】