Question: Aligning two genomes
gravatar for Nicolas Carter
4 months ago by
Nicolas Carter10 wrote:

Hello everyone, I'm trying to align two sequences with EMBOSS needle (global alignment). Specifically, they are the genomes of Sars-CoV-2, but one of the sequences (MT259240) doesn't have certain regions (the non-coding regions like 5'UTR and 3'UTR), so, I suspect that the alignment may be not quite right. Namely, the first non-coding region of the first sequence (NC_045512 the referent sequence) is being aligned with the first gene of the other sequence (ORF1ab gene). Is this incorrect and if so how would I correct the mistake? Thanks in advance.

alignment genome • 127 views
ADD COMMENTlink written 4 months ago by Nicolas Carter10

I don't think needle is up for the task. Try minimap2 or mummer

ADD REPLYlink written 4 months ago by Asaf8.4k

What format is your data in, fasta? You would need an aligner like LASTZ in addition to the two options mentioned by @Asaf.

ADD REPLYlink written 4 months ago by genomax92k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1565 users visited in the last hour