niccw
3/6/2019 - 3:15 PM

biostring

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)