Question: bowtie2 rename dataset from paired list
0
gravatar for brigida.rusconi
3.1 years ago by
United States
brigida.rusconi40 wrote:

Hi,

I have a paired reads dataset list that I use for a bowtie2 workflow. I try to rename the output file with the #{input_1} option but I get an empty string for the output file. I tried to use the #{own_file} that renames according to the reference sequence and that works. I have tried different combinations of input with or without 1 or even using _2 since it is paired end fastq, but nothing worked. I looked at the wrapper and in the input section it is called input_1 so I don't understand why it is not propagating the name. Is there another way to get the identifier from the paired list to rename the output file according to the initial fastq?

Best,

Brigida

 

bowtie galaxy • 1.0k views
ADD COMMENTlink modified 2.7 years ago • written 3.1 years ago by brigida.rusconi40

I tried the renaming with the paired and there it works. I think that in Bowtie2 it has an issue with recognizing the paired_collection. I can do the Bowtie2 with a paired_collection, but the workflow manager shows me this error when I am selecting the option "History does not include a dataset collection of the correct type or containing the correct types of datasets". Despite the error it will produce bam files for all the pairs in the collection, but I cannot get the name of each pair on the list. All the files will have an identical name e.g. Bowtie2 on data 11, data 4, and others: aligned reads (sorted BAM). This is why I need to rename them.

ADD REPLYlink written 3.1 years ago by brigida.rusconi40
0
gravatar for brigida.rusconi
3.1 years ago by
United States
brigida.rusconi40 wrote:

Hi,

So we changed the xml file of the Bowtie to

<data format="bam" name="output" label="${library.input_1.name}">


This renames the bam files to the same name as the pairs. I still see the error in the workflow regarding the absence of a paired list in my history, but it still runs. In case anyone else wants to have the output of the list named like the input for traceability I recommend doing the change in the output section.

 

Brigida

 

ADD COMMENTlink written 3.1 years ago by brigida.rusconi40

Does this work on a single BAM file as input.

ADD REPLYlink written 2.8 years ago by Bjoern Gruening5.1k

I have not checked in a while but I recall that Bowtie is one of the tools that will not rename. There are others and it has been that way for a while. A: Without 'Available inputs are:' it is not possible to rename datasets in workflo

ADD REPLYlink written 2.7 years ago by Guy Reeves1.0k
0
gravatar for brigida.rusconi
2.7 years ago by
United States
brigida.rusconi40 wrote:

I am not sure what you mean by single bam file input. There will be only one bam file generated for each paired end read dataset that you put into Bowtie2. If you are using a single reads file it should work too, because we take the name from the first input in the list.

ADD COMMENTlink written 2.7 years ago by brigida.rusconi40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 169 users visited in the last hour