Question: Aligning Against Multiple Reference Sequences
1
gravatar for John David Osborne
7.5 years ago by
John David Osborne • 160 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 COMMENT • link • modified 7.5 years ago by David Matthews • 630 • written 7.5 years ago by John David Osborne • 160
0
gravatar for David Matthews
7.5 years ago by
United Kingdom
David Matthews • 630 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 COMMENT • link written 7.5 years ago by David Matthews • 630
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