diffseq problems
Obi Griffith
obig at bcgsc.ca
Fri Jun 13 19:52:22 UTC 2003
I have been having two problems with your diffseq program. I am using
command line diffseq from the EMBOSS-2.5.0 package.
1) using the -columns option doesn't have any effect on the output. ie.
I can't get the nicely formatted tabular data promised:
# Report of diffseq of: AF129756 and AP000504
# AF129756 overlap starts at 6036
# AP000504 overlap starts at 1
# (AF129756) start end length sequence (AP000504) start end length sequence
6882 6882 1 't' 847 847 1 'a'
7830 7830 1 'a' 1795 1795 1 'g'
8307 8307 0 '' 2273 2273 1 't'
8500 8500 1 'a' 2466 2466 1 'g'
8688 8688 0 '' 2655 2658 4 'tgtg'
2) I get the following errors:
Find differences (SNPs) between nearly identical sequences
Input sequence:AF129756
Second sequence:AP000504
Word size [10]: 2
Output report [1.diffseq]:
Warning: bad /replace value 'T'
Warning: bad /replace value 'C'
Warning: bad /replace value 'C'
Warning: bad /replace value 'T'
Warning: bad /replace value 'C'
Warning: bad /replace value 'GCCTT'
Any idea what's going on here?
Obi.
More information about the EMBOSS
mailing list