This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Splitting fastq files

Hi there,

I was just wondering if anyone could offer any advice on splitting two merged fastq files (R1 and R2) into one per-sample fastq files? I've downloaded several biosamples from SRA via ftp, but they are merged into one file and I am unsure how to split them. Thanks!

fastq ncbi sra
• 2,402 views • link
updated by GenoMax GenoMax • written by nicole.kavanagh • nicole.kavanagh

How can one distinguish the two samples? Do they have individual indexes?

Reading your post again it seems like you are asking about interleaved data files where each R1 and R2 reads are present next to each other e.g. R1_1,R2_1,R1_2,R2_2,R1_3,R2_3 etc. If that is the case you can separate the interleaved reads using reformat.sh from BBMap suite.

reformat.sh in=interleaved.fq out1=R1.fq out2=R2.fq

If the two read files are "merged" by copying them end to end (R1_file_followed_by_R2_file) then you may need to use split by counting number of records (*4 lines).

• 0 views • link
written by GenoMax GenoMax

Thanks so much for your advice. Apologies, I wasn't very clear in my question. I have two files, the first file contains the R1_1, R1_2, R1_3, R1_4..etc and the second file contains R2_1, R2_2, R2_3, R2_4, so I don't think they fall into interleaved data or merged files. How would I go about splitting them into individual files? An example of the text contained within the fastq is as follows:

@SRR3138122.104397 104397/1
CAGCATGAGTGGCTATATTCATATATTCATTGACTCGATCTTCATTACTGTCAATAGCCAGTACTTCTTGTCCAGCTTCGATCAGTGTACGGCAGATGCTTCCGCCAAATCGTCCAAGGC
+
CCCCCGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGFGGGGGGGGGGGGGGGGGGGGGGGCFGGGGGGGGDFGGGGFGGGGGGGGGGGGGGGGGGGGGGG
@SRR3138122.104398 104398/1
TAGACATGCCGACAAAGAACCCTAAAACACGTAAGAAACCAAAAGTCAAAATCGGTGACATGGTTCGCTGCGAAGCAGAAGGGTTCATCTATCCGTTTCGTGGATATGTAGAACATGTTTATGATCACTCAGCAATCATTCGCATTGAAAACACGATGGAATGCGACAAGTGGTTAGCGAAAAGCAAAGATAATTTAGCAGTGGCTCGATTGGTGGATATGGAACTAATC
+
A6A@@FFFGFFGGGGFEGGFFEFGFGGGGGGGGGGGGFGGGGGGGGDFGGGF@EEG:FGGGF<CE@F@FGGGGGGGGCFDGGC<EFFGGGGGFGDFGGGGGFGGFFGGGGFAFFFFFFFFGCFFGFGGGGGCFF,5AFGDFG>FEGGGFGGEGGGGCGGGGGGGGGGGGGEDDEC=<FGGGGGA<CFFGGGGGGGEGGGGDBECE@7CECCF,?EEFCGGG9BECE9FCB
@SRR3138122.104399 104399/1
CTTTGTGGCAAAAGCGGAGAATCGCCATTGCTCCCTTGAAAACGTCTTTCCTGATTCCATTGTGTCTTCCCATGTCGGGTAAAATATAGTTCCATACGCAGCAAACTCCTGTCTTTCTTTTATTTTACTCTATAATAAACTGTCTAATCCAGCTTGTCATTTTCATAGTGATGCGTCAATTCTAAATTAGGGAATCCTTGCTGGCGAAGTGCTTCATACACTATCATTGCAGCTGTATTCGAT
+
CCCCCGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGGFFGGGGGGGGGGGGGGGGFGFGGGGGGGGGGGGGGGGGGFFGGGGGGGGGGGGGGGEGCFEFCFGGCGCGGGGGGGGGGGGGFGGGGGGGGGGGGGCFGGGGGGGGGGGGGGGGFGGG@FGFGGGGGGGGGGGGGEC>CCEFGFGGGFGGGFGGGGGFFGGGGGCACDE62,
@SRR3138122.104400 104400/1
CGTTGTTTCTGCTTTCTTGGAAAGACACCCTGATTTTGAGAGGATTGATGT
+
CCCCCGGGFGFGFGGAFG<FCFGGFGGGGG8FDFGGGD@FECF@FG?@FFF
@SRR3138122.104401 104401/1
GCTTTGAAAAATATTATTGCTATTGGTGCTGGCGCTATTCACGGCTTAGGGTTTGGCGATAATGCAAAGGCAGCAATCATGACGCGCGGATTAGCTGAGATCAGCCGTTTAGGTGTAGCAATGGATGCCAACCCGCTGACCTTTATCGGATTAAGCGGAGTAGGCGATCTAGTCGTTACTTGCACAAGTGTTCATTCTAGGAACTGGCGCGCAGGCAAATTGTTAGGACAAGGGCAGCCCCTAGAAG
+
A-6CCFGGF9AFEFGE9F9C,<<,<CEF,C,CFGGGFDF9@FFGGG+FC,CFGGC7@7CBFEGGCFC<AFCDEEEFCFFGFAGG@:F:E+CEEGGFGGGFF<F9FGGG+BFB8<=?BFF<FFC8<EF,?E=<+BC>:B@@8F<FB9AB7++33=>,3+5CFCACCF@*>:,>FFFFD:*;FCEBAFF7,?FCGFC;;;@,,6;<EC8C:838*3:55**<+2?;9>6++=88/;CE*;*;*:8:0<F

Thanks for your help!

• 0 views • link
updated by GenoMax GenoMax • written by nicole.kavanagh • nicole.kavanagh

There is no need to split anything. It looks like you have standard paired-end files (2 files per) for one sample. You should have 1 pair for SRR3138122 , another pair for SRR#### (next accession) etc.

• 0 views • link
written by GenoMax GenoMax

0 answers

No answers yet.

Log in to answer this question.

More posts like this