schatzlab / crossstitch

Code for phasing SVs with SNPs
52 stars 3 forks source link

MINQ error in extractHAIRS step #10

Open JYLeeBioinfo opened 5 years ago

JYLeeBioinfo commented 5 years ago

Hi

I am trying to phase SVs using Nanopore data

I encountered the following error.

MINQ must be at least 4

So I added --ONT 1 to the following command in the src/crossstitch.sh

(original)
$EXTRACTHAIRS --mbq 0 --bam $LONGREADSBAM --VCF $PHASEDSNPS.prehairs --out $OUTPREFIX.hairs

(changed)
$EXTRACTHAIRS --mbq 0 --bam $LONGREADSBAM --VCF $PHASEDSNPS.prehairs --out $OUTPREFIX.hairs --ONT 1

All processes finished well after changed it.

Is it okay to add ONT option to extractHAIRs step?

mschatz commented 5 years ago

We haven’t tested it but should be okay

Good luck

Mike

On Jun 22, 2019, at 11:52 AM, hd00ljy notifications@github.com wrote:

Hi

I am trying to phase SVs using Nanopore data

I encountered the following error.

MINQ must be at least 4 So I added --ONT 1 to the following command in the src/crossstitch.sh

(original) $EXTRACTHAIRS --mbq 0 --bam $LONGREADSBAM --VCF $PHASEDSNPS.prehairs --out $OUTPREFIX.hairs

(changed) $EXTRACTHAIRS --mbq 0 --bam $LONGREADSBAM --VCF $PHASEDSNPS.prehairs --out $OUTPREFIX.hairs --ONT 1 All processes finished well after changed it.

Is it okay to add ONT option to extractHAIRs step?

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub, or mute the thread.

JYLeeBioinfo commented 5 years ago

Thank you Mike!

JYLeeBioinfo commented 5 years ago

It seems I have to add reference if I add --ONT 1

image