<p><em>Answer from @cephbirk, converted from comment:</em></p>
<p>What I ended up doing was to choose the read coverage for the lowest covered nucleotide in each codon as each residue's read coverage.</p>

CC BY-SA 4.0
2021-06-25T11:17:26.477Z
2021-06-25T11:17:26.477Z
gringer
<ol>
<li>I used <code>samtools depth</code> to find the per base-pair read coverage over a number of isoform contigs from my <a href="https://github.com/trinityrnaseq/trinityrnaseq/wiki" rel="nofollow noreferrer">Trinity assembly</a>.</li>
<li>I have also conducted a multiple sequence alignment of those isoforms using Clustal Omega DNA aligner.</li>
</ol>
<p>In my final step, I want to assess the per-residue read coverage of the consensus amino acid sequence to demonstrate an idea of confidence in that sequence.</p>
<p>Does anyone have any tips on how to go about doing this and whether there are existing tools for such things?</p>

<phylogenetics><read-mapping>
RNASeq read coverage in protein space?
# RNASeq read coverage in protein space?

(Asked by: username_0 on 2019-06-07)

1. I used `samtools depth` to find the per base-pair read coverage over a number of isoform contigs from my [Trinity assembly](https://github.com/trinityrnaseq/trinityrnaseq/wiki).
2. I have also conducted a multiple sequence alignment of those isoforms using Clustal Omega DNA aligner.

In my final step, I want to assess the per-residue read coverage of the consensus amino acid sequence to demonstrate an idea of confidence in that sequence.

Does anyone have any tips on how to go about doing this and whether there are existing tools for such things?

**Tags:** phylogenetics, read-mapping

---

## 1 answers

### Answer by username_1 on 2021-06-25. Score: 1 | Accepted answer

*Answer from @cephbirk, converted from comment:*

What I ended up doing was to choose the read coverage for the lowest covered nucleotide in each codon as each residue's read coverage.


---

