Thanks for your answers. I understand that permit to load or create
"unsorted bam" (or sorted by name) can create bug with lot of tools
and it's difficult to restrict this manipulation just for few tools.
But I don't understand one thing. Why can I download from the toolshed
a tool which permit me to sort bam file by name (from devteam) while
it's not supported by galaxy ?!
Other question. Anyone can help me to write a new entry for an
‘unsorted bam’ (or namesorted) file type in the file which describe
the format file uploaded ?
Or do you think that I can run this sorting during the processing of
the tool which require this type of sorting : so run the tool of
interest and internally to the process run the sorting by name before
treating the bam file. In this way do you think that Galaxy will not
create an index during the process ?
2015-07-17 21:51 GMT+02:00 Frederik Coppens <frcop(a)psb.vib-ugent.be>:
We also encountered this issue.
As I understood it: if Galaxy gets a BAM file (either through upload or
generated by a tool), it will generate metadata and indexing is part of
that. This only works if the BAM file is coordinate sorted, otherwise it
fails (although there is nothing wrong with the BAM file).
As a lot of tools (and visualisation) require sorted BAM files, this avoids
a lot of problems for users (but does create this issue).
An ‘unsorted bam’ (or namesorted) file type could be a solution although
this would cause the illogic situation where a ‘bam' is hierarchically not
above ‘unsorted bam’ as is the case for e.g. the types of fastq. I guess
changing what is now ‘bam’ to ‘sorted bam’ will cause problems for a lot of
On 17 Jul 2015, at 21:29, Alexander Vowinkel <vowinkel.alexander(a)gmail.com>
Indeed, this would be the right approach.
But so far the current strategy is quite imputed also in tools.
E.g. BWA from devteam automatically sorts the bam afterwards - wanted or
2015-07-17 11:59 GMT-05:00 Ryan G <ngsbioinformatics(a)gmail.com>:
> I would think the correct approach would be to allow any bam file and if a
> tool requires a coordinate sorted bam and one is not provided it would make
> sense for the tool to fail with an error.
> Does galaxy itself require a coordinate sorted bam or this this so the
> majority of tools can function properly?
> Sent from my iPhone
> > On Jul 17, 2015, at 11:08 AM, John Chilton <jmchilton(a)gmail.com> wrote:
> > I'm not the right person to respond to this but since no one else has
> > I will explain my limited (probably incorrect) understanding of the
> > problem and what needs to be done. I believe Galaxy assumes all bam
> > files are coordinated sorted by default - so Galaxy's bam datatype
> > would be better thought of as a coordinate-sorted bam. Certain people
> > do not want to allow non-coordinate sorted bam files because these
> > would be allowed by Galaxy as valid as inputs to tools that expect
> > sorted bams - and things would break. For the same reason I imagine
> > these people would argue the bai file needs to be generated because
> > certain tools will depend on its existence.
> > I believe the proposed solution is to add an unsorted bam datatype to
> > Galaxy. I don't know if there is any work in progress on this but it
> > seems to be needed by many researchers.
> > Hope this clarifies things somewhat,
> > -John
> >> On Mon, Jun 29, 2015 at 11:06 AM, julie dubois <dubjulie(a)gmail.com>
> >> wrote:
> >> Hi all,
> >> My galaxy is the last update and I work with samtools 1.2.
> >> I have a bug with the sorting and the upload of bam file.
> >> This bug appears only when I try to upload bam sorted by name or when
> >> I try to sort by name.
> >> When I work with bam sorted by coordinate, there is no problem.
> >> So, when I try to sort a bam file by name, this error appears under
> >> the "bug icon" :
> >> Tool execution generated the following error message:
> >> [bam_sort_core] merging from 20 files...
> >> *** Error in `python': double free or corruption (!prev):
> >> 0x0000000002897990 ***
> >> Aborted (core dumped)
> >> The tool produced the following additional output:
> >> [E::hts_idx_push] chromosome blocks not continuous
> >> Galaxy attempted to build the BAM index with samtools 1.0+ but failed:
> >> [Errno 2] No such file or directory:
> >> And the same type of error appears when I try to upload bam sorted by
> >> name.
> >> What I understand is that the bam can be sorted properly (The command
> >> line to sort is ok in a terminal with no bug and I can follow this
> >> step in the job_working_directory during the run) but Galaxy seems to
> >> fail to create a bai.
> >> For our analysis pipeline, I think the bai is not crucial, is it
> >> possible to avoid this creation of bai during uploading and sorting of
> >> bam file by name ?
> >> Is there another problem than this creation of bai ?
> >> Thanks.
> >> Julie
> >> ___________________________________________________________
> >> Please keep all replies on the list by using "reply all"
> >> in your mail client. To manage your subscriptions to this
> >> and other Galaxy lists, please use the interface at:
> >> https://lists.galaxyproject.org/
> >> To search Galaxy mailing lists use the unified search at:
> >> http://galaxyproject.org/search/mailinglists/
> > ___________________________________________________________
> > Please keep all replies on the list by using "reply all"
> > in your mail client. To manage your subscriptions to this
> > and other Galaxy lists, please use the interface at:
> > https://lists.galaxyproject.org/
> > To search Galaxy mailing lists use the unified search at:
> > http://galaxyproject.org/search/mailinglists/
> Please keep all replies on the list by using "reply all"
> in your mail client. To manage your subscriptions to this
> and other Galaxy lists, please use the interface at:
> To search Galaxy mailing lists use the unified search at:
Please keep all replies on the list by using "reply all"
in your mail client. To manage your subscriptions to this
and other Galaxy lists, please use the interface at:
To search Galaxy mailing lists use the unified search at:
Frederik Coppens, Ph.D.
Project Leader Applied Bioinformatics & Biostatistics
Tel: 32 (0)9 3313692 Fax:32 (0)9 3313809
VIB DEPARTMENT OF PLANT SYSTEMS BIOLOGY, UGent
Technologiepark 927, 9052 Gent, BELGIUM