3：数据处理:sra转成fq文件

生信技能树已经很贴心地把sra格式数据全部下载。共5.3个T。存放在

代码语言：javascript复制

/data/fudan_TNBC/

查看共多少样本

代码语言：javascript复制

(base) pc@pc-System-Product-Name:/data/fudan_TNBC$ ls -l |grep "^-"|wc -l
727

把/data/fudan_TNBC/下的sra文件转成fastq文件，并存放在/project/raw_fq/下

代码语言：javascript复制

cd /project/raw_fq/
for id in `seq 8223 8454`; do    nohup sudo  fastq-dump --gzip --split-3 /data/fudan_TNBC/SRR851${id}.sra . & done

代码语言：javascript复制

nohup for id in `seq 854 999`; do  sudo fastq-dump --gzip --split-3 /data/fudan_TNBC/SRR8517{id}.sra -O .; done &

代码语言：javascript复制

for ((i=854;i<=999;i  ));do sudo fastq-dump --gzip --split-3 -A SRR35899$i.sra -O .;done &

raw_fq共5.3T，分批进行处理。 实际是，因为机房停电，很多文件没有转换完整。然后最后直接分两批转完了，共占用5.8T空间。

这一步前后用了一个周，除去开始转换太多爆内存和后期机房停电，实际转换大概需要60小时即可。

data project

0 人点赞