Skip to content

Commit

Permalink
unicycler: make gfa output really gfa (#3213)
Browse files Browse the repository at this point in the history
* unicycler: make gfa output really gfa

* generalize gfa output assertions

seems that the sequence in the gfa is reverse complememted (and maybe
also shifted). which is because its circular. so the old assertion did
not work.

* test unicycler
  • Loading branch information
bernt-matthias authored Nov 18, 2020
1 parent 0a6ad31 commit d7f4eb8
Show file tree
Hide file tree
Showing 2 changed files with 9 additions and 10 deletions.
1 change: 0 additions & 1 deletion .tt_skip
Original file line number Diff line number Diff line change
Expand Up @@ -7,5 +7,4 @@ tools/migmap/
tools/ncbi_entrez_eutils/einfo.xml
tools/sarscov2formatter/sarscov2formatter.xml
tools/trinotate
tools/unicycler
tools/valet
18 changes: 9 additions & 9 deletions tools/unicycler/unicycler.xml
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
<tool id="unicycler" name="Create assemblies with Unicycler" version="@[email protected]">
<tool id="unicycler" name="Create assemblies with Unicycler" version="@[email protected]" profile="20.09">
<macros>
<token name="@VERSION@">0.4.8</token>
</macros>
Expand Down Expand Up @@ -204,7 +204,7 @@ $pilon.no_pilon
</section>
</inputs>
<outputs>
<data name="assembly_graph" format="tabular" from_work_dir="assembly.gfa" label="${tool.name} on ${on_string}: Final Assembly Graph" />
<data name="assembly_graph" format="gfa1" from_work_dir="assembly.gfa" label="${tool.name} on ${on_string}: Final Assembly Graph" />
<data name="assembly" format="fasta" from_work_dir="assembly.fasta" label="${tool.name} on ${on_string}: Final Assembly"/>
</outputs>
<tests>
Expand Down Expand Up @@ -240,9 +240,9 @@ $pilon.no_pilon
<section name="lr_align">
<param name="scores" value="3,-6,-5,-2"/>
</section>
<output name="assembly_graph" ftype="tabular">
<output name="assembly_graph" ftype="gfa1">
<assert_contents>
<has_text text="TATCTGTTACTGAGAAGTTAATGGATGAATTGGCAC"/>
<has_line_matching expression="S\t1\t[ATCG]{5386,5386}\tLN:i:5386\tdp:f:1.0"/>
</assert_contents>
</output>
<output name="assembly" ftype="fasta">
Expand Down Expand Up @@ -295,9 +295,9 @@ $pilon.no_pilon
<section name="lr_align">
<param name="scores" value="3,-6,-5,-2"/>
</section>
<output name="assembly_graph" ftype="tabular">
<output name="assembly_graph" ftype="gfa1">
<assert_contents>
<has_text text="TATCTGTTACTGAGAAGTTAATGGATGAATTGGCAC" />
<has_line_matching expression="S\t1\t[ATCG]{5386,5386}\tLN:i:5386\tdp:f:1.0"/>
</assert_contents>
</output>
<output name="assembly" ftype="fasta">
Expand Down Expand Up @@ -342,9 +342,9 @@ $pilon.no_pilon
<section name="lr_align">
<param name="scores" value="3,-6,-5,-2"/>
</section>
<output name="assembly_graph" ftype="tabular">
<output name="assembly_graph" ftype="gfa1">
<assert_contents>
<has_text text="TATCTGTTACTGAGAAGTTAATGGATGAATTGGCAC" />
<has_line_matching expression="S\t1\t[ATCG]{5386,5386}\tLN:i:5386\tdp:f:1.0"/>
</assert_contents>
</output>
<output name="assembly" ftype="fasta">
Expand All @@ -362,7 +362,7 @@ $pilon.no_pilon
<param name="kmers" value="21,23"/>
</section>
<param name="long" value="only_long.fasta" ftype="fasta" />
<output name="assembly_graph" ftype="tabular">
<output name="assembly_graph" ftype="gfa1">
<assert_contents>
<has_text text="S" />
</assert_contents>
Expand Down

0 comments on commit d7f4eb8

Please sign in to comment.