I have two folders and each of them contain the same file names.
ls -1
10_S0_L001_R1_001.fastq.gz
10_S0_L001_R2_001.fastq.gz
11_S0_L001_R1_001.fastq.gz
11_S0_L001_R2_001.fastq.gz
Is there a way to extract from each dataset the flowcell name and use it as unique filename?
Thank you in advance.
1 answer
If you look inside the files then you should see headers that look something like
@HWI-EAS209_0006_FC706VJ:5:58:5894:21141
The flowcell serial is embedded in the header (e.g. FC706VJhere). Problem with using that for file names is all samples would have the same name unless you cat it to sample id with something like Sample_01_FC706VJ
You can use standard unix tools (such as cut, awk, tr, grep, etc) or non-standard ones (e.g.,bioawk) to extract metadata from your fastqs. Like genomax said, you'd need to extract some other identifier(s) in order to make your filenames unique.
Log in to answer this question.
-
Rename filenames based on list
written by genomes_and_MGEs •
-
How to merge folders that have similar (not the same) names, in a folder keeping as its name the co…
written by Mania •
-
How to move multiple files to multiple folders in linux/mac?
written by anikcropscience •
-
How to align multiple fastq files with Hisat2 ?
written by sarai •
-
Finding common and different genes between species
written by Ric
-
Print based on character strings
written by selplat21 •
-
Use PLINK2 to extract the specific variants
written by wangdp123
-
Renaming Bcftools lsec for multiple runs
written by Seigfried
-
Merging two fastq.gz files
written by tcf.hcdg
-
Parsing flowcell from hiseq 2000 read name
written by steven.davis •
