Skip to content
Snippets Groups Projects
Commit 8c4b480a authored by joan's avatar joan
Browse files

input args

parent abe46090
No related branches found
No related tags found
No related merge requests found
from RotTable import *
from Traj3D import *
import argparse
parser = argparse.ArgumentParser()
parser.add_argument("--filename", help="input filename of DNA sequence")
parser.parse_args()
args = parser.parse_args()
def main():
rot_table = RotTable()
traj = Traj3D()
if args.filename:
# Read file
lineList = [line.rstrip('\n') for line in open("plasmid_8k.fasta")]
# Formatting
seq = ''.join(lineList[1:])
traj.compute(seq, rot_table)
else:
traj.compute("AAAGGATCTTCTTGAGATCCTTTTTTTCTGCGCGTAATCTGCTGCCAGTAAACGAAAAAACCGCCTGGGGAGGCGGTTTAGTCGAAGGTTAAGTCAG", rot_table)
print(traj.getTraj())
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment