Local alignment, linear gap model on DNA (Pairwise Alignment)

• Category: Pairwise Alignment
• Validator: Pairwise Alignment Validator

Description

Given 2000 reads of length 100, compute the local alignment using the linear gap model of all pairs.

Match: 2, mismatch: -3, gap: -2.

Input

• 2000 reads of length 100 (i.e., data/genome.pairwise_alignment.reads.length.100.fa)

Output

For each pair of sequences, output the score of the alignment.

The output must be written into a file.

Example

For simplicity we assume 4 reads of length 10.

Reads:

> 1
TTTTATGGTA
> 2
CACTTATGTA
> 3
TTAACTGTTT
> 4
CAGTGGCTCA

Resulting alignments:

Sequence 0 and 1:
  TTATGGTA
  |||| |||
  TTAT-GTA

Sequence 0 and 2:
  TTT
  |||
  TTT

Sequence 0 and 3:
  TGG
  |||
  TGG

Sequence 1 and 2:
  TT-A-TGT
  || | |||
  TTAACTGT

Sequence 1 and 3:
  CA
  ||
  CA

Sequence 2 and 3:
  CT
  ||
  CT

Output:

0, 1: 12
0, 2: 6
0, 3: 6
1, 2: 8
1, 3: 4
2, 3: 4