随笔分类 - perl
摘要:perl 模块Getopt::Long使用示例 use strict; use warnings; use Getopt::Long; my @ARGV2 = @ARGV; my data ="x.txt"; my help; die("error
阅读全文
摘要:今天尝试用perl来画一个简单的图,根据数据我们这里画线图。使用的模块是SVG::TT::Graph,两年前更新过,相比perl模块动不动十几年没有新版本,这个模块算不错了,出图格式为SVG格式。 cpan仓库地址: https://metacpan.org/pod/SVG::TT::Graph::
阅读全文
摘要:perl脚本加密 这里只针对单个文件的perl程序做展示。 pp:https://metacpan.org/pod/pp PAR:https://metacpan.org/pod/PAR 模块安装: cpanm pp cpanm PAR::Filter::Crypto 编译加密: linux: pp
阅读全文
摘要:这里是指数组作为hash的value,即一个key对应多个值 这里利用perl中的特殊句柄DATA做示例用以备忘{perl 数据结构一旦复杂点,可读性急剧下降,坑🕳}。 demo: use strict; use warnings; my %ha; while(<DATA>){ chomp; my
阅读全文
摘要:
bam文件的第二列的flag值,代表了当前read比对到基因组序列的基本情况。flag值由2的0次方到2的11次方这11个数字组成,即1,2,4,8,16,32,64,128,256,512,1024,2048。每个数字表示的比对情况这里不赘述,当获得当前read的flag值时,如何拆解成flag最
阅读全文

摘要:这里就用perl自带的DATA句柄做示例。 1 #!/usr/bin/perl 2 use strict; 3 use warnings; 4 5 my %hash; 6 my @arr; 7 while(<DATA>){ 8 chomp; 9 my @line = split /\s+/,$_;
阅读全文
摘要:linux 下fasta序列按碱基显色 1 use strict; 2 use warnings; 3 4 my dary_gray = "\033[1;30m"; 6 my $light_gray = "\033[0;37m"; 7 my
阅读全文
摘要:打开多个文件时必备啊 1 use FileHandle; 2 3 my (%index,%fh); 4 # 创建句柄,这里利用gzip 压缩下 5 foreach my $k(keys %index){ 6 {k}_1.tx
阅读全文
摘要:substr EXPR,OFFSET,LENGTH,REPLACEMENT substr EXPR,OFFSET,LENGTH substr EXPR,OFFSET Extracts a substring out of EXPR and returns it. First character is
阅读全文
摘要:perl中的引用和C中的指针一样,用"\"标识,引用后可使用符号“->”取值。解引用则在对应的数据类型前加$,@ 或%。 这里这里用两数组求和做示例,引用传递的子函数的用法。 1 use strict; 2 3 my @s1=(1,13,4,25,6,27,8,22); 4 my @s2=(11,2
阅读全文
摘要:1 #!/usr/bin/perl -w 2 use strict; 3 use FindBin qw(Script); 4 5 my Bin; 6 print "the script dir:Script\n"; 7 8 use Cwd q
阅读全文
摘要:perl : while(defined(my f2=<FQ2>)){ condition } python3 : for line1,line2 in zip(fileHandle1,fileHandle2): condition
阅读全文
摘要:去除基因组序列中的未定位的scaffold、Contig序列和线粒体序,主要针对NCBI提供refseq基因组序列,组装到染色体级别的物种基本都通用。将所有碱基统一成大写字母,并计算每条染色体长度,每80个字符换行。处理脚本如下: 1 use strict; 2 open A,"ARGV[1]"; 4 open C,">$ARGV[2]...
阅读全文
摘要:目标如题,有多个fasta文件和一个文件名列表,将文件名列表中包含的文件匹配出来并提取第一条序列合并成一个fa文件。这个采用perl实现,用法和代码如下: 1 #!/usr/bin/perl -w 2 use strict; 3 4 sub usage{ 5 die "usage:perl $0 \n",unless(@ARGV==3); 6 } 7 usage(); 8 ...
阅读全文
摘要:这里记录两种perl数组去重的办法,一种利用哈希(hash),一种直接利用perl自带的模块List::MoreUtils内部的函数uniq。一、利用hash去重示例代码如下: 1 #!/usr/bin/perl -w 2 use strict; 3 4 my @list=qw /1 2 3 2 1 4 aa a bb c b bb d/; 5 foreach (@list){pri...
阅读全文