從原始的FASTQ生成TPM

比對使用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";

}

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。

推薦閱讀更多精彩內(nèi)容