git://source.jalview.org
/
jalview.git
/ commitdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
| commitdiff |
tree
raw
|
patch
|
inline
| side by side (parent:
51a7dcc
)
gracefully cope with non-standard nucleotide or amino acid symbols in pairwise alignment.
author
jprocter
<Jim Procter>
Mon, 2 Nov 2009 13:25:48 +0000
(13:25 +0000)
committer
jprocter
<Jim Procter>
Mon, 2 Nov 2009 13:25:48 +0000
(13:25 +0000)
src/jalview/analysis/AlignSeq.java
patch
|
blob
|
history
diff --git
a/src/jalview/analysis/AlignSeq.java
b/src/jalview/analysis/AlignSeq.java
index
49be7f9
..
999dd07
100755
(executable)
--- a/
src/jalview/analysis/AlignSeq.java
+++ b/
src/jalview/analysis/AlignSeq.java
@@
-113,6
+113,8
@@
public class AlignSeq
String type;
String type;
+ private int[] charToInt;
+
/**
* Creates a new AlignSeq object.
*
/**
* Creates a new AlignSeq object.
*
@@
-380,11
+382,13
@@
public class AlignSeq
if (type.equals(AlignSeq.PEP))
{
intToStr = pep;
if (type.equals(AlignSeq.PEP))
{
intToStr = pep;
+ charToInt = ResidueProperties.aaIndex;
defInt = 23;
}
else if (type.equals(AlignSeq.DNA))
{
intToStr = dna;
defInt = 23;
}
else if (type.equals(AlignSeq.DNA))
{
intToStr = dna;
+ charToInt = ResidueProperties.nucleotideIndex;
defInt = 4;
}
else
defInt = 4;
}
else
@@
-821,7
+825,6
@@
public class AlignSeq
public int[] stringToInt(String s, String type)
{
int[] seq1 = new int[s.length()];
public int[] stringToInt(String s, String type)
{
int[] seq1 = new int[s.length()];
-
for (int i = 0; i < s.length(); i++)
{
// String ss = s.substring(i, i + 1).toUpperCase();
for (int i = 0; i < s.length(); i++)
{
// String ss = s.substring(i, i + 1).toUpperCase();
@@
-834,33
+837,15
@@
public class AlignSeq
try
{
try
{
- if (type.equals("pep"))
- {
- seq1[i] = ResidueProperties.aaIndex[c];
- if (seq1[i] > 23)
- {
- seq1[i] = 23;
- }
- }
- else if (type.equals("dna"))
+ seq1[i] = charToInt[c]; // set accordingly from setType
+ if (seq1[i]<0 || seq1[i] > defInt) // set from setType: 23 for peptides, or 4 for NA.
{
{
- seq1[i] = ResidueProperties.nucleotideIndex[c];
- if (seq1[i] > 4)
- {
- seq1[i] = 4;
- }
+ seq1[i] = defInt;
}
} catch (Exception e)
{
}
} catch (Exception e)
{
- if (type.equals("dna"))
- {
- seq1[i] = 4;
- }
- else
- {
- seq1[i] = 23;
- }
+ seq1[i] = defInt;
}
}
}
}