Hi Bastien,
I would like to know how I can distinguish real singlet and identical pseudo
singlets?
Best regards, Lyu
----- Original Message -----
From: Bastien Chevreux <bach@xxxxxxxxxxxx>
To: mira_talk@xxxxxxxxxxxxx
Date: 2016/9/4, Sun 11:51
Subject: [mira_talk] Re: Same singlets,
On 03 Sep 2016, at 23:13 , wakamoto5959@xxxxxxxxxxx wrote:
I have thought about mapping before. However, most (or all?) of mapping
software use genome sequence as reference, not RNA-seq contigs.
Ummm, does it make a difference whether you feed genome sequences to a mapper
or RNASeq contigs?
If in doubt … try MIRA. I use it regularly for that.
I can count singlets from debris files, but again, they are "identical"
UNSAVED_SINGLETS.
This is the difficulty I have.
Yes, and if you save the singlets in the regular output, you will have
thousands of single read contigs. Where’s the difference from a counting point
of view?
B.
--
You have received this mail because you are subscribed to the mira_talk
mailing list. For information on how to subscribe or unsubscribe, please visit
http://www.chevreux.org/mira_mailinglists.html