I used samtools depth to find the per base-pair read coverage over a number of isoform contigs from my Trinity assembly. 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?

RNASeq read coverage in protein space?

Bioinformatics Asked by CephBirk on August 30, 2021

I used samtools depth to find the per base-pair read coverage over a number of isoform contigs from my Trinity assembly.
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?

alignment read mapping

One 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.

Correct answer by gringer on August 30, 2021

RNASeq read coverage in protein space?

One Answer

Add your own answers!

Related Questions

Genome scaffolding

If a gene is expressed at a level of 1/1200 compared to the average gene, how is probability 50:50 that we have a read mapped to it?

Discrepany between AutoDock Vina and AutoDock Tools

How many single copy gene are needed to construct the phylogenetic tree?

Time for running ADMIXTURE analysis

Output the linear predictor from a stratified cox model?

FastQC and Trimmomatic in Galaxy?

Measure the purifying selection for certain taxa along a phylogeny

What type of Protein ID is this?

How to find all WGS assemblies accessions of a species

What is ‘k’ in sequencing?

plot transcript expression vs length

How to translate amino acid sequences to Nucleotide sequences

How to find/build the evolutionary history of a protein from its sequence?

Practical use of RNA structure

Is there a Python program that is able to show convergent and divergent evolution on a phylogenetic tree?

Plot information stored in dataframe within a tree (ggtree)

Filtering pileup from site lists

How to modify DNA evolution model to fit actual data?

ConsensusClusterPlus for sample classification

Ask a Question