News forum

Counting genes in PSET 08

Picture of Eric Franzosa
Counting genes in PSET 08
by Eric Franzosa - Wednesday, 30 March 2016, 6:52 PM

Hi Folks,

One of your classmates pointed out that the MEME website (used in PSET 08) now strips the ">" character from FASTA headers in reports. In general this is a helpful feature, but it also breaks one of the suggested tricks for counting genes in MEME output. I have revised the command in the current posting of the homework (dated 30 March 2016). The new command looks like...

grep "^Y" meme.txt | grep "\[" | sed "s/\s.*//g" | sort | uniq | wc -l