reverse numbering

Stefanie Lager stefanielager at fastmail.ca
Fri Dec 6 16:02:21 UTC 2002


Hi,

I saw that GCG and EMBOSS use different numbering system when
reversing sequences. If I align 2 sequences, using water from EMBOOS
and bestfit from GCG, and ask to use the reverse strand of sequence 2.
All EMBOSS programs reverse and complements sequence 2, giving it a
new numbering. 
While bestfit (and other GCG programs) just compares the opposite
strand but keeps the numbering according to the ORIGINAL strand. 

Is it possible to keep the numbering according to the original strand
also in EMBOSS programs, when comparing the reverse strand?

Stefanie

########################################
# Program:  water
# Rundate:  Fri Dec 06 16:50:52 2002
# Report_file: A_B.water
########################################
#=======================================
#
# Aligned_sequences: 2
# 1: A.seq
# 2: B.seq
# Matrix: EDNAFULL
# Gap_penalty: 10.0
# Extend_penalty: 0.5
#
# Length: 33
# Identity:      33/33 (100.0%)
# Similarity:    33/33 (100.0%)
# Gaps:           0/33 ( 0.0%)
# Score: 165.0
# 
#
#=======================================

A.seq            268 gcccaggagagccaaaaagccgagctggaaggc    300
                     |||||||||||||||||||||||||||||||||
B.seq              1 gcccaggagagccaaaaagccgagctggaaggc     33


#---------------------------------------
#---------------------------------------


 BESTFIT of: A.gcg  check: 5728  from: 1  to: 420
 
A.seq
 
 to reverse of: B.gcg2  check: 1476  from: 1  to: 33
 
B.seq

 Symbol comparison table: /disk4/gcg/gcgcore/data/rundata/swgapdna.cmp
 CompCheck: 2335

         Gap Weight:     50      Average Match: 10.000
      Length Weight:      3   Average Mismatch: -9.000

            Quality:    330             Length:     33
              Ratio: 10.000               Gaps:      0
 Percent Similarity: 100.000   Percent Identity: 100.000

        Match display thresholds for the alignment(s):
                    | = IDENTITY
                    : =   5
                    . =   1

 A.gcg x B.gcg2            December 6, 2002 16:51  ..

                  .         .         .   
     268 gcccaggagagccaaaaagccgagctggaaggc 300
         |||||||||||||||||||||||||||||||||
      33 gcccaggagagccaaaaagccgagctggaaggc 1
_________________________________________________________________
    http://fastmail.ca/ - Fast Secure Web Email for Canadians


More information about the EMBOSS mailing list