View Single Post
Old 01-26-2014, 09:45 PM   #1
foolishbrat
Member
 
Location: South East Asia

Join Date: Nov 2008
Posts: 44
Default Adapter trimming NEBNext Library / MiSeq

I have a 51 single-end reads generated with MiSeq using NEBNext Multiplex Oligos for Illumina.

The sample sheet looks like this:

Code:
IEMFileVersion,4
Investigator Name,FB
Experiment Name,WT10104
Date,11/27/2013
Workflow,GenerateFASTQ
Application,FASTQ Only
Assay,TruSeq Small RNA
Description,
Chemistry,Default

[Reads]
51

[Settings]
ReverseComplement,0

[Data]
Sample_ID,Sample_Name,Sample_Plate,Sample_Well,I7_Index_ID,index,Sample_Project,Description
HS130333-1,,,,RPI3,TTAGGC,,
HS130333-2,,,,RPI4,TGACCA,,
HS130333-3,,,,RPI5,ACAGTG,,
The Primer index manual can be found here.

Sor for HS130333-1 file, according to the manual above the primer/adapter with index is:
5 ́-CAAGCAGAAGACGGCATACGAGATGCCTAAGTGACTGGAGTTCAGACGTGTGCTCTTCCGATC-s-T-3 ́

The document indicated that the expected index primer sequence read is TTAGGC which is the reverse complement of GCCTAA.


My question is if I use `trim_galore` or `cutadapt` to trim the data, what is the parameter -a I should use?

Is it the whole sequence above? Or first 5 ́-CAAGCAGAAGACGGCATACGAGAT?
Or GTGACTGGAGTTCAGACGTGTGCTCTTCCGATC-s-T-3 ́? (and what is 's' means here)

Or the reverse complement of the each of above?
foolishbrat is offline   Reply With Quote