GMAP error: invalid literal for int() with base 10

7 months ago by

United States

Jennifer Hillman Jackson ♦ 25k wrote:

Hello,

There are two versions for GMAP/MiGMAP in the Tool Shed: https://usegalaxy.org/toolshed

Both are a bit older and I am not sure of the exact usage or if there are any known issues. However, I can let you know that this error usually comes up when a tool is given an input that has formatting outside of what the 3rd party wrapped tool is expecting.

Example: MACS2 >> If unmapped reads are retained in a SAM input, there can be sequence identifiers that contain spaces, triggering this same error. For that case, there are two options: Remove the unmapped reads or convert the SAM (plain text) to BAM (binary) format to resolve the problem.

You'll need to review your inputs, this tool, and adjust to match what the tool is expecting. If you can't figure this out, do one of these (the second is much faster, as it is for most IUC tools):

Both of the tool versions in the ToolShed have authors that can be reached for usage advice. How to: log in or create an account, navigate to the tool, and a method for contact is under "Repository Actions" in the top right corner ).
Both authors can also be reached at the IUC Gitter channel for more direct communication plus help from the wider community who may be familiar with these older tool wrappers. https://gitter.im/galaxy-iuc/iuc

Thanks! Jen, Galaxy team

ADD COMMENT • link written 7 months ago by Jennifer Hillman Jackson ♦ 25k

I decided to ask for help at Gitter to get this started. Please join in on the conversation and supply more details, or they might reply here and ask about that first unless this is a known (similar to the MACS2 issue): https://gitter.im/galaxy-iuc/iuc?at=5acbc6fa27c509a774f48743

ADD REPLY • link written 7 months ago by Jennifer Hillman Jackson ♦ 25k

Thanks a lot ! and I add a custom GMAP to my local galaxy, eventually.

Another question about GMAP. I want to align genes to different reference sequences. but I'm not sure how many refererence sequences exactly. So, I hope to make a judge. It means if users have one reference, he just can align genes to the reference sequences one time. If he have two reference sequences, he can align them to the two sequences.

Can it come true by edit XML file?

I searched some answers. Someone said <if id="if-1"> <time from="5pm" to="9pm"/> </if> <then id="if-1"> <...some actions defined.../> </then> <else id="if-1"> <...other set of actions defined here.../> </else>

other said <xsl:choose>

<xsl:when test="expression">

    actions...

</xsl:when>

<xsl:when test="expression">

    actions...

</xsl:when>

<xsl:otherwise>

    actions that occur if none of the above are detected

</xsl:otherwise>

</xsl:choose>

but they seems no effect.

Thanks, Siyuan

ADD REPLY • link written 7 months ago by chenzhuod • 30

Most if not all of the existing Galaxy tools allow the selection of a single target reference genome at runtime. If there are multiple genome targets, the tool could be put into a workflow multiple times and the jobs will run together.

That said, if you still want to explore this option, try using Planemo for development. If you need specific help with any part during/after the development process, please consider joining the development community: galaxy-dev Mailing list, galaxyproject/iuc Gitter, publishing tools to the ToolShed, etc.

ADD REPLY • link written 7 months ago by Jennifer Hillman Jackson ♦ 25k

diff -r 93911bac43da gmap.xml --- a/gmap.xml Thu Jan 05 14:31:24 2012 -0600 +++ b/gmap.xml Wed May 09 23:42:32 2018 +0100 @@ -113,7 +113,7 @@ --wraplength=$advanced.wraplength #end if #end if - #if $split_output == True + #if $split_output $split_output #end if #if len($quality_protocol.__str__) > 0: @@ -123,7 +123,7 @@ #for $i in $inputs: ${i.added_input} #end for - #if $split_output == True + #if $split_output 2> $gmap_stderr #else 2> $gmap_stderr > $output

Similar posts • Search »