Hi Romain,

If it is just the conservation score you are interested in, then other software 
can give this to you too.  You can get this by uploading your alignment to the 
AMAS server or running AMAS locally, alternatively  ALSCRIPT also has a 
function to generate the conservation score.


Still, Mungo's solution does it nicely in Jalview, so probably best to stick 
with that!


Geoff.

On 27/11/2018 13:33, Mungo Carstairs (Staff) wrote:

Hi Romain,


Try running the attached script in the Jalview Groovy console - it should 
output residue Blosum62 scores against the alignment consensus.


Regards,


Mungo





[University of Dundee shield logo]<http://uod.ac.uk/sig-home>

Mungo Carstairs
Jalview Computational Scientist

The Barton Group
Division of Computational Biology

School of Life Sciences

University of Dundee, Dundee, Scotland, UK

www.jalview.org<http://www.jalview.org>

www.compbio.dundee.ac.uk<http://www.compbio.dundee.ac.uk>
[email protected]<mailto:[email protected]>



[University of                            Dundee 
Facebook]<http://uod.ac.uk/sig-fb> [University of                            
Dundee Twitter] <http://uod.ac.uk/sig-tw>  [University of                       
     Dundee LinkedIn] <http://uod.ac.uk/sig-li>  [University of                 
           Dundee YouTube] <http://uod.ac.uk/sig-yt>  [University of            
                Dundee Instagram] <http://uod.ac.uk/sig-ig>  [University of     
                       Dundee Snapchat] <http://uod.ac.uk/sig-sc>
We're Scottish University of the Year again!<http://uod.ac.uk/sig-strapline>
The Times / Sunday Times Good University Guide 2016 and 2017
________________________________
From: 
[email protected]<mailto:[email protected]> 
<[email protected]><mailto:[email protected]>
 on behalf of Romain Studer 
<[email protected]><mailto:[email protected]>
Sent: 27 November 2018 11:35:24
To: [email protected]<mailto:[email protected]>
Subject: [Jalview-discuss] Retrieve score for conservation at the individual 
amino acid level

Hi,

I am wondering if we can extract the score to do the colour map from an 
alignment?

For example, if we colour the alignment with Blosum62 (for example) and we 
select Colour by conservation, with a threshold of 30,
we will get column that are white, and some columns with various level of 
blue/purple.
And some outlier AA in those column will be kept still white or less bright. I 
guess this is based from the AMAS method?

Can we retrieve this information directly from Jalview and output the scores as 
a 2d matrix (genes x positions)?
Or is there a tool to compute such conservation score at the individual amino 
acid level?

Thanks,
Romain

The University of Dundee is a registered Scottish Charity, No: SC015096


_______________________________________________
Jalview-discuss mailing list
[email protected]<mailto:[email protected]>
http://www.compbio.dundee.ac.uk/mailman/listinfo/jalview-discuss



--
Geoff Barton | Professor of Bioinformatics | Head of Division of Computational 
Biology
School of Life Sciences | University of Dundee, Scotland, UK | 
[email protected]<mailto:[email protected]>
Tel: +44 1382 385860 | 
www.compbio.dundee.ac.uk<http://www.compbio.dundee.ac.uk> | twitter: @gjbarton


The University of Dundee is registered Scottish charity: No.SC015096


The University of Dundee is a registered Scottish Charity, No: SC015096
_______________________________________________
Jalview-discuss mailing list
[email protected]
http://www.compbio.dundee.ac.uk/mailman/listinfo/jalview-discuss

Reply via email to