/infernal

Permalink
Browse files

o src/Makefile.in & configure.ac: modified to correctly handle the 

  --with-mpi option to configure.

o scripts/cmsearch_pp.pl: revamped to potentially parse cmsearch output with
  hits from multiple CMs against same target db. Options to remove overlaps.

o matrices/*: all .mat files were rewritten to include the background
  (g) distribution. see ~nawrockie/notebook/7_0201_inf_rsearch/00LOG
  for details.

o benchmarks/cmsearch-rmark/infernal2glbf.pl: can handle cmsearch
  output that includes hits from multiple CMs against same target.



git-svn-id: https://svn.janelia.org/eddylab/eddys/src/infernal/trunk@1925 db3e70e4-12e1-0310-b23f-81a6e3d09ba7
  • Loading branch information...
1 parent 9eb81a2 commit efd5750d23142e288c0d120a08d6d67c38d19a81 @nawrockie nawrockie committed Apr 5, 2007
Unified Split
Showing with 4,650 additions and 4,071 deletions.
  1. +78 −61 benchmarks/cmsearch-rmark/infernal2glbf.pl
  2. +33 −16 configure.ac
  3. +25 −22 matrices/RIBOSUM100-25.mat
  4. +25 −22 matrices/RIBOSUM100-30.mat
  5. +25 −22 matrices/RIBOSUM100-35.mat
  6. +25 −22 matrices/RIBOSUM100-40.mat
  7. +25 −22 matrices/RIBOSUM100-45.mat
  8. +29 −26 matrices/RIBOSUM100-50.mat
  9. +26 −23 matrices/RIBOSUM100-55.mat
  10. +28 −25 matrices/RIBOSUM100-60.mat
  11. +25 −22 matrices/RIBOSUM100-65.mat
  12. +25 −22 matrices/RIBOSUM100-70.mat
  13. +29 −26 matrices/RIBOSUM100-75.mat
  14. +29 −26 matrices/RIBOSUM100-80.mat
  15. +28 −25 matrices/RIBOSUM100-85.mat
  16. +25 −22 matrices/RIBOSUM100-90.mat
  17. +25 −22 matrices/RIBOSUM100-95.mat
  18. +25 −22 matrices/RIBOSUM100.mat
  19. +25 −22 matrices/RIBOSUM45-30.mat
  20. +25 −22 matrices/RIBOSUM45-35.mat
  21. +25 −22 matrices/RIBOSUM45-40.mat
  22. +25 −22 matrices/RIBOSUM45-45.mat
  23. +29 −26 matrices/RIBOSUM45-50.mat
  24. +26 −23 matrices/RIBOSUM45-55.mat
  25. +29 −26 matrices/RIBOSUM45-60.mat
  26. +25 −22 matrices/RIBOSUM45-65.mat
  27. +25 −22 matrices/RIBOSUM45-70.mat
  28. +28 −25 matrices/RIBOSUM45-75.mat
  29. +29 −26 matrices/RIBOSUM45-80.mat
  30. +25 −22 matrices/RIBOSUM45.mat
  31. +27 −24 matrices/RIBOSUM55-25.mat
  32. +27 −24 matrices/RIBOSUM55-30.mat
  33. +25 −22 matrices/RIBOSUM55-35.mat
  34. +25 −22 matrices/RIBOSUM55-40.mat
  35. +25 −22 matrices/RIBOSUM55-45.mat
  36. +29 −26 matrices/RIBOSUM55-50.mat
  37. +25 −22 matrices/RIBOSUM55-55.mat
  38. +28 −25 matrices/RIBOSUM55-60.mat
  39. +25 −22 matrices/RIBOSUM55-65.mat
  40. +25 −22 matrices/RIBOSUM55-70.mat
  41. +29 −26 matrices/RIBOSUM55-75.mat
  42. +29 −26 matrices/RIBOSUM55-80.mat
  43. +29 −26 matrices/RIBOSUM55-85.mat
  44. +25 −22 matrices/RIBOSUM55-90.mat
  45. +25 −22 matrices/RIBOSUM55-95.mat
  46. +25 −22 matrices/RIBOSUM55.mat
  47. +27 −24 matrices/RIBOSUM60-25.mat
  48. +26 −23 matrices/RIBOSUM60-30.mat
  49. +25 −22 matrices/RIBOSUM60-35.mat
  50. +27 −24 matrices/RIBOSUM60-40.mat
  51. +25 −22 matrices/RIBOSUM60-45.mat
  52. +29 −26 matrices/RIBOSUM60-50.mat
  53. +29 −26 matrices/RIBOSUM60-55.mat
  54. +28 −25 matrices/RIBOSUM60-60.mat
  55. +25 −22 matrices/RIBOSUM60-65.mat
  56. +25 −22 matrices/RIBOSUM60-70.mat
  57. +29 −26 matrices/RIBOSUM60-75.mat
  58. +29 −26 matrices/RIBOSUM60-80.mat
  59. +28 −25 matrices/RIBOSUM60-85.mat
  60. +25 −22 matrices/RIBOSUM60-90.mat
  61. +25 −22 matrices/RIBOSUM60-95.mat
  62. +25 −22 matrices/RIBOSUM60.mat
  63. +26 −23 matrices/RIBOSUM65-25.mat
  64. +25 −22 matrices/RIBOSUM65-30.mat
  65. +25 −22 matrices/RIBOSUM65-35.mat
  66. +27 −24 matrices/RIBOSUM65-40.mat
  67. +25 −22 matrices/RIBOSUM65-45.mat
  68. +29 −26 matrices/RIBOSUM65-50.mat
  69. +29 −26 matrices/RIBOSUM65-55.mat
  70. +29 −26 matrices/RIBOSUM65-60.mat
  71. +25 −22 matrices/RIBOSUM65-65.mat
  72. +25 −22 matrices/RIBOSUM65-70.mat
  73. +29 −26 matrices/RIBOSUM65-75.mat
  74. +29 −26 matrices/RIBOSUM65-80.mat
  75. +29 −26 matrices/RIBOSUM65-85.mat
  76. +25 −22 matrices/RIBOSUM65-90.mat
  77. +25 −22 matrices/RIBOSUM65-95.mat
  78. +25 −22 matrices/RIBOSUM65.mat
  79. +25 −22 matrices/RIBOSUM70-25.mat
  80. +26 −23 matrices/RIBOSUM70-30.mat
  81. +25 −22 matrices/RIBOSUM70-35.mat
  82. +27 −24 matrices/RIBOSUM70-40.mat
  83. +25 −22 matrices/RIBOSUM70-45.mat
  84. +28 −25 matrices/RIBOSUM70-50.mat
  85. +27 −24 matrices/RIBOSUM70-55.mat
  86. +29 −26 matrices/RIBOSUM70-60.mat
  87. +25 −22 matrices/RIBOSUM70-65.mat
  88. +25 −22 matrices/RIBOSUM70-70.mat
  89. +29 −26 matrices/RIBOSUM70-75.mat
  90. +29 −26 matrices/RIBOSUM70-80.mat
  91. +29 −26 matrices/RIBOSUM70-85.mat
  92. +25 −22 matrices/RIBOSUM70-90.mat
  93. +25 −22 matrices/RIBOSUM70-95.mat
  94. +25 −22 matrices/RIBOSUM70.mat
  95. +25 −22 matrices/RIBOSUM75-25.mat
  96. +25 −22 matrices/RIBOSUM75-30.mat
  97. +25 −22 matrices/RIBOSUM75-35.mat
  98. +26 −23 matrices/RIBOSUM75-40.mat
  99. +25 −22 matrices/RIBOSUM75-45.mat
  100. +28 −25 matrices/RIBOSUM75-50.mat
  101. +28 −25 matrices/RIBOSUM75-55.mat
  102. +29 −26 matrices/RIBOSUM75-60.mat
  103. +25 −22 matrices/RIBOSUM75-65.mat
  104. +25 −22 matrices/RIBOSUM75-70.mat
  105. +28 −25 matrices/RIBOSUM75-75.mat
  106. +29 −26 matrices/RIBOSUM75-80.mat
  107. +28 −25 matrices/RIBOSUM75-85.mat
  108. +25 −22 matrices/RIBOSUM75-90.mat
  109. +25 −22 matrices/RIBOSUM75-95.mat
  110. +25 −22 matrices/RIBOSUM75.mat
  111. +25 −22 matrices/RIBOSUM80-25.mat
  112. +25 −22 matrices/RIBOSUM80-30.mat
  113. +25 −22 matrices/RIBOSUM80-35.mat
  114. +25 −22 matrices/RIBOSUM80-40.mat
  115. +25 −22 matrices/RIBOSUM80-45.mat
  116. +29 −26 matrices/RIBOSUM80-50.mat
  117. +27 −24 matrices/RIBOSUM80-55.mat
  118. +29 −26 matrices/RIBOSUM80-60.mat
  119. +25 −22 matrices/RIBOSUM80-65.mat
  120. +25 −22 matrices/RIBOSUM80-70.mat
  121. +29 −26 matrices/RIBOSUM80-75.mat
  122. +29 −26 matrices/RIBOSUM80-80.mat
  123. +28 −25 matrices/RIBOSUM80-85.mat
  124. +25 −22 matrices/RIBOSUM80-90.mat
  125. +25 −22 matrices/RIBOSUM80-95.mat
  126. +25 −22 matrices/RIBOSUM80.mat
  127. +25 −22 matrices/RIBOSUM85-25.mat
  128. +26 −23 matrices/RIBOSUM85-30.mat
  129. +25 −22 matrices/RIBOSUM85-35.mat
  130. +26 −23 matrices/RIBOSUM85-40.mat
  131. +25 −22 matrices/RIBOSUM85-45.mat
  132. +29 −26 matrices/RIBOSUM85-50.mat
  133. +27 −24 matrices/RIBOSUM85-55.mat
  134. +28 −25 matrices/RIBOSUM85-60.mat
  135. +25 −22 matrices/RIBOSUM85-65.mat
  136. +25 −22 matrices/RIBOSUM85-70.mat
  137. +29 −26 matrices/RIBOSUM85-75.mat
  138. +29 −26 matrices/RIBOSUM85-80.mat
  139. +28 −25 matrices/RIBOSUM85-85.mat
  140. +25 −22 matrices/RIBOSUM85-90.mat
  141. +25 −22 matrices/RIBOSUM85-95.mat
  142. +25 −22 matrices/RIBOSUM85.mat
  143. +25 −22 matrices/RIBOSUM90-30.mat
  144. +25 −22 matrices/RIBOSUM90-35.mat
  145. +25 −22 matrices/RIBOSUM90-40.mat
  146. +25 −22 matrices/RIBOSUM90-45.mat
  147. +29 −26 matrices/RIBOSUM90-50.mat
  148. +26 −23 matrices/RIBOSUM90-55.mat
  149. +29 −26 matrices/RIBOSUM90-60.mat
  150. +25 −22 matrices/RIBOSUM90-65.mat
  151. +25 −22 matrices/RIBOSUM90-70.mat
  152. +28 −25 matrices/RIBOSUM90-75.mat
  153. +29 −26 matrices/RIBOSUM90-80.mat
  154. +29 −26 matrices/RIBOSUM90-85.mat
  155. +25 −22 matrices/RIBOSUM90-90.mat
  156. +25 −22 matrices/RIBOSUM90-95.mat
  157. +25 −22 matrices/RIBOSUM90.mat
  158. +25 −22 matrices/RIBOSUM95-30.mat
  159. +25 −22 matrices/RIBOSUM95-35.mat
  160. +26 −23 matrices/RIBOSUM95-40.mat
  161. +25 −22 matrices/RIBOSUM95-45.mat
  162. +29 −26 matrices/RIBOSUM95-50.mat
  163. +26 −23 matrices/RIBOSUM95-55.mat
  164. +28 −25 matrices/RIBOSUM95-60.mat
  165. +25 −22 matrices/RIBOSUM95-65.mat
  166. +25 −22 matrices/RIBOSUM95-70.mat
  167. +28 −25 matrices/RIBOSUM95-75.mat
  168. +29 −26 matrices/RIBOSUM95-80.mat
  169. +29 −26 matrices/RIBOSUM95-85.mat
  170. +25 −22 matrices/RIBOSUM95-90.mat
  171. +25 −22 matrices/RIBOSUM95-95.mat
  172. +25 −22 matrices/RIBOSUM95.mat
  173. +46 −12 scripts/cmsearch_pp.pl
  174. +3 −2 src/Makefile.in
View
139 benchmarks/cmsearch-rmark/infernal2glbf.pl
@@ -40,96 +40,113 @@
&infernal::ParseINFERNAL($output);
# First determine if infernal was run with or without E-values
-if($infernal::nhit > 0)
+$at_least_one_hit = 0;
+for($c = 0; $c < $infernal::ncm; $c++)
{
- if (exists($infernal::hitevalue[0]))
+ if($infernal::nhit[$c] > 0)
{
- $has_evalues = 1;
- }
- else
- {
- $has_evalues = 0;
+ $at_least_one_hit = 1;
+ if (exists($infernal::hitevalue[$c][0]))
+ {
+ $has_evalues = 1;
+ }
+ else
+ {
+ $has_evalues = 0;
+ }
+ last;
}
}
-
-if(($infernal::nhit > 0) && ($use_evalues && (!$has_evalues)))
+if(!($at_least_one_hit))
+{
+ die("No hits found. Exiting.\n");
+}
+if($use_evalues && (!$has_evalues))
{
- die("ERROR in infernal2glbf.pl, trying to use E-values but none reported.\n");
+ die("ERROR, trying to use E-values but none reported.\n");
}
# if we're not sorting the scores, print them out in the order cmsearch reported them
if(!($sort_scores))
{
- for ($i = 0; $i < $infernal::nhit; $i++)
+ for ($c = 0; $c < $infernal::ncm; $c++)
{
- if ((($use_bitscores) && $infernal::hitbitscore[$i] > $b_cutoff) ||
- (($use_evalues) && $infernal::hitevalue[$i] < $e_cutoff))
+ for ($i = 0; $i < $infernal::nhit[$c]; $i++)
{
- #printf("%-24s %-6f\n", $infernal::targname[$i], $infernal::seqbitscore{$infernal::targname[$i]});
- if($infernal::hitsqfrom[$i] > $infernal::hitsqto[$i])
+ if ((($use_bitscores) && $infernal::hitbitscore[$c][$i] > $b_cutoff) ||
+ (($use_evalues) && $infernal::hitevalue[$c][$i] < $e_cutoff))
{
- #hit to reverse strand of query
- $orient = 1;
- }
- else
- {
- $orient = 0;
+ #printf("%-24s %-6f\n", $infernal::targname[$i], $infernal::seqbitscore{$infernal::targname[$i]});
+ if($infernal::hitsqfrom[$c][$i] > $infernal::hitsqto[$c][$i])
+ {
+ #hit to reverse strand of query
+ $orient = 1;
+ }
+ else
+ {
+ $orient = 0;
+ }
+ if($use_evalues)
+ {
+ printf("%-24s %-6f %d %d %d\n", $infernal::targname_byhit[$c][$i], $infernal::hitevalue[$c][$i], $infernal::hitsqfrom[$c][$i], $infernal::hitsqto[$c][$i], $orient);
+ }
+ else
+ {
+ printf("%-24s %-6f %d %d %d\n", $infernal::targname_byhit[$c][$i], $infernal::hitbitscore[$c][$i], $infernal::hitsqfrom[$c][$i], $infernal::hitsqto[$c][$i], $orient);
+ }
}
+ }
+ }
+}
+else # sort scores
+{
+ for ($c = 0; $c < $infernal::ncm; $c++)
+ {
+ for ($i = 0; $i < $infernal::nhit; $i++)
+ {
+ $key = "$c:$i";
if($use_evalues)
{
- printf("%-24s %-6f %d %d %d\n", $infernal::targname_byhit[$i], $infernal::hitevalue[$i], $infernal::hitsqfrom[$i], $infernal::hitsqto[$i], $orient);
+ $sc_H{$key} = $infernal::hitevalue[$c][$i];
}
else
{
- printf("%-24s %-6f %d %d %d\n", $infernal::targname_byhit[$i], $infernal::hitbitscore[$i], $infernal::hitsqfrom[$i], $infernal::hitsqto[$i], $orient);
+ $sc_H{$key} = $infernal::hitbitscore[$c][$i];
}
}
- }
-}
-else # sort scores
-{
- for ($i = 0; $i < $infernal::nhit; $i++)
- {
if($use_evalues)
{
- $sc_H{$i} = $infernal::hitevalue[$i];
+ @sorted_i_A = sort { $sc_H{$a} <=> $sc_H{$b} } (keys (%sc_H));
}
else
{
- $sc_H{$i} = $infernal::hitbitscore[$i];
+ @sorted_i_A = sort { $sc_H{$b} <=> $sc_H{$a} } (keys (%sc_H));
}
- }
- if($use_evalues)
- {
- @sorted_i_A = sort { $sc_H{$a} <=> $sc_H{$b} } (keys (%sc_H));
- }
- else
- {
- @sorted_i_A = sort { $sc_H{$b} <=> $sc_H{$a} } (keys (%sc_H));
- }
- for ($j = 0; $j < scalar(@sorted_i_A); $j++)
- {
- $i = $sorted_i_A[$j];
- if ((($use_bitscores) && $infernal::hitbitscore[$i] > $b_cutoff) ||
- (($use_evalues) && $infernal::hitevalue[$i] < $e_cutoff))
+ for ($j = 0; $j < scalar(@sorted_i_A); $j++)
{
- #printf("%-24s %-6f\n", $infernal::targname[$i], $infernal::seqbitscore{$infernal::targname[$i]});
- if($infernal::hitsqfrom[$i] > $infernal::hitsqto[$i])
- {
- #hit to reverse strand of query
- $orient = 1;
- }
- else
- {
- $orient = 0;
- }
- if($use_evalues)
- {
- printf("%-24s %-6f %d %d %d\n", $infernal::targname_byhit[$i], $infernal::hitevalue[$i], $infernal::hitsqfrom[$i], $infernal::hitsqto[$i], $orient);
- }
- else
+ $key = $sorted_i_A[$j];
+ ($c, $i) = split(":", $key);
+ if ((($use_bitscores) && $infernal::hitbitscore[$c][$i] > $b_cutoff) ||
+ (($use_evalues) && $infernal::hitevalue[$c][$i] < $e_cutoff))
{
- printf("%-24s %-6f %d %d %d\n", $infernal::targname_byhit[$i], $infernal::hitbitscore[$i], $infernal::hitsqfrom[$i], $infernal::hitsqto[$i], $orient);
+ #printf("%-24s %-6f\n", $infernal::targname[$c][$i], $infernal::seqbitscore{$infernal::targname[$c][$i]});
+ if($infernal::hitsqfrom[$c][$i] > $infernal::hitsqto[$c][$i])
+ {
+ #hit to reverse strand of query
+ $orient = 1;
+ }
+ else
+ {
+ $orient = 0;
+ }
+ if($use_evalues)
+ {
+ printf("%-24s %-6f %d %d %d\n", $infernal::targname_byhit[$c][$i], $infernal::hitevalue[$c][$i], $infernal::hitsqfrom[$c][$i], $infernal::hitsqto[$c][$i], $orient);
+ }
+ else
+ {
+ printf("%-24s %-6f %d %d %d\n", $infernal::targname_byhit[$c][$i], $infernal::hitbitscore[$c][$i], $infernal::hitsqfrom[$c][$i], $infernal::hitsqto[$c][$i], $orient);
+ }
}
}
}
View
49 configure.ac
@@ -58,12 +58,43 @@ AC_DEFINE_UNQUOTED(PACKAGE_LICENSE, "$PACKAGE_LICENSE")
if ((! test -d binaries) && test -d src); then
AC_MSG_NOTICE([This is a full source distribution - using full config])
+cwd=`pwd`
+AC_SUBST(TOPDIR,$cwd)
+# EPN, Thu Apr 5 07:58:07 2007
+# First, check for --with-mpi=<x> to build MPI cmsearch and cmalign progs.
+# if enabled, set compiler as <x>, if no "=<x>" provided, set compiler as mpicc
+# currently PROGS are not affected, though eventually we may want
+# to create 'mpicmsearch' and 'mpicmalign'
+usempi=0
+AC_ARG_WITH(mpi,
+[ --with-mpi enable MPI for search and alignment],
+[case $with_mpi in
+ yes) AC_MSG_NOTICE([configured cmalign and cmsearch for MPI])
+ mpi="mpicc"
+ usempi=1
+ AC_DEFINE(USE_MPI)
+ ;;
+ no) mpi=""
+ ;;
+ *) mpi=$with_mpi
+ usempi=1
+ ;;
+ esac])
+AC_SUBST(USEMPI,$usempi)
+
+# if $mpi == "", use default C compiler, set CFLAGS by default
+if test -z "$mpi"; then
+ AC_PROG_CC
+# if $mpi != "", use compiler $mpi, manually set CFLAGS
+else
+ AC_PROG_CC($mpi)
+ CFLAGS="-g -O2"
+fi
# Checks for programs.
-# Sets output variables CC, CFLAGS, and RANLIB.
+# Sets output variables RANLIB.
# AR might be hidden away on Solaris; provide some hints in the PATH.
#
-AC_PROG_CC
AC_PROG_CXX
AC_PROG_RANLIB
AC_PATH_PROG([AR], [ar], [:], [$PATH:/usr/ccs/bin:/usr/xpg4/bin])
@@ -250,20 +281,6 @@ AC_ARG_WITH(rigfilters,
esac])
#
-# --with-mpi MPI support for cmalign and cmsearch
-AC_ARG_WITH(mpi,
-[ --with-mpi enable MPI for search and alignment],
-[case $with_mpi in
- yes) AC_MSG_NOTICE([configured cmalign and cmsearch for MPI])
- AC_PROG_CC(mpicc)
- ;;
- no) ;;
- *) echo "Ignoring unknown argument to --with-mpi: $with_mpi"
- ;;
- esac])
-
-
-#
# end of optionally configured features
################################################################
View
47 matrices/RIBOSUM100-25.mat
@@ -1,30 +1,33 @@
RIBOSUM100-25
- A C G U
-A 1.97
-C -1.34 1.00
-G -0.92 -1.70 0.85
-U -1.01 -0.61 -1.11 1.38
+ A C G U
+ 0.258308 0.219504 0.297310 0.224878
+
+ A C G U
+A 1.971671
+C -1.343217 1.001064
+G -0.921788 -1.699044 0.854147
+U -1.005464 -0.613591 -1.111982 1.382532
H: 0.5617
E: -0.0957
- AA AC AG AU CA CC CG CU GA GC GG GU UA UC UG UU
-AA -2.60
-AC -8.47 -1.85
-AG -8.82 -9.59 -1.33
-AU -5.08 -1.65 -6.07 4.25
-CA -10.06 -10.55 -10.65 -5.77 -6.21
-CC -13.44 -8.85 -15.14 -2.37 -11.69 -3.46
-CG -3.87 -5.44 -2.57 2.03 -3.24 -5.07 5.22
-CU -12.11 -10.68 -9.36 -4.18 -8.33 -6.86 -4.39 -2.56
-GA -6.63 -7.17 -9.39 -3.79 -7.76 -12.90 -6.03 -7.14 -1.62
-GC -5.22 -2.70 -6.29 2.88 -5.40 -3.75 2.80 -2.20 -2.57 5.49
-GG -9.69 -10.38 -4.22 -4.93 -11.48 -12.06 -4.71 -9.14 -9.91 -4.77 -2.31
-GU -6.50 -4.92 -7.13 0.75 -6.20 -7.10 0.25 -4.41 -6.19 1.42 -5.53 3.32
-UA -2.49 -5.59 -2.82 1.88 -2.91 -6.47 2.93 -4.33 -4.87 1.89 -6.32 0.06 4.72
-UC -12.26 -8.88 -8.01 -4.54 -8.03 -8.95 -4.35 -4.67 -6.98 -4.16 -10.57 -4.61 -3.53 -4.14
-UG -3.55 -7.17 -5.83 -0.39 -6.21 -8.30 1.19 -6.47 -7.67 0.33 -4.34 -1.26 1.04 -4.20 3.19
-UU -8.87 -8.39 -9.76 -2.46 -8.85 -6.64 -3.01 -5.43 -10.15 -3.03 -5.02 -2.58 -1.58 -5.90 -3.64 -0.72
+ AA AC AG AU CA CC CG CU GA GC GG GU UA UC UG UU
+AA -2.604645
+AC -8.470563 -1.853570
+AG -8.822189 -9.594620 -1.327118
+AU -5.075545 -1.654552 -6.067763 4.252767
+CA -10.056671 -10.552723 -10.646974 -5.774248 -6.205123
+CC -13.442368 -8.854252 -15.136720 -2.374419 -11.693921 -3.456602
+CG -3.871417 -5.438003 -2.572185 2.030624 -3.240493 -5.067212 5.222069
+CU -12.105264 -10.682421 -9.357424 -4.179064 -8.325364 -6.862516 -4.386954 -2.559931
+GA -6.625421 -7.170250 -9.390651 -3.792122 -7.764439 -12.895778 -6.029141 -7.141504 -1.623429
+GC -5.219950 -2.695725 -6.290820 2.876080 -5.395295 -3.754027 2.801056 -2.204690 -2.569771 5.489061
+GG -9.691680 -10.383497 -4.223349 -4.930131 -11.479506 -12.063889 -4.707853 -9.136017 -9.906732 -4.774241 -2.305253
+GU -6.496390 -4.924699 -7.125769 0.745252 -6.201344 -7.102632 0.251253 -4.410597 -6.189448 1.424712 -5.534224 3.320225
+UA -2.485010 -5.590836 -2.822798 1.878540 -2.907659 -6.465319 2.926758 -4.327115 -4.869013 1.889057 -6.320650 0.064069 4.723641
+UC -12.263343 -8.877730 -8.011171 -4.540962 -8.030163 -8.946997 -4.345131 -4.669281 -6.980707 -4.161804 -10.568347 -4.611892 -3.529918 -4.138184
+UG -3.546296 -7.166265 -5.829573 -0.392702 -6.206307 -8.303615 1.189462 -6.470759 -7.672941 0.327564 -4.341610 -1.261343 1.043982 -4.199517 3.189382
+UU -8.867475 -8.385429 -9.762040 -2.456123 -8.850458 -6.641057 -3.008435 -5.429174 -10.145796 -3.034762 -5.019912 -2.575165 -1.581614 -5.904415 -3.644739 -0.715714
H: 3.4639
E: -2.6040
View
47 matrices/RIBOSUM100-30.mat
@@ -1,30 +1,33 @@
RIBOSUM100-30
- A C G U
-A 1.97
-C -1.34 1.00
-G -0.92 -1.70 0.85
-U -1.01 -0.61 -1.11 1.38
+ A C G U
+ 0.258284 0.219532 0.297342 0.224842
+
+ A C G U
+A 1.972360
+C -1.344047 1.001351
+G -0.922482 -1.699386 0.854425
+U -1.006582 -0.614272 -1.113000 1.383155
H: 0.5622
E: -0.0959
- AA AC AG AU CA CC CG CU GA GC GG GU UA UC UG UU
-AA -2.60
-AC -8.47 -1.85
-AG -8.82 -9.59 -1.33
-AU -5.07 -1.65 -6.07 4.25
-CA -10.06 -10.55 -10.65 -5.78 -6.20
-CC -13.44 -8.85 -15.14 -2.37 -11.69 -3.46
-CG -3.87 -5.44 -2.57 2.03 -3.24 -5.07 5.22
-CU -12.10 -10.68 -9.36 -4.18 -8.33 -6.86 -4.39 -2.56
-GA -6.62 -7.17 -9.39 -3.79 -7.76 -12.90 -6.03 -7.14 -1.62
-GC -5.22 -2.70 -6.29 2.88 -5.40 -3.75 2.80 -2.20 -2.57 5.49
-GG -9.69 -10.38 -4.22 -4.93 -11.48 -12.06 -4.71 -9.14 -9.91 -4.77 -2.31
-GU -6.50 -4.92 -7.13 0.75 -6.20 -7.10 0.25 -4.41 -6.19 1.42 -5.53 3.32
-UA -2.49 -5.60 -2.82 1.88 -2.91 -6.47 2.93 -4.33 -4.87 1.89 -6.32 0.06 4.72
-UC -12.26 -8.88 -8.01 -4.54 -8.03 -8.95 -4.35 -4.67 -6.98 -4.16 -10.57 -4.61 -3.53 -4.14
-UG -3.55 -7.17 -5.83 -0.39 -6.21 -8.30 1.19 -6.47 -7.67 0.33 -4.34 -1.26 1.04 -4.20 3.19
-UU -8.87 -8.38 -9.76 -2.46 -8.85 -6.64 -3.01 -5.43 -10.15 -3.03 -5.02 -2.57 -1.58 -5.90 -3.64 -0.71
+ AA AC AG AU CA CC CG CU GA GC GG GU UA UC UG UU
+AA -2.603894
+AC -8.470122 -1.853437
+AG -8.821726 -9.594466 -1.326941
+AU -5.074954 -1.654023 -6.067754 4.253602
+CA -10.056229 -10.552590 -10.646820 -5.775143 -6.204990
+CC -13.442235 -8.854428 -15.136875 -2.374198 -11.694097 -3.457087
+CG -3.871265 -5.438174 -2.572318 2.030071 -3.240654 -5.067683 5.221619
+CU -12.104720 -10.682187 -9.357168 -4.178441 -8.325129 -6.862591 -4.387006 -2.559595
+GA -6.624957 -7.170095 -9.390474 -3.791566 -7.764284 -12.895932 -6.029273 -7.141247 -1.623253
+GC -5.220724 -2.695882 -6.291214 2.875845 -5.395536 -3.754496 2.800486 -2.204757 -2.569904 5.488602
+GG -9.691504 -10.383629 -4.223460 -4.929956 -11.479638 -12.064330 -4.708279 -9.136048 -9.906843 -4.774661 -2.305651
+GU -6.495910 -4.924468 -7.125567 0.745141 -6.201176 -7.102715 0.250943 -4.410244 -6.189170 1.424606 -5.534233 3.320576
+UA -2.485164 -5.596511 -2.824172 1.878742 -2.907259 -6.465591 2.926251 -4.326510 -4.868630 1.888583 -6.320774 0.063396 4.724416
+UC -12.262800 -8.877496 -8.010915 -4.540369 -8.029929 -8.947071 -4.345193 -4.668945 -6.980451 -4.161876 -10.568378 -4.611534 -3.529380 -4.137848
+UG -3.545733 -7.166062 -5.829295 -0.393684 -6.206051 -8.303667 1.189222 -6.470420 -7.672718 0.327374 -4.341619 -1.261042 1.044086 -4.199159 3.189760
+UU -8.866521 -8.384784 -9.761374 -2.455107 -8.849813 -6.640721 -3.008077 -5.428428 -10.145129 -3.034405 -5.019533 -2.574439 -1.580581 -5.903669 -3.643971 -0.714557
H: 3.4640
E: -2.6041
Oops, something went wrong.

0 comments on commit efd5750

Please sign in to comment.