ClustalW und Omega von EBI und Blast von NCBI richten beide Sequenzen global aus. Smith Waterman von EBI richtet Sequenzen lokal aus, arbeitet aber mit nur zwei Sequenzen. Wie kann ich mehr als 2 Sequenzen lokal ausrichten?
Sie müssen es paarweise tun (es gibt keine andere Option).
Sie können das mit BLAST
oder sogar SSEARCH
(Smith-Watermann) tun. Wenn Sie eine fasta
Datei Ihrer Sequenzen haben, gehen Sie folgendermaßen vor:
makeblastdb -dbtype nucl -in <yourfastafile> -out <yourdbindex>
BLAST
mit Ihrer <yourfastafile>
als Eingabe und <yourdbindex>
als Datenbank mit den gewünschten Parametern aus. Dadurch wird ein paarweises Alignment zwischen allen Sequenzen durchgeführtEine ähnliche Strategie kann für verwendet werden SSEARCH
.
Even Water (in EBI selbst) ist ein lokales Alignment-Tool für multiples Sequenz-Alignment. Auch in BLAST erhalten Sie bla2seq, das auch für die paarweise Ausrichtung verwendet wird. In jedem dieser Tools müssen Sie nur FASTA-Sequenzen (beliebige Nummer) oder eine beliebige Sequenz mit > Zeichen eingeben und gegen das Subjekt (Datenbank oder Ihre erforderliche Sequenz) ausführen. Wenn Sie die Ergebnisse für alle Ausrichtungen auf einem Bildschirm sehen möchten, wäre ein Tool wie BioEdit ideal. Ich hoffe, das ist hilfreich, ich bin kein Experte für Bioinformatik, also fühl dich frei, mich zu korrigieren :)
Ich benutze Jalview, das ist ein kleines Programm, das in Java geschrieben ist. Es ist einfach zu bedienen und kann viele, viele Dinge tun. Es basiert auf Webdiensten und verfügt über 7 verschiedene Dienste, die nur für die Sequenzausrichtung eingebaut sind. Sie können Ihre Sequenzen aus vielen Dateiformaten laden oder sie per Referenz-ID abrufen. Es ist wirklich ein kleines Schweizer Taschenmesser. Ich empfehle wirklich einen Blick darauf zu werfen. Es ist auf http://www.jalview.org/ verfügbar .
BEARBEITEN:
Ich habe den Teil zur lokalen Ausrichtung vielleicht früher übersehen, aber ich habe ein wenig recherchiert und festgestellt, dass MULAN (MUltiple sequence Local Alignment and Conservation Visualization Tool) möglicherweise das ist, wonach Sie suchen.
jalview
ist nicht für lokale Ausrichtungen.
Benutzer11667
WYSIWYG
blast
@StackUnderFlow Dies ist der Befehl, der in der Befehlszeile (UNIX-Terminal oder Windows-Eingabeaufforderung) ausgeführt werden muss .blast
hat keine GUI oder selbst wenn es eine gibt, habe ich sie nie benutzt. Es ist viel besser, Befehle zu schreiben (es spart RAM und ermöglicht Ihnen auch etwas Automatisierung)Terdon
WYSIWYG
Terdon