Skip to content

Commit

Permalink
update gap names so consensus handles them properly
Browse files Browse the repository at this point in the history
  • Loading branch information
skoren committed Aug 10, 2024
1 parent bb80b76 commit 2fd2479
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions src/Snakefiles/4-processONT.sm
Original file line number Diff line number Diff line change
Expand Up @@ -88,9 +88,9 @@ echo Step 2a

cat alns-ont-nogap-2.gaf alns-ont-gap1.gaf alns-ont-gap2.gaf | {PYTHON} {VERKKO}/scripts/maybe_trim_alignment.py gapped-twice-unitig-unrolled-hifi-resolved.gfa 100 selected-winnowmap.ids > alns-cut-trim.gaf

{PYTHON} {VERKKO}/scripts/insert_aln_gaps.py gapped-twice-unitig-unrolled-hifi-resolved.gfa alns-cut-trim.gaf 10 100 alns-ont-nogap-3.gaf alns-ont-gap3.gaf cutgap o > gapped-thrice-unitig-unrolled-hifi-resolved.gfa
{PYTHON} {VERKKO}/scripts/insert_aln_gaps.py gapped-twice-unitig-unrolled-hifi-resolved.gfa alns-cut-trim.gaf 10 100 alns-ont-nogap-3.gaf alns-ont-gap3.gaf gapcut o > gapped-thrice-unitig-unrolled-hifi-resolved.gfa

{PYTHON} {VERKKO}/scripts/insert_aln_gaps.py gapped-thrice-unitig-unrolled-hifi-resolved.gfa alns-ont-nogap-3.gaf 2 0 alns-ont-nogap.gaf alns-ont-gap4.gaf extragap1 o > gapped-unitig-unrolled-hifi-resolved.gfa
{PYTHON} {VERKKO}/scripts/insert_aln_gaps.py gapped-thrice-unitig-unrolled-hifi-resolved.gfa alns-ont-nogap-3.gaf 2 0 alns-ont-nogap.gaf alns-ont-gap4.gaf gapextra1 o > gapped-unitig-unrolled-hifi-resolved.gfa

cat alns-ont-gap1.gaf alns-ont-gap2.gaf alns-ont-gap3.gaf alns-ont-gap4.gaf > ../{output.ont_gap_align}

Expand Down

0 comments on commit 2fd2479

Please sign in to comment.