khuch123
commented
5 years ago awk '/^>/ {if (seqlen){print "0\t"seqlen-1}; gsub(/^>/,"",$1);printf("%s\t",$1) ;seqlen=0;next; } { seqlen += length($0)}END{print "0\t"seqlen-1}' ${cazy} > ${cazy}.bed
Open khuch123 opened 5 years ago
khuch123
commented
5 years ago awk '/^>/ {if (seqlen){print "0\t"seqlen-1}; gsub(/^>/,"",$1);printf("%s\t",$1) ;seqlen=0;next; } { seqlen += length($0)}END{print "0\t"seqlen-1}' ${cazy} > ${cazy}.bed
Dear Sir/Madam,
I was following your pipeline on my own database but it has been a day that the awk command mentioned in the script is not working.
Please help at the earliest