RetroSeq
RetroSeq copied to clipboard
Defaults in -call
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