Question: Aligning two genomes
1
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.

Help
Access

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