-
Notifications
You must be signed in to change notification settings - Fork 4
Description
Hi,
Thank you for developing such a tool.
I have an issue in reading param files. Before this issue, I had another issue in reading input files but I figured it out.
I followed steps below:
- BASE="${HOME}/margin_phasing"
- INPUT_DIR="${BASE}/input/data"
- REF="Mus_musculus..fa"
- BAM="Remora.bam"
- paramfile="allParams.haplotag.ont-r94g507.json"
- THREADS="10"
- OUTPUT_DIR="${BASE}/output"
- OUTPUT_PREFIX="Ep300"
- OUTPUT_VCF="Ep300.vcf.gz"
- TRUTH_VCF="FVB_NJ.vcf"
- mkdir -p "${OUTPUT_DIR}"
- mkdir -p "${INPUT_DIR}"
- paramdirectory="${BASE}/paramfiles"
- mkdir -p "${paramdirectory}"
- sudo docker pull jmcdani20/hap.py:v0.3.12
- sudo docker pull kishwars/pepper_deepvariant:r0.8
- sudo docker run -v "${INPUT_DIR}":"${INPUT_DIR}"
-v "${OUTPUT_DIR}":"${OUTPUT_DIR}"
kishwars/pepper_deepvariant:r0.8 margin phase
"${INPUT_DIR}/${BAM}"
"${INPUT_DIR}/${REF}"
"${INPUT_DIR}/${TRUTH_VCF}"
"${paramdirectory}/${paramfile}"
-o "${OUTPUT_DIR}"
-p "${OUTPUT_PREFIX}"
-t ${THREADS}
and then
I got this error: "Could not read from params file: /home/ko/margin_phasing/paramfiles/allParams.haplotag.ont-r94g507.json".
I also created a directory where paramfile is existed. But, when my run read inputs, it did not read paramfile with the same way. In addition to that, I put a path where I have all param files default (I downloaded them from GitHub). Despite that, I got the same error again. How can I do my run can read param files? Could you help me with this?
I addition to that, another question in my mind is that is margin consistent with R.10.4.1 flowcell output of ONT?
Kind regards.