bash備忘錄

1.文件內容提取

  1. 去掉表頭,取第一列,存入新文件
    tail -n +2 sig-genes-deseq2.txt | awk '{print $1}' > sig-genes-list.txt
  2. 取第一行,將制表符換成換行符
    head -n1 111.txt | tr '\t' '\n\'
  3. echo "aaa" | xargs -I {} echo {}
    -I參數指定替換字符竄 大括號
  4. 刪除前三行 sed -i.bak '1,3d’ sed.ttt

2.批量操作

for fn in data/DRR0161{25..40};
dosamp=`basename ${fn}`
echo "Processing sample ${samp}"
salmon quant -i athal_index -l A \
       -1 ${fn}/${samp}_1.fastq.gz \
       -2 ${fn}/${samp}_2.fastq.gz \
       -p 8 -o quants/${samp}_quant
done 
for ((i=1;i<=19;i++)) 
for i in `seq 25 40`
for name in *.fastq; do echo $name; done
for id in {14..19}
do 
  mkdir DRR0161${i}; 
  cd DRR0161${i}; 
  wget ftp://ftp.sra.ebi.ac.uk/vol1/fastq/DRR016/DRR0161${i}/DRR0161${i}_1.fastq.gz; 
  wget ftp://ftp.sra.ebi.ac.uk/vol1/fastq/DRR016/DRR0161${i}/DRR0161${i}_2.fastq.gz; 
  cd ..; 
done

3.文件操作

  1. 搜索文件 find / -name httpd.conf
  2. 查看 ls -alth
  3. 查看目錄大小 du -sh
  4. 程序運行
echo 'java -jar ~/src/Trimmomatic-0.32/trimmomatic-0.32.jar $@' >> ~/bin/trimmomatic
chmod +x ~/bin/trimmomatic
  1. 除去不匹配的文件ls | grep -v gz
  2. 顯示目前后臺tophat進程,并kill
ps aux | grep tophat | awk ‘{print $2}’ | xargs kill -9
screen -r tmp_tophat_x_ten

try

  1. 刪除文件第一行:
    sed -i '1d' filename
  2. 刪除文件最后一列:
    awk '{print $NF}' filename
  3. 比較文件的方法:
    1)comm -3 --nocheck-order file1 file2
  1. grep -v -f file1 file2: 輸出file2中有file1中沒有的行
    diff file1 file2
最后編輯于
?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發布,文章內容僅代表作者本人觀點,簡書系信息發布平臺,僅提供信息存儲服務。

推薦閱讀更多精彩內容