Skip to content

Commit 61fd315

Browse files
committed
added -5SP to bwa mem
1 parent 4583536 commit 61fd315

1 file changed

Lines changed: 1 addition & 1 deletion

File tree

README.md

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -5,7 +5,7 @@ curl -L https:// | tar -jxf -
55
cd hickit-0.1_x64-linux
66

77
# Map Dip-C reads and extract contacts (skip if you use your own pipeline)
8-
./seqtk mergepe read1.fq.gz read2.fq.gz | ./pre-dip-c - | bwa mem -p hs37d5.fa - | gzip > aln.sam.gz
8+
./seqtk mergepe read1.fq.gz read2.fq.gz | ./pre-dip-c - | bwa mem -5SP -p hs37d5.fa - | gzip > aln.sam.gz
99
./k8 hickit.js vcf2tsv phased.vcf > phased_SNP.tsv # extract phased SNPs from VCF
1010
./k8 hickit.js sam2seg -v phased_SNP.tsv aln.sam.gz | ./k8 hickit.js chronly - | gzip > contacts.seg.gz
1111
./hickit -i contacts.seg.gz -o - | bgzip > contacts.pairs.gz # optional

0 commit comments

Comments
 (0)