scripts/codon_distribution_scer.pl at master · rhilt13/scripts · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
#! /usr/bin/perl

use Bio::SeqIO;
use Data::Dumper;

open(IN,$ARGV[0]);		# gly sites list file
$first_line=<IN>;
if ($_=~/^#/){
	$first_line=<IN>;
}

@a=split(/\t/,$first_line);
$id=$a[0];
$num=$a[1]*3 - 3;
#$num=$a[3]*3 - 3;
push @pos,$num;
while(<IN>){
	chomp $_;
	# next if ($_=~/^#/);
	@a=split(/\t/,$_);
	$no_of_sites{$a[0]}=$a[4];
	if ($a[0] eq $id){
		$num=$a[1]*3 - 3;
		#$num=$a[3]*3 - 3;
		push @pos, $num;
	}else{
		$id_hash{$id}= [@pos];
		$id=$a[0];
		@pos=();
		$num=$a[1]*3 - 3;
		#$num=$a[3]*3 - 3;
		push @pos,$num;
	}
}
# print Dumper(\%id_hash);
# print Dumper(\%no_of_sites);

open(IN2,$ARGV[1]);			# codon scores file
while(<IN2>){
	chomp;
	@b=split(/\t/,$_);
	# @c=split(/\(/,$b[1]);
	# $codon_hash{$b[0]}=$c[0];
	$codon_hash{$b[0]}=$b[1];
}
# print Dumper(\%codon_hash);

# =ccc
$new=Bio::SeqIO->new(-file=>$ARGV[2], -format=>"fasta");
print "+CDS_ID\tsites/protein\tsite\tloc\n";

$z=0;
$y=0;
while($seq=$new->next_seq){
	$id = $seq->id;
	$seq = $seq->seq;
	if (defined $id_hash{$id}){
		# print "aa\n";
		foreach $gly(@{$id_hash{$id}}){
			# print "$gly\n";
			@score=();
			@out=();
			if (($gly+300)<length($seq)){
				$str=substr $seq, $gly, 300;
			}else{
				$str=substr $seq, $gly;
			}
			# print "$gly\t$str\n";
			for ($i=0;$i<length($str);$i+=3){
				$codon=substr $str,$i,3;
				$codon=uc($codon);
				$score[$y][$z]=$codon_hash{$codon};
				# print "$gly\t$codon\n";
				# push(@score,$codon_hash{$codon});
				$y++;
			}
			# foreach $z(@score){
			# 	print "$z\n";
			# }
			# print "-----------------\n";
			# print "$id\t$gly\t@score\n";
			# @out=slide_win(\@score);
			# get_win(\@score,0);
			# $aa=$gly/3 +1;
			# foreach $result(@out){
			# 	print "$id\t$no_of_sites{$id}\t$aa\t$result\n";
			# }
			$z++;
		}
	}
}

foreach $site(0..@score-1){
	foreach $loc(0..@{$score[$site]}-1){
		print ""
	}
}