Attention:
The frontline support team will be unavailable to answer questions on April 15th and 17th 2019. We will be back soon after. Thank you for your patience and we apologize for any inconvenience!

I need some tips to run GenomeSTRiP - CNVDiscovery

Hi,
I'm trying to call CNVs by GenomeSTRiP software.
As I know, I should run SVpreprocess first and CNVDiscovery second.
So I tried SVpreprocess but It won't work.
The errer said they need more option "-ploidyMapFile" even I only put autosome.
Is it essential to call CNVs from only autosome?
then how can i get the file? I'm working on pig data..

And if my process was wrong, or whatever, I want to know the exact process(or cord) to call CNVs..

(+)

!/bin/bash

export SV_DIR=/directory/to/svtoolkit
classpath="${SV_DIR}/lib/SVToolkit.jar:${SV_DIR}/lib/gatk/GenomeAnalysisTK.jar:${SV_DIR}/lib/gatk/Queue.jar"

SVpreprocess

java -Xmx40g -cp ${classpath} \
org.broadinstitute.gatk.queue.QCommandLine \
-S ${SV_DIR}/qscript/SVPreprocess.q \
-S ${SV_DIR}/qscript/SVQScript.q \
-cp ${classpath} \
-gatk ${SV_DIR}/lib/gatk/GenomeAnalysisTK.jar \
-configFile ${SV_DIR}/conf/genstrip_parameters.txt \
-R Ref/Sus_scrofa.Sscrofa11.1.dna.toplevel.fa \
-I DURvsKWB.bam.list \
-md mdtest/ \
-bamFilesAreDisjoint true \
-jobLogDir CNVtest/logDir2 \
-run

CNVDiscovery

java -Xmx40g -cp ${classpath} \
org.broadinstitute.gatk.queue.QCommandLine \
-S ${SV_DIR}/qscript/discovery/cnv/CNVDiscoveryPipeline.q \
-S ${SV_DIR}/qscript/SVQScript.q \
-cp ${classpath} \
-gatk ${SV_DIR}/lib/gatk/GenomeAnalysisTK.jar \
-configFile ${SV_DIR}/conf/genstrip_parameters.txt \
-R Ref/Sus_scrofa.Sscrofa11.1.dna.toplevel.fa \
-I DURvsKWB.bam.list \
-genderMapFile DURvsKWB.gender.txt \
-md mdtest/ \
-runDirectory CNVtest/ \
-jobLogDir CNVtest/log \
-intervalList chr.list \
-tilingWindowSize 1000 \
-tilingWindowOverlap 500 \
-maximumReferenceGapLength 1000 \
-boundaryPrecision 100 \
-minimumRefinedLength 500 \
-run

I attached the cord i used.

Answers

Sign In or Register to comment.