RNASeq read coverage in protein space?

Bioinformatics Asked by CephBirk on August 30, 2021

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

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

Add your own answers!

Related Questions

Genome scaffolding

3  Asked on December 22, 2020


Time for running ADMIXTURE analysis

1  Asked on December 17, 2020 by iriel


Output the linear predictor from a stratified cox model?

1  Asked on December 17, 2020 by nienke


FastQC and Trimmomatic in Galaxy?

0  Asked on December 16, 2020


What type of Protein ID is this?

1  Asked on December 14, 2020 by firingam


How to find all WGS assemblies accessions of a species

1  Asked on December 12, 2020 by oren-milman


plot transcript expression vs length

0  Asked on December 10, 2020 by user3377241


Filtering pileup from site lists

0  Asked on December 8, 2020 by eliran-turgeman


How to modify DNA evolution model to fit actual data?

2  Asked on December 8, 2020 by anthony-guterres


Ask a Question

Get help from others!

© 2023 All rights reserved. Sites we Love: PCI Database, MenuIva, UKBizDB, Menu Kuliner, Sharing RPP, SolveDir