sh 提取所有FASTQ标头

Posted 2021-05-06

IFS=$'\n'
for line in `cat fastq_list.tsv`; do
    zcat $line | head -n 1 | awk -v fname=$line -v OFS="\t" '{ print $0, fname}' >> fastq_headers.tsv
done