Create a program that takes two sequences (genetic sequences) as input and can determine if there is a 100% match between the sequences. It should be able to conduct a pairwise comparison of two sequences. The code should output start and end positions of each occurrence of a match.

Also a FASTA sequence will be used as input.

Solution PreviewSolution Preview

This material may consist of step-by-step explanations on how to solve a problem or examples of proper writing, including the use of citations, references, bibliographies, and formatting. This material is made available for the sole purpose of studying and learning - misuse is strictly forbidden.

def readgene(filename):
    infile = open(filename, "r")

    seq = "";
    while True:
       line = infile.readline()
       if (line==""):

       if (line[0] != '>'):
            seq += line[:-1]    # to avoid \n character at the end of the line

    return seq

#gene1_filename = "titin_conserved.txt"
gene1_filename = input("Enter the query filename:")
#gene2_filename = "titin_full.txt"
gene2_filename = input("Enter the ref. filename:")...

By purchasing this solution you'll be able to access the following files:

for this solution

or FREE if you
register a new account!

PayPal, G Pay, ApplePay, Amazon Pay, and all major credit cards accepted.

Find A Tutor

View available Bioinformatics Tutors

Get College Homework Help.

Are you sure you don't want to upload any files?

Fast tutor response requires as much info as possible.

Upload a file
Continue without uploading

We couldn't find that subject.
Please select the best match from the list below.

We'll send you an email right away. If it's not in your inbox, check your spam folder.

  • 1
  • 2
  • 3
Live Chats