Skip to content

Commit

Permalink
Updating arcs-long with streaming reads (#111)
Browse files Browse the repository at this point in the history
* Stream reads from long-to-linked-pe to arcs --arks in long mode

* Update min molecule size for arcs-long

* Remove -g parameter for long-to-linked-pe
  • Loading branch information
lcoombe authored Apr 27, 2021
1 parent 2a81cc9 commit 491ff8f
Showing 1 changed file with 10 additions and 4 deletions.
14 changes: 10 additions & 4 deletions Examples/arcs-make
Original file line number Diff line number Diff line change
Expand Up @@ -255,14 +255,20 @@ else
endif

# long kmer
%_c$c_m$m_cut$(cut)_k$k_r$r_e$e_z$z_original.gv: %.renamed.fa $(reads).cut$(cut).fq.gz
$(reads).barcode-multiplicity.tsv: $(long_reads)
$(gtime) $(bin)/../src/long-to-linked-pe -l $(cut) -t $t -m2000 --bx-only -b $@ $<

%_c$c_m$m_cut$(cut)_k$k_r$r_e$e_z$z_original.gv: %.renamed.fa $(long_reads) $(reads).barcode-multiplicity.tsv
ifneq ($D, true)
$(gtime) arcs --arks -v -f $< -c $c -m $m -r $r -e $e -z $z -j $j -k $k -t $t -d $d --gap $(gap) -b $(patsubst %_original.gv,%,$@) $(word 2,$^) --barcode-counts $(barcode_counts).tsv
$(gtime) $(bin)/../src/long-to-linked-pe -l $(cut) -t $t -m2000 $(long_reads) |\
arcs --arks -v -f $< -c $c -m $m -r $r -e $e -z $z -j $j -k $k -t $t -d $d --gap $(gap) -b $(patsubst %_original.gv,%,$@) -u $(reads).barcode-multiplicity.tsv /dev/stdin

else ifneq ($(dist_upper), true)
$(gtime) arcs --arks -v -D -B $B -f $< -c $c -m $m -r $r -e $e -z $z -j $j -k $k -t $t -d $d --gap $(gap) -b $(patsubst %_original.gv,%,$@) $(word 2,$^) --barcode-counts $(barcode_counts).tsv
$(gtime) $(bin)/../src/long-to-linked-pe -l $(cut) -t $t -m2000 $(long_reads) |\
arcs --arks -v -D -B $B -f $< -c $c -m $m -r $r -e $e -z $z -j $j -k $k -t $t -d $d --gap $(gap) -b $(patsubst %_original.gv,%,$@) -u $(reads).barcode-multiplicity.tsv /dev/stdin
else
$(gtime) arcs --arks -v -D -B $B --dist_upper -f $< -c $c -m $m -r $r -e $e -z $z -j $j -k $k -t $t -d $d --gap $(gap) -b $(patsubst %_original.gv,%,$@) $(word 2,$^) --barcode-counts $(barcode_counts).tsv
$(gtime) $(bin)/../src/long-to-linked-pe -l $(cut) -t $t -m2000 $(long_reads) |\
arcs --arks -v -D -B $B --dist_upper -f $< -c $c -m $m -r $r -e $e -z $z -j $j -k $k -t $t -d $d --gap $(gap) -b $(patsubst %_original.gv,%,$@) -u $(reads).barcode-multiplicity.tsv /dev/stdin
endif

# Generate TSV from ARCS
Expand Down

0 comments on commit 491ff8f

Please sign in to comment.