<p>All transcripts should always be returned in the output, whether there's evidence supporting their expression or not. This is for the sake of simple convenience, since typically one would be merging these values across samples into a matrix. Missing values will then tend to become either <code>NA</code> or cause an error, when what would be appropriate would be to have values of 0. This also facilitates calculations such as "what's the median number of reads per transcript?" or "what percent of transcripts have an expected reads count of at least X".</p>

CC BY-SA 3.0
2018-03-13T07:50:56.173Z
Ok thanks. Only the transcripts with reads support should be then consider as 'expressed' right? And what about read support value between 0 and 1? Shall I consider them to be present or not?
aechchiki
CC BY-SA 3.0
2018-03-13T07:53:10.067Z
Normally one would use something like zFPKM to come up with some "expressed" threshold. Note that values above 1 don't necessarily mean that there is even a single read uniquely assignable to a given transcript.
Devon Ryan
CC BY-SA 3.0
2018-03-12T14:06:41.533Z
2018-03-12T14:06:41.533Z
Devon Ryan
<p>I quantified some samples using <a href="http://salmon.readthedocs.io/en/latest/salmon.html" rel="noreferrer">Salmon</a>. According to the documentation of <a href="http://salmon.readthedocs.io/en/latest/file_formats.html#fileformats" rel="noreferrer">output format</a>, the last column of the 'Quantification File' represents the number of reads that supported the given transcript, namely:</p>

<blockquote>
  <p>NumReads — This is salmon’s estimate of the number of reads mapping to each transcript that was quantified. It is an “estimate” insofar as it is the expected number of reads that have originated from each transcript given the structure of the uniquely mapping and multi-mapping reads and the relative abundance estimates for each transcript.</p>
</blockquote>

<p>In some case, I had in the output some transcripts that had <code>NumReads=0.000000</code>. Does this mean that I should discard them from the list of transcripts obtained in that sample? How comes that those transcripts end up in the output file in first place, since they don't have any reads support and the quantification was done using the reads themselves?</p>

<rna-seq><quantification><salmon>
Salmon output: transcripts quantified with zero reads support
# Salmon output: transcripts quantified with zero reads support

(Asked by: username_0 on 2018-03-12)

I quantified some samples using [Salmon](http://salmon.readthedocs.io/en/latest/salmon.html). According to the documentation of [output format](http://salmon.readthedocs.io/en/latest/file_formats.html#fileformats), the last column of the 'Quantification File' represents the number of reads that supported the given transcript, namely:

> NumReads — This is salmon’s estimate of the number of reads mapping to each transcript that was quantified. It is an “estimate” insofar as it is the expected number of reads that have originated from each transcript given the structure of the uniquely mapping and multi-mapping reads and the relative abundance estimates for each transcript.

In some case, I had in the output some transcripts that had `NumReads=0.000000`. Does this mean that I should discard them from the list of transcripts obtained in that sample? How comes that those transcripts end up in the output file in first place, since they don't have any reads support and the quantification was done using the reads themselves?

**Tags:** rna-seq, quantification, salmon

---

## 1 answers

### Answer by username_1 on 2018-03-12. Score: 4 | Accepted answer

All transcripts should always be returned in the output, whether there's evidence supporting their expression or not. This is for the sake of simple convenience, since typically one would be merging these values across samples into a matrix. Missing values will then tend to become either `NA` or cause an error, when what would be appropriate would be to have values of 0. This also facilitates calculations such as "what's the median number of reads per transcript?" or "what percent of transcripts have an expected reads count of at least X".

> **Comments:**
> **username_0**: Ok thanks. Only the transcripts with reads support should be then consider as 'expressed' right? And what about read support value between 0 and 1? Shall I consider them to be present or not?  
> **username_1**: Normally one would use something like zFPKM to come up with some "expressed" threshold. Note that values above 1 don't necessarily mean that there is even a single read uniquely assignable to a given transcript.  

---

