Question: How To Use Bwa For Multiple Reads Together With Reference Genome:
0
gravatar for HG
5.3 years ago by
HG1.1k
Germany
HG1.1k wrote:

HI,i have 50 ecoli whole genome sequence reads. Now i want to align whole the data with reference genome. Can any one suggest some idea how to do it using bwa??

bwa • 3.0k views
ADD COMMENTlink modified 5.3 years ago by Devon Ryan89k • written 5.3 years ago by HG1.1k
2

Please don't shout (all caps). Please read the bwa manual.

ADD REPLYlink written 5.3 years ago by Michael Dondrup46k

Sorry for caps ....but i already gone through bwa manual . I find i can upload single genome but if i would like to multiple genome together then how i will do ?

ADD REPLYlink written 5.3 years ago by HG1.1k
1

sorry, I don't get what you want to do, and what have you tried? what do you mean by 'upload single genome'? what do you mean by 'to multiple genome together'?

ADD REPLYlink written 5.3 years ago by Michael Dondrup46k
1

In case you have 50 different strains or genomes of Ecoli and you wish to align your reads against them, then you need to do it individually. I mean first align your reads against first genome. Then second and so on so forth.

ADD REPLYlink written 5.3 years ago by Ashutosh Pandey11k

Or it is possible to use tools like GSNAP which take into consideration of minor allele frequency while aligning. One of the colleague in my lab have tried to perform similar analysis. He first generated the gtf containing the difference between different strain of the same virus, then he aligned the reads using GSNAP. This allow him to use one single reference to align all our samples. The problem with this approach however is that it will require some more downstream processing to extract the strains information.

ADD REPLYlink written 5.3 years ago by Sam2.3k

As Michael mentioned, your question, as posed, is ambiguous. Do you want to align reads against a genome, concatenate multiple genomes together and align reads against that, align multiple genomes to each other, or something completely different? The easiest thing would be to just tell us what your biological goal is and what you've done experimentally so far. Odds are good that we'll know better than you how to proceed given that information.

ADD REPLYlink modified 5.3 years ago • written 5.3 years ago by Devon Ryan89k
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 2.3.0
Traffic: 1679 users visited in the last hour