anchor.diff:  
     sortedbam = per-fastq(
       fa:     fa,
  -    idx:    idx
  +    idx:    idx,
       fastq1: fastq1,
       fastq2: fastq2 );
   
     mergedbam = samtools-merge( bam: sortedbam );
   
     mpileup = samtools-mpileup(
  -    bam: mergedbam
  -    fa:  fa
  +    bam: mergedbam,
  +    fa:  fa,
       fai: fai );
   
     vcf = varscan( mpileup: mpileup );
  
positive.code:   fastq2: fastq2 );
    
  fun exonicfun = annovar(
    vcf:      vcf
    db:       db
    buildver: 'hg38' );
  
  
  
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";
  
  
  %% ============================================================
  