All pastes #2133465 Raw Edit

Someone

public python v1 · immutable
#2133465 ·published 2012-03-29 01:55 UTC
rendered paste body
#!/usr/bin/python2.7import stringfrom sys import *from Bio import SeqIOfrom Bio.Seq import Seqfrom Bio.SeqRecord import SeqRecorddef stripPunct(astr):    #For now, only output letters to satisfy EMBOSS    tmpstr = astr.translate(None, string.digits + string.punctuation + string.whitespace)    #tmpstr = astr.translate(None, string.whitespace)    return tmpstrdef main():    txtfile1 = open(argv[1],"r")    txtfile2 = open(argv[2],"r")    txt1 = txtfile1.read()    txt2 = txtfile2.read()    txtfile1.close()    txtfile2.close()    txtstrip1 = stripPunct(txt1)    txtstrip2 = stripPunct(txt2)    seq1 = Seq(txtstrip1)    seq2 = Seq(txtstrip2)    rec1 = SeqRecord(seq1,id="FILE_ONE")    rec2 = SeqRecord(seq2,id="FILE_TWO")    out1 = open(".tmp1.fasta", "w")    out2 = open(".tmp2.fasta", "w")    SeqIO.write(rec1, out1, "fasta")    SeqIO.write(rec2, out2, "fasta")            if __name__ == "__main__":    main()