Loading multiple sequence alignment fasta file for ggmsa
0
0
Entering edit mode
18 months ago
snowpin • 0

Hello,

I am trying to use ggmsa to visualize my amino acid multiple sequence alignment (MSA) but the ggmsa() function doesn't recognize my MSA. Any thoughts how I can correct this?

  1. I have been loading my MSA with the following...

align <- read.alignment("my_protein_alignment.fasta",format = "fasta")

  1. And then use the following to create a figure...

ggmsa(align, 1, 100, color = "Clustal", font = "DroidSansMono", char_width = 0.5, seq_name = TRUE )

  1. However, I get this error...

Error in prepare_msa(msa) : multiple sequence alignment object no supported...

How can I fix this? Thanks in advance.

ggmsa • 948 views
ADD COMMENT
0
Entering edit mode

Hi! I guess that align variable contains the output of msa, right? Can you show the command you used for running msa ?

ADD REPLY
0
Entering edit mode

Yes, the align variable contains my msa, which I made with mafft (mafft --auto input_sequences.fasta).

ADD REPLY

Login before adding your answer.

Traffic: 1853 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6