1
0
Entering edit mode
3 days ago
eozcan • 0

Hi,

I have 9 genomes, I would like to produce a metagenome like distribution using randomreads.sh. I concatenated genome fasta files in one reference file. Then, ran as below.

However, I would like to know if there is any way I can define the abundances for each genome based on qPCR results and then produce the reads accordingly. Would there be a way to produce reads from each genome separately with the absolute abundances?

Thank you!

randomreads.sh metagenome bbmap • 271 views
0
Entering edit mode

Can you not generate the reads independently and mix them as needed? You can then use shuffle.sh to mix the reads randomly giving you a representative metagenome.

0
Entering edit mode

I can generate the reads independently. But does shuffle.sh have a function of indicating the abundances? I didnt see any!

0
Entering edit mode

I was thinking that you would add known amounts of reads together based on your needs and then simply shuffle them so they represent a mixed metagenome.

0
Entering edit mode

Thats my question! How do I add known amount of reads? I will produce randomreads.sh from each genome, then shuffle them , right? But in which step exactly I am adding the known amount of reads for each genome.

0
Entering edit mode
2 days ago
GenoMax 107k

After you generate a certain number of reads for each genome using randomreads.sh. You can then use reformat.sh

Sampling parameters:

samplerate=1            Randomly output only this fraction of reads; 1 means sampling is disabled.
sampleseed=-1           Set to a positive number to use that prng seed for sampling (allowing deterministic sampling).
samplebasestarget=0     (sbt) Exact number of OUTPUT bases desired.


to select desired number of reads from each genome (e.g. 1 M from genome_1, 1.2 M from genome_2 etc). cat the sampled genome files together and then shuffle.sh that file.

0
Entering edit mode

Ah thank you! That is what I was looking for!

0
Entering edit mode

If you generate the reads individually remember to turn the metagenome option off.