anchor.diff:    fastq2: fastq2 );
     
   fun exonicfun = annovar(
  -  vcf:      vcf
  -  db:       db
  +  vcf:      vcf,
  +  db:       db,
     buildver: 'hg38' );
   
   
  
positive.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 );
  
negative.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";
  
  
  %% ============================================================
  