From 9088babd6a23a32508f4d6754dc48174bd2dd77b Mon Sep 17 00:00:00 2001 From: Cosmin SAVEANU <cosmin.saveanu@pasteur.fr> Date: Mon, 25 Feb 2019 17:04:57 +0100 Subject: [PATCH] Update README.md --- README.md | 19 ++++++++++++++++++- 1 file changed, 18 insertions(+), 1 deletion(-) diff --git a/README.md b/README.md index 135ff36..0eb3e6a 100644 --- a/README.md +++ b/README.md @@ -3,4 +3,21 @@ Two scripts to isolate two sequences from a multiple alignment and show them side by side. The result can be useful to: 1. find corresponding residue numbers between the aligned sequences; -2. compress information about similarity in a graph. \ No newline at end of file +2. compress information about similarity in a graph. + +Example of use: + +``` +python multiplalign_to_matrix.py -a testdata/DCP2_MAFFT.mfa -u "YEASX" -w "SCHPM" -o testdata/DCP2_MAFFT_matrix.txt +python matrix_to_graph_alignment.py -m testdata/DCP2_MAFFT_matrix.txt -o testdata/Scer_vs_Spom_DCP2.pdf +``` +From a multiple alignment (in jalview): + +<!--  --> +<img src="images/screenshot_jalview_dcp2.png" alt="jalview region of alignment" width="600"> + +one can get the following image, with only two of the sequences highlighted, lower is S. cerevisiae, upper is S. pombe: + + + +More information in the pymultialign_description.md file. \ No newline at end of file -- GitLab