[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