projectTitle: Detection of RNA Outlier Pipeline 1015
root: /vcu_gpfs2/home/gaynora/DROP/udx_1015/Output
htmlOutputPath: /vcu_gpfs2/home/gaynora/DROP/udx_1015/Output/html
indexWithFolderName: true

hpoFile: null
sampleAnnotation: /vcu_gpfs2/home/gaynora/DROP/udx_1015/sample_annotation_1015.tsv

geneAnnotation:
    v104: /vcu_gpfs2/home/gaynora/genomedirectory/collapsed_Homo_sapiens.GRCh38.104.gtf
genomeAssembly: hg38
genome:
    GRCh38: /vcu_gpfs2/home/gaynora/genomedirectory/Homo_sapiens.GRCh38.dna.primary_assembly.fa

#exportCounts:
    #geneAnnotations:
        #- v104
    #excludeGroups:
        #- mae
        #- outrider_external
        #- fraser_external

aberrantExpression:
    run: true
    groups:
        - outrider
    minlds: 1
    fpkmCutoff: 1
    implementation: autoencoder
    padjCutoff: 1.0
    zScoreCutoff: 0
    genesToTest: null
    maxTestedDimensionProportion: 3
    yieldSize: 2000000

aberrantSplicing:
    run: true
    groups:
        - fraser
    recount: false
    longRead: false
    keepNonStandardChrs: false
    filter: true
    minExpressionInOneSample: 20
    quantileMinExpression: 10
    minDeltaPsi: 0.05
    implementation: PCA
    padjCutoff: 1.0
    maxTestedDimensionProportion: 6
    genesToTest: null
    ### For FRASER2, use the follwing parameters instead of the 3 lines above:
    FRASER_version: "FRASER2"
    deltaPsiCutoff : 0.1
    quantileForFiltering: 0.75

mae:
    run: false
    groups:
        - group1
        - group2
        - group3
    gatkIgnoreHeaderCheck: true
    padjCutoff: 0.05
    allelicRatioCutoff: 0.8
    addAF: true
    maxAF: 0.001

rnaVariantCalling:
    run: false
    groups:
        - batch_0
    highQualityVCFs:
        - Data/Mills_and_1000G_gold_standard.indels.hg19.sites.chrPrefix.vcf.gz
        - Data/1000G_phase1.snps.high_confidence.hg19.sites.chrPrefix.vcf.gz
    dbSNP: Data/00-All.vcf.gz
    repeat_mask: Data/hg19_repeatMasker_sorted.chrPrefix.bed
    createSingleVCF: true
    addAF: true
    maxAF: 0.001
    maxVarFreqCohort: 0.05
    hcArgs: ""
    minAlt: 3
    yieldSize: 100000

tools:
    gatkCmd: gatk
    bcftoolsCmd: bcftools
    samtoolsCmd: samtools