标签: r bioinformatics fasta
我有一个平坦的.txt文件,其中包含FASTA format中的5个基因。我使用read.fasta()在R中将此文件读作FASTA文件。然后我使用getTrans()将每个基因的DNA序列翻译成多肽,有效翻译每个基因序列。
.txt
read.fasta()
getTrans()
然而,一旦序列被翻译,我就失去了每个基因的所有标题或注释,因此我不能轻易分辨哪个基因是哪个。
有没有办法在保存注释的同时在fasta文件中翻译多个基因?