thanks GenoMax
• 0 views
•
link
written
by
serene.s •
Subsampling metagenomics reads via reformat.sh
Hi Guys! I am using reformat.sh for subsampling my metagenome samples but I wanna know what does the parameter "sampleseed" do to the reads? Is it alright if I just use the samplereadstarget without using the seed value?
command used: reformat.sh in=file_R1.fastq.gz in2=file_R2.fastq.gz out=file_subsampled_R1.fastq.gz out2=file_subsampled_R2.fastq.gz samplereadstarget=10889622
Thanks in Advance Saraswati
• 1,467 views
•
link
written
by
serene.s •
1 answer
sampleseed is used as a seed if you want to get deterministic i.e. identical datasets for multiple runs of reformat.sh. Otherwise you may get a different set of reads everytime. If you are not interested in getting an identical dataset (which is why you may be sampling) then do not set an identical seed value for multiple runs.
Log in to answer this question.
More posts like this
-
bbduk paired-end RNA-seq trimming command
written by predeus
-
FastQC report of the metagenomic data
written by serene.s •
-
bbduk trimmed out too many sequences
written by AfinaM •
-
extract pathogens from metaphlan3 results
written by serene.s •
-
Dynamic string replacement issue in parallel for repair.sh bbmap
written by arunprasanna83 •
-
bbduk - trimming ITS primer sequence
written by AfinaM •
-
Subsampling fastq.gz reads using reformat.sh changed all the quality scores?
written by yinjiani900129 •
-
How to parallelize a rule in Snakemake ?
written by vin.darb
-
randomreads.sh simulated reads dont map to genome
written by grant.hovhannisyan
-
Tadpole error correction error killed:9
written by djackso •
