linux合并fa文件_与你分享生信好用的单行命令

刘小泽记录于 2019.5.10

将Turner写的一些好用的单行命令与大家分享，原文还有许多可以去看

https://github.com/stephenturner/oneliners

About Fastq/fasta

fastq sequences length distribution => 得到fq文件中序列长度的分布

$ zcat file.fastq.gz | awk 'NR%4 == 2 {lengths[length($0)]++} END {for (l in lengths) {print l, lengths[l]}}'

reverse complement => 反向互补

$ echo 'ATTGCTATGCTNNNT' | rev | tr 'ACTG' 'TGAC'

fastq2fasta

$ zcat file.fastq.gz | paste - - - - | perl -ane 'print ">$F[0]\n$F[1]\n";' | gzip -c > file.fasta.gz

split a multifasta file into single ones with csplit => fasta按>拆分

# * refers to the number of files 可以选择拆分的文件数量

$ csplit -z -q -n 4 -f test test.fa /\>/ {*}

# OR use awk 一次性全部按>拆分

$ awk '/^>/{s=++d".fa"} {print > s}' multi.fa

single line fasta to multi-line of 50 characters in each line => 单行fa变多行

$ awk -v FS= '/^>/{print;next}{for (i=0;i<=NF/50;i++) {for (j=1;j<=50;j++) printf "%s", $(i*50 +j); print ""}}' file

# fold -w 50 file

multi-line fasta to one-line => 一个多行fa文件变单行

# 方法一：

$ awk 'BEGIN{RS=">"}NR&