提取Rice RNA-seq Database（南方科大）数据

1. 数据来源于这个网站，下载水稻所有RNA-seq数据

image.png image.png

2. 准备自己的基因IDs文件, 注意基因是一行，而且空格隔开。

image.png

3. 利用脚本提取所有目标基因的所有样本的表达谱数据

perl extract_SRA.pl ALl_clusterI_genes.txt   ../epigenetic_regulators/ALL_PFKM.txt >ALL_clusterI_FPKM.txt

4. 处理提取的数据，把每一列的列名之间从\t换成空格。

image.png

5. 选择自己感兴趣的处理组(来自SRAIDs.xls)，可以选择一个基因在网站进行搜索 下载所有SRA号和对应处理，以方便后续选择。

image.png

6.准备SRA IDs，每个SRA ID一行

image.png

7.提取，注意目前Perl脚本有点bug最后一个SRA提取不到，需要写两遍（仅仅最后一个写两遍）。

perl Ectrac_SRA.pl  Title.txt  SRAIDs.txt   ALL_clusterI_FPKM.txt  >Submerge_stress_FPKM.txt