library(Biostrings) s <- readDNAStringSet("nm.fasta","fasta") subs <- subseq(s["seqA"], start=c(1, 2, 3), end=c(3, 6, 5) # convert it to string subschar <- as.character(subs)