#!/usr/bin/env python import sys #function that will count the nucleotides def count_nucleotide_types(seq): result = [] totalA = seq.count('A') totalC = seq.count('C') totalG = seq.count('G') totalT = seq.count('T') result.append(totalA) result.append(totalC) result.append(totalG) result.append(totalT) return result #open the sequence file determined in the command line sequencefile = open(sys.argv[1], 'r').readlines() #join the list elements sequence = ''.join(sequencefile) #replace the carriage returns sequence = sequence.replace('\n', '') #call the function and count the nucleotides values = count_nucleotide_types(sequence) print values