博文
使用awk实现：二代测序文件fastq转换为fasta格式
|||
方法1
awk '{if(NR%4 == 1){print ">" substr($0, 2)}}{if(NR%4 == 2){print}}' fastq > fasta
方法2：
awk 'BEGIN{P=1}{if(P==1||P==2){gsub(/^[@]/,">");print}; if(P==4)P=0; P++}' fastq > fasta
推荐使用方法1,方法2可能会有缺陷，有些数据的值也是@，如果恰好是第一个碱基，就会出错
https://blog.sciencenet.cn/blog-803390-819177.html
上一篇：perl printf打印数字时，如果位数不足，前面用0补齐
下一篇：fasta格式转化-fasta_formatter
