anchor.diff:    fastq2: fastq2 );
     
   fun exonicfun = annovar(
  -  vcf:      vcf
  -  db:       db
  +  vcf:      vcf,
  +  db:       db,
     buildver: 'hg38' );
   
   
  
positive.code: %% Input data
  %% ============================================================
  
  hg38-tar = 'hg38/hg38.tar';
  
  fastq1-gz = 'kgenomes/SRR359188_1.filt.fastq.gz'
              'kgenomes/SRR359195_1.filt.fastq.gz';
  
  fastq2-gz = 'kgenomes/SRR359188_2.filt.fastq.gz'
              'kgenomes/SRR359195_2.filt.fastq.gz';
  
  db = '/opt/data/annodb_hg38';
  
  
  %% ============================================================
  
negative.code: 
    sortedbam = per-fastq(
      fa:     fa,
      idx:    idx,
      fastq1: fastq1,
      fastq2: fastq2 );
  
    mergedbam = samtools-merge( bam: sortedbam );
  
    mpileup = samtools-mpileup(
      bam: mergedbam,
      fa:  fa,
      fai: fai );
  
    vcf = varscan( mpileup: mpileup );
  