| Example Program Global Alignments Computing an optimal global alignment between two sequences. A tutorial about global alignments.
1 | //#include <iostream>
| 2 |
| 3 | #include <seqan/basic.h>
| 4 | #include <seqan/file.h>
| 5 | #include <seqan/align.h>
| 6 |
| 7 | int main()
| 8 | {
| 9 | using namespace seqan;
| 10 | typedef Value<Gaps<Dna5String, ArrayGaps> >::Type TValue;
| 11 | /*
| 12 | using namespace seqan;
| 13 |
|
14 | typedef String<Dna> TSequence;
| 15 | TSequence seq1 = "atcgaatgcgga";
| 16 | TSequence seq2 = "actcgttgca";
|
17 | Score<int> scoringScheme(0, -1, -1, -2);
|
18 | Align<TSequence, ArrayGaps> align;
| 19 | resize(rows(align), 2);
| 20 | assignSource(row(align, 0), seq1);
| 21 | assignSource(row(align, 1), seq2);
| 22 |
| 23 | int score = globalAlignment(align, scoringScheme);
| 24 | std::cout << "Score = " << score << std::endl;
| 25 | std::cout << align << std::endl;
|
26 | score = globalAlignment(align, MyersHirschberg());
| 27 | std::cout << "Score = " << score << std::endl;
| 28 | std::cout << align << std::endl;
|
29 | typedef StringSet<TSequence, Dependent<> > TStringSet;
| 30 | typedef Graph<Alignment<TStringSet, void> > TAlignmentGraph;
| 31 |
| 32 | TStringSet string_set;
| 33 | appendValue(string_set, seq1);
| 34 | appendValue(string_set, seq2);
| 35 | TAlignmentGraph alignment_graph(string_set);
| 36 |
| 37 | score = globalAlignment(alignment_graph, scoringScheme, Gotoh());
| 38 | std::cout << "Score = " << score << std::endl;
| 39 | std::cout << alignment_graph << std::endl;
| 40 | */
| 41 | return 0;
| 42 | }
|
|