Short sequence tissue specific BLAST
Posted 15 January 2013 - 12:33 AM
I'm looking at a short motif sequence of 8 bases and I would like to find all genes and transscripts with this sequence in cardiac tissue. I tried blasting, but the short sequence seems to cloud results or I'm not supplying the correct parameters to the algorithm. I managed to get some hits from refseq RNA after playing with the "expected" and "word" settings, but surprisingly the transscript which the sequence is actually in (ubiquitin) was nowhere to be found, so I'm pretty sure I didn't get all the results (even though my max results were set at 20.000). For the tissue specific part, I have no idea where to begin. Any ideas how to solve this problem would be greatly appreciated!
Posted 15 January 2013 - 09:52 AM
As for the algorithm, I don't know much about the background theory of BLAST, but since I assume you expect only 100% match, you can probably increase the penalty for missmatches to limit it further. But that's just an idea.
Our country has a serious deficiency in lighthouses. I assume the main reason is that we have no sea.
I never trust anything that can't be doubted.
'Normal' is a dryer setting. - Elizabeth Moon