Re: forcing a full sequence comparison in blast
"Michael Nuhn" <[email protected]> Tue, 10 Apr 2007 09:55:06 +0200
| Newsgroups | gmane.science.biology.informatics.devel |
|---|---|
| Message-ID | <01da01c77b45$8dadb580$80baf683@nbc1> |
Hi, Phil! > Is there a way I can force blast to use the full 347 amino acids for > comparison. No. It does not make sense to put sequences into an alignment that do not belong together. > The researchers in my lab are most interested in places > with low similarities, since they are trying to find the portions > which make this organism virulent. Blast can create output in form of a table. Or you can choose the output option XML and convert that into a table. You can filter this table for genes which made only a bad hit. Where "bad" would have to be defined carefully. Then determine which genes did not make a blast hit at all. That way you would get a set of candidate genes causing virulence. (As far as this is possible with this kind of analysis.) Greetings from Germany, Michael.