SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Removing primers, adaptors, how to know if it's good? SS Santos Bioinformatics 25 12-12-2013 05:00 PM
Process to remove primers, adapters, etc. from Illumina data LizBent Bioinformatics 6 05-14-2012 04:08 AM
SE adapters and PE primers PaulineF Illumina/Solexa 2 02-09-2010 11:10 AM
Removing primers Khanjan 454 Pyrosequencing 1 02-05-2010 11:09 AM
removing adapters sequences from ChIPseq data? johannes.rainer Illumina/Solexa 0 02-05-2010 06:50 AM

Reply
 
Thread Tools
Old 07-01-2013, 10:32 PM   #1
xhuister
Member
 
Location: China

Join Date: Apr 2010
Posts: 41
Default Removing adapters, primers and polyA/T tails using Perl?

Hi all,

I've got many short sequences like the following:
ADAPTER1 (or 2) + insert + Poly(A) + GGGG(n) + ADAPTER 1 ( or 2)

I'd like to get the insert part for mapping. Also I need to record the adapter type, the length of poly(A) tail and length of GGGG(n) part.

Because it is more than just adpater trimming, I think it is better to write a script to record and trim the adapter, poly(A), GGGG(n) by myself. I can only write some Perl scripts.

But there should be some insertions or deletions in the adapter/polyA/GGG. Is there any Perl function or BioPerl function or some tool/package can search substrings with insertion/deletions/substitution?

Thank you!
xhuister is offline   Reply With Quote
Old 07-01-2013, 11:49 PM   #2
sphil
Senior Member
 
Location: Stuttgart, Germany

Join Date: Apr 2010
Posts: 192
Default

Maybe you can have a look at http://hannonlab.cshl.edu/fastx_toolkit/. There is trimming as well as clipping of reads in a fasta/q file. Using these and manually (of course using s script) monitor the length should do the job...
sphil is offline   Reply With Quote
Old 07-02-2013, 05:16 PM   #3
xhuister
Member
 
Location: China

Join Date: Apr 2010
Posts: 41
Default

Quote:
Originally Posted by sphil View Post
Maybe you can have a look at http://hannonlab.cshl.edu/fastx_toolkit/. There is trimming as well as clipping of reads in a fasta/q file. Using these and manually (of course using s script) monitor the length should do the job...
Thank you, sphil. The FASTX can only trim the adapter, but it cannot trim the polyA tail (A-stretch). I need to get the A-stretch and G-stretch length, is there any function (or regEx) to get this information? Sometimes, the A-stretch also has some deletions or insertions, like AAAAAAGAAAATAAA.
xhuister is offline   Reply With Quote
Old 07-03-2013, 04:04 AM   #4
sphil
Senior Member
 
Location: Stuttgart, Germany

Join Date: Apr 2010
Posts: 192
Default

hmmm in python there is a regex moduk named "re" since there are mismatches it will be quite tricky to filter thos except you have defined length, which I highly doubt...
sphil is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:46 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO