cosmo for particular sequences
0
0
Entering edit mode
Guest User ★ 13k
@guest-user-4897
Last seen 10.2 years ago
Hi! I'm trying to find motifs in a file of 50 sequences in FASTA format, they are particular in that they have no "G" so when I run this: >library(cosmo) > seqal <- system.file("Exfiles/short702k.FASTA", package="cosmo") > res <- cosmo(seqs=seqal, constraints="None", minW=10, maxW=10, models="TCM") I get the message: new Irange 202 is greater than pre-calculated Irange 0. numSites = 2. Exiting. If I add only one ficticious G in one of the sequences then it works, but I get some G's in the seqlogo plot Ana -- output of sessionInfo(): new Irange 202 is greater than pre-calculated Irange 0. numSites = 2. Exiting. -- Sent via the guest posting facility at bioconductor.org.
seqLogo seqLogo • 920 views
ADD COMMENT

Login before adding your answer.

Traffic: 521 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