Question: How to visualize multiple sequence alignment and read coverage simultaneously
0
gravatar for callamartyn
6 weeks ago by
callamartyn10
callamartyn10 wrote:

Hello,

I have rna-seq psuedoalignments for a genome with a lot of paralogs. I suspect there are some errors/complications with multi-mapping reads. Focusing on one transcript, I have a multiple sequence alignment of all its paralogs and can easily see their overlaps and homology. Separately, I have a bam file for my transcript and can see where the reads are aligning. I'm trying to get a sense for which of the reads are in a conserved region across the paralogs and which (if any) are in a region unique to the transcript of interest. But to do so I have to keep flipping back and forth between the two and counting positions. Does anyone know of a way to view both in one place? Like is there a way to convert the MSA file (currently a clustal .aln file) to a bam file with one of the sequences as a reference? Or a bunch of pairwise alignments? Or do I need to convert the information I want from the MSA alignment to a bed file and add that? Any suggestions would be greatly appreciated!

igv alignment • 92 views
ADD COMMENTlink written 6 weeks ago by callamartyn10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1941 users visited in the last hour