Untitled
unknown
plain_text
4 years ago
1.2 kB
7
Indexable
for file in ls *.temp
do
echo "Now Processing $file"
# Find lines that begin with numbers
awk '{ if ($1 ~ /^[0-9]/) print $1}' $file > temp_holding.tmp
# if lines found, truncate $file to ID and add entries to gff
if [ ! -s temp_holding.tmp ];
then
echo "No predicted genes found";
else
echo "Predicted Genes being processed";
prefix=$(echo $file | sed 's/_length.*//')
while read line
do
# $line contains the ID number we should parse out of $file
# need to truncate $file
ID=$(echo $file | sed 's/_length.*//')
grep "^$line" $file | awk '{print $0}' | awk -v ID="$ID" \
'{print ID, "glimer", "CDS", $4, $5, $9, $2, ".", "."}' >> final.gff
done < temp_holding.tmp
fi
done
Editor is loading...