Scrape Sequences

Tutorial

library(ape)

seq1 <- paste("JQ", seq(983161, 983255), sep = "")

sequences <- read.GenBank(seq1,
                          seq.names = seq1,
                          species.names = TRUE,
                          as.character = TRUE)

write.dna(sequences, "fish.fasta", format = "fasta")

Assigned Organism

library(ape)

seq1 <- read.table("sequence.seq")
seq1 <- seq1[,1]

sequences <- read.GenBank(seq1,
                          seq.names = seq1,
                          species.names = TRUE,
                          as.character = TRUE)

write.dna(sequences, "pan_paniscus.fasta", format = "fasta")