@@ -51,7 +51,7 @@ def invertSeq(sseq):
51
51
rnumber = 0
52
52
shuffleperseq = 1
53
53
if len (sys .argv ) < 6 :
54
- print "\n Usage: shuffleGenome <FastA input file> <output file> <Rearrangement size> <Number of rearrangements> <Shuffled sequences per seq>"
54
+ print ( "\n Usage: shuffleGenome <FastA input file> <output file> <Rearrangement size> <Number of rearrangements> <Shuffled sequences per seq>" )
55
55
sys .exit (1 )
56
56
else :
57
57
infileName = sys .argv [1 ]
@@ -77,22 +77,22 @@ def invertSeq(sseq):
77
77
header = infile .readline ()
78
78
infiled = infile .read ()
79
79
filelen = len (infiled )
80
- print "\n Input sequence:%s" % seqfile
81
- print "Input sequence length: %d\n " % filelen
80
+ print ( "\n Input sequence:%s" % seqfile )
81
+ print ( "Input sequence length: %d\n " % filelen )
82
82
83
83
#split genome into filelen/rsize parts
84
84
partlist = []
85
85
pos = 0
86
86
if rsize > filelen :
87
87
rsize = filelen
88
- for part in xrange (0 ,filelen / rsize ):
88
+ for part in range (0 ,filelen / rsize ):
89
89
partlist .append (infiled [pos :(part + 1 )* rsize ])
90
90
pos += rsize
91
91
92
92
partlist .append (infiled [pos :])
93
93
94
94
95
- for shuffleit in xrange (0 ,shuffleperseq ):
95
+ for shuffleit in range (0 ,shuffleperseq ):
96
96
count = 0
97
97
seq = ""
98
98
parttemp = ""
@@ -114,13 +114,13 @@ def invertSeq(sseq):
114
114
parttemp = partlist [part1 ]
115
115
partlist [part1 ] = partlist [part2 ]
116
116
partlist [part2 ] = parttemp
117
- print "Transposition"
118
- print " Positions %d and %d swapped" % (part1 ,part2 )
117
+ print ( "Transposition" )
118
+ print ( " Positions %d and %d swapped" % (part1 ,part2 ) )
119
119
elif operation == 2 :
120
120
#inversion
121
121
partlist [part1 ] = invertSeq (partlist [part1 ])
122
- print "Inversion"
123
- print " Position %d" % part1
122
+ print ( "Inversion" )
123
+ print ( " Position %d" % part1 )
124
124
125
125
126
126
@@ -137,7 +137,7 @@ def invertSeq(sseq):
137
137
#write output to file
138
138
fname = seqfile [:- 4 ]
139
139
fname += "_Shuffled_%d.fna" % (shuffleit + 1 )
140
- print "output: %s" % (fname )
140
+ print ( "output: %s" % (fname ) )
141
141
142
142
fout = open (fname ,'w' )
143
143
fout .write (output )
0 commit comments