Dear all,
Basically I want to trim DNA sequences using forward and reverse primer set. The trimmed sequences still contain the forward and reverse primer but what before and after would have been trimmed off
Two strategies would be: 1) "Perfect match then trim": i.e. find the location in the sequences that are perfect match to forward and reverse primers, then do the trim; I wrote a script can do it very well; 2) "Fuzzy match then trim", i.e. when searching the primer matching strings in the sequences, allow up to 2 mismatches, and then trim the sequence; I have a difficult time to do so; Bio::Grep can perform the fuzzy search, but how should I trim the seq and print the trimmed sequence out?
Thank you all for any suggestions.
-
Are you posting in the right place? Check out Where do I post X? to know for sure.
-
Posts may use any of the Perl Monks Approved HTML tags. Currently these include the following:
<code> <a> <b> <big>
<blockquote> <br /> <dd>
<dl> <dt> <em> <font>
<h1> <h2> <h3> <h4>
<h5> <h6> <hr /> <i>
<li> <nbsp> <ol> <p>
<small> <strike> <strong>
<sub> <sup> <table>
<td> <th> <tr> <tt>
<u> <ul>
-
Snippets of code should be wrapped in
<code> tags not
<pre> tags. In fact, <pre>
tags should generally be avoided. If they must
be used, extreme care should be
taken to ensure that their contents do not
have long lines (<70 chars), in order to prevent
horizontal scrolling (and possible janitor
intervention).
-
Want more info? How to link
or How to display code and escape characters
are good places to start.
|