How to not penalize a specific mismatch when aligning?
0
0
Entering edit mode
4.7 years ago
ervaulsy • 0

Hi, I'm trying to align bisulfite sequencing read file to a reference genome. I'm duplicating data analysis method in a paper. In the paper they are saying "... Following this three-letter alignment, we refined each alignment using the Smith–Waterman algorithm. We used a score matrix that does not penalize a mismatch between A in the read and G in the top strand of reference genome sequence" They used SOAP2 Aligner, I looked for options but I couldn't find one for this. How and where can i create a score matrix that doesn't penalize a specific mismatch (A and G)? Alignment rates were really low when I skipped this step.

alignment genome sequence • 526 views
ADD COMMENT

Login before adding your answer.

Traffic: 2710 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6