... | @@ -53,6 +53,10 @@ And unzip it: |
... | @@ -53,6 +53,10 @@ And unzip it: |
|
|
|
|
|
obi alignpairedend -R wolf/reads2 wolf/reads1 - > aligned_reads.fastq
|
|
obi alignpairedend -R wolf/reads2 wolf/reads1 - > aligned_reads.fastq
|
|
|
|
|
|
|
|
We can check the average alignment score (corresponding approximately to the similarity between the sequences in the overlap):
|
|
|
|
|
|
|
|
obi stats -a score_norm wolf/aligned_reads
|
|
|
|
|
|
### 3. Remove unaligned sequence records
|
|
### 3. Remove unaligned sequence records
|
|
|
|
|
|
obi grep -a mode:alignment wolf/aligned_reads wolf/good_sequences
|
|
obi grep -a mode:alignment wolf/aligned_reads wolf/good_sequences
|
... | | ... | |