Advertisement
Guest User

Untitled

a guest
Jun 25th, 2017
69
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.72 KB | None | 0 0
  1. from Bio import SeqIO
  2. from Bio.Seq import Seq, SeqRecord
  3. from Bio.Blast import NCBIXML
  4. from Bio.Alphabet import IUPAC
  5.  
  6.  
  7. fh = open("/home/local/GLOBANT/sebastian.bassi/Downloads/blastn_elin.xml")
  8.  
  9. br = NCBIXML.parse(fh)
  10.  
  11. tofile = []
  12.  
  13. for record in br:
  14.     for ali in record.alignments:
  15.         for hsp in ali.hsps:
  16.             seq_q = SeqRecord(seq=Seq(hsp.query, IUPAC.unambiguous_dna), id=record.query, name="", description="")
  17.             seq_s = SeqRecord(seq=Seq(hsp.sbjct, IUPAC.unambiguous_dna), id=ali.title, name="", description="")
  18.             tofile.append(seq_q)
  19.             tofile.append(seq_s)
  20.             break
  21.  
  22. # grabar a fasta
  23. oh = open("output.fasta", "w")
  24. SeqIO.write(tofile, oh, "fasta")
  25. oh.close()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement