The results can be found in Supplementary Data 9. Generally, the BLAST searches found very weak matches with poor e-values and/or low identity percentages (most under 30%), further validating the Novel Core HG obtained by combining BLAST, MCL, and our script.
Pattern ranking of the metabolite changes induced by the low dose Lisinopril treatment against the established specific metabolite patterns present in MetaMap Tox yielded weak matches with reduced feed consumption (in line with the clinical findings) in addition to that for kidney toxicity (ACE inhibitor like) (Supplemental Figure 8).
(Percentages do not add up to 100% because of a very small number of probe sets with weak matching to multiple clones53).
Even though expression in some cell types within the brain was strong, overall expression of variant Il4i1_2 in the brain was weak, matching the results obtained with pooled brain tissue in the Northern analysis.
Splice sites with weak matches to the consensus sequences are common, though it is not clear how such sites are efficiently utilized.
A script written with the BioPerl toolkit [ 70] was used to parse the BLAST output files for sequences that did not have any matches, or sequences with only weak matches using the criteria: (E-value ≥ 0.01), or (Percent Identity < 50%), or (<50% of the query length was used in the BLAST alignment).
Alignments were filtered with probcalc to eliminate weak matches (Pchance > 0.05), and short (<35 bp aligned, or < 32 matching bp) or ambiguous alignments removed with custom Perl scripts.
Known adaptor sequences were compared against the raw reads using standalone blast (blastn) with permissive settings, to allow detection of weak matches as well as perfect matches (-W 4, -F F, -e 6000).
The E-values thus obtained were converted to scores between 0 and 1, with 0 representing a strong match and 1 representing a weak match.
It's a weak match-up from the Darwinian perspective but I'm going to have to go with the idea that Jayhawks could at least peck you to death.
