[mira_talk] Re: virus genome assembly

  • From: Bastien Chevreux <bach@xxxxxxxxxxxx>
  • To: mira_talk@xxxxxxxxxxxxx
  • Date: Thu, 22 May 2014 21:40:18 +0200

On 22 May 2014, at 21:28 , Laurent MANCHON <lmanchon@xxxxxxxxxxxxxx> wrote:
>> That being said, maybe your problem can be reduced a bit. You wrote detected 
>> … does that mean “counting” of known variants or “discovery” of new variants?
> 
> just counting
> 
>> 
>> For counting of known variants, what I’d do would be a simple test by simply 
>> putting all splicing variants
> 
> this is here my question, how to select or isolate this bag of splicing 
> variants from mira results ?

Ummm … if these are known variants you want to count, you have the sequence of 
those right? Then make a fasta file with these, map your reads against these 
(strict settings etc.), then look up in the statistics file of MIRA how many 
reads mapped. Or am I missing something here?

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

Other related posts: