User Tools
2018pipeline
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision | |||
|
2018pipeline [2018/09/11 15:34] david [Pipeline script] |
2018pipeline [2018/09/11 15:44] david [Pipeline script] |
||
|---|---|---|---|
| Line 471: | Line 471: | ||
| # pipeline.sbatch - Perform grep on a GFF file, then run faFrag to extract DNA for | # pipeline.sbatch - Perform grep on a GFF file, then run faFrag to extract DNA for | ||
| # each line returned by grep. | # each line returned by grep. | ||
| - | + | ||
| + | annotations=chr_wormbase.gff | ||
| + | |||
| + | if [ ! -e "$annotations" ] | ||
| + | then | ||
| + | echo "Need to supply an annotation file" >&2 | ||
| + | exit 1 | ||
| + | fi | ||
| # Specify Gene list file | # Specify Gene list file | ||
| Line 479: | Line 486: | ||
| # Perform grep -f on Gene list file to get annotation lines from the GFF | # Perform grep -f on Gene list file to get annotation lines from the GFF | ||
| - | grep -f $gene_lst chr_wormbase.gff | grep five_prime_utr > tmpfile.gff | + | grep -f $gene_lst $annotations | grep five_prime_utr > tmpfile.gff |
| echo "temporary file contains $(wc -l tmpfile.gff) lines" | echo "temporary file contains $(wc -l tmpfile.gff) lines" | ||
| # | # | ||
2018pipeline.txt · Last modified: 2018/09/11 15:44 by david
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
