User:Lindenb/Notebook/UMR915/20101125

Dindel
Running Dindel on individual1. script so far: DINDEL_DIR=/usr/local/package/dindel-1.01 HG18=/GENOTYPAGE/data/pubdb/ucsc/hg18/chromosomes/hg18.fa BAM=XXXXXXXXXXXXXXXXXXXXXXXXXXXX.bam

stage1: ${DINDEL_DIR}/binaries/dindel-1.01-linux-64bit --analysis getCIGARindels --bamFile ${BAM} --outputFile sample.dindel_output --ref ${HG18} stage2: python ${DINDEL_DIR}/dindel-1.01-python/makeWindows.py --inputVarFile sample.dindel_output.variants.txt --windowFilePrefix sample.realign_windows --numWindowsPerFile 1000

stage3: for I in `ls sample.realign_windows.*.txt | cut -d '.' -f 3`; do ${DINDEL_DIR}/binaries/dindel-1.01-linux-64bit --analysis indels --doDiploid --bamFile ${BAM} --outputFile sample.dindel_stage2_output_windows.$$I  --ref ${HG18} --varFile sample.realign_windows.$$I.txt --libFile sample.dindel_output.libraries.txt ; done

Exome explorer
working on sources...