Entering edit mode
7.2 years ago
ram.mainali
▴
20
I produced list of kmers from WGS data which have coverage above a certain threshold or cut-off by using GenomeTester4 glistcompare. The command is: GenomeTester4/bin/glistcompare inputfile.list inputfile.list -i -c 15 -o output
Besides setting coverage above a certain threshold (cut-off), here in this case 15, I would like to set maximum coverage or cut-off for example 40 so that I will have list of k-mers with coverage 15 to 40. Could you please suggest me what I should do in this regard. Thank you
