This is a test version of Biostars. For the public version, visit https://www.biostars.org.
When to use Genome assembly and when to use short read alignment for consensus genration??

Let's say I have Illumina short-reads for a virus and its reference sequence is also available in the Database.

In this case what I should use to get a full-length genome from my Illumina short-reads?

  • Genome assembly
  • Or, short-read alignment for consensus generation.

Many times I get confused about when I should use which option.

Thanks in advance.

short-read assembly consensus alignment
• 2,280 views • link
updated by young_bioinformatician young_bioinformatician • written by Anisur Rahman • Anisur Rahman

I would say this boils down mostly to what you intend to do downstream, and how good the reference is.

If the reference is a complete genome with lots of validation, and you only care about variant identification, I'd say just mapping the reads is the way to go.

If the reference isn't that good, or you intend to do further analysis, there's no harm in assembling the genome.

Of course, you can do both pretty quickly.

• 0 views • link
written by Joe Joe

1 answer

If you want a full-length genome, you should do a de novo assembly. The reason is that if there are long indels or structural mutations between your genome and the reference, your reads may not align to those regions.

• 0 views • link
written by shelkmike shelkmike

long indels or structural mutations

Unless one is dealing with really large viruses those should not be applicable?

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this