Question

如何在bash下求一基因序列的反向互补序列？

例如：ATTCAG 变为 CTGAAT

Answer 1

啊，这个需要生物学专业知识，虽然我不懂，但我通过Google搜索解决了你的问题。

 

Bash代码：

echo "ATTCAG" | tr "[ATGCatgcNn]" "[TACGtacgNn]" | rev

Python代码：

print ''.join(["ATCG"["TAGC".index(n)] for n in "ATTCAG"[::-1]])

Answer 2

这个应该不行，要么自己写一个脚本。最省事儿的还是用一些软件啊~~oligo或者primer premier都可以

追问

谢谢 其实我就是问怎么写script的。