RetroSeq icon indicating copy to clipboard operation
RetroSeq copied to clipboard

Defaults in -call

Open MaximilianStammnitz opened this issue 5 years ago • 0 comments

Hi @tk2,

Many thanks for maintaining RetroSeq - great tool!

After some performance tests, I just realised that the default parameter settings in your retroseq.pl -call raw description are actually not in line with the values set (currently: -reads = 10 and -q = 20). Maybe you could quickly replace the source script's lines 199-201 accordingly?

`[-depth Max average depth of a region to be considered for calling. Default is $DEFAULT_MAX_DEPTH.]

[-reads It is the minimum number of reads required to make a call. Default is $DEFAULT_READS.]

[-q Minimum mapping quality for a read mate that anchors the insertion call. Default is $DEFAULT_ANCHORQ.]`

Cheers, Max

MaximilianStammnitz avatar Oct 26 '20 16:10 MaximilianStammnitz