Extract Simple Alignment From Sam/Bam
Hi
I have annotated a number of RNA reads with the originating gene, and keep them in a bam file.
Now I want to extract a simple alignment for each gene from this sam file. I don't care about genomic location, quality of read or any other data. Only the aligned sequences corresponding to that gene. One alignment file (something like .aln) for each gene.
Any ideas on how to do that?
1 answer
I would actually suggest collecting the the reads and re-aligning against that gene using a Smith-Waterman implementation.
The Bio* packages have mature tools to display pairwise alignments - CIGAR is too goofy a format for this task.
• 1 views
•
link
written
by
Jeremy Leipzig
Log in to answer this question.
More posts like this
-
Matching read IDs to mapped transcripts/genes
written by joshcylee •
-
no alignment from plasma derived small RNA
written by memento1984_1 •
-
Does TCGA rna-seq data include fusion gene expression?
written by 9521ljh •
-
Calculate Mapping Rate for an alignment to each sequence of a multi-fasta reference
written by kspata
-
How do I extract aligned reads from BAM file?
written by MAPK
-
Nucletotide distribution, at each position, in a .sam/.bam file ?
written by remi.maglione •
-
SAM File target extraction
written by geneart$$ •
-
Extract Raw Counts From Sam File
written by mad.cichlids
-
Bam Viewers That Show Alignment Location On Query As Well As Reference?
written by jnfass •
-
Reconstructing Alignments From Cigar Strings
written by Revh •
what about using an alignment viewer like IGV ?