Question

Function 2. Find Best Match (40 points) We will use the term genome to refer to the string that represents the complete set of genes in an organism, and sequence to refer to some substring or sub-sequence in the genome. Write a function called findBestMatch that takes a genome and a sequence and returns the similarity score of the best match found in the genome as a float. HINT: Problem 3 from Recitation 5 is very similar to this function. We encourage you to write that function before attempting to write this one. Your function should take two parameters in this order: · a string parameter for the genome a string parameter for the sequence o Your function should return the highest similarity score as a float. Your function should not print anything. Your function MUST be named findBestMatch. · o Edge Cases: sequence is empty genome is empty sequence is longer than genome → return-1 → return-1 -return-2 · · Examples: "ATACGC", "ACT" → 0.66 ·

In C++ please!

Answer 1

#include <iostream>

#include <string>

#include <iostream>

#include <vector>

using namespace std;

float findBestMatch(string sequence1 , string sequence2){

  

if(sequence1.size() ==0 || sequence2.size()==0)

return -1;

if(sequence1.size() < sequence2.size())

return 0;

float count = 0.0;

int j = 0, max = -1;

for(int i = 0;i<sequence1.size();++i ){

count = 0;

for(j = 0;j<sequence2.size();++j, ++i ){

if(sequence1.at(i) != sequence2.at(j))

break;

  

++count;

}

if(j ==sequence2.size()){

count = sequence2.size();

}

if(max < count)

max = count;

}

return 1.0*max/sequence2.size();

  

}

int main()

{

string str1, str2;

cout<<"Enter string 1: "<<endl;

getline(cin,str1);

cout<<"Enter string 2: "<<endl;

getline(cin,str2);

  

cout<<"Similarity Score is: "<<findBestMatch(str1, str2)<<endl;

return 0;

}

=========================================================================
See Output

Enter string 1: ATACGC main.cpp Enter string 2: ACT Similarity Score is: 0.666667 9 float findBestMatch(string sequence1, string sequence2)t 10 if(sequence1.size) 0I sequence2.sizeO-0) return -1; 12 13 14 15 16 if(sequence1.sizeO) < sequence2.sizeO) return 0; float count 0.0; int j = 0, max = -1; for(int 1 = 0;i<sequence! . size();++i )( count0 for(j 0;jesequence2.size();++j, 18 19- 20 21 ++1 )( if(sequence1.at(i)!-sequence2.at(j)) break; ++count; 23 24 25- 26 27 28 29 30 31 return 1.0*max/sequence2.sizeO; 32 if(jsequence2.size))I count sequence2.sizeO if(max < count) max-count; 34 int mainO 35- 36 37 38 39 40 41 string str1, str2; cout<<"Enter string 1: "<<endl; getline(cin,str1); cout<<"Enter string 2: "<<endl; getline(cin,str2);



Thanks, PLEASE RATE