Sigcleave and report

Gary Williams, Tel 01223 494522 gwilliam at hgmp.mrc.ac.uk
Thu May 20 15:38:15 UTC 2004


It looks as if the same problem affects restrict:
The text in the tail of the report is not cleared.
For example the second sequence in the list has no hits, but the
fragments from the first sequence are reported again.

Gary

% restrict @list -enz acci -site 4 stdout -frag 
Finds restriction enzyme cleavage sites
########################################
# Program: restrict
# Rundate: Thu May 20 16:34:01 2004
# Report_format: table
# Report_file: stdout
########################################

#=======================================
#
# Sequence: AY582533     from: 1   to: 3846
# HitCount: 2
#
# Minimum cuts per enzyme: 1
# Maximum cuts per enzyme: 2000000000
# Minimum length of recognition site: 4
# Blunt ends allowed
# Sticky ends allowed
# DNA is linear
# Ambiguities allowed
#
#=======================================

  Start     End Enzyme_name Restriction_site 5prime 3prime 5primerev
3primerev
   1773    1778 AccI        GTMKAC             1774   1776        
.         .
   2018    2023 AccI        GTMKAC             2019   2021        
.         .

#---------------------------------------
#
# Fragment lengths:
#     1827
#     1774
#     245                  
#
#---------------------------------------
#=======================================
#
# Sequence: AY582532     from: 1   to: 4635
# HitCount: 0
#
# Minimum cuts per enzyme: 1
# Maximum cuts per enzyme: 2000000000
# Minimum length of recognition site: 4
# Blunt ends allowed
# Sticky ends allowed
# DNA is linear
# Ambiguities allowed
#
#=======================================

  Start     End Enzyme_name Restriction_site 5prime 3prime 5primerev
3primerev

#---------------------------------------
#
# Fragment lengths:
#     1827
#     1774
#     245
#
#---------------------------------------  



Henrikki Almusa wrote:
> 
> Hello,
> 
> There was a bug in sigcleave. The tail text "No scores over " is left in tail
> after printing report without any hits. If there is multiple sequences in
> input file and sequence later have hits, then same text appears in report
> with hits as well.
> 
> However it seems that there is no function like "ajReportClearTail" so it's
> fixed by putting empty AjPString in there after writing a report. Would it be
> useful functions to clear header and tail? If not then, I'll put this to
> emboss-submit as is.
> 
> --
> Henrikki Almusa
> 
>   ------------------------------------------------------------------------
>                                Name: sigcleave_tail_fix.patch
>    sigcleave_tail_fix.patch    Type: text/x-diff
>                            Encoding: 7bit

-- 
Gary Williams
MRC Rosalind Franklin Centre for Genomics Research
Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SB, UK
Tel: +44 1223 494522			Fax: +44 1223 494512
E-mail: gwilliam at rfcgr.mrc.ac.uk	Web: http://www.rfcgr.mrc.ac.uk



More information about the emboss-dev mailing list