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 on December 23, 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.

Answered by gringer on December 23, 2021

RNASeq read coverage in protein space?

One Answer

Add your own answers!

Related Questions

Downloading SRA Files from AWS

NCT search in PubMed via Entrez (python)

Comparing groups of 16S rRNA sequences – how to?

Relationship between the instantaneous rate of change of a continuous time Markov model and the rate of a global clock model?

load the phenotype data for ballgown

using snakemake shadow rules to store temp files on local nodes

Need help with Pairwise Alignment module in iterating over the alignment

Software for taxonomic assignment?

Any user friendly way to find rare mutations in whole genome raw?

marge the matrix from computematrix of deeptools

How to get rows with similar values in two different columns using command line?

How to quantify similarity of genomes and find differences in set of S aureus genomes?

Python equivalent to R GRanges and IRanges

Detecting Allelic Imbalance

Mapping statistics from bam file using bbtools and sambamba

Call NCBI Multiple Sequence Alignment Viewer with target URL?

annotation using ChIPseeker package error

Preprocessing Affymetrix SNP Array Matrix

Get protein names corresponding to PDB ID

Details of DESeq2 modeling a batch effect

Ask a Question