grep "NZ_CP027786.1" ~/term1/genome/GCF_003795145.1_ASM379514v1_feature_table.txt | cut -f 1,2,8,9,10 > table grep -v "^#" sd_direct.txt > res_d sed 's/^ *\([^ ].*\)/\1/' res_d > fres_d sed -E 's/ +/\t/g' fres_d | sed 's/\(.*\)\t\(.*\)\t\(.*\)\t\(.*\)\t\(.*\)\t\(.*\)/\4\t\6\t\1\t\2\t\3/' > nres_d grep -v "^#" sd_complement.txt > res_c sed 's/^ *\([^ ].*\)/\1/' res_c > fres_c sed -E 's/ +/\t/g' fres_c | sed 's/\(.*\)\t\(.*\)\t\(.*\)\t\(.*\)\t\(.*\)\t\(.*\)/\4\t\6\t\1\t\2\t\3/' > nres_c cat nres_c nres_d table | sort -k3n > table_with_sd.txt