Advertisement
Guest User

Sequence 2 fasta converter

a guest
Aug 23rd, 2010
124
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Perl 2.04 KB | None | 0 0
  1. seq2fasta.pl script:
  2.  
  3. use Bio::SeqIO;
  4. my $in = Bio::SeqIO->new (-fh => \*STDIN);
  5.  
  6. # output in Fasta format
  7. my $out = Bio::SeqIO->new (-format => 'Fasta');
  8.  
  9. while (my $seq = $in->next_seq) {
  10.    $out->write_seq($seq);
  11. }
  12.  
  13.  
  14. sequence file:
  15. >Contig1244/1236-2383 2383 Weight: 1.02
  16. ACCCAGGCCAATGTATAGTGCACATACGCAAATAAGACATTATAAAAAGGTTCATCGCATACTAACATTCGGACTTGCAT
  17. ACAAACCTATAGCAGTGGCTAGAACAGGACCACCAGCATACTCAGATGTGTTTATGATGACTTCACTTGCACATATTCCT
  18. ATAGAATATCCATTCATGTATTTGAATCCATCAGAATTTGCTTTATTACCACCAGGCTCTAAAGTTAAACGTGTTTTAAG
  19. CAGAGTAAAAGCTGAAAAGGTTCGTGTTGCATTTCCAACTAACTCTACAGATAATAATTTGGCTACATTACATACAAACA
  20. GCTTCTTGAGAGTTGGAATTAATTTAAATCAAAAGATTCAATCAGTCAACGCTCAACCAGGTGGATTTGCTGCTGCTCAA
  21. CCTATGATAGCCACTTCTATTACACCATTTGATGCAAAAAATTATCAATTTTGGATTGATAATTTTTATGGTGTTCCAAA
  22. CAATGATAATAAAGATCCAGATCAATTTACAACTCAAACTCCAAGACATCAATTCGGTATACCATGGGTAGCTCAACATT
  23. ACTGTTGTCCAGTTTCTCAAACTAACGACCCTACTAATAGTGGATGGGAAGATTTTCAAGCTCACTTGGAAGAAATTAAA
  24. GTAGATGGACCTAGCGGATATATCGCTGAAATAGAATATTCACCAGCTATGGGATTATTAAAAGCACCAATTCAAGGTAT
  25. ATGGACTGGATTGCCTAGTGTAGCTACCCCTGGAACAGCTAGAAGCATTCCCGTAAATACTGGAGCAGGAAACACTCAAC
  26. ATAGAAGAATTGAAGCATCCGTGACAAATTCTCTTAATACTGGCATAACAGAAGTAACATCAAACTTTAACAGACAACCA
  27. ACAGCAGAAAACTTTACATTGTTAACTCCTATTGAAAAAAATCAGAGAATTTGTGCAGGAATAACTCCTGAATATTTTGC
  28. AAAGTCACAACCTACATTACATGTTGGAGTAATGCCAGTTCCAGCTTTAACAACATCAGCAATTGATAACGTTACAAACA
  29. ACAGTTCATTTACTGACGCTCAAGCATATTTTGAAATAACATGTGAAATGGAAGTAGAATGTGCTTTCCCTACACGCAGA
  30. GCATTAGCTACAATTGCAAAAACAACTG
  31.  
  32.  
  33.  
  34.  
  35. testing:
  36.  
  37. OK
  38. $ seq2fasta.pl < SEQFILE
  39.  
  40. not OK
  41. $ cat SEQFILE | seq2fasta.pl
  42.  
  43. ------------- EXCEPTION: Bio::Root::Exception -------------
  44. MSG: Failed resetting the filehandle; IO error occurred
  45. STACK: Error::throw
  46. STACK: Bio::Root::Root::throw /home/lupey/downloads/bioperl/bioperl-live/Bio/Roo
  47. t/Root.pm:472
  48. STACK: Bio::Tools::GuessSeqFormat::guess /home/lupey/downloads/bioperl/bioperl-l
  49. ive/Bio/Tools/GuessSeqFormat.pm:507
  50. STACK: Bio::SeqIO::new /home/lupey/downloads/bioperl/bioperl-live/Bio/SeqIO.pm:3
  51. 75
  52. STACK: ./seq2fasta.pl:5
  53. -----------------------------------------------------------
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement