SEQanswers (
-   Bioinformatics (
-   -   Amplicon sequencing remove primer & adapters (

wzhang1001 07-17-2012 11:03 AM

Amplicon sequencing remove primer & adapters

I have amplicon sequencing experiment with amplicon length between 45 ~124 bases. The sequencing length is 150, which completely covers the amplicon.

I wanted to remove the primers and adapters in the fastq file, which seems to be a little challenging becase we don't know the sequence of primers (but we do know the adapter sequence).

But we know the sequence of amplicons, so I'm thinking if there is any way that I can retrieve the sequence of amplicons from the the fastQ file, which will automatically get rid of both primer and adapters. Of course I need also retrieve the correspnding Q score for the amplicon sequences.

Is there any tool I can use for this purpose?

Thanks in advance for your help!

JackieBadger 07-17-2012 12:38 PM

You can use PRINSEQ
As far as I remember you can tell the program to trim a certain amount of bases from either end.
You can also do this with Galaxy:

Prinseq may be better as you will retain your Q data.

If you wanted to sort your amplicons while you remove primers and adapters you could use jMHC
For this you need to know your primer sequence but you should be able to figure out what your primers were from your sequence.

jlli 07-18-2012 05:17 AM

You can try

cutadapt :
fastx clipper :

Tha cutadapt has more options.

Krish_143 07-18-2012 05:23 AM

You can also try

Trimmomatic :

wzhang1001 07-19-2012 10:32 AM

Hi all,

Thank youf or all your help!

All times are GMT -8. The time now is 08:51 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.