[EMBOSS] Oddcomp behaves oddly ...
Marc Logghe
Marc.Logghe at DEVGEN.com
Wed Mar 8 09:00:14 UTC 2006
... Or rather, how should I use it properly ?
OK, suppose your run compseq to obtain the frequency for individual
residues:
compseq tsw:Q62671 -word 1
Apparently this example protein sequence is rather rich in leucine (106
L out of 889).
In order to detect this leucine bias, a little file was created
(leu.comp) that had the following content:
<file leu.comp>
Word size 1
Total count 0
# bias should be detected as 106 > 100
L 100
</file leu.comp>
Oddcomp was run like this:
oddcomp tsw:Q62671 -infile leu.comp -window 889
But the sequece is not reported.
When I change the L count to 10 in leu.comp it does not work neither.
Strangely enough, when the default window is taken (30) the sequence is
reported.
What is happening here ?
Regards,
Marc
More information about the EMBOSS
mailing list