DECIPHER maximum allowable size
0
0
Entering edit mode
@d27aed9d
Last seen 2.3 years ago
Germany

Hi there I want to align up to 9 sequences using the package DECIPHER. I have managed to read the fasta files and apply readDNAStringSet. When I then run AlignSeqs I get following error

Aligning Sequences:
  |======================                                                                                         |  20% 

Error in f(p.profile[, anchors[2, n - 1]:anchors[1, n], drop = FALSE],  : 
  Alignment larger (3,062,288,908) than the maximum allowable size (2,147,483,647).

Can you help me understand where the issue is and how to solve this?

Thank you!

DECIPHER Alignment alig • 1.8k views
ADD COMMENT
0
Entering edit mode

This error typically occurs when aligning genomes that are not collinear. How long are the sequences you are aligning?

ADD REPLY
0
Entering edit mode

Approximately 104073951 bp. I want to compare several genes. Is DECIPHER limited to short sequences?

ADD REPLY
0
Entering edit mode

AlignSeqs() can handle longish sequences if they are collinear. If not, you can use FindSynteny() followed by AlignSynteny().

ADD REPLY
0
Entering edit mode

I run FindSynteny() and I got following Error:

'dbFile' must be a character string or SQLiteConnection.

I used a Large DNAStringSet as input for the FindSynteny()

ADD REPLY
0
Entering edit mode

FindSynteny() takes a database as input rather than a DNAStringSet. I recommend reading the documentation for alignment.

ADD REPLY

Login before adding your answer.

Traffic: 686 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6