Question: Assembly of DNA from .bam file
0
gravatar for Triple Nipple
14 months ago by
Triple Nipple10 wrote:

From a .bam file of RNA sequence reads that are aligned to a mouse reference genome, how can I get a new bam file that contains this same genome assembled in reference to mm10 reference genome. If I made things clear I need to get one long read sequence for each location and not many reads like it's the case with this .bam file.

For bases that in some reads (and not all) they seem to be mutated compared to reference, I need to have the most "probable" base, and in bases that aren't available, I also need them to be specified.

I would really appreciate your help and advice on doing this.

Thanks

ADD COMMENTlink modified 14 months ago by swbarnes29.2k • written 14 months ago by Triple Nipple10
2
gravatar for WouterDeCoster
14 months ago by
Belgium
WouterDeCoster44k wrote:

The term you are looking for is a consensus sequence. Take a look at this post: Generating consensus sequence from bam file
By chance, I also saw this today using Alfred for consensus.

ADD COMMENTlink modified 14 months ago • written 14 months ago by WouterDeCoster44k
0
gravatar for swbarnes2
14 months ago by
swbarnes29.2k
United States
swbarnes29.2k wrote:

From a .bam file of RNA sequence reads that are aligned to a mouse reference genome, how can I get a new bam file that contains this same genome assembled in reference to mm10 reference genome.

It's probably simplest to just realign to the genome.

I need to get one long read sequence for each location

You are going to have gaps between genes and between introns. Wouldn't you be better off making a consensus from your transcriptome alignment? (For those RNA sequences where you have enough coverage)

ADD COMMENTlink written 14 months ago by swbarnes29.2k
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: 1203 users visited in the last hour