标签：arr gff 对照 start length key cds gene my

 1 #!/usr/bin/perl
 2 
 3 use strict;
 4 use warnings;
 5 
 6 my $gff = $ARGV[0];my %cut = &cut($gff);
 7 
 8 my $gene_number = $ARGV[1];my $gene_id = $ARGV[2];my $gene_name = $ARGV[3];my $start = $ARGV[4];my $end = $ARGV[5];
 9 
10 foreach my $key_1(keys %cut)
11 {
12     my $hash_2 = $cut{$key_1};
13 
14     foreach my $key_2(keys %$hash_2)
15     {
16         my @arr = $hash_2 -> {$key_2};
17         
18 #        print "$key_1,$key_2,$arr[0][0],$arr[0][1]\n";
19         
20         if($key_1 eq $gene_name)
21         {
22 
23             
24             if($arr[0][0] <= $start && $arr[0][1] >= $start )
25             {
26                 if($start+30>=$arr[0][1])
27                 {
28                     print "$gene_number\t$gene_id\t$gene_name\t$start".'-'."$end\tstart:exon$key_2\n";
29                 }
30             }
31             if($arr[0][0] <= $end && $arr[0][1] >= $end )
32             {
33                 if($end-30>=$arr[0][0])
34                 {
35                     print "$gene_number\t$gene_id\t$gene_name\t$start".'-'."$end\tend:exon$key_2\n";
36                 }
37             }
38         }
39     }
40 }
41 
42 
43 #############sub################
44 
45 sub cut
46 {
47     my %gene;
48 
49     my $exon = 1;my $start = 0;my $length = 0;
50 
51     my $gff = shift;open GFF,"$gff"; 
52     
53     while(my $line = <GFF>)
54     {
55         chomp $line;
56 
57         my @q = split /\s/,$line;
58 
59         if($q[2] =~ /mRNA/)
60         {
61             $exon = 1;$start = $q[3];$length = 0;
62         }
63         elsif($q[2] =~ /CDS/)
64         {
65             $q[8] =~ /Parent=(.*);S/m;my $key = $1;  
66 
67             my $new_length = $q[4]-$q[3]; 
68 
69             $exon =~/(\w*)/; $gene{$key}{$1} = [$length+1,$length+$new_length];
70 
71             $exon++; $length += $new_length;
72         }    
73     }   
74     return %gene;
75 }

标签：arr,gff,对照,start,length,key,cds,gene,my
来源： https://www.cnblogs.com/yuanjingnan/p/11087860.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。