Я делаю простое попарное выравнивание последовательностей ДНК с помощьюpairwiseAlignment из пакета Biostrings в Bioconductor:
library('Biostrings')
seq1 = 'ATGCTA'
seq2 = 'ATGTA'
pairwiseAlignment(pattern = seq1, subject = seq2)
Вывод выглядит следующим образом:
Global PairwiseAlignedFixedSubject (1 of 1)
pattern: [1] ATGCTA
subject: [1] ATG-TA
score: -4.091219
Для очень длинных последовательностей вывод усекается и отображается только одна строка:
Global PairwiseAlignedFixedSubject (1 of 1)
pattern: [1] AT-------------------------------------------------...----------------TGTCTTCCAKATCTGGCGCGCCTGGGTTGATATC
subject: [1] ATTGGCGGCCGCGCCACCATGCCAGAGCCAGCGAAGTCTGCTCCCGCCCCG...GAAGGCTGTATGCTGTTGTCTTCAAGATCTGGTACCGCTGGGTTGATATC
score: -29418.8
Как я могу вывести полное выравнивание в текстовый файл?
as.character(pattern())
даст вам что-то полезное. - person Martin Morgan   schedule 13.04.2012