Question: Demultiplexing And Splitting Miseq Fastq
2
gravatar for peris
6.4 years ago by
peris120
United States
peris120 wrote:

I got the following files from a service provider while doing MiSeq sequencing. Can someone help me in demultiplexing and splitting this files? I have used 24 independent indexes for indexing.

Read1:

@130110_A5T23:1:1101:15646:1320:1#0/1
NAGTAAGATACAGTCTATCGGGTTTAAGTTATACAACATAGTACAGTACATTCATACCTACCTCTGCAATTAAATTTGGCGGTGTCATAATGTCTTTCAGCACAACTAAAAGAAAAGTTTAAAAGTGATATAAAGTTATCTTTTGCACTT
+
#>>AABFFFFFFGGGGGGGGGGCFGHHHHHHHHHHHHHHHHGFHHHHHHHHGHGHHHHHHGHHGHHHGHHHHHHHHHEHGGGEEFGHHHHFHHHGHHHHFHHHHGGGHHHHGGHHGHHHHHHHGHBGHHHHHHHHHFHFHHHHHGHFHHH
@130110_A5T23:1:1101:17282:1320:0#0/1
NGGGAGGATGGCTTGAGCCTAGAGGGTTGAGGCTGCAGTGAGCTGTGAGCATGCCATTGCACTCCAGCCTGGGCAACAGAGTGAGGCTCTGTCTCAAAAAAAAAAAAAAAAAAAAAAGAAATAAAATATAAGGGACATAAAATAAAGCAA
+
#>>AAA1A1B>F1A1FAFGCG0GBFEAAEHGBFCFFHHA1FEGGFFGGFFHHHGHHHHGHHHHHHHHHGHGGGHHGHHHFHGHGFGGHHHHHHHHHFGHHGGGEGGGGGGGGGGGCC.<<G0DGGFC:CCC0CC.::00:<00:9C####

Index

@130110_A5T23:1:1101:15646:1320:1#0/2
NTTGTATC
+
#>>AAAFF
@130110_A5T23:1:1101:17282:1320:0#0/2
NTTAGTCT
+
#1>>AAFF
miseq • 8.8k views
ADD COMMENTlink modified 5.2 years ago by yinghua10 • written 6.4 years ago by peris120
2

Have you just tried using FastqMultx from ea-utils (http://code.google.com/p/ea-utils/wiki/FastqMultx)? You'll have to make a file specifying index:sample relationship, but it seems otherwise straightforward and supports paired-end demultiplexing.

ADD REPLYlink written 6.4 years ago by Devon Ryan94k

+1 I didn't know about this tool. This looks like the best way to go.

ADD REPLYlink written 6.4 years ago by Dan D7.0k
1

So all you have available is two FASTQ files, one for the index read and one for the single read, is that correct? Because if you have the basecall data we can run CASAVA. If these two files are all you have I can code up something for you Friday morning.

ADD REPLYlink modified 6.4 years ago • written 6.4 years ago by Dan D7.0k

Hi Deedee, Thanks for you reply. Yah you are right, I have only the fastq files, not the raw files. And this is actually a paired end run. I have Read1, Read2 and Index ;these three fastq files. It will be of great help if you can write a code for me as I dont know to write codes for such conditions. Thanks and regards

ADD REPLYlink written 6.4 years ago by peris120
1
gravatar for Matt Shirley
6.4 years ago by
Matt Shirley9.2k
Cambridge, MA
Matt Shirley9.2k wrote:

I came up with a solution for this problem a while back, although I didn't know about FastqMultx, and might have done something different now

ADD COMMENTlink written 6.4 years ago by Matt Shirley9.2k
1
gravatar for yinghua
5.2 years ago by
yinghua10
United States
yinghua10 wrote:

I recently wrote a small code for our miseq machine to demultiplex I1, R1, R2 reads assuming barcode reads are in the I1 file. It does a little more than exact matching. For barcodes reads not exactly matched with designed barcodes, which normally amount about 10% of total reads, edit distances with designed barcodes are calculated and a cut off is set by both user and the minimum edit distance among designed barcodes.

https://github.com/yhwu/idemp

ADD COMMENTlink written 5.2 years ago by yinghua10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2260 users visited in the last hour