How to extract FASTA headers in R
1 answer
Here's one example using R:
library(seqinr)
fa <- read.fasta(file = "samples.fasta", seqtype = "DNA")
CON=file("output.csv", "w")
writeLines("genes", CON)
writeLines(names(fa), CON)
close(CON)
However, is R necessary in this case? Can you just use simpler UNIX tools?
echo "genes" > output.csv
grep ">" samples.fasta | cut -c2- >> output.csv
edit: if you want specific parts of the fasta headers, you will need to be more specific in your question
• 0 views
•
link
written
by
dr_cfoster
Log in to answer this question.
More posts like this
-
unable to get feature count results
written by Ravita •
-
Uniprot Database for download
written by SHN •
-
Creating Gene ID and Gene name colum from FASTA IDs
written by WUSCHEL
-
Extract information from prokka gff files
written by Space_Life •
-
Exacting all the headers in a fasta file using R
written by xuenanwang •
-
Subset multiple columns in R or unix
written by julie.sawitzke •
-
Extract fasta-headers that have the same word
written by genomes_and_MGEs •
-
How to extract the last 1000 nt from a group of sequences in a FASTA file?
written by Andrei •
-
How to find common data in replicates?
written by WUSCHEL
-
How to replace Accession Number for GI in a fasta header
written by marcelavillegasp •
I came across this handy tool. CSV⬌FASTA https://cdcgov.github.io/CSV2FASTA/
This does not answer your question - it does not use R. Please accept the answer that addresses your question, do not add answers that explicitly break a requirement without specifying that the requirement is broken.
I'm moving this answer to a comment now.