... | ... | @@ -55,11 +55,11 @@ And unzip it: |
|
|
|
|
|
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
|
|
|
obi stats -a score_norm wolf/aligned_reads
|
|
|
|
|
|
### 3. Remove unaligned sequence records
|
|
|
### 3. Remove sequence records with a low overlap alignment score:
|
|
|
|
|
|
obi grep -a mode:alignment wolf/aligned_reads wolf/good_sequences
|
|
|
obi grep -p "sequence['score_norm'] > 0.8" wolf/aligned_reads wolf/good_sequences
|
|
|
|
|
|
### 4. Assign each sequence record to the corresponding sample/marker combination
|
|
|
|
... | ... | |