awk '{if($0 ~ />/) {print $0} else {printf("%s", $0)}}' [fasta file] | sed 's/>/\n>/g' | sed '$ s/$/\n/' | sed 1d | awk '{if($0 ~ /^>/) {printf("%s\t", $0)} else {print length($0)}}'
统计基因组fasta文件中每条scaffold的碱基数量
最新推荐文章于 2024-09-30 08:15:30 发布