matcher output bug ?
David Bauer
bauer at genprofile.com
Fri Jan 11 15:48:29 UTC 2002
Hi,
it seems that matcher somtimes forgets to label the matches with the
usual ":" in the output if the matching region is very short.
Example1 shows no label at all and Example2 the labeling starts somwhere
in the middle of the matching region.
Example1:
########################################
# Program: matcher
# Rundate: Fri Jan 11 16:40:50 2002
# Report_file:
########################################
#=======================================
#
# Aligned_sequences: 2
# 1: c1
# 2: AW370204
# Matrix: EDNAFULL
# Gap_penalty: 12
# Extend_penalty: 4
#
# Length: 166
# Identity: 165/166 (99.4%)
# Similarity: 1/166 ( 0.6%)
# Gaps: 165/166 (99.4%)
# Score: 809
#
#
#=======================================
200 210 220 230 240
c1 GTGTCCGAGTTGGAGAGTGCTCTGCAGATGGAGCCAGCTGCCTTCCAGGC
AW3702 GTGTCCGAGTTGGAGAGTGCTCTGCAGATGGAGCCAGCTGCCTTCCAGGC
40 50 60 70 80
250 260 270 280 290
c1 TTTATATTCTGCTGAGAAGCCAAAGCTGGAAGATGAGCATCTCGTTTTCT
AW3702 TTTATATTCTGCTGAGAAGCCAAAGCTGGAAGATGAGCATCTCGTTTTCT
90 100 110 120 130
300 310 320 330 340
c1 TCTGTCAGATGGGCAAGCGGGGCCTCCA-GGCCACGCAGCTGGCCCGGAG
AW3702 TCTGTCAGATGGGCAAGCGGGGCCTCCAAGGCCACGCAGCTGGCCCGGAG
140 150 160 170 180
350 360
c1 TCTTGGATACACTGGG
AW3702 TCTTGGATACACTGGG
190 200
#---------------------------------------
#---------------------------------------
Example 2:
########################################
# Program: matcher
# Rundate: Fri Jan 11 16:43:55 2002
# Report_file:
########################################
#=======================================
#
# Aligned_sequences: 2
# 1: c1
# 2: BE695392
# Matrix: EDNAFULL
# Gap_penalty: 12
# Extend_penalty: 4
#
# Length: 230
# Identity: 222/230 (96.5%)
# Similarity: 4/230 ( 1.7%)
# Gaps: 222/230 (96.5%)
# Score: 1030
#
#
#=======================================
140 150 160 170 180
c1 TTCGACGTGCGCTCTCGCGAGGAGGCGGCAGCTGGGACCATCCCAGGGGC
BE6953 TGCGATGTGCGCTCTCGCGAGGAG-CGGCAGCTGGGACCATCCCAGGG-C
20 30 40 50 60
190 200 210 220 230
c1 GCTCAACATCCCGGTGTCCGAGTTGGAGAGTGCTCTGCAGATGGAGCCAG
BE6953 GCTCAACATCCCGGTGTCCGAGTTGGAGAGTGCTCTGCAGACTGAGCCAG
70 80 90 100 110
240 250 260 270 280
c1 CTGCCTTCCAGGCTTTATATTCT-GCTGAGAAGCCAAAGCTGGAAGATGA
:::::::::::::
BE6953 CTGCCTTCCAGGCTTTATATTCTCGCTGAGAAGCCAAAGCTGGAAGATGA
120 130 140 150 160
290 300 310 320 330
c1 GCATCTCGTTTTCTTCTGTCAGATGGGCAAGCGGGGCCTCCAGGCCACGC
::::::::::::::::::::::::::::::::::::::::::::::::::
BE6953 GCATCTCGTTTTCTTCTGTCAGATGGGCAAGCGGGGCCTCCAGGCCACGC
170 180 190 200 210
340 350 360
c1 AGCTGGCCCGGA-GTCTTGGATACACTGGG
:::::::::::: :::::::::::::::::
BE6953 AGCTGGCCCGGAAGTCTTGGATACACTGGG
220 230 240
#---------------------------------------
#---------------------------------------
--
Dr. David Bauer
GenProfile AG, Max-Delbrueck-Center, Erwin-Negelein-Haus
Robert-Roessle-Str. 10, D-13125 Berlin, Germany
bauer at genprofile.com, Tel:49-30-94892165, FAX:49-30-94892151
More information about the emboss-dev
mailing list