Question: Aligning Against Multiple Reference Sequences
1
gravatar for John David Osborne
7.5 years ago by
John David Osborne160 wrote:
Are there any tools in Galaxy to align short reads against multiple reference sequences? I have a dozen microbial genomes sequenced for which there are 2 reference genomes already sequenced. We have tried aligning each of these individually against either of the reference genomes - some align better against the first reference genome, some align better against the second reference genome. Ideally though I would like to be able to align against both at the same time. Is this possible? I have found a tool called GenomeMapper and hints of 2 other tools in development that do something like this, but nothing for Galaxy yet. How do others proceed with this type of problem? Workflows appreciated! :) -John
galaxy • 1.4k views
ADD COMMENTlink modified 7.5 years ago by David Matthews630 • written 7.5 years ago by John David Osborne160
0
gravatar for David Matthews
7.5 years ago by
United Kingdom
David Matthews630 wrote:
Hi John, Probably the simplest thing for you to do would be to concatenate the two genomes one after the other using the concatenate tool under "text manipulation". This will generate a new organism with apparently two chromosomes one from bacteria A and one from bacteria B. When you run tophat or bowtie the sam file will indicate which "chromosome" (i.e. which bacteria) it assigned the read to. Hope this helps. Best Wishes, David. __________________________________ Dr David A. Matthews Senior Lecturer in Virology Room E49 Department of Cellular and Molecular Medicine, School of Medical Sciences University Walk, University of Bristol Bristol. BS8 1TD U.K. Tel. +44 117 3312058 Fax. +44 117 3312091 D.A.Matthews@bristol.ac.uk
ADD COMMENTlink written 7.5 years ago by David Matthews630
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: 180 users visited in the last hour