See approx. line 165:
N<-min(max(data$A+data$C+data$G+data$T+data$N),50) ## Total number of reads
Where it now says 50 should be changed to a parameter.
Continue using 50 as the default value.
Discussion:
... coverage >50 should almost never be seen in a normal MSG MseI experiment with a Drosophila sized
genome.
This said, it would be nice to relax this for applications to reduced representation libraries or smaller
genomes ...
See approx. line 165: N<-min(max(data$A+data$C+data$G+data$T+data$N),50) ## Total number of reads
Where it now says 50 should be changed to a parameter.
Continue using 50 as the default value.
Discussion:
... coverage >50 should almost never be seen in a normal MSG MseI experiment with a Drosophila sized genome. This said, it would be nice to relax this for applications to reduced representation libraries or smaller genomes ...