[Bioperl-l] standaloneblastplus

Smithies, Russell Russell.Smithies at agresearch.co.nz
Thu Sep 9 22:54:43 UTC 2010


Is that a typo in your email or are some of your fasta headers in your db incorrect?
Eg.
>301333052|GU725063|Xiphinema adenohystherum  internal transcribed
>301333052|GU725063|spacer 1
AAAGACGAAAAATATATACTTTCTCACCGAAATATCAGACTATTCGGTTCCGATTTGATTCGCGGTAAGGCGCGCTAGCGCGGTAGCCCGCTATCGGTTCTGACTGTCTGGACCCCGAAAAATAGTAAGAAAGACGGTGTCGTCGGATCTCGGGTTAGTATTGTATATATCGGGGACTTATCGGTCTGT

Shouldn't that be:
>301333052|GU725063|Xiphinema adenohystherum  internal transcribed spacer 1
AAAGACGAAAAATATATACTTTCTCACCGAAATATCAGACTATTCGGTTCCGATTTGATTCGCGGTAAGGCGCGCTAGCGCGGTAGCCCGCTATCGGTTCTGACTGTCTGGACCCCGAAAAATAGTAAGAAAGACGGTGTCGTCGGATCTCGGGTTAGTATTGTATATATCGGGGACTTATCGGTCTGT

Maybe the invalid fasta headers are breaking the db formatter?


Russell Smithies

Technical Support
T +64 3 489 9085
E  russell.smithies at agresearch.co.nz
Invermay  Research Centre
Puddle Alley,
Mosgiel,
New Zealand
T  +64 3 489 3809
F  +64 3 489 9174
www.agresearch.co.nz

> -----Original Message-----
> From: bioperl-l-bounces at lists.open-bio.org [mailto:bioperl-l-
> bounces at lists.open-bio.org] On Behalf Of sally roberts
> Sent: Friday, 10 September 2010 4:10 a.m.
> To: bioperl-l at lists.open-bio.org
> Subject: [Bioperl-l] standaloneblastplus
>
> I am running a test for standaloneblastplus but getting data back that
> does not exist in my query or my local database. Below is a outline of my
> script small database, query list, and erroneous results. As you will
> notice the query list is comprised of the first four sequences found in
> the database. The results say it can not find the first two and then the
> mathces for the last two do not exist!
>
> Thanks for any help!
>
>
>
> Program
>
>
> #!/usr/bin/perl
>
> use Bio::Tools::Run::StandAloneBlastPlus;
>
>
> $fac = Bio::Tools::Run::StandAloneBlastPlus->new(
>   -db_name => 'ITS',
>   -db_data => 'smallDB.fas',
>   -create => 1
> );
>
> $result = $fac->blastn( -query => , 'sequences.fasta',
>                         -outfile => 'ITStest2.bls');
>
>
> smallDB.fas Data
>
> >302585252|HM807352|Waitea circinata  internal transcribed spacer 1
> ATGATTGGTGGCTGTTGCTGGCTAGTGTTTCTAGTATGTGCACGCCACACCTTCAATCCCACTTACACCTGTGC
> ACCTTTTGTAGTATTACTTGTGGATATCGAGAGAAAGTTTAGTCTTTCACTCTGTTGAAACCGGTTTACTACGT
> TTTTTTATACACACACAATAGTCATTGAATGTATTTTTTATTTCTTATGATAAAAA
>
> >302585252|HM807352|Waitea circinata  internal transcribed spacer 2
> GAATCTCTCAAATACAATAATTTTTTTTGATTGTTGTATTTGGACTTGGAAGCTGTTGGCGCAAGTCGACTCTT
> CTCAAATGTATTAGCTGGGGTTTATATAGTTGGATCCTTGGTGTGATAATTATCTACGCCTTGAAGTCCCTGTA
> GACTCTGCTTCAAATCGTCTCTTCATGAGACAATATTTGAATCA
>
> >302585250|HM802273|Fusarium oxysporum  contains 18S ribosomal RNA,
> internal transcribed spacer 1, 5.8S ribosomal RNA, internal transcribed
> spacer 2, and 28S ribosomal RNA"
> CTCGAGGTCTTACCGAGTTTCACTCCCAACCCCTGTGAACATCACCAATTGTTGCCTCGGCGGATCAGCCCGCT
> CCCGGTAAAACGGGACGGCCCGCCAGAGGACCCCTAAACTCTGTTTCTATATGTAACTTCTGAGTAAAACCATA
> AATAAATCAAAACTTTCAACAACGGATCTCTTGGTTCTGGCATCGATGAAGAACGCAGCAAAATGCGATAAGTA
> ATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACATTGCGCCCGCCAGTATTCTGGCGGGCAT
> GCCTGTTCGAGCGTCATTTCAACCCTCAAGCCCCCGGGTTTGGTGTTGGGGATCGGCGAGCCCTTGCGGCAAGC
> CGGCCCCGAAATCTAGTGGCGGTCTCGCTGCAGCTTCCCCTGCATAGAAAACCCGCGGGGGGGGGAC
>
> >302585249|HM802272|Fusarium oxysporum  contains 18S ribosomal RNA,
> internal transcribed spacer 1, 5.8S ribosomal RNA, internal transcribed
> spacer 2, and 28S ribosomal RNA"
> GAGGTCTTACCGAGTTTCACTCCCAACCCCTGTGAACATACCACTTGTTGCCTCGGCGGATCAGCCCGCTCCCG
> GTAAAACGGGACGGCCCGCCAGAGGACCCCTAAACTCTGTTTCTATATGTAACTTCTGAGTAAAACCATAAATA
> AATCAAAACTTTCAACAACGGATCTCTTGGTTCTGGCATCGATGAAGAACGCAGCCAAATGCGATAAGTAATGT
> GAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACATTGCGCCCGCCAGTATTCTGGCGGGCATGCCT
> GTTCGAGCGTCATTTCAACCCTCAAGCACAGCTTGGTGTTGGGACTCGCGTTAATTCGCGTTCCTCAAATTGAT
> TGTCGTTCACGTCGAGCTTCCATAGCGTAGTAGTAAAACCCTCGTTACTGGTAATCGTCGCGGCCACGCCGTTA
> AACCCCAACTTCTGAATGTTGACCTCGGATCAGTAAGGAATACCCGCTGAACTTAAGCATATCATTAAGCGGAG
> GAA
>
> >302585248|HM802271|Fusarium oxysporum  contains 18S ribosomal RNA,
> internal transcribed spacer 1, 5.8S ribosomal RNA, internal transcribed
> spacer 2, and 28S ribosomal RNA"
> CCGAGGTCTTACCGAGTTTCACTCCCAACCCCTGTGAACATACCAATTGTTGCCTCGGCGGATCAGCCCGCTCC
> CGGTAAAACGGGACGGCCCGCCAGAGGACCCCTAAACTCTGTTTCTATATGTAACTTCTGAGTAAAACCATAAA
> TAAATCAAAACTTTCAACAACGGATCTCTTGGTTCTGGCATCGATGAAGAACGCAGCAAAATGCGATAAGTAAT
> GTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACATTGCGCCCGCCAGTATTCTGGCGGGCATGC
> CTGTTCGAGCGTCATTTCAACCCTCAAGCCCCCGGGTTTGGTGTTGGGGATCGGCGAGCCCTTGCGGCAAGCCG
> GCCCCGAAATCTAGTGGCGGTCTCGCTGCAGCTTCCATTGCGTAGTAGTAAAACCCTCGCAACTGGTACGCGGC
> GCGGCCAAGCCGTTAAACCCCCAACTTCTGAATGTTGACCTCGGATCAGGTAGGAATACCCGCTGAACTTAAGC
> ATATCATTAAAGCGGAGGAA
>
> >301333053|GU725064|Xiphinema turcicum  internal transcribed spacer 1
> GGAGAGATTATATCTTTCTCGAAAAGAGAAAAAATATCCGAGCCGAGCGAACCGACCGAAAAACGCGGTGAGGC
> GCCTTTTGCGCAAAGTCCGTACGTCGGTTCTTAGCGAATATAGCCTCGGCCTGGGACCCGAAAGATGTTTCCTA
> TATGTATCTCGAGACCGACCGTTTAAGACGGTAGCCGGAAAAAAGATTATACCGTGGGTGAAGGTGTCGAAAAG
> AATAATGTAGGTAAAAAAGAAAGACAGACAGAGGAGAGAAAGAACGAAAGTAGAACTCGAACGTAGTTTGAGCT
> ACGCAGTAACGGTATCCGTCGTGGGACATCGCGGTGCGTCGGTTGTAGGGAGTTAAGATTACCTACCCGACACC
> TCGATATTAATCCCGCGCGAATAAATGCGGATTACCGTGAATGTACGCTCTGCTTCGATATCGGGCTTCTTTTG
> ACACCGAAAATATATATATGAATAAAAATAAAGTCACCCTCGTTGCAACGGTATATATCAAAGCGGTTTTCCGT
> GAAAAGAAAGAAGGCGGCTTCGGTTCTCGTTATATTAGGAATAATCTAAGTAATTTCAGACGTCCCGGGAATCG
> TTACTATAGATAGAGAGCGATAGTAACGGTTTCTCCTTCGGGTACTTATCGAACGTTAACACTGCGGTAATCCG
> TCTGGCCGCAAGGAGAGAGGTGTTACGTTCGGCAGCCCTAAATTTCGACCCGTTCGACTAATGCGACGGCCCTA
> CCGAGAAAATGTAGGGCCTATGTACATAGTCCGAAAGAAATACGATCGGAATATTAAGGGTTAGGTTTAAAGAG
> TCATCGGTTCCGAGTACGCGTTCGTTCGGCACGATGCGTGTGTGTATATATCGTAGAGGAGTATTGACGATATA
> TATGTATGCGTATTCGCCCTTACGATAAGAGAATATCGCGTAATTCGGAGCGGCCGTTCTTCGCGAGAGAGAGA
> ACGCA
> CGCGTTAGAAGCTTACGAGTCGGTGTTAAGTTCGAAGGAGAGAGGTTCGAACCGAAGCCGGCGAGTACGCGTTA
> AGTCGTTTCGCGAGAGACGGTCCGGGACGAAAAGGAGAGAGTATCGTCCGGGTGTCCGCCCGAAATAGATATCT
> TATCGAGAATATTTTTATATAGTTCGTTAGAAAGAATGCGAACTTTAAA
>
> >301333052|GU725063|Xiphinema adenohystherum  internal transcribed spacer
> 1
> AAAGACGAAAAATATATACTTTCTCACCGAAATATCAGACTATTCGGTTCCGATTTGATTCGCGGTAAGGCGCG
> CTAGCGCGGTAGCCCGCTATCGGTTCTGACTGTCTGGACCCCGAAAAATAGTAAGAAAGACGGTGTCGTCGGAT
> CTCGGGTTAGTATTGTATATATCGGGGACTTATCGGTCTGTCGAGTTTCTTTCCGGGGTTCTTTGAGTTTATTG
> GGACAGTGTCGGTTGTAGCGGAGTTTGGATACCTACCCGACTGTCCTTATAGTCTCGTGAACACGAGCCCGGGA
> ATAGAAGAGACTCGGCTGATAACGACCGACTATATCTCGTTATATACTCAGAGTTGAATAACTGAGTGGCTCGA
> AACGGCGACATTGTACTTACTATTTTATGTAGACTCTGGAAATATCAGACGTCCCGGGGAATCGTTACAGAGGA
> AATATAGGGTACCTGGAAAAAGAATGGTACCCGTTCCTGTAATGATTCCTTATTCGGGTACCTATCGAATACTA
> ACGGCGCGGATCCCCCGTCTGGCCGCGACGGAATAAGCGTTAGATTCGGTATCCCTATATTCGCGAGTATTCGA
> CTAGTCATGAAATAGAGCCCTTATCGGGGTATCGACTGTCGATCGGATAGAAAGCGAATTAGGGTTAGGTTTAA
> AGAGTCATTGGTTCCGTATATATGGGTGGAACGTACCCGTAAAGGAACAGCCGTAGACGCGAGTTCGGAAATAA
> GTATATTCTCGCGAGAAAGAGGGTCCGTGTACCTTCAAGGTACTTGAATTTAGACCCAGTCTCGTGAATATACG
> TAACTCGTCGAATGGCTCGGGACATGTAGAATACTATGTCCGGGTGACCGCCCGAAATAAGAATATTCATCAGA
> AACTTTTATATATAGTTCGCCGAATAATAGCGAAC
>
> >301333051|GU725062|Xiphinema sphaerocephalum  internal transcribed spacer
> 1
> AAAGTCGAAAAAATATACTTTCTCGCGGAGAAATAATACGGACCGTTCAGTCCGACTCTATACGCGGTAAGGCG
> CTCTTGCGCGAGAGCCCGCTGTCGGTTCTGACGGTCCGGACCCCGAAAAGTAGTAAGTACGACTACGATATATC
> GTGGTCGAGTATCGGTTAGTAATAGTATATCGGGACTGACCGATCGGTCGGTCGAGTTTCTACCGGCTTCTTTG
> AGTCTATTCGGGCAGCGTCGGTTGTAGCGGAGTTTGGATACCTACCCGACTGTCCTTGTAGAACTCGTGAATTC
> GAGCTCGGTAACCGGGAACTCGGCTGAGAACGACCGATTACTTCTCGATACGCTCGAACGTATATATCTAACCG
> AGAAAAGGCGACGTTGTACTTACTATTTATATCAGACGTCCCGAGAGTCGTTACGGTCGGAAATATTGGGTACC
> GGTATCGGACCCGTTTCCGTATCGGCTCTTTATTCGGGTACCTATCGAATACTAACGCCGCGGTTCACCGTCTG
> GCCGCGACGGAATACGCGTTAGATTCGGCACCCCCTATATTCGTATATATATCGACTAGTCTCGAAATAGAGCC
> CTTACTAGGGTGAAGACTATGTCGATCGGAAAGAATCGGATTAGGGGTAGGTTTAAAGAGTCATCGGTTCCGTG
> TATCCGGGCGAAATATATACCCGTAACGGAACGACCGTTGACGCGAGTTTGAAGATATATACATGTACGTATAT
> GAGACAAAAAAACGAGGGTCTGTACCGTGAATTTTTTAGGTACCGAAAAGAGGACCCCCGGTCTCGTGAATATG
> TATTACTCGCCGAACGGTTCGGGACATGGAGAATATTATGTCCGGGTGACCGCCCGAAATAGAAATTTTTTTCT
> ATAAAGTTTTGATATACGTATAGTTCGTCGAATAAAAGC
>
> >301333050|GU725061|Xiphinema hispanum  internal transcribed spacer 1
> AAAGCCGAAAAATATATACTTTCTCAGAGAAATACTAGACTAGTCGATTCCGACTTGATTCGCGGTAAGGCGCT
> TTCGCGCGATAGCCCGCTGTCGGTTCCGACCGTCTGGACCCCGAAAAATAGTAAGAACGACGGTGTCGTCGATC
> TCGGTTAGAAATTGTATATATGTCGGGACGGATCGGTCGGTCGAGTTCCTTTCGGTGTTCTTAGAGTTTATTCG
> GGCAGTGTCGGTTGTAGCGGAGTTTGGATACCTACCCGACTGTCCTTATAATCTCGTGAACTTAGAACCCGGAA
> TAGAGGGAACTCGGCTGATAACGACCGACTTATGTCTCGCCGTATACCGTGAGTTATTTGACCGAGTGGCTCGA
> AACGGCGGTATTGTACTTACTATTTATCTAGTCTCTGGAAATATCAGACGTCCCGGGAATCGTTACAGCGGAAA
> TATAGGGTACCCGAAAAACTGGTACCCGTTTCTGAAACGACTCCTTATTCGGGTACCTATCGAATACTAACGCC
> GCAGTTTCCCGTCTGGCTGCGATGGAAAAAGCGTTAGATTCGGGATCTCTATATTCGCGGGTGTTCGATTAGTC
> GTGAAATACAGCCCTTACGCGGGTGACGACGGTCGATCGGAAAGAAAGCGAATTAGGGTTAGGTTTAAAGAGTC
> ATTGGTTCCGTGTACGGGCGAAAAAGTACCCGTTACGGAACGGCCGTCGACGCGAGTGTGGAAATAAGTATATA
> GTTACGAGAAAGAGGGTCTGTACCTCGGAGTTTTTTGAAGGTACCGTAATCAGGACCCTGTCTCGTGAATATAC
> AAGTTACTCGCCGAACGGTTCGGCCAATGTAGAATTTTATGTCCGGGTGACCGCCCGAAATAGAAATATTTCAT
> AAAAAGCTTTTATATATAGTTTGCCGAATAATAGCAAACG
>
> >301333049|GU725060|Xiphinema pyrenaicum  internal transcribed spacer 1
> AAAGCGGAAAAATTACTTTCTCACCCGGAAAAAACAGACCGTTTATCGGTCCGACTTGAAACGCGGTAAGGCGC
> TCTTGCGCGATAGCCCGCCGTCGGTTCCGATGGTCTGGACCCCGAAAAATAGTAAGAACGACGGTGTCGTCGAT
> TCTCGGTTAGTAGTATATCCGGTCGGATCGATATATATCGGTCGGTCGAGTTTCTATCGGGTTCTTTGAGTTTC
> TTCGGACAGCGTCGGTTGTAGTGGAGTTTGGATACCTACCCGACTGTCCTTATAATCTCGTGAACTCTAGCCCG
> ATAATAATACGGAACTCGGCTGAGAACGACCGACTTAGGTCTGAGTAGATATACTGAGAATATTACCTAGCCGA
> GATGAACGAAACGGCGACATTGGAGTTTTACTATTTACTCGTATCAGACGTCCCGGGAATCGTTGCAGTTGAAT
> TACATATATACGGGTACCTGTAATTGGACTCGTTTCTGTAACGGTTCTTTAGTCGGGTACCTATCGAATACTAA
> CGCCGCGGTTATCCGTCTGGCCGCGATGGAATAAGCGTTAGATTCGGCATCCCTTTATTCGTATACGTTCGAGT
> AGTCGTGAATTAGAACCCTTTAACCGGGGTGAAGACTATCGACGGGAGATAAGCGAATTAGGGGTAGGTTTAAA
> GAGTCATCGGTTCCGGATACGGAGAGAAAAATGCCCGTAATGGAACGACCATTGAAGCGGGATCTATATATATA
> TATATATGATTCGCCCGATGGTTCGGGACATGGAGAATTTTATGTCCGGGTGACCGCCCGAAATAGAAATATTT
> ACTTCAAAGTTATTTATATATAGTTCGCCTTATAAGAGCGAACG
>
>
>
> sequences.fasta data
>
> >Test1
> ATGATTGGTGGCTGTTGCTGGCTAGTGTTTCTAGTATGTGCACGCCACACCTTCAATCCCACTTACACCTGTGC
> ACCTTTTGTAGTATTACTTGTGGATATCGAGAGAAAGTTTAGTCTTTCACTCTGTTGAAACCGGTTTACTACGT
> TTTTTTATACACACACAATAGTCATTGAATGTATTTTTTATTTCTTATGATAAAAA
>
> >Test2
> GAATCTCTCAAATACAATAATTTTTTTTGATTGTTGTATTTGGACTTGGAAGCTGTTGGCGCAAGTCGACTCTT
> CTCAAATGTATTAGCTGGGGTTTATATAGTTGGATCCTTGGTGTGATAATTATCTACGCCTTGAAGTCCCTGTA
> GACTCTGCTTCAAATCGTCTCTTCATGAGACAATATTTGAATCA
>
> >Test3
> CTCGAGGTCTTACCGAGTTTCACTCCCAACCCCTGTGAACATCACCAATTGTTGCCTCGGCGGATCAGCCCGCT
> CCCGGTAAAACGGGACGGCCCGCCAGAGGACCCCTAAACTCTGTTTCTATATGTAACTTCTGAGTAAAACCATA
> AATAAATCAAAACTTTCAACAACGGATCTCTTGGTTCTGGCATCGATGAAGAACGCAGCAAAATGCGATAAGTA
> ATGTGAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACATTGCGCCCGCCAGTATTCTGGCGGGCAT
> GCCTGTTCGAGCGTCATTTCAACCCTCAAGCCCCCGGGTTTGGTGTTGGGGATCGGCGAGCCCTTGCGGCAAGC
> CGGCCCCGAAATCTAGTGGCGGTCTCGCTGCAGCTTCCCCTGCATAGAAAACCCGCGGGGGGGGGAC
>
> >Test4
> GAGGTCTTACCGAGTTTCACTCCCAACCCCTGTGAACATACCACTTGTTGCCTCGGCGGATCAGCCCGCTCCCG
> GTAAAACGGGACGGCCCGCCAGAGGACCCCTAAACTCTGTTTCTATATGTAACTTCTGAGTAAAACCATAAATA
> AATCAAAACTTTCAACAACGGATCTCTTGGTTCTGGCATCGATGAAGAACGCAGCCAAATGCGATAAGTAATGT
> GAATTGCAGAATTCAGTGAATCATCGAATCTTTGAACGCACATTGCGCCCGCCAGTATTCTGGCGGGCATGCCT
> GTTCGAGCGTCATTTCAACCCTCAAGCACAGCTTGGTGTTGGGACTCGCGTTAATTCGCGTTCCTCAAATTGAT
> TGTCGTTCACGTCGAGCTTCCATAGCGTAGTAGTAAAACCCTCGTTACTGGTAATCGTCGCGGCCACGCCGTTA
> AACCCCAACTTCTGAATGTTGACCTCGGATCAGTAAGGAATACCCGCTGAACTTAAGCATATCATTAAGCGGAG
> GAA
>
>
>
>
> Results
>
> BLASTN 2.2.24+
>
>
> Reference: Zheng Zhang, Scott Schwartz, Lukas Wagner, and Webb
> Miller (2000), "A greedy algorithm for aligning DNA sequences", J
> Comput Biol 2000; 7(1-2):203-14.
>
>
>
> Database: ITS
>            5 sequences; 1,102 total letters
>
>
>
> Query=  Test1
> Length=204
>
>
> ***** No hits found *****
>
>
>
> Lambda     K      H
>     1.33    0.621     1.12
>
> Gapped
> Lambda     K      H
>     1.28    0.460    0.850
>
> Effective search space used: 202071
>
>
> Query=  Test2
> Length=192
>
>
> ***** No hits found *****
>
>
>
> Lambda     K      H
>     1.33    0.621     1.12
>
> Gapped
> Lambda     K      H
>     1.28    0.460    0.850
>
> Effective search space used: 189507
>
>
> Query=  Test3
> Length=437
>
> Score     E
> Sequences producing significant alignments:
> (Bits)  Value
>
> dbj|AB581518.1|  Uncultured fungus genes for 18S rRNA, ITS1 and 5...
> 300    2e-085
> dbj|AB581521.1|  Uncultured fungus genes for 18S rRNA, ITS1 and 5...
> 69.4    6e-016
> dbj|AB581519.1|  Uncultured fungus genes for 18S rRNA, ITS1 and 5...
> 58.4    1e-012
> dbj|AB581522.1|  Uncultured fungus genes for 18S rRNA, ITS1 and 5...
> 56.5    4e-012
>
>
> >dbj|AB581518.1| Uncultured fungus genes for 18S rRNA, ITS1 and 5.8S rRNA,
> partial
> sequence, clone: G59F
> Length=203
>
>  Score =  300 bits (162),  Expect = 2e-085
>  Identities = 176/182 (96%), Gaps = 4/182 (2%)
>  Strand=Plus/Plus
>
> Query  10   TTACCGAGTTT-C-ACTCCC-AACCCCTGTGAACATCACCAATTGTTGCCTCGGCGGATC
> 66
>             ||||||||||| | |||||| |||||| |||||||| |||| ||||||||||||||||||
> Sbjct  23   TTACCGAGTTTACAACTCCCAAACCCCAGTGAACAT-ACCACTTGTTGCCTCGGCGGATC
> 81
>
> Query  67   AGCCCGCTCCCGGTAAAACGGGACGGCCCGCCAGAGGACCCCTAAACTCTGTTTCTATAT
> 126
>             ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
> Sbjct  82   AGCCCGCTCCCGGTAAAACGGGACGGCCCGCCAGAGGACCCCTAAACTCTGTTTCTATAT
> 141
>
> Query  127  GTAACTTCTGAGTAAAACCATAAATAAATCAAAACTTTCAACAACGGATCTCTTGGTTCT
> 186
>             ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
> Sbjct  142  GTAACTTCTGAGTAAAACCATAAATAAATCAAAACTTTCAACAACGGATCTCTTGGTTCT
> 201
>
> Query  187  GG  188
>             ||
> Sbjct  202  GG  203
>
>
> >dbj|AB581521.1| Uncultured fungus genes for 18S rRNA, ITS1 and 5.8S rRNA,
> partial
> sequence, clone: G64F
> Length=217
>
>  Score = 69.4 bits (37),  Expect = 6e-016
>  Identities = 39/40 (97%), Gaps = 0/40 (0%)
>  Strand=Plus/Plus
>
> Query  149  AATAAATCAAAACTTTCAACAACGGATCTCTTGGTTCTGG  188
>             ||||| ||||||||||||||||||||||||||||||||||
> Sbjct  178  AATAAGTCAAAACTTTCAACAACGGATCTCTTGGTTCTGG  217
>
>
> >dbj|AB581519.1| Uncultured fungus genes for 18S rRNA, ITS1 and 5.8S rRNA,
> partial
> sequence, clone: G60F
> Length=206
>
>  Score = 58.4 bits (31),  Expect = 1e-012
>  Identities = 39/42 (92%), Gaps = 3/42 (7%)
>  Strand=Plus/Plus
>
> Query  146  ATAA-ATAAATCAAAACTTTCAACAACGGATCTCTTGGTTCT  186
>             |||| || ||| ||||||||||||||||||||||||||||||
> Sbjct  165  ATAACAT-AAT-AAAACTTTCAACAACGGATCTCTTGGTTCT  204
>
>
> >dbj|AB581522.1| Uncultured fungus genes for 18S rRNA, ITS1 and 5.8S rRNA,
> partial
> sequence, clone: G65F
> Length=256
>
>  Score = 56.5 bits (30),  Expect = 4e-012
>  Identities = 30/30 (100%), Gaps = 0/30 (0%)
>  Strand=Plus/Plus
>
> Query  157  AAAACTTTCAACAACGGATCTCTTGGTTCT  186
>             ||||||||||||||||||||||||||||||
> Sbjct  225  AAAACTTTCAACAACGGATCTCTTGGTTCT  254
>
>
>
> Lambda     K      H
>     1.33    0.621     1.12
>
> Gapped
> Lambda     K      H
>     1.28    0.460    0.850
>
> Effective search space used: 442850
>
>
> Query=  Test4
> Length=521
>
> Score     E
> Sequences producing significant alignments:
> (Bits)  Value
>
> dbj|AB581518.1|  Uncultured fungus genes for 18S rRNA, ITS1 and 5...
> 309    4e-088
> dbj|AB581521.1|  Uncultured fungus genes for 18S rRNA, ITS1 and 5...
> 69.4    7e-016
> dbj|AB581519.1|  Uncultured fungus genes for 18S rRNA, ITS1 and 5...
> 58.4    1e-012
> dbj|AB581522.1|  Uncultured fungus genes for 18S rRNA, ITS1 and 5...
> 56.5    5e-012
>
>
> >dbj|AB581518.1| Uncultured fungus genes for 18S rRNA, ITS1 and 5.8S rRNA,
> partial
> sequence, clone: G59F
> Length=203
>
>  Score =  309 bits (167),  Expect = 4e-088
>  Identities = 177/181 (97%), Gaps = 3/181 (1%)
>  Strand=Plus/Plus
>
> Query  7    TTACCGAGTTT-C-ACTCCC-AACCCCTGTGAACATACCACTTGTTGCCTCGGCGGATCA
> 63
>             ||||||||||| | |||||| |||||| ||||||||||||||||||||||||||||||||
> Sbjct  23   TTACCGAGTTTACAACTCCCAAACCCCAGTGAACATACCACTTGTTGCCTCGGCGGATCA
> 82
>
> Query  64   GCCCGCTCCCGGTAAAACGGGACGGCCCGCCAGAGGACCCCTAAACTCTGTTTCTATATG
> 123
>             ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
> Sbjct  83   GCCCGCTCCCGGTAAAACGGGACGGCCCGCCAGAGGACCCCTAAACTCTGTTTCTATATG
> 142
>
> Query  124  TAACTTCTGAGTAAAACCATAAATAAATCAAAACTTTCAACAACGGATCTCTTGGTTCTG
> 183
>             ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
> Sbjct  143  TAACTTCTGAGTAAAACCATAAATAAATCAAAACTTTCAACAACGGATCTCTTGGTTCTG
> 202
>
> Query  184  G  184
>             |
> Sbjct  203  G  203
>
>
> >dbj|AB581521.1| Uncultured fungus genes for 18S rRNA, ITS1 and 5.8S rRNA,
> partial
> sequence, clone: G64F
> Length=217
>
>  Score = 69.4 bits (37),  Expect = 7e-016
>  Identities = 39/40 (97%), Gaps = 0/40 (0%)
>  Strand=Plus/Plus
>
> Query  145  AATAAATCAAAACTTTCAACAACGGATCTCTTGGTTCTGG  184
>             ||||| ||||||||||||||||||||||||||||||||||
> Sbjct  178  AATAAGTCAAAACTTTCAACAACGGATCTCTTGGTTCTGG  217
>
>
> >dbj|AB581519.1| Uncultured fungus genes for 18S rRNA, ITS1 and 5.8S rRNA,
> partial
> sequence, clone: G60F
> Length=206
>
>  Score = 58.4 bits (31),  Expect = 1e-012
>  Identities = 39/42 (92%), Gaps = 3/42 (7%)
>  Strand=Plus/Plus
>
> Query  142  ATAA-ATAAATCAAAACTTTCAACAACGGATCTCTTGGTTCT  182
>             |||| || ||| ||||||||||||||||||||||||||||||
> Sbjct  165  ATAACAT-AAT-AAAACTTTCAACAACGGATCTCTTGGTTCT  204
>
>
> >dbj|AB581522.1| Uncultured fungus genes for 18S rRNA, ITS1 and 5.8S rRNA,
> partial
> sequence, clone: G65F
> Length=256
>
>  Score = 56.5 bits (30),  Expect = 5e-012
>  Identities = 30/30 (100%), Gaps = 0/30 (0%)
>  Strand=Plus/Plus
>
> Query  153  AAAACTTTCAACAACGGATCTCTTGGTTCT  182
>             ||||||||||||||||||||||||||||||
> Sbjct  225  AAAACTTTCAACAACGGATCTCTTGGTTCT  254
>
>
>
> Lambda     K      H
>     1.33    0.621     1.12
>
> Gapped
> Lambda     K      H
>     1.28    0.460    0.850
>
> Effective search space used: 530378
>
>
>   Database: ITS
>     Posted date:  Aug 27, 2010  9:43 AM
>   Number of letters in database: 1,102
>   Number of sequences in database:  5
>
>
>
> Matrix: blastn matrix 1 -2
> Gap Penalties: Existence: 0, Extension: 2.5
>
>
>
>
>
>
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
=======================================================================
Attention: The information contained in this message and/or attachments
from AgResearch Limited is intended only for the persons or entities
to which it is addressed and may contain confidential and/or privileged
material. Any review, retransmission, dissemination or other use of, or
taking of any action in reliance upon, this information by persons or
entities other than the intended recipients is prohibited by AgResearch
Limited. If you have received this message in error, please notify the
sender immediately.
=======================================================================




More information about the Bioperl-l mailing list