3
$\begingroup$

I am aligning a large number of ESTs. It seems poly-A tails show in many different ways. In addition to occurring at the very end, they can be flanked by the cloning sequence one one end, or have mismatches/errors. What is a good rule or available tools that will handle the usual cases?

A few examples of the non-trivial cases I found, with their Genbank Accs:

>EE409337
... AAAAAAAAAAAAAAAAAAAAAAAAAGGAAAAAAAAAAAAAAAAAAAAAAAAAAAACCTTGTC
>EE409340
... TTTCTACAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACTTGTC
>EE409361
... TTGTTAAACTGAAAAAAAAAAAAAAAAAAAAAAAAAAAACCATGTCGGC
TTACTGAATTGAA
>EE420306
.... AAAAAAAGTTATGTTAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAGGGAAAAAAA
AAAAAAAAAAAAAAAAA

Cross-posted on SeqAnswers,Biostars.

$\endgroup$
  • $\begingroup$ Did the answer on Biostars not work for you? (If so, what is it that it doesn't do that you would like it to do?) $\endgroup$ – blep Jul 10 '13 at 20:22
  • $\begingroup$ @dd3: trimest has an option for trimming the poly-A regions, but we have to specify the rules. I have already done the basic filtering, but what are the "good" rules? Is there something better than min-length and max-mismatches that will handle real sequences I have given examples of? $\endgroup$ – highBandWidth Jul 10 '13 at 20:34
  • 1
    $\begingroup$ You could try the empirical approach. Take a random sample of N ESTs and calculate the distributions for min-length and max-mismatches. Use the distributions to pick appropriate cutoffs (some sort of statistics thing). $\endgroup$ – blep Jul 10 '13 at 20:55
  • $\begingroup$ @dd3, that definitely makes sense and seems like a good idea. (although there's some judgement to be made for things like ...GAGAACTGGGCAAAA: where does the tail start?). Since EST processing is well established, shouldn't something like this already exist in the tools or the literature? How does everyone else do it? $\endgroup$ – highBandWidth Jul 11 '13 at 3:09
  • 1
    $\begingroup$ Use dust and remove all low complexity sequences.. Or simply use a regex search to filter out poly-As $\endgroup$ – WYSIWYG Aug 22 '14 at 6:07

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.