This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Splitting a large FASTA file into smaller FASTA files by Scaffolds

I am working with a large FASTA file that is subdivided into "Scaffolds", similar to this sample:

>Scaffold1
AAGTCCTGACTNNNNNNNNNNNNCGTACGTATCGATCG
>Scaffold2
ACGTAGCTAATCGAGATCGATCNNNNNNNNNNNNNCGATTACGTACGATGTG
>Scaffold3
ACGTATCGTACGTGTANNNNNNNNNNNNCGTACGATGCTACTCACGGATTACAAA

I was wondering if there was a way in R to split the one large FASTA file into smaller FASTA files each containing one Scaffold. Thank you!

alignment fasta r
• 9,066 views • link
updated by Ram Ram • written by cameron.gudobba • cameron.gudobba

Any particular reason you want to do this in R? You should always try and use the right tool for the job; in this case something like csplit, perhaps? Certainly not R.

• 0 views • link
written by Alexander Skates Alexander Skates

4 answers 

$ (sudo) pip install pyfaidx
$ faidx --split-files scaffolds.fa

I know it's not R but I can't recommend an R solution for this with a straight face.

pyfaidx: https://github.com/mdshw5/pyfaidx/

• 0 views • link
updated by Ram Ram • written by Matt Shirley Matt Shirley

why R?

awk 'BEGIN {O="";} /^>/ { O=sprintf("%s.fa",substr($0,2));} {if(O!="") print >> O;}' input.fa
• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

While I agree that using a base package is much faster and more convenient, this question is for R. So in case someone searches for this in the future for R, here is a functional solution:

library(iterators)
f <- ireadLines('filepath/inputfile.fasta')
outpath <- 'output/directory/fileprefix' ## don't include the .fasta
count <- 1
while(TRUE) {
        d <- try(nextElem(f))
        if(class(d) == 'try-error') break
        if(grepl(">", d)){
                try(close(outfile), silent=TRUE)
                outfile <- file(paste(c(outpath, count, '.fasta'), collapse=''), open='a')
                write(d, outfile, sep='\n')
                count <- count + 1
        } else {
                write(d, outfile, append=TRUE, sep='\n')
        }
}
rm(f)
closeAllConnections()
• 0 views • link
updated by Ram Ram • written by Steven Lakin Steven Lakin

Here is my R solution, assuming that you want to split your fasta file by identifier:

library(Biostrings)
scaffolds = readAAStringSet(filepath = "scaffolds.fasta", use.names = T)
t2 = split(scaffolds, names(scaffolds)) 
sapply(names(t2), function (x) writeXStringSet(t2[[x]], filepath = paste0(x,".fasta")))
• 0 views • link
updated by Ram Ram • written by polepole40 • polepole40

Log in to answer this question.

More posts like this