您好,欢迎来到伴沃教育。
搜索
您的当前位置:首页XiaLab的Read2Counts.pl

XiaLab的Read2Counts.pl

来源:伴沃教育

写在前面

写了一个脚本,整理了一些命令。方便课题组所有人,具体功能是:

从Fastq格式的数据直接到Counts数据

下面介绍用法,

输入数据

1.样本数据

ls -1  S*.fastq
SRR3080055.sra_1.fastq
SRR3080055.sra_2.fastq
SRR3080056.sra_1.fastq
SRR3080056.sra_2.fastq
SRR3080057.sra_1.fastq
SRR3080057.sra_2.fastq
SRR3080058.sra_1.fastq
SRR3080058.sra_2.fastq
SRR3080059.sra_1.fastq
SRR3080059.sra_2.fastq

整理成一个表格

ls S*.fastq|paste - - > sample.list.txt
cat sample.list.txt

表格内容如下

SRR3080055.sra_1.fastq  SRR3080055.sra_2.fastq
SRR3080056.sra_1.fastq  SRR3080056.sra_2.fastq
SRR3080057.sra_1.fastq  SRR3080057.sra_2.fastq
SRR3080058.sra_1.fastq  SRR3080058.sra_2.fastq
SRR3080059.sra_1.fastq  SRR3080059.sra_2.fastq

2.参考基因组

A_subgenome.fa

3.基因结构注释信息

A_subgenome.gff3

运行流程

perl /tools/Read2Counts.pl A_subgenome.fa sample.list A_subgenome.gff3

然后就等着就可以了

Copyright © 2019- bangwoyixia.com 版权所有 湘ICP备2023022004号-2

违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务