Yes, you are right. Its a mixture of plant- microbe RNA. I want to map the reads to both the genomes and separate them out. My concern is that the output with the ones mapped to plant might as well by mistake contain some from the bacteria and vice versa. Anyway to compare the output files generated from the two alignments?
How to check if the mapping results from tophat2 are correct?
Hi, I have RNASeq data (paired-end) from a plant-microbe interaction study at different conditions. I mapped the reads to the plant genome using tophat2 and got 32.3% to 89.2% overall read alignment rate (30.6% to 86.4% concordant pair alignment rate). 1. How to ensure that the mapping results are good? 2. How to check whether in fact the mapped reads are correct, I mean how to check that all mapped reads are of the plant and some reads from the bacteria are not contaminating the result in the output file?
2 answers
You can always view the alignments with a genome viewer (IGV is easy to use) and confirm that the alignments look reasonable. Using a different aligner (BBMap/STAR are good alternates) will provide an independent confirmation.
As for the origin of the reads (plant/microbe) you can check a few using Blast at NCBI. Other than some chloroplast/mitochondrial DNA (which may have some similarities to microbes) other data should map specifically to plants.
The only way to check for contamination is to map the reads against likely contaminants and see how many align at least as well there compared to how they align to your plant. I'm guessing that you really want to align your reads against both genomes and then assign them according to how they align. I expect the bbsplit tool from BBMap would work for something like this.
• 1 views
•
link
written
by
Devon Ryan
BBSplit syntax for generating builds for the reference genome and how to call different builds.
Use appropriate references as input. Bacterial sequences should not be aligning to plants (with some exceptions noted above).
If you are sure that there should be only one plant (and a limited number of bacteria, do you know how many or is it general contamination from plant surfaces etc) then using the BBsplit method is the fastest way to get this done. lakhs of reads are no problem for BBMap, which can process a couple of thousand reads per second.
Yes, I am sure as it is a planned experiment- one plant and one bacteria. Is the link you have sent and this one (https://www.biostarhandbook.com/tools/bbmap/bbmap-help.html ), are sufficient or you can direct me to some tutorial, as I am new to linux and transcriptome analysis?
Log in to answer this question.
More posts like this
-
Adaptor trimming issue
written by 1769mkc
-
Optimization of RNA-Seq data mapping with tophat2
written by Denis
-
Problem converting bam to sam (samtools)
written by hinkel2 •
-
PolyA reads in RNA-seq and how do aligners handle them?
written by Ruolin Liu •
-
The comparison between HISAT2 and Tophat2
written by Peng Huang
-
TopHat alignment on Trinity Transcripts giving a lot of unmapped results.
written by kanika.151
-
Chip-Seq's mapping rate problem.(too low)
written by leenaehyeon •
-
What does it mean that: 100% of my alignment are discordant?
written by tiago211287
-
tophat output files containing the reads which mapped uniquely as a pair
written by trakhtenberg
-
Tophat Both Pairs Must Map Concordantly
written by AW
You could use a metagenomics tool, like Kraken, to quickly estimate what contaminants are present in your sample. It won't tell you about issues mapping reads, but you'll at least know the source of your problem.
thanks joe. I'll try it. Kraken will tell the source of contamination, if it finds such reads? This means I'll have to give my reference. Can you please direct me to some links for Kraken manual and tutorials?