rev_compl_V0.py

########################################################
# nimarafati@gmail.com	                               #
# Please cite the script by referencing to github      #
# repository 					       #
########################################################

from Bio import SeqIO
from Bio.SeqRecord import SeqRecord
import argparse 
import re

parser = argparse.ArgumentParser(description="This script generates reverse complementary sqeuence\n.")
parser.add_argument("-fa", "--fasta", help = "A fasta sequence is provided here.", action = "store")
parser.add_argument("-up", "--uppercase", help = "Convert sequence to uppercase.", action = "store_true")
args = parser.parse_args()

def make_rc_record(record):
#	"""Returns a new SeqRecord with the reverse complement sequence."""
	return SeqRecord(seq = record.seq.reverse_complement(), id = "rc_" + record.id, description = "reverse complement")

print(args.fasta)
records = map(make_rc_record, SeqIO.parse(args.fasta, "fasta"))
SeqIO.write(records, "rev_comp_" + args.fasta, "fasta")