Skip to content

Commit

Permalink
s
Browse files Browse the repository at this point in the history
  • Loading branch information
lixin4306ren committed Feb 10, 2017
1 parent f242bb2 commit 8722153
Show file tree
Hide file tree
Showing 2 changed files with 7 additions and 27 deletions.
32 changes: 6 additions & 26 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -3,45 +3,25 @@
此流程依赖于array suite,安装array suite参考此[网页](http://www.arrayserver.com/wiki/index.php?title=Oshell#Overview)

## 主要参数
主程序为`RNA_seq_align_all_module.pl``-h`显示所有参数如下
主程序为`Run_OSA_all_modules.pl`主要参数用法
```
RNA-Seq pipeline
Usage: RNA_seq_align_all_module.pl <options>
-sample_list sample list
-out_dir output result folder
-ref_dir reference OSA base path
-ref reference name
-gene_model gene_model name
-gtf gene gtf file
-gtf_exon exon gff file for DEXSeq
-step align or merge_rmdup or count or count_exon or tdf
-h or -help Show Help , have a choice
perl Run_OSA_all_modules.pl sample.list [output_dir] [database_dir] [genome_version] [gene_model_version]
```


`sample.list`为所有样品目录路径,目录下为fastq.gz文件,格式如下:
```
/home/others/xli/raw_data/2016_08_01_ZY/RNA-Seq/S1
/home/others/xli/raw_data/2016_08_01_ZY/RNA-Seq/S2
```
`-ref`默认基因组版本Human.B37.3
`-gene_model`默认版本为RefGene
`genome_version`默认基因组版本Human.B37.3
`gene_model_version`默认版本为RefGene
基因组和基因注释版本对应信息参考此[网页](http://www.arrayserver.com/wiki/index.php?title=A_list_of_compiled_genome_and_gene_model_from_OmicSoft)


## alignment
如何是human
如果是human
```
perl ~/scripts/rnaseq_pipeline/Run_OSA_all_modules.pl sample.list
```
如果是其他物种,需指明所用genome版本和gene model版本,例如mouse
```
Expand Down
2 changes: 1 addition & 1 deletion Run_OSA_all_modules.pl
Original file line number Diff line number Diff line change
Expand Up @@ -28,7 +28,7 @@
closedir DIR;


print "mono /home/others/xli/soft/oshell/16-8-11/oshell.exe --runscript $ref_dir $pwd/$sample_name.oscript /tmp/ /opt/apps/mono/bin/mono\n";
print "mono /home/others/xli/soft/oshell/16-8-11/oshell.exe --runscript $ref_dir $pwd/$sample_name.oscript ./ /opt/apps/mono/bin/mono\n";

open O,">$pwd/$sample_name.oscript"||die;

Expand Down

0 comments on commit 8722153

Please sign in to comment.