where should I found the meaning of results' header

55 views
Skip to first unread message

gnili...@gmail.com

unread,
Mar 31, 2020, 10:54:02 PM3/31/20
to CTK User Group

Dear author,
    I have been use CTK to my analysis, but parts of header in the results are hard to understand.
   Such as, results of `CIMS.pl`, what's the meaning of 'count(>m,k). Apparently, it's not the reads count of the mutation site.
   So, where should I found the meaning of results' header ?
   Thank you for your jobs.

Huijuan Feng

unread,
Apr 5, 2020, 10:56:57 AM4/5/20
to CTK User Group
Hi,

Thanks for choosing CTK.

The count(>m, k) is the number of clustered mutation sites with k tags and have >=m mutation frequency. We calculated to cumulative number to estimate the FDR. 
For example, if the site has k=15 and m=2, then count(>m, k) is the total sites number with m>=2 and k=15. You can get the same results using this command:

awk '{if($7>15&&$8>=2)print}' cims.output | wc -l


For more information, you can refer to the examples here:

Please let me know if you have more questions.
Thanks!
Reply all
Reply to author
Forward
0 new messages