Mapping DNAseq reads on cds ?
It may be a dumb question but, does anyone know a method / tool to map DNAseq reads to cds ? Reads falling in between exon and intro may have to be automatically trimmed. Thanks for help.
alignment
next-gen
• 1,879 views
•
link
written
by
thomasbersez
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How do I consider that my reads might map to the 3' UTR (3'-tag RNA seq)?
written by applepeanut •I used a 3'-tag library kit and sequenced RNA samples. I have read that because of the library preparation method, the reads will be biased …
-
Mapping reads to mitochondrial genome
written by GA •I need to map sequencing reads to the Drosophila mitochondrial genome. Does anyone know if I can accomplish this using the STAR aligner? Can the …
-
Filtering a GFF to retain only exons with CDS annotation
written by elcorteganoHi all, I'm trying to filter a large GFF file with many gene/exon entries to make it retain only those gene/exons entries that contain CDS …
-
Extract reads mapped within chosen intronic region from BAM file (no junctions)
written by gspirito •Hi everyone, here's my question: I have some BAM files resulting from mapping RNA-seq reads. I would like to know whether there is a way …
-
Racon polishing script that maps long reads and does racon polishing continuously until no changes …
written by zack.saud •Has anyone come across a script that does that maps long reads and does racon polishing continuously until no changes are seen? I have Nanopore …
-
MMR for variant calling ?
written by thomasbersezHi, Does anyone ever used MMR or others ambiguous reads resolution tools on bam files before to do variant calling? Does it improvement results accuracy? …
-
How to Remove Reads Found in Negative Controls from Experimental Samples
written by kdevakan •Hi I have paired-end reads from human DNA samples, where I am trying to determine the metagenomic viral profile for each sample. I also have …
-
detect SNPs using DNAseq fastq files which were used to assemble the genome
written by shl198Hello all, I find a problem that interest me. The chinese hamster genome was assembled using DNAseq. And when I tried to use the same …
-
Merging Data From Pfam Pdb And Uniprot
written by Aurobhima •<p>Hi,</p> <p>Does anyone know of a method to map between Pfam, PDB and UniProt. I have very specific criteria I want to select data with, …
-
Metabolic Network Drawing
written by Kwet •<p>Hi all,</p> <p>As part of a project I have a to draw a large metabolic map. I was wondering what was the most efficient way …
Reads in that region will fail to map if you are only aligning against CDS sequences in your reference. Aligners will generally soft clips reads (if the alignment ends partway in the read).
What exactly are you wanting to do?
I am developing a pipeline for variant calling on wild relative of cultivated plants. The overall idea is to find cds variants in wild specie. Because wild specie genome isn't always available, I may to map my DNAseq read onto cds of the gene(s) of interest.
Thanks for the answer!
It is doable with BWA or Bowtie2. The main issue is that you'll need to tweak the gap parameters to assign a very low penalty to gaps, since I guess that some of your CDS are entire mRNA without the introns.
In BWA you might be interested in the following parameters:
Source: http://bio-bwa.sourceforge.net/bwa.shtml
Starting from there, I think it is just trial and error to find the correct parameters