view seq_shuffle_dinucleotide.py @ 2:e72172a0a429 default tip

Uploaded
author xuebing
date Sat, 31 Mar 2012 21:49:00 -0400
parents
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()