Pode usar gsub
e uma expressão regular que insira o marcador de quebra de linha (\n
) a cada n caracteres. Por exemplo:
gsub("(.{1,4})", "\\1\n", "ABCDEFGHIJKLMN")
#> [1] "ABCD\nEFGH\nIJKL\nMN\n"
Pode escrever como uma função, aproveitando para incluir nela a criação da linha com o nome:
generateNucSeq <- function(n) {
v = sample(c("A", "C", "G", "T"), n, replace = TRUE)
paste(v, collapse = "")
}
makeFastaLine <- function(sequence, n = 60) {
seqName <- deparse(substitute(sequence))
seqWrpd <- gsub(paste0("(.{1,", n, "})"), "\\1\n", sequence)
paste0(">", seqName, "\n", seqWrpd)
}
sequence1 <- generateNucSeq(350)
sequence2 <- generateNucSeq(270)
multiFASTA <- file("multiFASTA.fasta")
writeLines(makeFastaLine(sequence1), makeFastaLine(sequence2), con = multiFASTA)
close(multiFASTA)
readLines("multiFASTA.fasta")
#> [1] ">sequence1"
#> [2] "GTTTCGGCTCGTGCGGGTGTGCAGTCCATCGTCTCAGTCCAAGTGGTATTCAACCTTGGT"
#> [3] "AACAGGAGGATCTCCGCCGTTATCTTCCACACTCATATGCATATGACATCGGACAAGGGG"
#> [4] "AATCTGTCATGGCTGCCTAGTAGGTCGCTATAAACCAGGCTGCTCTCCGGCAGACAGCAT"
#> [5] "TCCGCATCATCCGTCGTTGGGTTCCATTCGTCGCTCTATCATAATAGTTAATCAGGCAAC"
#> [6] "CCGGTGGTAGTCCCATCACGGACTTAGATGAGAACCTACGATACTGCCGTAACAGCGCGC"
#> [7] "GACTTACCACATCGTCGAGTTGGACGGGTCACGCTGCTAGGTCATTACGT"
#> [8] ">sequence2"
#> [9] "CCCGCTCAGTTAATCTTGATCACTGAAAGCCTATCCAGCTGCACACAGACGAACTGTGTA"
#> [10] "GTTATTATCTTTAATCGTAACCGTTACAGCGTTCCACCAAACTTTATCAATATTGGGAAG"
#> [11] "GAGTGATATACATAAAGCTTTTCCGGACCGGGGTTTGCATTGGGAACAGGAGCCTTTTAG"
#> [12] "TATCCAAAGATACGAGGTTAGCGCCGCCGTATCCAAGATCTATCTGCGTCGGTAACTTGA"
#> [13] "GTTGGAATGCGAGCACGAGTGTTAGCTCCG"
help('strwrap')
.strwrap
divide palavras; não sei se há alguma maneira de indicar número fixo de caracteres para quebrar uma string sem espaços brancos.