[EMBOSS] dreg and reverse strand

[Andres Pinzon]

Hi,
Im using dreg to find some patterns on a xanthomonas* genome reverse strand.
This is the command im using:

 dreg -sequence ./campestrisVesicatoria.gb -pattern
'TTC(G|T|C){14,17}TTC(G|A|T)' -outfile campestrisVes-rev.dreg.gb
-rformat3 genbank -sask1

But there are 2 problems:

1) the output file is, in fact, a genbank output file, for instance:

-------------------------------------
   misc_feature    450..471
                     /note="*pat regex1"
-------------------------------------

But it lacks the "complement" keyword, so "Artemis" the annotation
program im using, can not read it as a feature in the reverse strand.
So,  It should look like this:
----------------------------------------------
   misc_feature    complement(450..471)
                     /note="*pat regex1"
-----------------------------------------------

2) It seems that dreg is identifying a wrong pattern. I edited the
output file and put the complement keyword in one of the hits (the
example above).
But theres no pattern at that position.

This is a screenshot:
http://bioinf.ibun.unal.edu.co/~problem/art-dreg.png

Am i doing something wrong?

Thanks in advance,

-- 
Andrés Pinzón cPh.D.
http://groups.google.com/group/Bioinformatica_es
Bioinformatics Center, Colombia EMBnet node
http://bioinf.ibun.unal.edu.co
Tel +57 3165000 ext 16961 Fax +571 3165415
Micology and Phytopathology Laboratory - Los Andes University.
Tel +571 3394949 ext. 2768