Hmm, the last time I studied biology and genetics was 40 years ago, not only did I forget most of it, but the knowledge has changed tremendously in the interval. In short, I really don't know much on the subject.

However, I know that the human genome sequencing project has used very widely Perl, thanks in part to the very good tools that have been developed for that, including a number of good reusable modules. I would think that you should take a look at BioPerl and at related Perl modules, it is quite likely that the tools to do what you want already exist and are available on the CPAN or eslewhere.

Je suis Charlie.

In reply to Re: Filtering matches of near-perfect-matched DNA sequence pairs by Laurent_R
in thread Filtering matches of near-perfect-matched DNA sequence pairs by onlyIDleft

Title:
Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post, it's "PerlMonks-approved HTML":


  • Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
  • Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
  • Read Where should I post X? if you're not absolutely sure you're posting in the right place.
  • Please read these before you post! —
  • Posts may use any of the Perl Monks Approved HTML tags:
    a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
  • You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
            For:     Use:
    & &amp;
    < &lt;
    > &gt;
    [ &#91;
    ] &#93;
  • Link using PerlMonks shortcuts! What shortcuts can I use for linking?
  • See Writeup Formatting Tips and other pages linked from there for more info.