SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA unique mapping metheuse Bioinformatics 3 04-13-2015 07:03 AM
How does BWA deal with multi mapping reads? metheuse Bioinformatics 5 03-17-2015 12:08 PM
Unique mapping with Tophat2 sazz Bioinformatics 7 07-01-2014 08:47 AM
How to do unique mapping in bowtie2 metheuse Epigenetics 9 11-06-2013 06:10 AM
behaviour of bwa mapping for homologous regions. KevinLam Bioinformatics 1 07-28-2010 09:39 AM

Reply
 
Thread Tools
Old 08-26-2015, 05:00 AM   #1
harlequin
Junior Member
 
Location: Europe

Join Date: Oct 2012
Posts: 7
Default Mapping: how to deal not unique mapping to homologous genes in different contigs

Hello,

I'm trying to map some RNA-Seq reads to a fungus genome using the STAR aligner, and am subsequently using HTSeq to estimate the counts per gene. I am using an ENSEMBL genome fasta file and the corresponding GTF file. My problem is, the assembly consists of 413 super contigs, so I am unable to only use half of the homologous chromosomes (i.e. having a haploid reference) to avoid having my reads matched to multiple locations. Any suggestions for this? What is the most popular way to deal with this problem, namely the lack of fully assembled chromosomes and therefore the resulting inability to filter the homologous ones to avoid non-unique alignments?

Thanks in advance.
harlequin is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:21 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO