正在寻找java "A to T, T to A"找到这个建议:
String sequence = "AATTTCTCGGTTTCAAT";
sequence = sequence.replace("A", "t")
.replace("T", "a")
.replace("C", "g")
.replace("G", "c")
.toUpperCase();
System.out.println(sequence);
Run Code Online (Sandbox Code Playgroud)
这是一个简单而简洁的解决方案,适用于您的特定情况,如果您的DNA字符串相对较短,将具有可接受的性能.对于处理大量数据的更通用的解决方案,您应该逐个迭代字符并构建新字符串.或者作为多基因润滑剂指出 - 考虑一种存储格式,每个基数仅使用2位而不是16位.
小智 6
我会寻求一个更通用的解决方案:
public String tr(String original, String trFrom, String trTo) {
StringBuilder sb = new StringBuilder();
for (int i = 0; i < original.length(); ++i) {
int charIndex = trFrom.indexOf(original.charAt(i));
if (charIndex >= 0) {
sb.append(trTo.charAt(charIndex));
} else {
sb.append(original.charAt(i));
}
}
return sb.toString();
}
Run Code Online (Sandbox Code Playgroud)
调用这样的函数会得到你需要的结果:
tr("ACGTA", "ATCG", "TAGC")
Run Code Online (Sandbox Code Playgroud)
所以该函数与unix tr实用程序几乎相同:
echo ACGTA | tr ATCG TAGC
Run Code Online (Sandbox Code Playgroud)
就像我昨天解释的那样,字符串是不可变的,您无法更改字符串,您必须创建一个新字符串并替换旧字符串。
您可以这样解决您的问题:
String s = "ACGTA";
StringBuilder sb= new StringBuilder();
for (char c:s.toCharArray()) {
switch(c) {
case 'A': sb.append('T');break;
case 'T': sb.append('A');break;
case 'C': sb.append('G');break;
case 'G': sb.append('C');break;
default: //handle error here -> invalid char in String
}
}
s = sb.toString();
Run Code Online (Sandbox Code Playgroud)
该解决方案的优点是您不必创建太多的 String 对象(每个“替换”操作都会创建一个新的 String,如果您必须恢复大量DNA 序列,这可能会降低性能)
这是基于来自Polygenelubricants和rsp的非常有用的评论的更高性能版本:
String s = "ACGTA";
char[] reverse = new char[s.length()];
for (int i = 0; i < reverse.length; i++) {
switch(s.charAt(i)) {
case 'A': reverse[i] = 'T';break;
case 'T': reverse[i] = 'A';break;
case 'C': reverse[i] = 'G';break;
case 'G': reverse[i] = 'C';break;
default: //handle error here -> invalid char in String
}
}
s = new String(reverse);
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
2269 次 |
| 最近记录: |