palomo11
palomo11
I have the same question. @jacodela Did you figure out how to do it?
Hi Stephen, Thanks for the advice. That's what I've done and it is working so far.
Hi, Not sure if MIDAS now support *gff, otherwise you can get the .genes files as follows (I guess there must be simplest way of doing it, but this one...
I think there should be a space after -v: awk -v OFS='\t' '{$5 = "CDS"; print}' ${file%.faa}_OK3.txt > ${file%.faa}_OK3_OK.txt awk -v OFS='\t' 'NR==FNR {h[$1] = $2; next} {print $1,h[$1],$2,$3,$4,$5}' ${file%.faa}_OK4.txt...