Preferred isoschizomer (bis)
    Guy Bottu 
    gbottu at black.vub.ac.be
       
    Thu Apr 24 08:56:15 UTC 2003
    
    
  
	Dear colleagues,
I took a second look, and it is even worse : the file withrefm.304
contains as many as 149 enzymes with restriction site CTGCAG. The file
embossre.enz contains 2 sites CTGCAG (BspMAI and PstI) and 147 sites ctgcag.
When I run restrict with parameter -nolimit it finds the 2 sites and
when I run it in default mode it finds only BspMAI.
There is clearly a bug+misfeature in the programs rebaseextract+restrict. 
One would expect that restrict by default finds PstI and with -nolimit
finds all 149 enzymes (although this would give a monstrous output).
	Sincerely,
	Guy Bottu
    
    
More information about the EMBOSS
mailing list