This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Align Then Visualize The Alignment Of Rna Sequences With R Or Python

Hi,

I have two RNA-sequences along with there bracket annotation ..

i have 2 questions:

  1. is it possible to align then inside R??
  2. if i already have the alignment outside R, is there anyway to present it in a colrful way such to show the matched pairs, the gaps, the non-matched pairs in different colors ..

i can align the sequences inside python, but i'm not sure if there is a good way to perform this inside python, so the question would be in R or python...

Regards,

r alignment
• 5,965 views • link
updated by Leandro Lima Leandro Lima • written by Abdullah Abdullah

3 answers

Hi,

Answer 1: i would call an external command to compute such alignment and recover the result from R.

Answer 2: ok, this is neither R nor Python, but what about Ralee (an emacs plugin):http://personalpages.manchester.ac.uk/staff/sam.griffiths-jones/software/ralee/ ?

• 17 views • link
written by Fabrice Jossinet Fabrice Jossinet

seems like a good idea ... but what i want is something more dynamic that can be included into a script, because i want to do this for a lot of alignments ...

• 0 views • link
written by Abdullah Abdullah

You can align them using MAFFT or PRANK calling them from a R/python/perl/younameit script. To view alignments and generate figures I particularly like jalview:http://www.jalview.org/ It has a lot of different view options.

• 32 views • link
written by Biojl Biojl

Maybe this can help/answer.

http://www.bioconductor.org/help/workflows/high-throughput-sequencing/

• 135 views • link
written by Leandro Lima Leandro Lima

Log in to answer this question.

More posts like this