[Bioperl-l] something wrong with my code?

Ambrose Ng ng at cs.usask.ca
Thu Oct 2 15:50:53 EDT 2003


my code here:
		
my @s_ind = $hsp->seq_inds('hit', 'id', 1);
print "conserved=", my $s_ind = join (" ", at s_ind), "\n";

reports this:

conserved=7896-8225 8249-8949


but the numbers are different from the blast report.

Is it more likely to do with my code or is there something going on with the 
bioperl module compatibility with this blast report that I'm 
misunderstanding?

Ambrose

-------------- next part --------------
Query=  gnl|UG|Os#S16394561 Oryza sativa (japonica cultivar-group) cDNA clone:00
2-141-A12, full insert sequence /gb=AK108259 /gi=32993468 /ug=Os.4 /len=1031
        (1031 letters)

Database:  /fdf/coe05/gm0/0/rice_Chr1f.dna
           426 sequences; 63,195,892 total letters.
Searching.......................................................done.

                                                               Score     E
Sequences producing significant alignments:                    (bits)  Value

o_sativa|chr_1|OJ1116_H09|6985                                  1933  0.0

>o_sativa|chr_1|OJ1116_H09|6985
            Length = 110086

 Score = 1933 bits (975), Expect = 0.0
 Identities = 1031/1124 (91%)
 Strand = Plus / Plus


Query: 1    aagctagctcgcgcgcaaattaatcgccgatcgatcctacgatgaagcggaagcggccgg 60
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 7826 aagctagctcgcgcgcaaattaatcgccgatcgatcctacgatgaagcggaagcggccgg 7885


Query: 61   cggcgctgcgcggcggggaggaggcggcggcggcggcgctgaagcgtgggccgtggacgc 120
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 7886 cggcgctgcgcggcggggaggaggcggcggcggcggcgctgaagcgtgggccgtggacgc 7945


Query: 121  ccgaggaggacgaggtgctggcgcggttcgtggcgcgggaggggtgcgaccggtggcgca 180
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 7946 ccgaggaggacgaggtgctggcgcggttcgtggcgcgggaggggtgcgaccggtggcgca 8005


Query: 181  cgctgccgcggcgcgcgggcctgctgcgctgcggcaagagctgccgcctccggtggatga 240
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8006 cgctgccgcggcgcgcgggcctgctgcgctgcggcaagagctgccgcctccggtggatga 8065


Query: 241  actacctccgccccgacatcaagcgctgccccatcgccgacgacgaggaggacctcatcc 300
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8066 actacctccgccccgacatcaagcgctgccccatcgccgacgacgaggaggacctcatcc 8125


Query: 301  tccgcctccaccgcctcctcggcaaccggt---------------------------ggt 333
            ||||||||||||||||||||||||||||||                           |||
Sbjct: 8126 tccgcctccaccgcctcctcggcaaccggtaagca<--[8161:8243]-->taggtggt 8251


Query: 334  cgctgatcgccgggaggttgccggggcgcacggacaacgagatcaagaactactggaact 393
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8252 cgctgatcgccgggaggttgccggggcgcacggacaacgagatcaagaactactggaact 8311


Query: 394  cgcatctcagcaagaagctcatcgcgcagggcatcgacccgcggacgcacaagccgctga 453
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8312 cgcatctcagcaagaagctcatcgcgcagggcatcgacccgcggacgcacaagccgctga 8371


Query: 454  cggccgccgccgatcactccaacgccgccgctgccgtcgccgccacttcttacaagaagg 513
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8372 cggccgccgccgatcactccaacgccgccgctgccgtcgccgccacttcttacaagaagg 8431


Query: 514  cggtgccggccaagccgccgaggacggcatcctcgccggccgctggcattgagtgcagcg 573
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8432 cggtgccggccaagccgccgaggacggcatcctcgccggccgctggcattgagtgcagcg 8491


Query: 574  acgatcgtgcccgaccggccgacggtggcggtgacttcgcagcgatggtgagcgccgccg 633
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8492 acgatcgtgcccgaccggccgacggtggcggtgacttcgcagcgatggtgagcgccgccg 8551


Query: 634  atgccgagggattcgaaggaggatttggcgatcagttctgtgccgaggatgcagttcatg 693
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8552 atgccgagggattcgaaggaggatttggcgatcagttctgtgccgaggatgcagttcatg 8611


Query: 694  gtggcttcgacatgggttccgcttccgccatggtgggtgacgacgacttctcctcgtttc 753
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8612 gtggcttcgacatgggttccgcttccgccatggtgggtgacgacgacttctcctcgtttc 8671


Query: 754  ttgattctctgatcaacgacgagcagttaggcgatctcttcgtcgtcgagggcaacgatc 813
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8672 ttgattctctgatcaacgacgagcagttaggcgatctcttcgtcgtcgagggcaacgatc 8731


Query: 814  acgagcatggcaatggtgagattggtcatggagacgtcatggaatccaaacagtaatctt 873
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8732 acgagcatggcaatggtgagattggtcatggagacgtcatggaatccaaacagtaatctt 8791


Query: 874  tcgggagattgatcagggaataagttgaccatgagaaaacatgtaagaacagcttgtctg 933
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8792 tcgggagattgatcagggaataagttgaccatgagaaaacatgtaagaacagcttgtctg 8851


Query: 934  ctgagtcatgaccggtgatgtgtatgtgatggaagaaaacgtatgtacgacttaatgctt 993
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 8852 ctgagtcatgaccggtgatgtgtatgtgatggaagaaaacgtatgtacgacttaatgctt 8911


Query: 994  gcacctatactactgctgtgtgcattattcgagtgatc 1031
            ||||||||||||||||||||||||||||||||||||||
Sbjct: 8912 gcacctatactactgctgtgtgcattattcgagtgatc 8949


More information about the Bioperl-l mailing list