StephenFordham

import re
DNA = input('Please enter your DNA sequence...').upper()
while re.search('[^AGCT]', DNA):
     exit("DNA sequence entered is not valid")

StephenFordham

DNA = input('Please enter your DNA sequence...')

StephenFordham

import re
DNA = input('Please enter your DNA sequence...').upper()
while re.search('[^AGCT]', DNA):
    DNA = input('Please enter a valid DNA sequence...').upper()

StephenFordham

import sys
pattern = sys.argv[1].upper()
pattern2 = sys.argv[2].upper()
count_no = int(sys.argv[3])

StephenFordham

import re
import sys
# DNA = input('Please enter your DNA sequence...').upper()
# while re.search('[^AGCT]', DNA):
#     DNA = input('Please enter a valid DNA sequence...').upper()

DNA = open('E.coli_genome.fasta').read()

pattern = sys.argv[1].upper()
pattern2 = sys.argv[2].upper()

StephenFordham

import argparse
parser = argparse.ArgumentParser(description='Finding restriction sites in DNA sequences')

StephenFordham

parser.add_argument('input_filename')
args = parser.parse_args()

StephenFordham

DNA = open(args.input_filename).read().rstrip('\n')

StephenFordham

parser.add_argument('pattern',
                    help='Enter the pattern of the restriction enzyme')

parser.add_argument('count_no',
                    help='count number to be specified')

parser.add_argument('input-filename',
                    help='Enter the name of your DNA sequence')

args = parser.parse_args()

StephenFordham

parser.add_argument("-i",
                    "--input_filename",
                    help='Enter the name of your DNA sequence')