-
Notifications
You must be signed in to change notification settings - Fork 44
Cpg_ambiguous
Paul Lott edited this page Jan 23, 2013
·
2 revisions
##CpG Model
With added ambiguous character definitions
[Borodovsky, M., Ekisheva, S., Problems and Solutions in Biological Sequence Analysis. Cambride University Press, 2006, pg 80]
##Model File
#STOCHHMM MODEL FILE
<MODEL INFORMATION>
======================================================
MODEL_NAME: CpG
MODEL_DESCRIPTION: Taken from CH3 Durbin/Eddy 3.17
MODEL_CREATION_DATE: August 28,2009
<TRACK SYMBOL DEFINITIONS>
======================================================
TRACK1: A,C,G,T
<AMBIGUOUS SYMBOL DEFINITIONS>
======================================================
TRACK1: N[A,C,G,T], R[A,G], Y[C,T], S[C,G], W[A,T], K[G,T]
TRACK1: M[A,C], B[C,G,T], D[A,G,T], H[A,C,T], V[A,C,G]
<STATE DEFINITIONS>
######################################################
STATE:
NAME: INIT
TRANSITION: STANDARD: P(X)
HIGH: 0.5
LOW: 0.5
####################################################
STATE:
NAME: HIGH
GFF_DESC: HIGH
PATH_LABEL: H
TRANSITION: STANDARD: P(X)
HIGH: 0.5
LOW: 0.5
END: 1
EMISSION: TRACK1: P(X)
ORDER: 0 AMBIGUOUS: AVG
0.2 0.3 0.3 0.2
###################################################
STATE:
NAME: LOW
GFF_DESC: LOW
PATH_LABEL: L
TRANSITION: STANDARD: P(X)
HIGH: 0.4
LOW: 0.6
END: 1
EMISSION: TRACK1: P(X)
ORDER: 0 AMBIGUOUS: AVG
0.3 0.2 0.2 0.3
##################################################
//END
##Sequence File
>Example 3.16 Example
GGCACTGAA