Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Mapping: how to deal not unique mapping to homologous genes in different contigs harlequin RNA Sequencing 0 08-26-2015 05:00 AM
mapping contigs to reference genome sadiexiaoyu Bioinformatics 2 07-04-2014 06:19 AM
Predicted Ensembl genes, but not in RefSeq sindrle RNA Sequencing 7 04-25-2014 09:24 AM
RNA-seq Mapping to Many Contigs Dario1984 Bioinformatics 10 06-26-2013 06:00 PM
Too few reads mapping back to contigs cerebralrust Bioinformatics 8 02-26-2012 11:25 PM

Thread Tools
Old 09-16-2019, 07:58 AM   #1
Junior Member
Location: Europe

Join Date: Sep 2019
Posts: 1
Default Re-mapping contigs or predicted genes?

I'm quite new to bioinformatics and am currently running into a dilemma on how to approach my re-mapping.
I have successfully assembled my FASTQ-files (reads) and now have my contigs in FASTA files.
From these FASTA-files, I want to predict its genes (with Prodigal) and find the depth (abundance) of the predicted genes. However, I am unsure if I should re-map my contigs back to my reads first, or if I should gene predict first, and then re-map the predicted genes back to the reads. As I see it, by re-mapping the predicted genes, I risk losing some of the reads not entirely in the gene. However, by re-mapping the contigs, I get a depth for all the reads on the contig, but will then have to correlate that data to the predicted genes somehow. Do you have any suggestions or preferences? And if you do, why would you choose one over the other.

Thank you very much in advance!
Drengo is offline   Reply With Quote
Old 09-16-2019, 10:37 AM   #2
Senior Member
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,082

Cross-posted on biostars:
GenoMax is offline   Reply With Quote

contigs, fasta, fastq, mapping, prodigal

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 12:12 AM.

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