基礎(chǔ)命令:
cd 空格 #進(jìn)入家目錄
cd ~ #進(jìn)入家目錄
cd - #與上次所在目錄進(jìn)行切換
pwd #查看當(dāng)前所在路徑
tree #以樹狀結(jié)構(gòu)展現(xiàn)文件層級
history #查看所有運(yùn)行命令的歷史記錄
cat/more/less #查看文件內(nèi)容
ls #查看當(dāng)前目錄的文件
ls -lh #查看當(dāng)前文件的詳細(xì)信息
alias ll='ls -lh' #建立命令快捷方式
echo #打印字符串
rm #刪除文件
mv #移動文件或改名(剪切)
cp #拷貝
1、在任意文件夾下面創(chuàng)建形如 1/2/3/4/5/6/7/8/9 格式的文件夾。
>mkdir -p 1/2/3/4/5/6/7/8/9 #創(chuàng)建多級目錄要加參數(shù)-p
>tree 1
1
└── 2
└── 3
└── 4
└── 5
└── 6
└── 7
└── 8
└── 9
8 directories, 0 files
2、在創(chuàng)建好的文件夾下面,比如我的是 ~/1/2/3/4/5/6/7/8/9 ,里面創(chuàng)建文本文件 me.txt
>cd ~/1/2/3/4/5/6/7/8/9 #進(jìn)入文件夾
>touch me.txt #創(chuàng)建文件
>ls #查看目錄文件
me.txt
3、在文本文件 me.txt 里面輸入內(nèi)容:
I love bioinfomatics.
And you ?
>vim me.txt #esc : wq 保存
>cat me.txt #查看文件內(nèi)容
I love bioinfomatics.
And you ?
4、刪除上面創(chuàng)建的文件夾 1/2/3/4/5/6/7/8/9 及文本文件 me.txt
>cd ~
>rm -rf 1 #-rf表示逐級刪除目錄,刪除文件不需要
5、在任意文件夾下面創(chuàng)建 folder1~5這5個文件夾,然后每個文件夾下面繼續(xù)創(chuàng)建 folder1~5這5個文件夾
>for i in `seq 5`;do mkdir folder_$i;for n in `seq 5`;do mkdir folder_$i/folder_$n;done;done
>ls */
folder_1/:
folder_1 folder_2 folder_3 folder_4 folder_5
folder_2/:
folder_1 folder_2 folder_3 folder_4 folder_5
folder_3/:
folder_1 folder_2 folder_3 folder_4 folder_5
folder_4/:
folder_1 folder_2 folder_3 folder_4 folder_5
folder_5/:
folder_1 folder_2 folder_3 folder_4 folder_5
6、在第五題創(chuàng)建的每一個文件夾下面都創(chuàng)建第二題文本文件 me.txt ,內(nèi)容也要一樣。
>for i in `seq 5`;do for n in `seq 5`;do echo 'I love bioinfomatics.And you ?' > ~/folder_$i/folder_$n/me.txt;done;done
>tree folder_*
folder_1
├── folder_1
│ └── me.txt
├── folder_2
│ └── me.txt
├── folder_3
│ └── me.txt
├── folder_4
│ └── me.txt
└── folder_5
└── me.txt
folder_2
├── folder_1
│ └── me.txt
├── folder_2
│ └── me.txt
├── folder_3
│ └── me.txt
├── folder_4
│ └── me.txt
└── folder_5
└── me.txt
folder_3
├── folder_1
│ └── me.txt
├── folder_2
│ └── me.txt
├── folder_3
│ └── me.txt
├── folder_4
│ └── me.txt
└── folder_5
└── me.txt
folder_4
├── folder_1
│ └── me.txt
├── folder_2
│ └── me.txt
├── folder_3
│ └── me.txt
├── folder_4
│ └── me.txt
└── folder_5
└── me.txt
folder_5
├── folder_1
│ └── me.txt
├── folder_2
│ └── me.txt
├── folder_3
│ └── me.txt
├── folder_4
│ └── me.txt
└── folder_5
└── me.txt
25 directories, 25 files
7、再次刪除掉前面幾個步驟建立的文件夾及文件
>rm -rf folder_*
8、下載 http://www.biotrainee.com/jmzeng/igv/test.bed 文件,后在里面選擇含有 H3K4me3 的那一行是第幾行,該文件總共有幾行。
>wget http://www.biotrainee.com/jmzeng/igv/test.bed #下載
>grep -n 'H3K4me3' test.bed #-n顯示行號,如下所示有 H3K4me3 的那一行是第8行
8:chr1 9810 10438 ID=SRX387603;Name=H3K4me3%20(@%20HMLE);Title=GSM1280527:%20HMLE%20Twist3D%20H3K4me3%20rep2%3B%20Homo%20sapiens%3B%20ChIP-Seq;Cell%20group=Breast;<br>source_name=HMLE_Twist3D_H3K4me3;cell%20type=human%20mammary%20epithelial%20cells;transfected%20with=Twist1;culture%20type=sphere;chip%20antibody=H3K4me3;chip%20antibody%20vendor=Millipore; 222 . 9810 10438 0,226,255
>wc -l test.bed #查看文件行數(shù),該文件共10行
10 test.bed
9、下載 http://www.biotrainee.com/jmzeng/rmDuplicate.zip 文件,并且解壓,查看里面的文件夾結(jié)構(gòu)
>wget http://www.biotrainee.com/jmzeng/rmDuplicate.zip
>unzip rmDuplicate.zip
>tree rmDuplicate
rmDuplicate
├── picard
│ ├── paired
│ │ ├── readme.txt
│ │ ├── tmp.header
│ │ ├── tmp.MarkDuplicates.log
│ │ ├── tmp.metrics
│ │ ├── tmp.rmdup.bai
│ │ ├── tmp.rmdup.bam
│ │ ├── tmp.sam
│ │ └── tmp.sorted.bam
│ └── single
│ ├── readme.txt
│ ├── tmp.header
│ ├── tmp.MarkDuplicates.log
│ ├── tmp.metrics
│ ├── tmp.rmdup.bai
│ ├── tmp.rmdup.bam
│ ├── tmp.sam
│ └── tmp.sorted.bam
└── samtools
├── paired
│ ├── readme.txt
│ ├── tmp.header
│ ├── tmp.rmdup.bam
│ ├── tmp.rmdup.vcf.gz
│ ├── tmp.sam
│ ├── tmp.sorted.bam
│ └── tmp.sorted.vcf.gz
└── single
├── readme.txt
├── tmp.header
├── tmp.rmdup.bam
├── tmp.rmdup.vcf.gz
├── tmp.sam
├── tmp.sorted.bam
└── tmp.sorted.vcf.gz
6 directories, 30 files
10、打開第九題解壓的文件,進(jìn)入 rmDuplicate/samtools/single 文件夾里面,查看后綴為 .sam 的文件,搞清楚生物信息學(xué)里面的SAM/BAM 定義是什么。
>cd rmDuplicate/samtools/single
>ls
readme.txt tmp.header tmp.rmdup.bam tmp.rmdup.vcf.gz tmp.sam tmp.sorted.bam tmp.sorted.vcf.gz
>less -S tmp.sam #.bam文件是二進(jìn)制問件,查看要安裝Samtools,使用samtools view查看
SRR1042600.42157053 0 chr1 629895 42 51M * 0 0 ATAACCAATACTACCAATCANTACTCATCATTAATAATCATAATGGCTA
SRR1042600.42212881 0 chr1 629895 42 51M * 0 0 ATAACCAATACTACCAATCANTACTCATCATTAATAATCATAATGGCTA
SRR1042600.12010763 16 chr1 629895 24 51M * 0 0 ATAACCAATACTTCTAATCAAAACTCATCATTAATAATCATAATGGCTA
>