This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Demuxbyname.sh dual index problem

Hi, im trying to demultiplex a fastq file with demuxbyname.sh. Im using the following code:

demuxbyname.sh in=1-1_S1_L001_R1_001.fastq in2=1-1_S1_L001_R2_001.fastq out=out_%_R1.fq out2=out_%_R2.fq names=pool1.txt substringmode=t  delimiter=:

My samples are dual index (barcode in forward and barcode in reverse file). So in the txt file i have the 2 barcodes per ID in the first line (something like ACTCCCG+GTCGGG) and the other ones below like in a column. The fastq files are like the following:

@M03485:29:000000000-J4R9F:1:1101:16497:1353 1:N:0:NCCGGAGA+TATAGCCT
GCCACTAGATCCAAGAGATCCGNTGTTGNAAGTTGTATTTCCAAAATAACCCAAATAATACATTCCATTTAAGAGTTTATAAATAAAATGAAGATCAAGGGGACAATAATAAAAAAGAAAGAGAATACACTCTTCCCTCTTAAACCCTTCTTGACTTCACAAAGTTCACAGGTGTGTAAAGGATGAAAAGTTAATTA$
+
3AA33FCFBFFFFGGGGAGGGG#BBFFG#B2AFFHHHHFFHFBHHGFHHHEEFHHHHHHFFHEGGHFGHHHHHHF5FGHHHHHGHGHDF3@FHFGFHHBGGCEEGHHHHHHHFFFE?FHHHHHHHHFBGHHHHHGGCHCCFFHHGHGGHHGHFHFFHHHFHGCFFFGHHHHHFFDDGHFBFHFGHHHGHHDFHFGFG$
@M03485:29:000000000-J4R9F:1:1101:16985:1388 1:N:0:NCCGGAGA+TATAGCCT
ACTGACTGATGCAAGAGATCCGNTGTTGNAAGTGAACCTTGTTACGACTTTTACTTCCGTATAGTAAGTAGATCGGAAGAGCACACGGCTGAACTCCAGTCACTCCGTAGAATCTCGGATGCCGTCTTCTGCTTGAAAAAAATCAAATATTATGCTTCGCAGATTTCAAACCCGGACCCTTTGTTACGACTTTTCCT$

The program dont show any error, i just dont get any output from the process

Set INTERLEAVED to false
Input is being processed as paired
Time:               3.899 seconds.
Reads Processed:    6470118     1659.47k reads/sec
Bases Processed:    1623999618  416.53m bases/sec
Reads Out:          0
Bases Out:          0

IN another post, the issue was fixed with substringmode=t but not in this case. Any help is appreciated.

demultiplex demuxbyname.sh bbmap dual-index

0 answers

No answers yet.

Log in to answer this question.