Hi Everyone
I am using rna seq data with tophat2 and getting unmapped.bam file in addition to usual accepted.hits.bam file.
The unmapped.bam file is quite big in size as compared with alignment file. Is there anything i should worry about.
Regards V
•
written
by
Varun Gupta
It is larger because you have more unmapped reads than mapped reads. Whether or not that is something to worry about depends on many factors - the origins of samples, the genome you align against etc.
•
written
by
Istvan Albert
I am using the default parameters only...
•
written
by
Varun Gupta
I DID
Some of them have N's . So one reason is that for sure
But the ones which does not have N's how to deal with them..
Regards
•
written
by
Varun Gupta
Well, did you examine the ones without N's? Are they poor quality? Are they genomic contamination? Another species contamination?
There might be a very big problem. Maybe the sample prep people are screwing up. Maybe a kit was defective. Maybe your reference file is screwed up. Maybe you were totally misinformed as to what the project was about. You won't know unless you actually look at your raw data.
Perhaps you need to modify the arguments passed to tophat as the contraints (like the number of mismatches allowed, maximum size of intron and so on) to increase the number of accepted hits.
The other possibility is try to provide a better reference genome, if it's not a good one, or perhaps to see if your RNA-seq experiment is well done.
Did you check your "junctions.bed" file( It should be empty in your case)? I think thats because tophat is not able to find the splice junctions.
I was getting this problem and i think I changed the -r (innermatedist) parameter in tophat to fix this
Hi
This was single end data. So -r was not required. My junctions.bed file was not empty either..
Regards
varun
Log in to answer this question.
Hi Gupta
Did you find answers about this problem? now I met same situation, my
accepted-hits.bamandummapped.bamare about 1M and 1.2G respectively. If you have any ideas to solve it, please help me.Thank you!
Not really Chen. I moved to STAR aligner