Skip to content

Commit

Permalink
Updated our cutoff for an alignment to be 6 mismatches in a window of…
Browse files Browse the repository at this point in the history
… 8 aminos
  • Loading branch information
alexander-nord committed Oct 24, 2021
1 parent 9c4f43b commit f863924
Showing 1 changed file with 7 additions and 7 deletions.
14 changes: 7 additions & 7 deletions src/Diviner.pl
Original file line number Diff line number Diff line change
Expand Up @@ -1703,15 +1703,15 @@ sub FindGhostExons
$exon_str = $exon_str.' '.$end_exon;
}

my $target_info = "MSA Ali Region : $exon_str\n";
$target_info = $target_info." Target Genome : $target_species ($chr:$SearchRanges[0]..$SearchRanges[1])\n";
$target_info = $target_info." Source Species : $source_species\n";
my $target_info = "MSA Ali Region : $exon_str\n";
$target_info = $target_info." Target Genome : $target_species ($chr:$SearchRanges[0]..$SearchRanges[1])\n";
$target_info = $target_info." Source Species : $source_species\n";

# Is it an especially elusive ghost we're chasing?
if ($num_tbn_hits == 0) {
print $outf "[ ] Search failure (no tblastn hits)\n";
print $outf " $target_info";
print $outf " Search Sequence: $SearchSeqs[$q]\n\n";
print $outf " Search Sequence : $SearchSeqs[$q]\n\n";
next;
}

Expand Down Expand Up @@ -3090,9 +3090,9 @@ sub LocalMatchMismatchAli
}
}

# We'll set our condition for killing the alignment as being a window of 5
# aminos where 4 of the positions are gaps or mismatches.
my $window_size = 5;
# We'll set our condition for killing the alignment as being a window of 8
# aminos where 6 of the positions are gaps or mismatches.
my $window_size = 8;
my $min_matches = 2;
my $kill_trigger # Scores below this terminate our walk
= ($window_size-$min_matches)*Max($gap,$mismatch) + $min_matches*$match;
Expand Down

0 comments on commit f863924

Please sign in to comment.