Raw read counts in GMGC10.sample-abundance.tsv.xz is not integer?

35 views
Skip to first unread message

songcan chen

unread,
Feb 7, 2023, 4:38:55 AM2/7/23
to gmgc-users
Why does 'raw' column in GMGC10.sample-abundance.tsv.xz contain not-integer values? According to the explanation in the GMGC webpage (https://gmgc.embl.de/download.cgi), 'raw' column gives the number of reads mapped to the unigene, which should be an integer. Yet, this value for some unigenes is less than 1 or not integer (see the example below). Why is that?


Example: 

GMGC10.204_015_426.DMSA1 (unigene)        SAMEA1906561 (sample)    0.10801909817490765 (scaled)     0.5 (raw)     0.0 (raw_unique)     0.0 (normed10m)



Luis Pedro Coelho

unread,
Feb 7, 2023, 7:25:15 AM2/7/23
to songcan chen, gmgc-users
Dear Songcan et al.,

This column (raw) includes fractional assignment of reads when there are multiply-mapped reads, using the dist1 strategy from NGLess, see


So, it does not normalize by either the number of reads in the sample nor the length of reference gene, but it can still be fractional

HTH
Luis
--
You received this message because you are subscribed to the Google Groups "gmgc-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to gmgc-users+...@googlegroups.com.

Reply all
Reply to author
Forward
0 new messages