Question: Extracting alignment coordinates from SAM file
0
gravatar for dasakan
10 months ago by
dasakan0
dasakan0 wrote:

Hi, I have a SAM alignment file. Alignment was done using BWA MEM between reference genome (in scaffolds, not a finished genome) and assembled transcripts (as contigs). I am interested to extract the genome coordinates of each contig/split contigs (as these are transcripts, so while aligning to the genomic sequence, may be aligned to different regions as split sequence due to exon/intron boundaries). Your suggestion will be highly appreciated:)

Thank you

sequencing alignment assembly • 408 views
ADD COMMENTlink written 10 months ago by dasakan0
1

contig,split,transcripts... I don't understand, please provide an example.

ADD REPLYlink written 10 months ago by Pierre Lindenbaum124k

Hi,

The meaning goes like this:

Contig: assembled sequences of RNA-seq reads, where each read is considered as transcript. As these sequences were generated from RNA sequences, so there will not be intronic sequences. Hence, if we align a contig sequence (mRNA/DNA) to genomic sequence, it may align to different parts of the genome to corresponding exons, which I mean different fragments of same contig sequence as 'split'. (N.B. I have on hand only the assembled contig sequence, without the RNA-seq reads).

ADD REPLYlink modified 10 months ago by finswimmer13k • written 10 months ago by dasakan0
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: 2093 users visited in the last hour