anchor.diff:  : idx( File ) [fastq1( File ) fastq2( File )] )in bash *{
   
     tar xf $idx
  +  bam=alignment.bam
     bowtie2 -D 5 -R 1 -N 0 -L 22 -i S,0,2.50 --no-unal -x bt2idx \
     -1 $fastq1 -2 $fastq2 -S - | samtools view -b - > $bam
     rm bt2idx.*
  
positive.code: 
  % gunzip
  deftask gunzip( out( File ) : gz( File ) )in bash *{
    gzip -c -d $gz > $out
  }*
  
  
negative.code: fun exonicfun = annovar(
    vcf:      vcf,
    db:       db,
    buildver: "hg38" );
  
  
  %% ============================================================
  