Выравнивание последовательности печати в файл

Я делаю простое попарное выравнивание последовательностей ДНК с помощью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

Как я могу вывести полное выравнивание в текстовый файл?


person Martin Preusse    schedule 13.04.2012    source источник
comment
почему бы не задать вопрос в списке рассылки Bioconductor? Подписка не требуется. Я думаю, что ответ "вы не можете, напрямую", но что-то вроде as.character(pattern()) даст вам что-то полезное.   -  person Martin Morgan    schedule 13.04.2012
comment
только что отправил письмо в список ...   -  person Martin Preusse    schedule 13.04.2012
comment
Было ли когда-либо найдено решение для этого?   -  person user1357015    schedule 03.06.2012
comment
Он находится в разработке, проверьте список рассылки биокондукторов.   -  person Martin Preusse    schedule 14.06.2012


Ответы (2)


Вот ссылка на обсуждение в списке рассылки биокондукторов. Пока нет простого способа распечатать отформатированное выравнивание, но, возможно, его стоит реализовать.

https://stat.ethz.ch/pipermail/bioconductor/2012-April/044904.html

person Martin Preusse    schedule 20.04.2012

Я думаю, что функция R printPairwiseAlignment() из пакета Biostrings предназначена для того, чтобы сделать это так, как я думаю, вы ищете.

person user3767180    schedule 23.06.2014