[EMBOSS] Note about restrict -alphabetic
Guy Bottu
gbottu at ben.vub.ac.be
Mon Sep 22 15:18:30 UTC 2003
from : BEN
Dear colleagues,
I just tried something with restrict under EMBOSS 2.7.1. The program sorts by
default the cut sites from begin to end. With parameter -alphabetic it sorts
instead the output line according to the name of the restriction enzyme. All
well, but if there are several cut sites for one enzyme, these are apparently
output randomly. Cf. this example for the command
restrict vec:puc18 -enzymes=haeiii,mboi -alphabetic
-------------------
########################################
# Program: restrict
# Rundate: Mon Sep 22 17:11:15 2003
# Report_format: table
# Report_file: puc18.restrict
########################################
#=======================================
#
# Sequence: PUC18 from: 1 to: 2686
# HitCount: 26
#
# 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
2675 2678 HaeIII GGCC 2676 2676 . .
1741 1744 HaeIII GGCC 1742 1742 . .
1821 1824 HaeIII GGCC 1822 1822 . .
389 392 HaeIII GGCC 390 390 . .
287 290 HaeIII GGCC 288 288 . .
646 649 HaeIII GGCC 647 647 . .
2088 2091 HaeIII GGCC 2089 2089 . .
820 823 HaeIII GGCC 821 821 . .
831 834 HaeIII GGCC 832 832 . .
849 852 HaeIII GGCC 850 850 . .
1283 1286 HaeIII GGCC 1284 1284 . .
2067 2070 MboI GATC 2066 2070 . .
2325 2328 MboI GATC 2324 2328 . .
2342 2345 MboI GATC 2341 2345 . .
2378 2381 MboI GATC 2377 2381 . .
1545 1548 MboI GATC 1544 1548 . .
430 433 MboI GATC 429 433 . .
1373 1376 MboI GATC 1372 1376 . .
1448 1451 MboI GATC 1447 1451 . .
1459 1462 MboI GATC 1458 1462 . .
1467 1470 MboI GATC 1466 1470 . .
277 280 MboI GATC 276 280 . .
1557 1560 MboI GATC 1556 1560 . .
1662 1665 MboI GATC 1661 1665 . .
2003 2006 MboI GATC 2002 2006 . .
2021 2024 MboI GATC 2020 2024 . .
#---------------------------------------
#---------------------------------------
------------------------------------------------------------------------------
IMHO it would be better if the lines were sorted by enzyme name AND for each
enzyme by position from begin to end. Is this also your opinion ?
Sinceely,
Guy Bottu
More information about the EMBOSS
mailing list