[mira_talk] Re: New Solexa seqence identifier format - Paired end reads

  • From: Bastien Chevreux <bach@xxxxxxxxxxxx>
  • To: mira_talk@xxxxxxxxxxxxx
  • Date: Mon, 16 Apr 2012 00:28:15 +0200

On Apr 14, 2012, at 23:33 , Alexis Blanchet-Cohen wrote:
> Will MIRA understand the new format?

The 3.4 series does not. If you are feeling adventurous, give the 3.9.0 
development version a try: this one does.

> Do I have to convert the new format to the old format for MIRA to distinguish 
> the paired end reads?

Yes.

> If so, what is the sed command to do this?

I wrote earlier "this is left as an exercise to the reader." Oh well :-)

Try
  sed 's/^@\(.+\) \(.\)/@\1\/\2 \2/' infile >outfile

Should work.

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: