View Single Post
Old 09-22-2014, 02:07 PM   #7
kbradnam
Member
 
Location: Davis, CA

Join Date: May 2011
Posts: 53
Default

I don't see this:

$ cut -f 1 -d " " completeness_cutoff.tbl > 248_CEG_ids.txt
$ grep ">" kogs.fa | sed 's/.*___//' | sort -u > 458_CEG_ids.txt
$ grep -f 248_CEGs_ids.txt 458_CEG_ids.txt | wc -l

The last step reports a count of 248. I.e. all CEG IDs present in completeness_cutoff.tbl are present in logs.fa.
kbradnam is offline   Reply With Quote