Mercurial > repos > xuebing > sharplabtool
view seqshuffle.py @ 13:292186c14b08
Uploaded
author | xuebing |
---|---|
date | Sat, 10 Mar 2012 08:17:36 -0500 |
parents | b7f1d9f8f3bc |
children |
line wrap: on
line source
import sys from altschulEriksonDinuclShuffle import * def readFastaFile(infile,outfile): fin = open(infile) out = open(outfile,'w') currSeq='' currSeqname=None for line in fin: if '>' in line: if currSeqname !=None: out.write(currSeqname+dinuclShuffle(currSeq)+'\n') currSeqname=None currSeq='' currSeqname=line else: currSeq=currSeq+line.strip().upper() if currSeqname!=None: out.write(currSeqname+dinuclShuffle(currSeq)+'\n') fin.close() out.close() def readrawseq(infile,outfile): ''' each line is a sequence ''' fin = open(infile) out = open(outfile,'w') for line in fin: out.write(dinuclShuffle(line.strip().upper())+'\n') fin.close() out.close() def main(): seqfile = sys.argv[1] outfile = sys.argv[2] fasta = sys.argv[3] if fasta == 'fasta': readFastaFile(seqfile,outfile) else: readrawseq(seqfile,outfile) main()