Perl练习01之计算一条read的ATGC个数GC含量格式化输出

gccount_read.pl

#!/usr/bin/perl
use strict;
use warnings;
my $read = "AACAAACCCCTTTTCTCTATTAAAAAATACAAAATAGCTTAGCTGCGGCATAGTGGAGCACG";
my $G = ($read =~ s/G/G/g);#s为匹配表示匹配,($read =~ s/G/G/g)返回匹配次数个数
my $C = ($read =~ s/C/C/g);
my $A = ($read =~ s/A/A/g);
my $T = ($read =~ s/T/T/g);
my $total = $G + $C + $A +$T;
my $GC_P = ($G + $C)/$total;#格式化输出printf
printf "\$G:%d\n\$C:%d\n\$A:%d\n\$T:%d\n\$total:%d\n\$GC_P:%.2f\n",($G,$C,$A,$T,$total,$GC_P);#sprintf格式化不输出
my $result = sprintf("\$G:%d\n\$C:%d\n\$A:%d\n\$T:%d\n\$total:%d\n\$GC_P:%.2f\n",($G,$C,$A,$T,$total,$GC_P));#输出到文件GC.file
open OUT,">",&#

Perl练习01之计算一条read的ATGC个数GC含量格式化输出相关推荐

下载FATSQ，读取10条序列并计算每条序列的长度和GC含量
FASTQ 这是目前存储测序数据最普遍.最公认的一个数据格式,另一个是uBam格式,但这篇文章中不打算对其进行介绍.上面所讲的FASTA文件,它所存的都是已经排列好的序列(如参考序列),FASTQ存的 ...
perl语言学习第一篇：将文本perl.txt内容读入，再格式化输出（菜鸟小试牛刀）
实习较闲,自学perl,刚看完菜鸟教程,小试牛刀下面的程序主要是读入perl.txt文件 cat perl.txt 4 karsa 1 uzi 3 xiaohu 2 ming 5 langx 然后格 ...
JAVA计算两条直线的交点（判断两条线段相交）
//计算两条直线的交点 //输入四个点,前两个点确定线段1,后两个点确定线段2 //输出交点坐标 public class Main{public static void main(String ar ...
OpenCV计算两条平行线之间的距离
代码来自www.opencvchina.com #include "cv.h" #include "highgui.h" #include "cxco ...
使用C++面向对象思想计算两条直线交点
使用C++面向对象思想计算两条直线交点以下是使用C++面向对象思想计算两条直线交点的示例代码: #include <iostream>using namespace std;class ...
两条直线的交点 c语言,计算两条直线的交点(C#)
PS:从其他地方看到的源码是有问题的.下面是修正后的 /// /// 计算两条直线的交点 /// /// L1的点1坐标 /// L1的点2坐标 /// L2的点1坐标 /// L2的点2坐标 /// ...
第2关：计算二叉树的深度和节点个数
#include <iostream> #include "binary_tree.h"BiTreeNode* CreatBiTree(char* s, int &am ...
计算最长公共数字串个数
一．华为OJ题目(计算最长公共数字个数) 样例输出输出123058789,函数返回值9 输出54761,函数返回值5 接口说明函数原型: unsignedint Continumax(char** ...
计算n!中结尾零的个数
计算n!中结尾零的个数(C++实现) . [cpp] view plaincopy #include /*计算n!结尾零的个数,返回零的个数.*/ int CalZeroNum(int n) { in ...
python——变量的类型、不同类型变量的计算、变量的输入以及格式化输出
变量的类型在 python 中定义变量是不需要指定类型(在其他很多高级语言中都需要) 数据类型可以分为数字型.非数字型数字型:整型(int).浮点型(float). 布尔型(bool): 真 Tr ...

Perl练习01之计算一条read的ATGC个数GC含量格式化输出

Perl练习01之计算一条read的ATGC个数GC含量格式化输出相关推荐

最新文章

热门文章