prepare illumina fastq, extract barcodes

195 views
Skip to first unread message

heikieh

unread,
Feb 26, 2014, 1:41:14 PM2/26/14
to qiime...@googlegroups.com
Dear all,

I would like to analize illumina data (fastq) with the qiime pipeline.
However, in the tutorial for preparing the data I realized that for the split_libraries_fastq.py I need an extra fastq file for my barcodes.
I have found the extract_barcodes.py protocol. My problem with the protocol is that my sequences do not all start with the barcode.
Instead I have e.g. ACACTAG, NACACTGA, NNACAGATC or NNNACAGCTA and so on with N being the basepairs that are in front of the barcodes.
How could I best resolve this? How could I efficiently remove basepairs in front of the barcodes?
I would be really happy about any suggestions!

Heike

Greg Caporaso

unread,
Feb 27, 2014, 10:56:30 AM2/27/14
to qiime...@googlegroups.com
Hello,
We don't currently have support for this, but are planning on adding it. You can follow progress here. Sorry for the inconvenience! 

Greg

Tony Walters

unread,
Feb 27, 2014, 11:11:02 AM2/27/14
to qiime...@googlegroups.com
Heikieh,

Just to clarify the format of your sequences (for the above issue tracking that Greg opened), do your sequences look like this:
spacer of 0-X bases--Barcode Sequence--Primer Sequence--rest of amplicon read


--
 
---
You received this message because you are subscribed to the Google Groups "Qiime Forum" group.
To unsubscribe from this group and stop receiving emails from it, send an email to qiime-forum...@googlegroups.com.
For more options, visit https://groups.google.com/groups/opt_out.

Heike Staudacher

unread,
Feb 27, 2014, 11:35:22 AM2/27/14
to qiime...@googlegroups.com
Hello Tony and Greg,
yes they look like that, with 0-3 nucleotides before the barcode.
I will follow the progress with the link that you provided, thank you!
In think in my case it would be great if I could tell the program that my barcode are the 6 basepairs preceeding the primer sequence. This is not possible yet, right?

Thank you!
Heike

Greg Caporaso

unread,
Feb 28, 2014, 9:35:41 AM2/28/14
to qiime...@googlegroups.com
Correct, that is currently not possible. Sorry for the inconvenience. 

Greg

Sanjeev Sariya

unread,
Jun 9, 2015, 11:04:56 AM6/9/15
to qiime...@googlegroups.com
Hi Dr. Greg,

Pardon me to hack someone else's thread. 
May I know if there's any feature newly added by QIIME for the issue discussed ?

Best,
Sanjeev
--------------------------------------------------------------------------------------------------------------
Reply all
Reply to author
Forward
0 new messages