Commit c59ad42e authored by alain viari's avatar alain viari

corrected Annotation reordering for CDS

parent 30df24d0
......@@ -55,17 +55,17 @@ pushTmpDir ORG.organnot
loginfo "Done."
loginfo "Ordering annotations..."
awk '/^.....[^ ]/ { \
awk '/^.....(misc|repeat|rRNA|tRNA|gene)/ { \
match($3,"[0-9][0-9]*"); \
pos=substr($3,RSTART,RLENGTH)*1000 + 1; \
print pos,$0} \
/^..... / { \
pos++; \
print pos,$0; \
next} \
{ pos++; \
print pos,$0}' "${RESULTS}.annot" | \
sort -nk1 | \
sort -nk1 |\
awk '{ \
match($0,"^[0-9]* ");\
line=substr($0,RLENGTH+1,100000);\
line=substr($0,RLENGTH+1);\
print line}'
loginfo "Done."
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment