Can I use nwalign to align one sequence to the reverse-complement of a second sequence?

3 次查看（过去 30 天）
Dan Allison
Dan Allison 2017-3-9
评论： Dan Allison 2017-3-13
I want to see if single stranded DNA will fold back on itself. I know I can just reverse-complement the sequence first, then call nwalign, but I want the alignment to show complementary bases aligning, not matching bases. I tried creating a modified version of the nuc44 scoring matrix, giving AT, TA, GC, CG all weight 5 and other weight -4, but that didn't seem to work. I assume the order in nuc44 is ATGCSWRYKMBVHDN.

请先登录，再回答此问题。

回答（1 个）

Luuk van Oosten
Luuk van Oosten 2017-3-10
Dear Dan Allison,
Wouldn't it be easier if you use the oligoprop function and see if your oligo's form hairpins or dimerize?
Let us take the [personalized] example from MATLAB's oligoprop page:
Dan = oligoprop('ACGTAGAGGACGTN')
and inspect Dan:
Dan =
struct with fields:
GC: 53.5714
GCdelta: 3.5714
Hairpins: 'ACGTagaggACGTn'
Dimers: [3×14 char]
etc etc etc etc
So in this case, yes, there are hairpins.
  1 个评论
Dan Allison
Dan Allison 2017-3-13
Thanks for the response. I was not aware of the oligoprop function. I see that it does identify some hairpins, but not all the same ones that the mFold software does. I'll play around with some of the input settings since I'm just using the defaults for now.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Sequence Alignment 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by