比對使用STAR

計算counts使用featureCounts

轉(zhuǎn)化counts使用腳本 readsToTPM.pl

$usage = "perl feartureCounts featureCounts.summary";die? "\t\"the usage is? $usage \"\t " if @ARGV ==0;open FILE,$ARGV[0] or die $!;open FILE2,$ARGV[1] or die $1;open OUT1,">TPM.txt";open OUT2,">reads.txt";print OUT1 "geneName\tTPM\n";print OUT2 "geneName\treadsCount\n";while(){chomp;if ($_ =~ /Assigned\t(\d+)/){$totalReads = $1;last;}}while(){

chomp;

next if $_ =~ /[Program][Geneid]/;

@a = split;

printf OUT1 $a[0]."\t"."%.2f",$a[6]/$a[5]*1000/$totalReads*1000000;

print OUT1 "\n";

print OUT2 "$a[0]\t$a[6]\n";

}