The following entry did not occur in all.csv: NC_008463.1 RefSeq sequence_feature 759485 759542 . + . ID=id13 Note=23S ribosomal RNA rRNA prediction is too short gbkey=misc_feature # pipeline: Rscript diff-exp-deseq4.r for i in de*csv do awk -f get-padj-le0.05.awk $i >> $i"2" done