I'd like to simulate 10% sequencing error using
art_illumina. The simulator doesn't have a parameter that I can just give the 10%, but it has this:
-qs --qShift the amount to shift every first-read quality score by -qs2 --qShift2 the amount to shift every second-read quality score by NOTE: For -qs/-qs2 option, a positive number will shift up quality scores (the max is 93) that reduce substitution sequencing errors and a negative number will shift down quality scores that increase sequencing errors. If shifting scores by x, the error rate will be 1/(10^(x/10)) of the default profile.
wgsim has an option for base error rate. Can I do the same for
Q: Does the "error rate" in the description mean sequencing error? If I set both
qShift2 to 10,
1/(10^(x/10)) == 0.1, where x is 10, does that mean 10 is the value I should give to the simulator?