07 2012 档案

[笔记]reads_classifier
摘要：reads_classifier.pl 1 #!/usr/bin/perl 2 #=============================================================================== 3 # 4 # FILE: reads_classifier.pl 5 # 6 # USAGE: ./reads_classifier.pl 7 # 8 # DESCRIPTION: 9 # 10 # OPTIONS: --- 11 # REQUIREMENTS: --- 12 # ... 阅读全文
posted @ 2012-07-19 19:21 Puriney 阅读(199) 评论(0) 推荐(0) 编辑
[笔记]sam2bed
摘要：sam2bed.pl 1 #!/usr/bin/perl -w 2 # 3 ######################################## 4 # 5 # File Name: 6 # sam2bed.pl 7 # 8 # Description: 9 # 10 # 11 # Usage:12 # 13 # 14 # Author:15 # Xi Wang, wang-xi05@mails.thu.edu.cn16 # 17 # Date:18 # Mon Nov 2 23:08:26 CST 200919 #20 ##################... 阅读全文
posted @ 2012-07-16 18:43 Puriney 阅读(502) 评论(0) 推荐(0) 编辑
[笔记] sam文件挑取uniq reads
摘要：本来呢，samtools本来就可以一句话就能提取，只是没得网络没得下载samtools，只能保留自传脚本。输入文件为raw.sam输出为提取uniq reads后的uniq.sam（仍然保留着sam格式）终端里会打印出简单的报告 1 use strict; 2 use warnings; 3 4 my ($in, $out) = @ARGV; 5 die <DATA> unless (@ARGV == 2); 6 open IN, "$in"; 7 open OUT, ">$out"; 8 9 my ($sum, $uniq, $xnd 阅读全文
posted @ 2012-07-16 16:57 Puriney 阅读(741) 评论(0) 推荐(0) 编辑