Unix Bootcamp
Adapted by : Jason Walker, McDonnell Genome Institute
Additional adaptation by : Alex Wagner, McDonnell Genome Institute
Original author : Keith Bradnam, UC Davis Genome Center
Version 1.04 --- 2016-11-11
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Please send feedback, questions, money, or abuse to [email protected]
This 'bootcamp' is intended to provide the reader with a basic overview of essential Unix/Linux commands that will allow them to navigate a file system and move, copy, and edit files. It will also introduce a brief overview of some 'power' commands in Unix. It was orginally developed as part of a Bioinformatics Core Workshop taught at UC Davis (Using the Linux Command-Line for Analysis of High Throughput Sequence Data).
The Unix operating system has been around since 1969. Back then there was no such thing as a graphical user interface. You typed everything. It may seem archaic to use a keyboard to issue commands today, but it is much easier to automate keyboard tasks than mouse tasks. There are several variants of Unix (including Linux and OSX), though the differences do not matter much for most basic functions.
Increasingly, the raw output of biological research exists as in silico data, usually in the form of large text files. Unix is particularly suited to working with such files and has several powerful (and flexible) commands that can process your data for you. The real strength of learning Unix is that most of these commands can be combined in an almost unlimited fashion. So, if you can learn just five Unix commands, you will be able to do a lot more than just five things.
Command-line examples that you are meant to type into a terminal window will be shown indented in a constant-width font, e.g.
ls -lrh
Sometimes the accompanying text will include a reference to a Unix command. Any such text will also be in a constant-width, boxed font. E.g.
Type the
pwdcommand again.
From time to time this documentation will contain web links to pages that will help you find out more about certain Unix commands. Usually, the first mention of a command or function will be a hyperlink to Wikipedia. Important or critical points will be styled like so:
This is an important point!
The lessons from this point onwards will assume very little apart from the following:
- You have access to a Unix/Linux system
- You know how to launch a terminal program on that system
- You have a home directory where you can create/edit new files
In the following documentation, we will also assume that the logged in user has a username 'ubuntu' and the home directory is located at /home/ubuntu.
A terminal is the common name for the program that does two main things. It allows you to type input to the computer (i.e. run programs, move/view files etc.) and it allows you to see output from those programs. All Unix machines will have a terminal program available.
Open the terminal application. You should now see something that looks like the following:
There will be many situations where it will be useful to have multiple terminals open and it will be a matter of preference as to whether you want to have multiple windows, or one window with multiple tabs (there are typically keyboard shortcuts for switching between windows, or moving between tabs).
It's important to note that you will always be inside a single directory when using the terminal. The default behavior is that when you open a new terminal you start in your own home directory (containing files and directories that only you can modify). To see what files and directories are in our home directory, we need to use the ls command. This command lists the contents of a directory. If we run the ls command we should see something like:
ubuntu@:~$ ls
R bin tools workspace
ubuntu@:~$There are four things that you should note here:
- You will probably see different output to what is shown here, it depends on your computer setup. Don't worry about that for now.
- The
ubuntu@:~$text that you see is the Unix command prompt. In this case, it contains a user name ('ubuntu') and the name of the current directory ('~', more on that later). Note that the command prompt might not look the same on different Unix systems. In this case, the$sign marks the end of the prompt. - The output of the
lscommand lists two things. In this case, they are both directories, but they could also be files. We'll learn how to tell them apart later on. These directories were created as part of a specific course that used this bootcamp material. You will therefore probably see something very different on your own computer. - After the
lscommand finishes it produces a new command prompt, ready for you to type your next command.
The ls command is used to list the contents of any directory, not necessarily the one that you are currently in. Try the following:
ubuntu@:~$ ls workspace/
data lib lost+found
ubuntu@:~$ ls /etc/perl
CPAN NetLooking at directories from within a Unix terminal can often seem confusing. But bear in mind that these directories are exactly the same type of folders that you can see if you use any graphical file browser. From the root level (/) there are usually a dozen or so directories. You can treat the root directory like any other, e.g. you can list its contents:
ubuntu@:~$ ls /
bin dev home lib lost+found mnt proc run snap sys usr vmlinuz
boot etc initrd.img lib64 media opt root sbin srv tmp var workspaceYou might notice some of these names appearing in different colors. Many Unix systems will display files and directories differently by default. Other colors may be used for special types of files. When you log in to a computer you are working with your files in your home directory, and this is often inside a directory called 'users' or 'home'.
There may be many hundreds of directories on any Unix machine, so how do you know which one you are in? The command pwd will Print the Working Directory and that's pretty much all this command does:
ubuntu@:~$ pwd
/home/ubuntuWhen you log in to a Unix computer, you are typically placed into your home directory. In this example, after we log in, we are placed in a directory called 'ubuntu' which itself is a subdirectory of another directory called 'home'. Conversely, 'users' is the parent directory of 'clmuser'. The first forward slash that appears in a list of directory names always refers to the top level directory of the file system (known as the root directory). The remaining forward slash (between 'home' and 'ubuntu') delimits the various parts of the directory hierarchy. If you ever get 'lost' in Unix, remember the pwd command.
As you learn Unix you will frequently type commands that don't seem to work. Most of the time this will be because you are in the wrong directory, so it's a really good habit to get used to running the pwd command a lot.
If we want to make a new directory (e.g. to store some work related data), we can use the mkdir command:
ubuntu@:~$ mkdir workspace/Learning_unix
ubuntu@:~$ ls workspace
data Learning_unix lib lost+foundWe are in the home directory on the computer but we want to to work in the new Learning_unix directory. To change directories in Unix, we use the cd command:
cd workspace/Learning_unix
ubuntu@:~/workspace/Learning_unix$Notice that — on this system — the command prompt has expanded to include our current directory. This doesn't happen by default on all Unix systems, but you should know that you can configure what information appears as part of the command prompt.
Let's make two new subdirectories and navigate into them:
ubuntu@:~/workspace/Learning_unix$ mkdir Outer_directory
ubuntu@:~/workspace/Learning_unix$ cd Outer_directory
ubuntu@:~/workspace/Learning_unix/Outer_directory$
ubuntu@:~/workspace/Learning_unix/Outer_directory$ mkdir Inner_directory
ubuntu@:~/workspace/Learning_unix/Outer_directory$ cd Inner_directory/
ubuntu@:~/workspace/Learning_unix/Outer_directory/Inner_directory$Now our command prompt is getting quite long, but it reveals that we are four levels beneath the home directory. We created the two directories in separate steps, but it is possible to use the mkdir command in way to do this all in one step.
Like most Unix commands, mkdir supports command-line options which let you alter its behavior and functionality. Command-like options are — as the name suggests — optional arguments that are placed after the command name. They often take the form of single letters (following a dash). If we had used the -p option of the mkdir command we could have done this in one step. E.g.
mkdir -p Outer_directory/Inner_directoryNote the spaces either side the -p!
Sometimes options are entire words, usually preceded by a double-dash. For example, using the --parents option:
mkdir --parents Outer_directory/Inner_directoryis identical to using the -p option.
Many programs will provide information about the command being called by passing a -h or --help option.
ubuntu@:~$mkdir --help
Usage: mkdir [OPTION]... DIRECTORY...
Create the DIRECTORY(ies), if they do not already exist.
Mandatory arguments to long options are mandatory for short options too.
-m, --mode=MODE set file mode (as in chmod), not a=rwx - umask
-p, --parents no error if existing, make parent directories as needed
-v, --verbose print a message for each created directory
-Z set SELinux security context of each created directory
to the default type
--context[=CTX] like -Z, or if CTX is specified then set the SELinux
or SMACK security context to CTX
--help display this help and exit
--version output version information and exit
GNU coreutils online help: <http://www.gnu.org/software/coreutils/>
Full documentation at: <http://www.gnu.org/software/coreutils/mkdir>
or available locally via: info '(coreutils) mkdir invocation'In addition, many commands will have man (manual) pages, which are often more detailed descriptions of the program and its usage.
man ls
man cd
man man # yes even the man command has a manual pageWhen you are using the man command, press space to scroll down a page, b to go back a page, or q to quit. You can also use the up and down arrows to scroll a line at a time. The man command is actually using another Unix program, a text viewer called less, which we'll come to later on.
Let's change directory to the root directory, and then into our home directory
ubuntu@:~/workspace/Learning_unix/Outer_directory/Inner_directory$ cd /
ubuntu@:/$ cd home
ubuntu@:/home$ cd ubuntu
ubuntu@:~$In this case, we may as well have just changed directory in one go:
cd /home/ubuntu/The leading / is incredibly important. The following two commands are very different:
cd /home/ubuntu/
cd home/ubuntu/The first command says go the unbuntu directory that is beneath the home directory that is at the top level (the root) of the file system. There can only be one /home/ubuntu directory on any Unix system.
The second command says go to the unbuntu directory that is beneath the home directory that is located wherever I am right now. There can potentially be many home/ubuntu directories on a Unix system (though this is unlikely).
Learn and understand the difference between these two commands.
Frequently, you will find that you want to go 'upwards' one level in the directory tree. Two dots .. are used in Unix to refer to the parent directory of wherever you are. Every directory has a parent except the root level of the computer. Let's go into the Learning_unix directory and then navigate up three levels:
ubuntu@:~$ cd workspace/Learning_unix/
ubuntu@:~/workspace/Learning_unix$ cd ..
ubuntu@:~/workspace$ cd ..
ubuntu@:~$ cd ..
ubuntu@:/home$What if you wanted to navigate up two levels in the file system in one go? It's very simple, just use two sets of the .. operator, separated by a forward slash:
cd ../..Using cd .. allows us to change directory relative to where we are now. You can also always change to a directory based on its absolute location. E.g. if you are working in the /home/ubuntu/workspace/Learning_unix directory and you then want to change to the /tmp directory, then you could do either of the following:
$ cd ../../../../tmp
or...
$ cd /tmp
They both achieve the same thing, but the 2nd example requires that you know about the full path from the root level of the computer to your directory of interest (the 'path' is an important concept in Unix). Sometimes it is quicker to change directories using the relative path, and other times it will be quicker to use the absolute path.
Remember that the command prompt shows you the name of the directory that you are currently in, and that when you are in your home directory it shows you a tilde character (`~``) instead? This is because Unix uses the tilde character as a short-hand way of specifying a home directory.
See what happens when you try the following commands (use the pwd command after each one to confirm the results if necessary):
cd /
cd ~
cd
cd -Hopefully, you should find that cd and cd ~ do the same thing, i.e. they take you back to your home directory (from wherever you were). You will frequently want to jump straight back to your home directory, and typing cd is a very quick way to get there.
You can also use the ~ as a quick way of navigating into subdirectories of your home directory when your current directory is somewhere else. I.e. the quickest way of navigating from the root directory to your Learning_unix directory is as follows:
ubuntu@:~$ cd /
ubuntu@:/$ cd ~/workspace/Learning_unixThe .. operator that we saw earlier can also be used with the ls command, e.g. you can list directories that are 'above' you:
ubuntu@:~/workspace/Learning_unix$ cd ~/Learning_unix/Outer_directory/
ubuntu@:~/workspace/Learning_unix/Outer_directory$ ls ../../
command_line_course Learning_unix linux_bootcampTime to learn another useful command-line option. If you add the letter 'l' to the ls command it will give you a longer output compared to the default:
ubuntu@:~/workspace/Learning_unix$ ls -l /home
total 4
drwxr-xr-x 12 ubuntu ubuntu 4096 Nov 12 01:45 ubuntuFor each file or directory we now see more information (including file ownership and modification times). The 'd' at the start of each line indicates that these are directories. There are many, many different options for the ls command. Try out the following (against any directory of your choice) to see how the output changes.
ls -l
ls -R
ls -l -t -r
ls -lhNote that the last example combine multiple options but only use one dash. This is a very common way of specifying multiple command-line options.
We now have a few (empty) directories that we should remove. To do this use the rmdir command, this will only remove empty directories so it is quite safe to use. If you want to know more about this command (or any Unix command), then remember that you can just look at its man page.
ubuntu@:~$ cd ~/workspace/Learning_unix/Outer_directory/
ubuntu@:~/workspace/Learning_unix/Outer_directory$ rmdir Inner_directory/
ubuntu@:~/workspace/Learning_unix/Outer_directory$ cd ..
ubuntu@:~/workspace/Learning_unix$ rmdir Outer_directory/
ubuntu@:~/workspace/Learning_unix$ ls
ubuntu@:~/workspace/Learning_unix$EXERCISE:
Recreate the directories you just removed with a single command. Use the--helpoption if you need to!
Now remove them again, but this time in a single command.
*** Note, you have to be outside a directory before you can remove it with rmdir ***
Saving keystrokes may not seem important, but the longer that you spend typing in a terminal window, the happier you will be if you can reduce the time you spend at the keyboard. Especially, as prolonged typing is not good for your body. So the best Unix tip to learn early on is that you can tab complete the names of files and programs on most Unix systems. Type enough letters that uniquely identify the name of a file, directory or program and press tab...Unix will do the rest. E.g. if you type 'tou' and then press tab, Unix should autocomplete the word to 'touch' (this is a command which we will learn more about in a minute). In this case, tab completion will occur because there are no other Unix commands that start with 'tou'. If pressing tab doesn't do anything, then you have not have typed enough unique characters. In this case pressing tab twice will show you all possible completions. This trick can save you a LOT of typing!
Navigate to your home directory, and then use the cd command to change to the Learning_unix directory. Use tab completion to complete directory name. If there are no other directories starting with 'L' in your home directory, then you should only need to type 'cd' + 'L' + 'tab'.
Tab completion will make your life easier and make you more productive!
Another great time-saver is that Unix stores a list of all the commands that you have typed in each login session. You can access this list by using the history command or more simply by using the up and down arrows to access anything from your history. So if you type a long command but make a mistake, press the up arrow and then you can use the left and right arrows to move the cursor in order to make a change.
The following sections will deal with Unix commands that help us to work with files, i.e. copy files to/from places, move files, rename files, remove files, and most importantly, look at files. First, we need to have some files to play with. The Unix command touch will let us create a new, empty file. The touch command does other things too, but for now we just want a couple of files to work with.
ubuntu@:~$ cd workspace/Learning_unix/
ubuntu@:~/workspace/Learning_unix$ touch red_fish.txt
ubuntu@:~/workspace/Learning_unix$ touch blue_fish.txt
ubuntu@:~/workspace/Learning_unix$ ls
red_fish.txt blue_fish.txttouch also accepts multiple files as arguments.
ubuntu@:~/workspace/Learning_unix$ touch one_fish.txt two_fish.txt
ubuntu@:~/workspace/Learning_unix$ ls
blue_fish.txt one_fish.txt red_fish.txt two_fish.txtNow, let's assume that we want to move these files to a new directory ('colors'). We will do this using the Unix mv (move) command. Remember to use tab completion:
ubuntu@:~/workspace/Learning_unix$ mkdir colors
ubuntu@:~/workspace/Learning_unix$ mv red_fish.txt colors/
ubuntu@:~/workspace/Learning_unix$ mv blue_fish.txt colors/
ubuntu@:~/workspace/Learning_unix$ ls
colors one_fish.txt two_fish.txt
ubuntu@:~/workspace/Learning_unix$ ls colors/
blue_fish.txt red_fish.txtEXERCISE:
Make a new directory called 'counts', and move one_fish.txt and two_fish.txt into it. Can you move the two files with a singlemvcommand?
For the mv command, we always have to specify a source file (or directory) that we want to move, and then specify a target location. If we had wanted to we could have moved both files in one go by typing any of the following commands:
mv *.txt Temp/
mv *t Temp/
mv *ea* Temp/The asterisk * acts as a wild-card character, essentially meaning 'match anything'. The second example works because there are no other files or directories in the directory that end with the letters 't' (if there was, then they would be moved too). Likewise, the third example works because only those two files contain the letters 'ea' in their names. Using wild-card characters can save you a lot of typing.
The '?' character is also a wild-card but for only a single character.
In the earlier example, the destination for the mv command was a directory name (colors). So we moved a file from its source location to a target location, but note that the target could have also been a (different) file name, rather than a directory. E.g. let's make a new file and move it whilst renaming it at the same time:
ubuntu@:~/workspace/Learning_unix$ touch rags
ubuntu@:~/workspace/Learning_unix$ ls
colors counts rags
ubuntu@:~/workspace/Learning_unix$ mv rags counts/riches
ubuntu@:~/workspace/Learning_unix$ ls counts/
earth.txt heaven.txt richesIn this example we create a new file ('rags') and move it to a new location and in the process change the name (to 'riches'). So mv can rename a file as well as move it. The logical extension of this is using mv to rename a file without moving it (you have to use mv to do this as Unix does not have a separate 'rename' command):
ubuntu@:~/workspace/Learning_unix$ mv counts/riches counts/ragsIt is important to understand that as long as you have specified a 'source' and a 'target' location when you are moving a file, then it doesn't matter what your current directory is. You can move or copy things within the same directory or between different directories regardless of whether you are in any of those directories. Moving directories is just like moving files:
ubuntu@:~/workspace/Learning_unix$ mkdir fish
ubuntu@:~/workspace/Learning_unix$ mv counts fish
ubuntu@:~/workspace/Learning_unix$ ls -R .
.:
colors fish
./colors:
blue_fish.txt red_fish.txt
./fish:
counts
./fish/counts:
one_fish.txt rags two_fish.txtThis step moves the counts directory inside the fish directory.
EXERCISE:
Try creating a 'net' directory inside 'Learning_unix' and thencdto your home directory. Can you movefishinsidenetwithout usingcd?
You've seen how to remove a directory with the rmdir command, but rmdir won't remove directories if they contain any files. So how can we remove the files we have created (inside Learning_Unix/Temp)? In order to do this, we will have to use the rm (remove) command.
Please read the next section VERY carefully. Misuse of the
rmcommand can lead to needless death & destruction
Potentially, rm is a very dangerous command; if you delete something with rm, you will not get it back! It is possible to delete everything in your home directory (all directories and subdirectories) with rm, that is why it is such a dangerous command.
Let me repeat that last part again. It is possible to delete EVERY file you have ever created with the rm command. Are you scared yet? You should be. Luckily there is a way of making rm a little bit safer. We can use it with the -i command-line option which will ask for confirmation before deleting anything (remember to use tab-completion):
ubuntu@:~/workspace/Learning_unix$ cd net/fish/counts
ubuntu@:~/workspace/Learning_unix/net/fish/counts$ ls
one_fish.txt rags two_fish.txt
ubuntu@:~/workspace/Learning_unix/net/fish/counts$ rm -i one_fish.txt rags two_fish.txt
rm: remove regular empty file 'one_fish.txt'? y
rm: remove regular empty file 'rags'? y
rm: remove regular empty file 'two_fish.txt'? y
ubuntu@:~/workspace/Learning_unix/net/fish/counts$ lsWe could have simplified this step by using a wild-card (e.g. rm -i *.txt) or we could have made things more complex by removing each file with a separate rm command. Let's finish cleaning up:
ubuntu@:~/workspace/Learning_unix/net/fish/counts$ cd ~/workspace/Learning_unix/
ubuntu@:~/workspace/Learning_unix$ rmdir -p net/fish/counts/
ubuntu@:~/workspace/Learning_unix$ rm -ir colors/
rm: descend into directory 'colors/'? y
rm: remove regular empty file 'colors/red_fish.txt'? y
rm: remove regular empty file 'colors/blue_fish.txt'? y
rm: remove directory 'colors/'? yCopying files with the cp (copy) command is very similar to moving them. Remember to always specify a source and a target location. Let's create a new file and make a copy of it:
ubuntu@:~/workspace/Learning_unix$ touch file1
ubuntu@:~/workspace/Learning_unix$ cp file1 file2
ubuntu@:~/workspace/Learning_unix$ ls
file1 file2What if we wanted to copy files from a different directory to our current directory? Let's put a file in our home directory (specified by ~ remember) and copy it to the current directory (Learning_unix):
ubuntu@:~/workspace/Learning_unix$ touch ~/file3
ubuntu@:~/workspace/Learning_unix$ ls ~
file3 tools workspace
ubuntu@:~/workspace/Learning_unix$ cp ~/file3 .
ubuntu@:~/workspace/Learning_unix$ ls
file1 file2 file3This last step introduces another new concept. In Unix, the current directory can be represented by a . (dot) character. You will mostly use this only for copying files to the current directory that you are in. Compare the following:
ls
ls .
ls ./In this case, using the dot is somewhat pointless because ls will already list the contents of the current directory by default. Also note how the trailing slash is optional. You can use rm to remove the temporary files.
The cp command also allows us (with the use of a command-line option) to copy entire directories. Use cp --help to see how the -R or -r options let you copy a directory recursively.
EXERCISE:
Create a new directory called 'filing_cabinet' and move the files into it withmv. Make a second copy of this directory called 'trash_can'. Were all of the files from 'filing_cabinet' also in 'trash_can'? Remove both directories and their contents with a single command.
So far we have covered listing the contents of directories and moving/copying/deleting either files and/or directories. Now we will quickly cover how you can look at files. The less command lets you view (but not edit) text files. We will use the echo command to put some text in a file and then view it:
ubuntu@:~/workspace/Learning_unix$ echo "Call me Ishmael."
Call me Ishmael.
ubuntu@:~/workspace/Learning_unix$ echo "Call me Ishmael." > opening_lines.txt
ubuntu@:~/workspace/Learning_unix$ ls
opening_lines.txt
ubuntu@:~/workspace/Learning_unix$ less opening_lines.txtOn its own, echo isn't a very exciting Unix command. It just echoes text back to the screen. But we can redirect that text into an output file by using the > symbol. This allows for something called file redirection.
Careful when using file redirection (>), it will overwrite any existing file of the same name
When you are using less, you can bring up a page of help commands by pressing h, scroll forward a page by pressing space, or go forward or backwards one line at a time by pressing j or k. To exit less, press q (for quit). The less program also does about a million other useful things (including text searching).
Let's add another line to the file:
ubuntu@:~/workspace/Learning_unix$ echo "The primroses were over." >> opening_lines.txt
ubuntu@:~/workspace/Learning_unix$ cat opening_lines.txt
Call me Ishmael.
The primroses were over.Notice that we use >> and not just >. This operator will append to a file. If we only used >, we would end up overwriting the file. The cat command displays the contents of the file (or files) and then returns you to the command line. Unlike less you have no control on how you view that text (or what you do with it). It is a very simple, but sometimes useful, command. You can use cat to quickly combine multiple files or, if you wanted to, make a copy of an existing file:
cat opening_lines.txt > file_copy.txtubuntu@:~/workspace/Learning_unix$ ls
opening_lines.txt
ubuntu@:~/workspace/Learning_unix$ ls -l
total 4
-rw-rw-r-- 1 ubuntu ubuntu 42 Jun 15 04:13 opening_lines.txt
ubuntu@:~/workspace/Learning_unix$ wc opening_lines.txt
2 7 42 opening_lines.txt
ubuntu@:~/workspace/Learning_unix$ wc -l opening_lines.txt
2 opening_lines.txtThe ls -l option shows us a long listing, which includes the size of the file in bytes (in this case '42'). Another way of finding this out is by using Unix's wc command (word count). By default this tells you many lines, words, and characters are in a specified file (or files), but you can use command-line options to give you just one of those statistics (in this case we count lines with wc -l).
Nano is a lightweight editor installed on most Unix systems. There are many more powerful editors (such as 'emacs' and 'vi'), but these have steep learning curves. Nano is very simple. You can edit (or create) files by typing:
nano opening_lines.txtYou should see the following appear in your terminal:
The bottom of the nano window shows you a list of simple commands which are all accessible by typing 'Control' plus a letter. E.g. Control + X exits the program.
One other use of the echo command is for displaying the contents of something known as environment variables. These contain user-specific or system-wide values that either reflect simple pieces of information (your username), or lists of useful locations on the file system. Some examples:
ubuntu@:~/workspace/Learning_unix$ echo $USER
ubuntu
ubuntu@:~/workspace/Learning_unix$ echo $HOME
/home/ubuntu
ubuntu@:~/workspace/Learning_unix$ echo $PATH
/home/ubuntu/bin:/home/ubuntu/.local/bin:/home/ubuntu/tools/perl5/bin:/home/ubuntu/tools/bin:/home/ubuntu/workspace/data/anaconda/bin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/usr/local/games:/snap/bin:/home/ubuntu/tools/bowtie-1.1.2:/home/ubuntu/tools/bowtie2-2.2.9:/home/ubuntu/tools/trinityrnaseq-2.2.0:/home/ubuntu/tools/hisat2-2.0.4:/home/ubuntu/tools/sambamba_v0.6.4:/home/ubuntu/tools/stringtie-1.3.0.Linux_x86_64:/home/ubuntu/tools/gffcompare-0.9.8.Linux_x86_64:/home/ubuntu/tools/RSEM-1.2.31:/home/ubuntu/tools/cufflinks-2.2.1.Linux_x86_64:/home/ubuntu/tools/bedtools2/bin:/home/ubuntu/tools/MUMmer3.23:/home/ubuntu/tools/allpathslg-52488/bin:/home/ubuntu/tools/bin/Sniffles/bin/sniffles-core-1.0.0:/home/ubuntu/tools/ensembl-tools-release-86/scripts/variant_effect_predictor:/home/ubuntu/tools/VAAST_2.2.0/bin:/home/ubuntu/tools/speedseq/bin:/home/ubuntu/tools/hall_miscThe last one shows the content of the $PATH environment variable, which displays a — colon separated — list of directories that are expected to contain programs that you can run. This includes all of the Unix commands that you have seen so far. These are files that live in directories which are run like programs (e.g. ls is just a special type of file in the /bin directory).
Knowing how to change your $PATH to include custom directories can be necessary sometimes (e.g. if you install some new bioinformatics software in a non-standard location).
Use nano to add the following lines to opening_lines.txt:
Now is the winter of our discontent.
All children, except one, grow up.
The Galactic Empire was dying.
In a hole in the ground there lived a hobbit.
It was a pleasure to burn.
It was a bright, cold day in April, and the clocks were striking thirteen.
It was love at first sight.
I am an invisible man.
It was the day my grandmother exploded.
When he was nearly thirteen, my brother Jem got his arm badly broken at the elbow.
Marley was dead, to begin with.
You will often want to search files to find lines that match a certain pattern. The Unix command grep does this (and much more). The following examples show how you can use grep's command-line options to:
- show lines that match a specified pattern
- ignore case when matching (
-i) - only match whole words (
-w) - show lines that don't match a pattern (
-v) - Use wildcard characters and other patterns to allow for alternatives (
*,., and[])
grep was opening_lines.txt
The Galactic Empire was dying.
It was a pleasure to burn.
It was a bright, cold day in April, and the clocks were striking thirteen.
It was love at first sight.
It was the day my grandmother exploded.
When he was nearly thirteen, my brother Jem got his arm badly broken at the elbow.
Marley was dead, to begin with.
grep -v was opening_lines.txt
Call me Ishmael.
The primroses were over.
Now is the winter of our discontent.
All children, except one, grow up.
In a hole in the ground there lived a hobbit.
I am an invisible man.
grep all opening_lines.txt
Call me Ishmael.
grep -i all opening_lines.txt
Call me Ishmael.
All children, except one, grow up.
grep in opening_lines.txt
Now is the winter of our discontent.
The Galactic Empire was dying.
In a hole in the ground there lived a hobbit.
It was a bright, cold day in April, and the clocks were striking thirteen.
I am an invisible man.
Marley was dead, to begin with.
grep -w in opening_lines.txt
In a hole in the ground there lived a hobbit.
It was a bright, cold day in April, and the clocks were striking thirteen.
grep -w o.. opening_lines.txt
Now is the winter of our discontent.
All children, except one, grow up.
grep [aeiou]t opening_lines.txt
In a hole in the ground there lived a hobbit.
It was love at first sight.
It was the day my grandmother exploded.
When he was nearly thirteen, my brother Jem got his arm badly broken at the elbow.
Marley was dead, to begin with.
grep -w -i [aeiou]t opening_lines.txt
It was a pleasure to burn.
It was a bright, cold day in April, and the clocks were striking thirteen.
It was love at first sight.
It was the day my grandmother exploded.
When he was nearly thirteen, my brother Jem got his arm badly broken at the elbow.One of the most poweful features of Unix is that you can send the output from one command or program to any other command (as long as the second commmand accepts input of some sort). We do this by using what is known as a pipe. This is implemented using the '|' character (which is a character which always seems to be on different keys depending on the keyboard that you are using). Think of the pipe as simply connecting two Unix programs. Here's an example which introduces some new Unix commands:
ubuntu@:~/workspace/Learning_unix$ grep was opening_lines.txt | wc -c
316
ubuntu@:~/workspace/Learning_unix$
grep was opening_lines.txt | sort | head -n 3 | wc -c
130The first use of grep searches the specified file for lines matching 'was', it sends the lines that match through a pipe to the wc program. We use the -c option to just count characters in the matching lines (316).
The second example first sends the output of grep to the Unix sort command. This sorts a file alphanumerically by default. The sorted output is sent to the head command which by default shows the first 10 lines of a file. We use the -n option of this command to only show 3 lines. These 3 lines are then sent to the wc command as before.
Whenever making a long pipe, test each step as you build it!
Interactive Tutorial Cheat sheet Another Cheat Sheet
scp for transferring files over SSH
wget for downloading files from URL
watch for monitoring
xargs for complex inputs
ps for process status
top for running processes
kill for stopping processes
The following examples introduce some other Unix commands, and show how they could be used to work on a fictional file called file.txt. Remember, you can always learn more about these Unix commands from their respective man pages with the man command. These are not all real world cases, but rather show the diversity of Unix command-line tools:
tail -n 20 file.txt | head- Show lines of a file that begin with a start codon (ATG) (the
^matches patterns at the start of a line):
grep "^ATG" file.txt- Cut out the 3rd column of a tab-delimited text file and sort it to only show unique lines (i.e. remove duplicates):
cut -f 3 file.txt | sort -u- Count how many lines in a file contain the words 'cat' or 'bat' (
-coption ofgrepcounts lines):
grep -c '[bc]at' file.txt- Turn lower-case text into upper-case (using tr command to 'transliterate'):
cat file.txt | tr 'a-z' 'A-Z'- Change all occurrences of 'Chr1' to 'Chromosome 1' and write changed output to a new file (using sed command):
cat file.txt | sed 's/Chr1/Chromosome 1/' > file2.txt- 2015-06-14 - Version 1.0, adapted from Unix and Perl for Biologists Primer
- 2015-06-24 - Version 1.01: clarified that this material is assuming user name is 'ubuntu'and made other minor clarifications (such as what this material was first produced for).
- 2015-11-13 - Version 1.02: further adapted for CSHL 2015 Advanced Sequencing Technologies & Applications course.
- 2016-11-11 - Version 1.04: additional material for CSHL 2016 Advanced Sequencing Technologies & Applications course.
| Previous Section | This Section | Next Section |
|---|---|---|
| Logging-into-Amazon-Cloud | Unix Bootcamp | Environment |