2.7 Bonus : Finding Things

Teaching: 20 min || Exercises: 10 min

Overview

Questions:
  • How can I find files?
  • How can I find things in files?
Learning Objectives:
  • Use grep to select lines from text files that match simple patterns.
  • Use find to find files and directories whose names match simple patterns.
  • Use the output of one command as the command-line argument(s) to another command.
  • Explain what is meant by ‘text’ and ‘binary’ files, and why many common tools don’t handle the latter well.
Key Points:
  • find finds files with specific properties that match patterns.
  • grep selects lines in files that match patterns.
  • --help is an option supported by many bash commands, and programs that can be run from within Bash, to display more information on how to use these commands or programs.
  • man [command] displays the manual page for a given command.
  • $([command]) inserts a command’s output in place.

2.7.1 grep

In the same way that many of us now use ‘Google’ as a verb meaning ‘to find’, Unix programmers often use the word ‘grep’. ‘grep’ is a contraction of ‘global/regular expression/print’, a common sequence of operations in early Unix text editors. It is also the name of a very useful command-line program.

grep finds and prints lines in files that match a pattern. For our examples, we will use a file that contains three haiku taken from a 1998 competition in Salon magazine. For this set of examples, we are going to be working in the 02_unix_intro/writing subdirectory:

cd
cd Desktop/workshop_files_Bact_Genomics_2023/02_unix_intro/writing
cat haiku.txt
The Tao that is seen
Is not the true Tao, until
You bring fresh toner.

With searching comes loss
and the presence of absence:
"My Thesis" not found.

Yesterday it worked
Today it is not working
Software is like that.
Forever, or Five Years

We haven’t linked to the original haiku because they don’t appear to be on Salon’s site any longer. As Jeff Rothenberg said, ‘Digital information lasts forever — or five years, whichever comes first.’ Luckily, popular content often has backups.

Let’s find lines that contain the word ‘not’:

grep not haiku.txt
Is not the true Tao, until
"My Thesis" not found
Today it is not working

Here, not is the pattern we’re searching for. The grep command searches through the file, looking for matches to the pattern specified. To use it type grep, then the pattern we’re searching for and finally the name of the file (or files) we’re searching in.

The output is the three lines in the file that contain the letters ‘not’.

By default, grep searches for a pattern in a case-sensitive way. In addition, the search pattern we have selected does not have to form a complete word, as we will see in the next example.

Let’s search for the pattern: ‘The’.

grep The haiku.txt
The Tao that is seen
"My Thesis" not found.

This time, two lines that include the letters ‘The’ are outputted, one of which contained our search pattern within a larger word, ‘Thesis’.

To restrict matches to lines containing the word ‘The’ on its own, we can give grep with the -w option. This will limit matches to word boundaries.

Later in this lesson, we will also see how we can change the search behavior of grep with respect to its case sensitivity.

grep -w The haiku.txt
The Tao that is seen

Note that a ‘word boundary’ includes the start and end of a line, so not just letters surrounded by spaces. Sometimes we don’t want to search for a single word, but a phrase. This is also easy to do with grep by putting the phrase in quotes.

grep -w "is not" haiku.txt
Today it is not working

We’ve now seen that you don’t have to have quotes around single words, but it is useful to use quotes when searching for multiple words. It also helps to make it easier to distinguish between the search term or phrase and the file being searched. We will use quotes in the remaining examples.

Another useful option is -n, which numbers the lines that match:

grep -n "it" haiku.txt
5:With searching comes loss
9:Yesterday it worked
10:Today it is not working

Here, we can see that lines 5, 9, and 10 contain the letters ‘it’.

We can combine options (i.e. flags) as we do with other Unix commands. For example, let’s find the lines that contain the word ‘the’. We can combine the option -w to find the lines that contain the word ‘the’ and -n to number the lines that match:

grep -n -w "the" haiku.txt
2:Is not the true Tao, until
6:and the presence of absence:

Now we want to use the option -i to make our search case-insensitive:

grep -n -w -i "the" haiku.txt
1:The Tao that is seen
2:Is not the true Tao, until
6:and the presence of absence:

Now, we want to use the option -v to invert our search, i.e., we want to output the lines that do not contain the word ‘the’.

grep -n -w -v "the" haiku.txt
1:The Tao that is seen
3:You bring fresh toner.
4:
5:With searching comes loss
7:"My Thesis" not found.
8:
9:Yesterday it worked
10:Today it is not working
11:Software is like that.

If we use the -r (recursive) option, grep can search for a pattern recursively through a set of files in subdirectories.

Let’s search recursively for Yesterday in the workshop_files_Bact_Genomics_2023/02_unix_intro/writing directory:

grep -r Yesterday .
./LittleWomen.txt:"Yesterday, when Aunt was asleep and I was trying to be
 as still as a
./LittleWomen.txt:Yesterday at dinner, when an Austrian officer stared at
 us and then
./LittleWomen.txt:Yesterday was a quiet day spent in teaching, sewing, and
 writing in my
./haiku.txt:Yesterday it worked

grep has lots of other options. To find out what they are, we can type:

grep --help
Usage: grep [OPTION]... PATTERN [FILE]...
Search for PATTERN in each FILE or standard input.
PATTERN is, by default, a basic regular expression (BRE).
Example: grep -i 'hello world' menu.h main.c

Regexp selection and interpretation:
  -E, --extended-regexp     PATTERN is an extended regular expression (ERE)
  -F, --fixed-strings       PATTERN is a set of newline-separated fixed strings
  -G, --basic-regexp        PATTERN is a basic regular expression (BRE)
  -P, --perl-regexp         PATTERN is a Perl regular expression
  -e, --regexp=PATTERN      use PATTERN for matching
  -f, --file=FILE           obtain PATTERN from FILE
  -i, --ignore-case         ignore case distinctions
  -w, --word-regexp         force PATTERN to match only whole words
  -x, --line-regexp         force PATTERN to match only whole lines
  -z, --null-data           a data line ends in 0 byte, not newline

Miscellaneous:
...        ...        ...
Exercise 2.7.1.1: Using grep

Which command would result in the following output:

and the presence of absence:
  1. grep "of" haiku.txt
  2. grep -E "of" haiku.txt
  3. grep -w "of" haiku.txt
  4. grep -i "of" haiku.txt

The correct answer is 3, because the -w option looks only for whole-word matches. The other options will also match ‘of’ when part of another word.

Wildcards

grep‘s real power doesn’t come from its options, though; it comes from the fact that patterns can include wildcards. (The technical name for these is regular expressions, which is what the ’re’ in ‘grep’ stands for.) Regular expressions are both complex and powerful; if you want to do complex searches, please visit his link.

As taster, we can find lines that have an ‘o’ in the second position like this:

grep -E "^.o" haiku.txt
You bring fresh toner.
Today it is not working
Software is like that.

We use the -E option and put the pattern in quotes to prevent the shell from trying to interpret it. (If the pattern contained a *, for example, the shell would try to expand it before running grep.) The ^ in the pattern anchors the match to the start of the line. The . matches a single character (just like ? in the shell), while the o matches an actual ‘o’.

Exercise 2.7.1.2: Tracking a Species

Leah has several hundred data files saved in one directory, each of which is formatted like this:

24/01/2023,Escherichia_coli,5
24/01/2023,Mycobacteria_tuberculosis,22
24/01/2023,Salmonella_enterica,7
24/01/2023,Staphylococcus_aureus,19
25/01/2023,Streptococcus_pneumoniae,2
25/01/2023,Pseudomonas_aeruginosa,1
25/01/2023,Haemophilus_influenzae ,18
26/01/2023,Carsonella_ruddii,1
26/01/2023,Escherichia_coli,6

She wants to write a shell script that takes a species as the first command-line argument and a directory as the second argument. The script should return one file called <species>.txt containing a list of dates and the number of that species seen on each date. For example using the data shown above, Escherichia_coli.txt would contain:

24/01/2023,5
26/01/2023,6

Below, each line contains an individual command, or pipe. Arrange their sequence in one command in order to achieve Leah’s goal:

cut -d : -f 2
grep -w $1 -r $2
$1.txt
cut -d , -f 1,3

Hint: use man grep to look for how to grep text recursively in a directory and man cut to select more than one field in a line. An example of such a file is provided in workshop_files_Bact_Genomics_2023/02_unix_intro/bacteria_counts.csv

grep -w $1 -r $2 | cut -d : -f 2 | cut -d , -f 1,3 > $1.txt

Actually, you can swap the order of the two cut commands and it still works. At the command line, try changing the order of the cut commands, and have a look at the output from each step to see why this is the case.

You would call the script above like this:

bash count-species.sh Mycobacteria_tuberculosis .
Exercise 2.7.1.3: Little Women

You and your friend, having just finished reading Little Women by Louisa May Alcott, are in an argument. Of the four sisters in the book, Jo, Meg, Beth, and Amy, your friend thinks that Jo was the most mentioned. You, however, are certain it was Amy. Luckily, you have a file LittleWomen.txt containing the full text of the novel (workshop_files_Bact_Genomics_2023/02_unix_intro/writing/LittleWomen.txt). Using a for loop, how would you tabulate the number of times each of the four sisters is mentioned? If you haven’t encountered for loops yet, you may want to go review the next bonus lesson on loops.

Hint: One solution might employ the commands grep and wc and a |, while another might utilize grep options. There is often more than one way to solve a programming task, so a particular solution is usually chosen based on a combination of yielding the correct result, elegance, readability, and speed.
for sis in Jo Meg Beth Amy
do
    echo $sis:
    grep -ow $sis LittleWomen.txt | wc -l
done

Alternative, slightly inferior solution:

for sis in Jo Meg Beth Amy
do
    echo $sis:
    grep -ocw $sis LittleWomen.txt
done

This solution is inferior because grep -c only reports the number of lines matched. The total number of matches reported by this method will be lower if there is more than one match per line.

Perceptive observers may have noticed that character names sometimes appear in all-uppercase in chapter titles (e.g. ‘MEG GOES TO VANITY FAIR’). If you wanted to count these as well, you could add the -i option for case-insensitivity (though in this case, it doesn’t affect the answer to which sister is mentioned most frequently).

2.7.2 find

While grep finds lines in files, the find command finds files themselves. Again, it has a lot of options; to show how the simplest ones work, we’ll use the workshop_files_Bact_Genomics_2023/02_unix_intro/ directory tree shown below.

.
├── bacteria_counts.csv
├── bacteria_rpob/
│   ├── bacteria_truncated_rpob.fasta
|   ├── Mycobacterium_tuberculosis_H37Rv_rpob.fasta
|   ├── rpob_gene.fasta
|   ├── Salmonella_enterica.fasta
|   ├── Staphylococcus_aureus.fasta
|   ├── Streptococcus_pneumoniae.fasta
├── bacteria.txt
├── creatures/
│   ├── basilisk
|   ├── minotaur.dat
|   ├── unicorn.dat
├── exp/
│   ├── 2015-10-23-calibration.txt
|   ├── 2015-10-23-dataset_overview.txt
|   ├── 2015-10-23-dataset1.txt
     ...
├── G26832.gff3.gz
├── G26832.tsv
├── molecules/
│   ├── cubane.pdb
│   ├── ethane.pdb
│   ├── methane.pdb
│   ├── octane.pdb
│   ├── pentane.pdb
│   └── propane.pdb
...
more here
...
└── writing/
    ├── haiku.txt
    └── LittleWomen.txt

The 02_unix_intro directory contains thirteen files, [bacteria_counts.csv, bacteria.txt, G26832.gff3.gz, G26832.tsv etc.] and ten directories: bacteria_rpob, creatures, exp, molecules etc., containing various files.

For our first command, let’s run find . (remember to run this command from the workshop_files_Bact_Genomics_2023/02_unix_intro folder).

find .
.
.
./sequencing_run1
./sequencing_run1/sample2_run1.fastq
./sequencing_run1/sample5_run1.fastq
./sequencing_run1/sample3_run1.fastq
./sequencing_run1/sample1_run1.fastq
./sequencing_run1/sample4_run1.fastq
./molecules
./molecules/propane.pdb
./molecules/.DS_Store
./molecules/octane.pdb
./molecules/cubane.pdb
./molecules/ethane.pdb
./molecules/pentane.pdb
./molecules/methane.pdb
./TBNmA041_annotation_truncated_1.tsv
...

As always, the . on its own means the current working directory, which is where we want our search to start. find’s output is the names of every file and directory under the current working directory. This can seem useless at first but find has many options to filter the output and in this lesson we will discover some of them.

The first option in our list is -type d that means ‘things that are directories’. Sure enough, find’s output is the names of the nine directories (including .):

find . -type d
.
./sequencing_run1
./molecules
./sequencing
./sequencing/run2
./sequencing/run1
./bacteria_rpob
./sequencing_run2
./north-pacific-gyre
./north-pacific-gyre/2012-07-03
./creatures
./exp
./mycobacteria_rpob
./writing

Notice that the objects find finds are not listed in any particular order. If we change -type d to -type f, we get a listing of all the files instead:

find . -type f
./sequencing_run1/sample2_run1.fastq
./sequencing_run1/sample5_run1.fastq
./sequencing_run1/sample3_run1.fastq
./sequencing_run1/sample1_run1.fastq
./sequencing_run1/sample4_run1.fastq
./molecules/propane.pdb
./molecules/.DS_Store
./molecules/octane.pdb
./molecules/cubane.pdb
./molecules/ethane.pdb
./molecules/pentane.pdb
./molecules/methane.pdb
./TBNmA041_annotation_truncated_1.tsv
...

Now let’s try matching by name:

find . -name *.txt
find: morse.txt: unknown primary or operator

We expected it to find all the text files, but it only prints out ./morse.txt. The problem is that the shell expands wildcard characters like * before commands run. Since *.txt in the current directory expands to ./morse.txt, the command we actually ran was:

find . -name morse.txt

find did what we asked; we just asked for the wrong thing.

To get what we want, let’s do what we did with grep: put *.txt in quotes to prevent the shell from expanding the * wildcard. This way, find actually gets the pattern *.txt, not the expanded filename morse.txt:

find . -name "*.txt"
./bacteria.txt
./morse.txt
./north-pacific-gyre/2012-07-03/NENE01729B.txt
./north-pacific-gyre/2012-07-03/NENE02040B.txt
./north-pacific-gyre/2012-07-03/NENE01729A.txt
./north-pacific-gyre/2012-07-03/NENE02040A.txt
./north-pacific-gyre/2012-07-03/NENE01971Z.txt
...
Listing vs. Finding

ls and find can be made to do similar things given the right options, but under normal circumstances, ls lists everything it can, while find searches for things with certain properties and shows them.

As we said earlier, the command line’s power lies in combining tools. We’ve seen how to do that with pipes; let’s look at another technique. As we just saw, find . -name "*.txt" gives us a list of all text files in or below the current directory. How can we combine that with wc -l to count the lines in all those files?

The simplest way is to put the find command inside $():

wc -l $(find . -name "*.txt")
       8 ./bacteria.txt
      50 ./morse.txt
     300 ./north-pacific-gyre/2012-07-03/NENE01729B.txt
     300 ./north-pacific-gyre/2012-07-03/NENE02040B.txt
     300 ./north-pacific-gyre/2012-07-03/NENE01729A.txt
     300 ./north-pacific-gyre/2012-07-03/NENE02040A.txt
     300 ./north-pacific-gyre/2012-07-03/NENE01971Z.txt
...

When the shell executes this command, the first thing it does is run whatever is inside the $(). It then replaces the $() expression with that command’s output. Since the output of find are all the filenames ./bacteria.txt, ./morse.txt, ./north-pacific-gyre/2012-07-03/NENE01729B.txt, etc. the shell constructs the command:

wc -l ./bacteria.txt ./morse.txt ./north-pacific-gyre/2012-07-03/NENE01729B.txt etc.

which is what we wanted. This expansion is exactly what the shell does when it expands wildcards like * and ?, but lets us use any command we want as our own ‘wildcard’.

It’s very common to use find and grep together. The first finds files that match a pattern; the second looks for lines inside those files that match another pattern. Here, for example, we can find txt files that contain the word “searching” by looking for the string ‘searching’ in all the .txt files in the current directory:

grep "searching" $(find . -name "*.txt")
./writing/LittleWomen.txt:sitting on the top step, affected to be searching
 for her book, but was
./writing/haiku.txt:With searching comes loss
Exercise 2.7.2.1: Matching and Subtracting

The -v option to grep inverts pattern matching, so that only lines which do not match the pattern are printed. Given that, which of the following commands will find all .pdb files in molecules/ except octane.pdb? Once you have thought about your answer, you can test the commands in the workshop_files_Bact_Genomics_2023/02_unix_intro/molecules directory.

  1. find molecules -name "*.pdb" | grep -v octane
  2. find molecules -name *.pdb | grep -v octane
  3. grep -v "octane" $(find molecules -name "*.pdb")
  4. None of the above.

Option 1 is correct. Putting the match expression in quotes prevents the shell expanding it, so it gets passed to the find command.

Option 2 also works in this instance because the shell tries to expand *.pdb but there are no *.pdb files in the current directory, so the wildcard expression gets passed to find.

Option 3 is incorrect because it searches the contents of the files for lines which do not match ‘octane’, rather than searching the file names.

Binary Files

We have focused exclusively on finding patterns in text files. What if your data is stored as images, in databases, or in some other format?

A handful of tools extend grep to handle a few non text formats. But a more generalizable approach is to convert the data to text, or extract the text-like elements from the data. On the one hand, it makes simple things easy to do. On the other hand, complex things are usually impossible. For example, it’s easy enough to write a program that will extract X and Y dimensions from image files for grep to play with, but how would you write something to find values in a spreadsheet whose cells contained formulas?

A last option is to recognize that the shell and text processing have their limits, and to use another programming language. When the time comes to do this, don’t be too hard on the shell: many modern programming languages have borrowed a lot of ideas from it, and imitation is also the sincerest form of praise.

The Unix shell is older than most of the people who use it. It has survived so long because it is one of the most productive programming environments ever created — maybe even the most productive. Its syntax may be cryptic, but people who have mastered it can experiment with different commands interactively, then use what they have learned to automate their work. Graphical user interfaces may be easier to use at first, but once learned, the productivity in the shell is unbeatable. And as Alfred North Whitehead wrote in 1911, ‘Civilization advances by extending the number of important operations which we can perform without thinking about them.’

Exercise 2.7.2.2: find Pipeline Reading Comprehension

Write a short explanatory comment for the following shell script:

wc -l $(find . -name "*.pdb") | sort -n
  1. Find all files with a .pdb extension recursively from the current directory
  2. Count the number of lines each of these files contains
  3. Sort the output from step 2, numerically

2.7.3 Credit

Information on this page has been adapted and modified from the following source:

  • https://github.com/swcarpentry/shell-novice