This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Large Scale SARS CoV 2 multiple Sequence Alignment

Dear members,

I want to align all the SARS-CoV-2 complete genome sequences that are available till date (444424 in NCBI viruses). I have tried command line muscle but it is unable to handle such a large number of sequences. I would appreciate if anyone can let us know how to align these huge number of sequences.

multiple sars-cov-2 msa sequence genome alignment
• 2,845 views • link
written by s_bio • s_bio

A large number of those assemblies are identical with one another or subsequences of longer assemblies. Quite a few also include degenerate bases which you would have to deal with somehow. My suggestion is that you first remove all redundancy from the data set by clustering..

• 0 views • link
written by 5heikki 5heikki

in addition, there may be a polyA tail at the end, and that may have various lengths, for many genomes the sequences end like this:

>OU534909.1 Severe acute respiratory syndrome coronavirus 2 genome assembly, complete genome: monopartite [29797:29896]
TAATGTGTAAAATTAATTTTAGTAGTGCTATCCNNNNNNNNNNNNNNNNNNNNNNNNNNN
NNNNNNNAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
• 0 views • link
written by Istvan Albert Istvan Albert

Thank you for your reply. We intend to do that but after the alignment process is over..

• 0 views • link
written by s_bio • s_bio

use high blosum matrix for alignment

• 0 views • link
written by cpad0112 cpad0112

Thank you for the suggestion. Will surely try high blosum..

• 0 views • link
written by s_bio • s_bio

1 answer

nextalign is able to handle large numbers of sequences and produces a robust result: https://docs.nextstrain.org/projects/nextclade/en/stable/user/nextalign-cli.html

Another strategy involves aligning each individual sequence against the SARS-CoV-2 reference genome one-by-one, then combining the results; Check out the methods at https://github.com/roblanf/sarscov2phylo

You can also reverse-engineer the method used by pangolin. The pangolin method maps genomes against the SARS-CoV-2 reference genome using minimap2, then converts the sam-format mapped file to a multifasta file using gofasta. Check out the 'align_to_reference' rule here: https://github.com/cov-lineages/pangolin/blob/master/pangolin/scripts/pangolearn.smk

• 0 views • link
written by dr_cfoster dr_cfoster

Log in to answer this question.

More posts like this