[mira_talk] Re: PROJ_default.(un)padded.fasta vs PROJ_ReferenceStrain vs PROJ_AllStrains
- From: Bastien Chevreux <bach@xxxxxxxxxxxx>
- To: mira_talk@xxxxxxxxxxxxx
- Date: Mon, 17 Oct 2011 21:24:05 +0200
On Oct 17, 2011, at 21:13 , Bastien Chevreux wrote:
> So, assume you mapped a strain against a backbone, and you get in the
> multiple alignment:
>
> ReferenceStrain gtacg*tcagtcg
> default/StrainX read1 gtacgAtTag*cg
> default/StrainX read2 gtacgAtTag*cg
> default/StrainX read3 gtacgAtTag*cg
> etc.pp
>
> (note the differences)
>
> Now, if you extract a consensus, convert_project will create one consensus
> sequence for each strain and one for all strains together. In the case above
> it will be three sequences:
>
> ReferenceStrains: gtacgtcagtcg
> default/StrainX: gtacgAtTagcg
> AllStrains: gtacgWYWRKcg
Gah, wrong "AllStrains". Here's the correct one:
AllStrains: gtacgatyagtcg
B.
Other related posts:
