User Tools

Site Tools


2018pipeline

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
2018pipeline [2018/09/11 15:34]
david [Pipeline script]
2018pipeline [2018/09/11 15:44] (current)
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