Problemastic reads quality score distribution by dwgsim
Hi there,
I was using dwgsim to do illumina reads simulation for my study. The command I used is like this >>
dwgsim -e 0.0001-0.01 -E 0.0001-0.01 -d 100 -s 75 -1 100 -2 100 -C 2.5 -r 0 -R 0 -y 0.01 -c 0 -H hap1.fasta sim.hap1.dwgsim
So I should expect a decreasing quality score from 40 to 20 along the simulated reads on average. But the quality scores I got from the simulated reads revealed a very different picture>>
What might be the reason for this? Thank you very much!
Best,
Jia-Xing
I was using dwgsim to do illumina reads simulation for my study. The command I used is like this >>
dwgsim -e 0.0001-0.01 -E 0.0001-0.01 -d 100 -s 75 -1 100 -2 100 -C 2.5 -r 0 -R 0 -y 0.01 -c 0 -H hap1.fasta sim.hap1.dwgsim
So I should expect a decreasing quality score from 40 to 20 along the simulated reads on average. But the quality scores I got from the simulated reads revealed a very different picture>>
What might be the reason for this? Thank you very much!
Best,
Jia-Xing
Comment