split fastq by @SEQID
Hi all,
I have a couple of fastq files containing reads starting with different name like: @HWI-ST865:463:C7C8KACXX:2:2316:21016:100943 1:N:0:TAAGGCGA @HWI-ST1178:227:C7C95ACXX:7:1101:1581:2125 1:N:0:TAAGGCGA
My question is: how can I split them in two parts? I tried to use some tools like fastx_toolkit but I cannot create a proper barcode file Is there any easy way to do that such as a grep command, cause i also tried with grep but i got an output containing only the first line of the reads and missed the other three
Thank you in advance!
1 answer
Log in to answer this question.
More posts like this
-
Why extracting reads from BAM to fastq produces small amount of reads?
written by artemd •
-
BBDuk.sh and BBDuk2.sh barcodefilter does not filter the reads
written by chilifan
-
Separate the reverse sequence from R1 and corresponding forward seq from R2
written by Sajan Raju •
-
SolexaQA++ lengthsort - c output
written by gtho123
-
demultiplex a dataset when you have barcodes as a separate fastq
written by IP
-
[solved] Read a file once and get next lines that match a list of pattern with awk ?
written by Franck8413 •
-
Removing barcode and adapters?
written by manekineko
-
Adding barcodes to forward and reverse reads after sequencing
written by cmogren •
-
Java Error "Concurrent.Executionexception" While Running Trimmomatic V0.32
written by arnstrm
-
Parsing Fastq Files
written by deepthitheresa •
So you want to split on the "HWI-STXXX" bit? Or every unique ID should be a different output file?
Split them by "HWI-STXXX"
Probably a mixed data set. Of late some submitters have been merging data from multiple flowcells/machines into one file for SRA submission (beats me why they do it) and this could be a case of that sort.
yes, this is exactly the case, but it was accidentally done. 2 different persons sequenced the same sample in 2 different sequencers without being aware of and then they decided to merge the outputs
Hahahah :D Awesome. Was it the exact same biological sample? If so, is the data publicly avalible? If so, i'd be interested in looking at the QC data. See how much of an effect sequencing machine/etc really plays on the downstream statistics.
Yes exactly the same sample but different capture processes with the same kit (exome sequencing). Unfortunately data are not publicly available....my boss will kill me if i do that!sorry....hahaha... Anyway, I suppose this fact will affect the analysis anyhow, cause the capture process was different despite the fact they used the same kit and protocol. You know sometimes things are working almost 100% and sometimes not.
No worries man - there's more than enough data to go around :) And yeah, maybe a different capture process will highlight different exons better, who knows, it might not be a waste at all!