首先，要用到數(shù)組，請(qǐng)先看關(guān)于數(shù)組的基礎(chǔ)知識(shí)
要求：把電腦中的
SRR開頭的文件名改為容易識(shí)別的Library_Name
電腦中的SRR文件如下：
SRR文件
Run和Library_Name的對(duì)應(yīng)關(guān)系如下
也就是要以第二列代替第三列，所以簡(jiǎn)單的rename命令不行，因?yàn)檫@些名字之間?簡(jiǎn)單的替換
SAMN09837892 Lib_FUSCCTNBC158.TT_WES SRR7696207 FUSCCTNBC158_FrozenPrimaryTumorTissue
SAMN09838014 Lib_FUSCCTNBC337.TT_WES SRR8517853 FUSCCTNBC337_FrozenPrimaryTumorTissue
SAMN09837870 Lib_FUSCCTNBC123.TT_WES SRR8517854 FUSCCTNBC123_FrozenPrimaryTumorTissue
SAMN09837938 Lib_FUSCCTNBC228.TT_WES SRR8517855 FUSCCTNBC228_FrozenPrimaryTumorTissue
SAMN09837940 Lib_FUSCCTNBC230.TT_WES SRR8517856 FUSCCTNBC230_FrozenPrimaryTumorTissue
SAMN09837941 Lib_FUSCCTNBC233.TT_WES SRR8517857 FUSCCTNBC233_FrozenPrimaryTumorTissue
SAMN09837942 Lib_FUSCCTNBC234.TT_WES SRR8517858 FUSCCTNBC234_FrozenPrimaryTumorTissue
SAMN09837934 Lib_FUSCCTNBC220.TT_WES SRR8517859 FUSCCTNBC220_FrozenPrimaryTumorTissue
SAMN09837935 Lib_FUSCCTNBC223.TT_WES SRR8517860 FUSCCTNBC223_FrozenPrimaryTumorTissue
SAMN09837936 Lib_FUSCCTNBC225.TT_WES SRR8517861 FUSCCTNBC225_FrozenPrimaryTumorTissue
SAMN09837937 Lib_FUSCCTNBC227.TT_WES SRR8517862 FUSCCTNBC227_FrozenPrimaryTumorTissue
SAMN09837943 Lib_FUSCCTNBC235.TT_WES SRR8517863 FUSCCTNBC235_FrozenPrimaryTumorTissue
SAMN09837944 Lib_FUSCCTNBC236.TT_WES SRR8517864 FUSCCTNBC236_FrozenPrimaryTumorTissue
SAMN09838036 Lib_FUSCCTNBC387.TT_WES SRR8517865 FUSCCTNBC387_FrozenPrimaryTumorTissue
SAMN09838037 Lib_FUSCCTNBC389.TT_WES SRR8517866 FUSCCTNBC389_FrozenPrimaryTumorTissue
SAMN09838032 Lib_FUSCCTNBC380.TT_WES SRR8517867 FUSCCTNBC380_FrozenPrimaryTumorTissue
SAMN09838033 Lib_FUSCCTNBC383.TT_WES SRR8517868 FUSCCTNBC383_FrozenPrimaryTumorTissue
SAMN09838077 Lib_FUSCCTNBC440.TT_WES SRR8517869 FUSCCTNBC440_FrozenPrimaryTumorTissue
SAMN09838075 Lib_FUSCCTNBC438.TT_WES SRR8517870 FUSCCTNBC438_FrozenPrimaryTumorTissue
SAMN09837967 Lib_FUSCCTNBC262.TT_WES SRR8517871 FUSCCTNBC262_FrozenPrimaryTumorTissue
SAMN09837968 Lib_FUSCCTNBC263.TT_WES SRR8517872 FUSCCTNBC263_FrozenPrimaryTumorTissue
SAMN09837965 Lib_FUSCCTNBC260.TT_WES SRR8517873 FUSCCTNBC260_FrozenPrimaryTumorTissue
提示
- 1根據(jù)上面的對(duì)應(yīng)關(guān)系，首先創(chuàng)建( touch)這些gz文件出來
- 2 然后 shell腳本改名
比如SRR7696207改為L(zhǎng)ib_FUSCCTNBC158.TT_WES
1 創(chuàng)建gz文件
先寫入上述第二個(gè)表的文件
cat >tmp
SAMN09837892 Lib_FUSCCTNBC158.TT_WES SRR7696207 FUSCCTNBC158_FrozenPrimaryTumorTissue
SAMN09838014 Lib_FUSCCTNBC337.TT_WES SRR8517853 FUSCCTNBC337_FrozenPrimaryTumorTissue
SAMN09837870 Lib_FUSCCTNBC123.TT_WES SRR8517854 FUSCCTNBC123_FrozenPrimaryTumorTissue
SAMN09837938 Lib_FUSCCTNBC228.TT_WES SRR8517855 FUSCCTNBC228_FrozenPrimaryTumorTissue
SAMN09837940 Lib_FUSCCTNBC230.TT_WES SRR8517856 FUSCCTNBC230_FrozenPrimaryTumorTissue
.....
創(chuàng)建1.fastq.gz和2.fastq.gz文件
cat tmp|cut -f3|while read id;do touch ${id}_{1,2}.fastq.gz;done
或者
cat config|awk '{print$3}'|while read id;do touch ${id}_{1,2}.fastq.gz;done
ls
SRR6269851_1.fastq.gz SRR6269859_1.fastq.gz SRR6269867_1.fastq.gz SRR6269875_1.fastq.gz SRR8517854_1.fastq.gz SRR8517862_1.fastq.gz SRR8517870_1.fastq.gz
SRR6269851_2.fastq.gz SRR6269859_2.fastq.gz SRR6269867_2.fastq.gz SRR6269875_2.fastq.gz SRR8517854_2.fastq.gz SRR8517862_2.fastq.gz SRR8517870_2.fastq.gz
SRR6269852_1.fastq.gz SRR6269860_1.fastq.gz SRR6269868_1.fastq.gz SRR6269876_1.fastq.gz SRR8517855_1.fastq.gz SRR8517863_1.fastq.gz SRR8517871_1.fastq.gz
SRR6269852_2.fastq.gz SRR6269860_2.fastq.gz SRR6269868_2.fastq.gz SRR6269876_2.fastq.gz
2重命名
cat tmp|while read id
do
arr=($id)
Lib_name=${arr[1]}
SRR_name=${arr[2]}
echo $Lib_name echo $SRR_name
mv ${SRR_name}_1.fastq.gz ${Lib_name}_1.fastq.gz
mv ${SRR_name}_2.fastq.gz ${Lib_name}_2.fastq.gz
done
解釋如下
- 1 arr=($id)
把read到的id也就是tmp的所有數(shù)據(jù)放在了名為arr數(shù)組里，()的作用就是建立數(shù)組- 2 arr[1]讀取arr第2列，命名為L(zhǎng)ib_name
arr[12]讀取arr第3列，命名為SRR_name
課程分享
