|
Lines 7-14
Link Here
|
| 7 |
# license. Please see the LICENSE file that should have been included |
7 |
# license. Please see the LICENSE file that should have been included |
| 8 |
# as part of this package. |
8 |
# as part of this package. |
| 9 |
|
9 |
|
| 10 |
import os, sys, getopt, re, time |
10 |
import re, time |
| 11 |
from string import maketrans |
|
|
| 12 |
from Bio import SeqIO |
11 |
from Bio import SeqIO |
| 13 |
from Bio import Translate |
12 |
from Bio import Translate |
| 14 |
from Bio.Seq import Seq |
13 |
from Bio.Seq import Seq |
|
Lines 277-284
Link Here
|
| 277 |
nice looking 6 frame translation with GC content - code from xbbtools |
276 |
nice looking 6 frame translation with GC content - code from xbbtools |
| 278 |
similar to DNA Striders six-frame translation |
277 |
similar to DNA Striders six-frame translation |
| 279 |
""" |
278 |
""" |
| 280 |
comp = complement(seq) |
279 |
from Bio.Seq import reverse_complement |
| 281 |
anti = reverse(comp) |
280 |
anti = reverse_complement(seq) |
|
|
281 |
comp = anti[::-1] |
| 282 |
length = len(seq) |
282 |
length = len(seq) |
| 283 |
frames = {} |
283 |
frames = {} |
| 284 |
for i in range(0,3): |
284 |
for i in range(0,3): |
|
Lines 399-404
Link Here
|
| 399 |
# {{{ |
399 |
# {{{ |
| 400 |
|
400 |
|
| 401 |
if __name__ == '__main__': |
401 |
if __name__ == '__main__': |
|
|
402 |
import sys, getopt |
| 402 |
# crude command line options to use most functions directly on a FASTA file |
403 |
# crude command line options to use most functions directly on a FASTA file |
| 403 |
options = {'apply_on_multi_fasta':0, |
404 |
options = {'apply_on_multi_fasta':0, |
| 404 |
'quick':0, |
405 |
'quick':0, |
