Question: Extracting alignment coordinates from SAM file
gravatar for dasakan
21 months ago by
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 • 711 views
ADD COMMENTlink written 21 months ago by dasakan0

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

ADD REPLYlink written 21 months ago by Pierre Lindenbaum131k


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 21 months ago by finswimmer14k • written 21 months ago by dasakan0
Please log in to add an answer.


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