This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Generating a consensus sequence of a large number of amino acid sequences in python

Hi,

I want to generate a consensus sequence of a large number of amino acid sequences using python.

For example,

If there are three different amino acid sequences:

   MAHEDYIG ... HLY
   MAAEDYLP ... HLY
   MAHEDYIG ... TLY

Then MAHEDYIG ... HLY has to be generated as a consensus sequence.

Could anyone create python code for this algorithm? I want to put FASTA file, which is aligned with Multiple Sequence Alignment tool, in the code and want to put a specific weight for each conserved residue.

consensus msa sequence

Yes! I've used MUSCLE/Clustalw/Clustalomega packages for running MSA but I am struggling with generating consensus sequence since I cannot find any packages for it in python.

0 answers

No answers yet.

Log in to answer this question.