How can I create a new fasta file after modifying an existing fasta file?
I have a fasta file with 10 protein sequences. Using R I generate a new protein sequence and using write.fasta I insert the new sequence in the original file. Is there a way to create a new fasta file with a new name containing the 11 sequences? Basically I need to "save as".
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
R - How to convert a sequence in a fasta file into a matrix
written by kbaitsi •
-
create sequence similarity network directly from FASTA file
written by taojincs •
-
Augustus output, fasta with original contig name
written by belliardocaro •
-
Extract several parts from fasta header
written by rah •
-
How to download all complete protein sequences from BLAST search
written by milnerk3 •
-
How to get script for inserting information in a table to work?
written by Mimmi Ahlmén •
-
Question: Converting Blast Results/Output to FASTA format
written by apulunuj •
-
SignalP results - extract full sequence for proteins with a signal peptide
written by JMac •
-
how to convert a protein sequence containing fasta file into DNA sequence containing file?
written by sansritisinha •
-
Creating A Dna Sequence Database Locally For Blastplus
written by Eric Normandeau
What have you tried? Show us your code please.
The original fasta file is called histone4.fa. After I have generated the new sequence(new hist) I write
write.fasta(newhist, names = names("histone4.fa"), file.out = "histone4.fa", open = "a")and a new line is added to my open fasta file. If I writewrite.fasta(newhist, names = names("histone4.fa"), file.out = "histone4new.fa", open = "a")a new file is created but only with the new sequence. If I writeit works but I was wondering if there is a more compact way to do it.
*histone = readAAStringSet("histone4.fa")
You could concatenate
histoneandnewhistand write it in a single line.