Quick Start
From PileLine
| Line 23: | Line 23: | ||
* Variants (vs reference genome) pileup file. | * Variants (vs reference genome) pileup file. | ||
| - | '''2.''' You may compare Case1 vs Control1 at variant level using this command line: | + | '''2.''' You may compare Case1 vs Control1 at variant level using pileline-2smc functionality. Try this command line: |
$ cd DOWNLOADED_FILES_DIRECTORY | $ cd DOWNLOADED_FILES_DIRECTORY | ||
$ sh YOUR_PATH_TO_PILELINE/cmd/pileline-2smc.sh | $ sh YOUR_PATH_TO_PILELINE/cmd/pileline-2smc.sh | ||
–a ./Control1.pileup –b ./Case1.pileup | –a ./Control1.pileup –b ./Case1.pileup | ||
–v ./Control1varfilter.pileup –w ./Case1varfilter.pileup | –v ./Control1varfilter.pileup –w ./Case1varfilter.pileup | ||
| - | –o ./ | + | –o ./myoutput1.txt |
Executing this code you will obtain 3 output files: | Executing this code you will obtain 3 output files: | ||
| Line 35: | Line 35: | ||
* [[FileNAME3]]: Variants found in Case1 and in Control1 (i.e. germinal mutations or SNPs) | * [[FileNAME3]]: Variants found in Case1 and in Control1 (i.e. germinal mutations or SNPs) | ||
| - | '''3.''' You can also compare several | + | Run pileline-2smc to compare Case2 vs Control2: |
| + | $ sh YOUR_PATH_TO_PILELINE/cmd/pileline-2smc.sh | ||
| + | –a ./Control2.pileup –b ./Case2.pileup | ||
| + | –v ./Control2varfilter.pileup –w ./Case2varfilter.pileup | ||
| + | –o ./myoutput2.txt | ||
| + | |||
| + | '''3.''' You can also compare several samples at variant level: | ||
$ sh YOUR_PATH_TO_PILELINE/cmd/pileline-nsmc.sh | $ sh YOUR_PATH_TO_PILELINE/cmd/pileline-nsmc.sh | ||
| - | -a | + | -a -a |
| - | -b | + | -b -b |
'''4.''' At this point it could be useful to annotate SNPs. Execute this command line to this end : | '''4.''' At this point it could be useful to annotate SNPs. Execute this command line to this end : | ||
Revision as of 15:41, 18 June 2010
PipeLine input files
PileLine is capable to handle, filter and compare genomic position files (GP) including standard pileup, BED,GFF, or VCF files.
Basically, GP are tabular files where the two first columns contain chromosome name and position coordinate respectively. Additional optional fields are accepted in PileLine, see an example of GP input file below:
10 118829 optional1 optional2 optional3 ... 10 121207 optional1 optional2 optional3 ... 10 121337 optional1 optional2 optional3 ... 10 121636 optional1 optional2 optional3 ...
PipeLine Guided example
1. Download GP example files (pileup format) to your working directory:
Each .zip file contains 2 pileup files:
- Whole pileup file
- Variants (vs reference genome) pileup file.
2. You may compare Case1 vs Control1 at variant level using pileline-2smc functionality. Try this command line:
$ cd DOWNLOADED_FILES_DIRECTORY $ sh YOUR_PATH_TO_PILELINE/cmd/pileline-2smc.sh –a ./Control1.pileup –b ./Case1.pileup –v ./Control1varfilter.pileup –w ./Case1varfilter.pileup –o ./myoutput1.txt
Executing this code you will obtain 3 output files:
- FileNAME1: Variants found in Case1 but not in Control1 (i.e. somatic point mutations or SNPs)
- FileNAME2: Variants found in Control1 but not in Case1 (i.e. somatic point mutations or SNPs)
- FileNAME3: Variants found in Case1 and in Control1 (i.e. germinal mutations or SNPs)
Run pileline-2smc to compare Case2 vs Control2:
$ sh YOUR_PATH_TO_PILELINE/cmd/pileline-2smc.sh –a ./Control2.pileup –b ./Case2.pileup –v ./Control2varfilter.pileup –w ./Case2varfilter.pileup –o ./myoutput2.txt
3. You can also compare several samples at variant level:
$ sh YOUR_PATH_TO_PILELINE/cmd/pileline-nsmc.sh -a -a -b -b
4. At this point it could be useful to annotate SNPs. Execute this command line to this end :
