Hello!
I have fastq files from HiSeq2500 that all the reads are tagged with our additional custom indexes along with TruSeq indexes. So, my reads sort of look like below.
TruSeq adapter w/ index NNNN our index NNNN read
After split reads by TruSeq index, I want to split again with our indexes.
Is there a tool that can separate reads by a certain index sequence at certain position?
(i.e. if a sequence from 5th base to 10th base is 'AGTGC', then it is index #1)
I hope my explanation is enough for you. Thank you!
1 answer
Hello!
You can try out the Checkout routine from our MiGEC package, which was specifically designed for such kind of tasks. It can also output those "NNN.." sequences to output read header, in case those are the part of UMI barcode (see this overview article for details).
PS In your case you should provide barcodes.txt file containing
S1 NNNNatgtgcATTGatgcNNNatgc
S2 NNNNatgtgcGATTatgcNNNatgc
S3 NNNNatgtgcTGATatgcNNNatgc
...
where ATTG, GATT, TGAT,... are your barcodes, lower case characters indicate primer sequence (those regions allow mismatches) and NNNN is the degenerate region. IUPAC ambiguous could also be used.
•
written
by
mikhail.shugay
Log in to answer this question.
-
Removing overrepresented sequences in paired end RNA-seq
written by LayneSadler
-
FastQ files: Double demultiplexing
written by compuTE
-
How to demultiplex SmartSeq2 run where settings were set to one index instead of two?
written by RhiRhiO •
-
Trim NEBnext Adapters with trimmomatic - Same as TruSeq adapters?
written by alexgoro13 •
-
Help Troubleshooting FASTQC
written by mac03pat •
-
bcl2fastq2: how to correctly use the --use-bases-mask for different sequencing methods by Illumina ?
written by 乙
-
Which truseq trimmomatic adapters file to use when removing truseq adapters?
written by salamandra
-
Please advise me on PE sequencing reads prepared by NEBNex kit and Hiseq2000
written by seta
-
NIPT data analysis - TruSeq Indexed Adapter sequence reads and GC bias
written by Ning •
-
Simple shell command for extracting indexed reads from fastq?
written by mbk0asis
Have you discovered this answer.
Demultiplex Illumina With Barcodes On Identifier Line
is that what you were meaning.