#!/bin/sh
awk -f /data/results/tools/rnaseq/collect-locus-coords-from-merged-gtf.awk /data/images/proton/run59/merged_asm/merged.gtf > merged-xloc-coords.csv
awk -v f="merged-xloc-coords.csv" -f /data/results/tools/rnaseq/filter-cuffdiff-genes1.awk gene_exp.diff  > gene_exp.diff.minlen100
# keep only status=="OK" and recalc q_values
Rscript /data/results/tools/rnaseq/re_adjust_q_values.r gene_exp.diff.minlen100