primer3-org / primer3

Primer3 is a command line tool to select primers for polymerase chain reaction (PCR).
GNU General Public License v2.0
222 stars 64 forks source link

A question about ntthal to cacualte "ANY" bounds between two oligos #64

Open zenghp88 opened 2 years ago

zenghp88 commented 2 years ago

When I run following command line with ntthal "ANY" mode, the result seems unreasonable. ./ntthal -a ANY -mv 50 -dv 1.5 -n 0.6 -d 50 -s1 ATTAGGAACCCACTCCCAAGATAAC -s2 GTTATCTAGGGAGTGGGTTTCTGAT Calculated thermodynamical parameters for dimer: dS = -285.288 dH = -102700 dG = -14218 t = 46.3398

SEQ     ATT  G           AAGATAAC
SEQ        AG AACCCACTCCC
STR        TC TTGGGTGAGGG
STR     TAG  T           ATCTATTG

but in fact, it shoud be like this:

SEQ     ATT  G           A
SEQ        AG AACCCACTCCC AGATAAC
STR        TC TTGGGTGAGGG TCTATTG
STR     TAG  T           A

And I found that the value of "stackint2Enthalpies[3][0][0][0]" is "inf" which should be illogical. and it may lead to the unreasonable result above. Both v2.5.0 and the latest version have been tested.

Best regards, Huaping Zeng