8月31日 perl 学习笔记   基因序列练习

任务: 把测序得到结果构成哈希,基因名字为哈希的键值,基因序列为哈希的元素。 然后去除重复的基因。

实战操作

#! /usr/bin/perl

open (OUT,’>/Users/mac/Desktop/liunew3.txt’); #利用文件句柄创建输出文件

open(OUTTTT,’>/Users/mac/Desktop/liunew666666.txt’); #利用文件句柄创建输出文件

open (IN,'< /Users/mac/Desktop/name.txt’)|| die “can’t open”; #利用文件句柄调取文件

open (OUTTT,’>/Users/mac/Desktop/liunew666.txt’);

#先将数据生成一行

my %hash;

$n=2;

while (){

if(/\A>\s*/){

$id{[$n – 1 ]} = “$_”;  #利用哈希

$n++;

}else{

$hash{$id{[$n – 1 ]}} = $_;    #利用哈希

}

print OUT ;

}

my %new_hash = reverse %hash;

print %hash;

#对哈希排序

foreach $key(sort keys %new_hash){

$value = $new_hash{$key};    #本质上是一种输出

print OUTTT $key;

print OUTTT “\n”;

print OUTTT $value;

}

my %new_hash2 = reverse %new_hash;

print OUTTTT %new_hash2;

    原文作者:小郑的学习笔记
    原文地址: https://www.jianshu.com/p/2f7996ae62e7
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞