rfilefile-writingfastadna-sequence

Converting multiple text sequence files into Fasta files in R


I have 100s of DNA sequence files that are text files, and i want to convert them in Fasta format. I have tried with cat but its not giving expected output. How can I convert these files into fasta format in R?

example :

file1.txt

ATCTACGTACGTGCATG


file2.txt

CGTAGCATTGCATGATC

Expected output

file1.fa

>file1
ATCTACGTACGTGCATG


file2.fa

>file2
CGTAGCATTGCATGATC

Solution

  • I can do using the for loop, but I didn't find a package or more simplest way.

    This approach is working.

      files1<-list.files(pattern = "*.txt")
      for (i in 1:length(files1))
      {
      logFile = read.table(paste0(files1[i]))
      
      write.table(rbind(paste0(">",files1[i]),logFile),paste0(files1[i],".fa"),row.names = FALSE,col.names = FALSE,quote = FALSE)
      
      }