User talk:Darek Kedra/sandbox 28: Difference between revisions
Darek Kedra (talk | contribs) |
Darek Kedra (talk | contribs) |
||
Line 125: | Line 125: | ||
</pre> | </pre> | ||
=== Absolute vs relative directory naming== | === Absolute vs relative directory naming=== | ||
We can do all the operations on directories or files using two conventions: | We can do all the operations on directories or files using two conventions: | ||
<pre> | <pre> | ||
Line 140: | Line 140: | ||
This is useful shortcut for saving typing and avoiding errors when i.e. working on different systems. | This is useful shortcut for saving typing and avoiding errors when i.e. working on different systems. | ||
===Permissions=== | |||
In order to restrict actions associated with a given file or directory, each entity has 3 flags (r = read, w = write, e = execute) for 3 groups of users (file owner, group to which file owner belongs, i.e. students, and "all" for all the remaining users on that computer). So we have 9 fields describing what each of these groups can do with the file. On the top of it (or rather in the front of the string) we have another flag to tell us about what kind of thing it is ("-" = just a regular file, d = directory, l = link, etc.) | In order to restrict actions associated with a given file or directory, each entity has 3 flags (r = read, w = write, e = execute) for 3 groups of users (file owner, group to which file owner belongs, i.e. students, and "all" for all the remaining users on that computer). So we have 9 fields describing what each of these groups can do with the file. On the top of it (or rather in the front of the string) we have another flag to tell us about what kind of thing it is ("-" = just a regular file, d = directory, l = link, etc.) | ||
Line 170: | Line 170: | ||
This is often useful for sharing data / results with other users on the same machine or when writing scripts/executing some programs downloaded from the net. | This is often useful for sharing data / results with other users on the same machine or when writing scripts/executing some programs downloaded from the net. | ||
===copy, rename/move files, create symbolic links=== | |||
Often we need to organize files in directories from the existing ones, which requires moving thing around, making copies, renaming directories and files. To save space or to avoid using long path names, we can create links to other existing files or folders. | |||
CAVEAT: both cp and mv commands will perform the task without warning if the destination file exist. It is one of the ways of destroying your data. | |||
To avoid this error use "cp -i" or "mv -i" which will ask you for confirmation when you are overwriting existing file. | |||
<pre> | |||
#copy move rename pattern: | |||
cp/mv source destination | |||
cp file1.txt file2.txt | |||
cp file1.txt /some/directory/of/your/choice/ | |||
cp /home/linus/projecs/linus/file1.txt . #notice the space and the dot "." at the end of the command. The command is copy the file to the current directory | |||
</pre> | |||
Instead of doing copy (original file stays intact) one can do mv (move/rename) where only the new copy remains and the old one is destroyed. | |||
<pre> | |||
mv file1.txt file2.txt | |||
mv file1.txt /some/directory/of/your/choice/ | |||
mv /home/linus/projecs/linus/file1.txt . #notice the space and the dot "." at the end of the command. The command is copy the file to the current directory | |||
</pre> | |||
Linking is a way of organizing your data. I.e. you have one big genome fasta file, and want to map your next generation sequencing reads using different mappers, which will require using this one file as an input to several programs. You can store your fasta file in one directory (i.e. /home/vagrant/projects/chicken/genome1.fa) and create links to this file in several subdirectories /home/vagrant/projects/chicken/bwa_mapping/, /home/vagrant/projects/chicken/bowtie2_mapping/, etc.). The command: | |||
<pre> | |||
cd /home/vagrant/projects/chicken/bwa_mapping/ | |||
ln -s /home/vagrant/projects/chicken/genome1.fa . | |||
cd /home/vagrant/projects/chicken/bowtie2_mapping/ | |||
ln -s /home/vagrant/projects/chicken/genome1.fa . | |||
</pre> | |||
===view files (more/less, head, tail), count (wc)=== | |||
There are two basic types of files, text files and binary files. A lot of file formats in bioinformatics belong to simple text files, so the basic viewing/processing them is essential. Never ever do not try to open big text files (say few GB FASTQ file) in an editor, be it an editor on Linux or Microsoft Word... | |||
Basic checking of the file content: | |||
<pre> | |||
more myfile04.txt #displays content of the file screen by screen | |||
less myfile04.txt #same as above but you can scroll back | |||
head myfile04.txt #displays just the first 10 top lines of the file | |||
head -100 myfile04.txt #displays the first 100 lines | |||
tail myfile04.txt #displays just the last 10 lines of the file | |||
tail -100 myfile04.txt #displays the last 100 lines of the file | |||
wc myfile04.txt #newline, word, and byte counts for the file | |||
wc -l myfile04.txt # just count the lines of the file | |||
</pre> | |||
Using these commands is very useful for basic sanity check, so we | |||
#search for strings / replace strings (grep & sed) | #search for strings / replace strings (grep & sed) | ||
#compressing / uncompressing files (gzip, bzip2, tar) | #compressing / uncompressing files (gzip, bzip2, tar) |
Revision as of 08:28, 7 November 2013
Winterschool program
Software list
Basics
- linux Ubuntu 12.04.3 vs Debian 7.1 (think about 32 vs 64 bit versions)
- java http://www.java.com/en/download/linux_manual.jsp?locale=en
Specific tools 1
- TagDust: http://genome.gsc.riken.jp/osc/english/software/src/tagdust.tgz
- fastareformat from fastareformat exonerate-2.2.0 [1]
- fixing fasta headers (gff fields) with python? small script
- GEM [2]
- CAVEAT: (problem with cores on different laptops...)
http://sourceforge.net/projects/gemlibrary/files/gem-library/Binary%20pre-release%203/
- BWA http://sourceforge.net/projects/bio-bwa/files/
- Stampy http://www.well.ox.ac.uk/~gerton/software/Stampy/stampy-1.0.22r1848.tgz
- last http://last.cbrc.jp/ (the 362 versiona has split and splice-mapping options)
- bowtie http://bowtie-bio.sourceforge.net/bowtie2/index.shtml (bowtie2)
- samtools http://sourceforge.net/projects/samtools/files/
- picard http://sourceforge.net/projects/picard/files/
- IGV/ IGVtools http://www.broadinstitute.org/software/igv/download
- bamtools https://github.com/pezmaster31/bamtools
- requires cmake: http://www.cmake.org/files/v2.8/cmake-2.8.12.tar.gz (or apt get)
- bedtools http://code.google.com/p/bedtools/downloads/list
- GATK http://www.broadinstitute.org/gatk/auth?package=GATK (download yourself: license!)
- vcftools http://sourceforge.net/projects/vcftools/files/
Specific tools 2/RNA-Seq
- tophat http://tophat.cbcb.umd.edu/
- cufflinks http://cufflinks.cbcb.umd.edu/ (may require Boost libs!)
- GEMtools https://github.com/gemtools/gemtools
Vagrant fixes
For X11 forwarding the Vagrantfile has to contain
config.ssh.forward_x11 = true
Introduction to Linux and the command line
- why Linux?
- runs on everything from cell phones to supercomputers
- long history of stable tools
- most of the bioinformatics software was written and intended to run on Linux
logging in, connecting to other servers with ssh / sftp
As with other computers, one requires username password combination to connect to a specific computer. This combination can be specific to each of the computers or shared between i.e. all workstations at a given location.
SSH is a name for secure, encrypted connections between computers. It consist of two components, ssh server running on a remote machine and a ssh client on your laptop / workstation. The client is included in the default installations of recent Linux and OS X (Mac), but on Windows one has to install it ( http://www.chiark.greenend.org.uk/~sgtatham/putty/download.html ). It is not required you have it for the course, but it is a good idea to have it on your computer if you want to access remote servers on the command line.
With properly configured ssh connection (on Linux and Mac) you can run not just command line programs but also graphical user interface. This is also possible on Windows, but it is way more complicated to set up.
XXX Emilio: start vagrant or Linux from inside your virtual box machine, log in XXX
File and directories naming
Linux is case sensitive, do MyFile.txt is different from MYFILE.TXT or myfile.txt. Try to use some consistent naming schemes for your project directories, input data or result files. You can use very long, descriptive names, like these result files from ENCODE project:
wgEncodeUwTfbsNhdfneoCtcfStdAlnRep0.bam_VS_wgEncodeUwTfbsNhdfneoInputStdAlnRep1.bam.regionPeak.gz
Things to keep in mind:
- never ever use space/tabs in your file/directory names
- avoid Unix special characters in file names (!?"'%&^~*$|/\{}[]()<>:)
As with Windows and Mac systems one can imagine the systems of directories (synonym of folders) as a giant tree, where each of the branches has an upper directory and may contain lower level directories (sub-directories). At the top of this branched hierarchy it is just one place from where everything starts (no C:, D: X: drives), and in Unix speak it is called root directory. Some examples of directory naming:
/home/linus/ /usr/local/bin /home/linus/bioinf/programs /home/linus/projects/chicken
Individual users have separate directories in which they store their data, results etc.
Some shortcuts to remember (what is after # sign is a description):
/ #root directory ~ #your home directory .. #one directory above . #current directory
There are few commands to move between directories, create new ones, and list what is located there:
pwd #print working directory = show where you are ls #list what is located in the current directory ls -l #list what is located in the current directory with details cd /home/linus/projects/chicken #go to this directory cd - # special shortcut to go back to the directory you have been before cd / #go to the root directory cd ~ #go to your home directory spec cd .. #go todirectory one level above the current one mkdir mynew_directory #create new directory rmdir myold_directory #remove some directory (it must be empty!)
Absolute vs relative directory naming
We can do all the operations on directories or files using two conventions:
# absolute directory path starting with root "/" ls /home/linus/projects/chicken/genome.fasta #relative directory path (lets assume we are in /home/linus/projects/banana/) ls ../chicken/genome.fasta
This is useful shortcut for saving typing and avoiding errors when i.e. working on different systems.
Permissions
In order to restrict actions associated with a given file or directory, each entity has 3 flags (r = read, w = write, e = execute) for 3 groups of users (file owner, group to which file owner belongs, i.e. students, and "all" for all the remaining users on that computer). So we have 9 fields describing what each of these groups can do with the file. On the top of it (or rather in the front of the string) we have another flag to tell us about what kind of thing it is ("-" = just a regular file, d = directory, l = link, etc.)
By default, owner has read+write permissions, his group members can read but not write (modify or delete) his files, and the rest of the word should have no right to see content of his files. These permissions are visible when listing content of a directory with "ls -l":
-rwxr-xr-x 1 vagrant users 49 Nov 7 14:43 my_program.py -rw-r----- 1 vagrant users 1812 Nov 7 14:39 myfile01.txt ---------- 1 vagrant users 9045 Nov 7 14:41 myfile02.txt -r--r--r-- 1 vagrant users 2016 Nov 7 14:39 myfile03.txt -rw------- 1 vagrant users 67863 Nov 7 14:40 myfile04.txt -rw-rw-rw- 1 vagrant users 8125 Nov 7 14:40 myfile05.txt -rw-rw---- 1 vagrant users 7233 Nov 7 14:40 myfile06.txt lrwxrwxrwx 1 vagrant users 12 Nov 7 14:50 myfile07.txt -> myfile06.txt drwxr-xr-x 2 vagrant users 0 Nov 7 14:50 test_dir1
To change permissions, we have to specify first the group we want to modify, action (add permissions or remove them) then the permission themselves, finaly the name of the file(s):
chmod a+x my_new_script.py #add execute permission to for all chmod a-w # remove write permission for owner, group and the rest
This is often useful for sharing data / results with other users on the same machine or when writing scripts/executing some programs downloaded from the net.
copy, rename/move files, create symbolic links
Often we need to organize files in directories from the existing ones, which requires moving thing around, making copies, renaming directories and files. To save space or to avoid using long path names, we can create links to other existing files or folders.
CAVEAT: both cp and mv commands will perform the task without warning if the destination file exist. It is one of the ways of destroying your data. To avoid this error use "cp -i" or "mv -i" which will ask you for confirmation when you are overwriting existing file.
#copy move rename pattern: cp/mv source destination cp file1.txt file2.txt cp file1.txt /some/directory/of/your/choice/ cp /home/linus/projecs/linus/file1.txt . #notice the space and the dot "." at the end of the command. The command is copy the file to the current directory
Instead of doing copy (original file stays intact) one can do mv (move/rename) where only the new copy remains and the old one is destroyed.
mv file1.txt file2.txt mv file1.txt /some/directory/of/your/choice/ mv /home/linus/projecs/linus/file1.txt . #notice the space and the dot "." at the end of the command. The command is copy the file to the current directory
Linking is a way of organizing your data. I.e. you have one big genome fasta file, and want to map your next generation sequencing reads using different mappers, which will require using this one file as an input to several programs. You can store your fasta file in one directory (i.e. /home/vagrant/projects/chicken/genome1.fa) and create links to this file in several subdirectories /home/vagrant/projects/chicken/bwa_mapping/, /home/vagrant/projects/chicken/bowtie2_mapping/, etc.). The command:
cd /home/vagrant/projects/chicken/bwa_mapping/ ln -s /home/vagrant/projects/chicken/genome1.fa . cd /home/vagrant/projects/chicken/bowtie2_mapping/ ln -s /home/vagrant/projects/chicken/genome1.fa .
view files (more/less, head, tail), count (wc)
There are two basic types of files, text files and binary files. A lot of file formats in bioinformatics belong to simple text files, so the basic viewing/processing them is essential. Never ever do not try to open big text files (say few GB FASTQ file) in an editor, be it an editor on Linux or Microsoft Word...
Basic checking of the file content:
more myfile04.txt #displays content of the file screen by screen less myfile04.txt #same as above but you can scroll back head myfile04.txt #displays just the first 10 top lines of the file head -100 myfile04.txt #displays the first 100 lines tail myfile04.txt #displays just the last 10 lines of the file tail -100 myfile04.txt #displays the last 100 lines of the file wc myfile04.txt #newline, word, and byte counts for the file wc -l myfile04.txt # just count the lines of the file
Using these commands is very useful for basic sanity check, so we
- search for strings / replace strings (grep & sed)
- compressing / uncompressing files (gzip, bzip2, tar)
- pipelines and redirection
- awk in 5 minutes
- where to go from there (clusters, python)
FASTQ
- Illumina file formats (quality encodings)
- paired / unpaired reads
- quality checking (fastqc)
- trimming & filtering (TagDust)
- source of published FASTQ data: Short Read Archive vs ENA
Genomic fasta and gtf/gff gene annotation
- resources at ENSEMBL
- basic checks and reformatting
- grepping fasta headers
- fasta reformat from exonerate??
Mapping genomic reads
- overview of mappers
- GEM
- bwa +/- stampy
- last / bowtie
- mapping steps (for each mapper)
- genome indexing
- mapping
- +/- postprocessing
SAM and BAM file formats
- Analyzing BAM files
- sorting / indexing
- viewing the mappings in IGV
tools for processing BAM files
- samtools
- picard
- bamtools
getting mapping stats
- extracting reads mapping to regions
- getting coverage info for selected regions
Detecting SNPs
- general procedure
- GATK pipeline
- other SNP calling programs [tba]
Working with VCF files
- VCF file format
- viewing VCFs in IGV
- filtering SNPs by quality
- set operations on VCF files (common SNPs, unique SNPs)
RNASeq
- caveats (ribosomal RNA contamination)
- mapping RNASeq
- tophat
- GRAPE
- creating gene models from RNASeq (cufflinks)