On Mon, Jan 16, 2017 at 06:56:39PM +0800, Colin Hercus wrote:
> I always thought "*" was for sequence not stored rather than that the
> sequence is zero length. But as it works I'll stick with it.
You are correct that "*" is unspecified rather than zero size, but the
specification states the sequence has to be "*" or 1+ bases so the
empty string is not legal.
It's one of the odd corner cases, but perhaps no more odd than wanting
to store a zero length sequence! I assume you want it there simply as
evidence that it existed? I wonder what other tools will break though.
Another oddity is the one base long sequence with quality value 9. :)
James
--
James Bonfield (j...@sanger.ac.uk) | Hora aderat briligi. Nunc et Slythia Tova
| Plurima gyrabant gymbolitare vabo;
A Staden Package developer: | Et Borogovorum mimzebant undique formae,
https://sf.net/projects/staden/ | Momiferique omnes exgrabure Rathi.
--
The Wellcome Trust Sanger Institute is operated by Genome Research
Limited, a charity registered in England with number 1021457 and a
company registered in England with number 2742969, whose registered
office is 215 Euston Road, London, NW1 2BE.
------------------------------------------------------------------------------
Developer Access Program for Intel Xeon Phi Processors
Access to Intel Xeon Phi processor-based developer platforms.
With one year of Intel Parallel Studio XE.
Training and support from Colfax.
Order your platform today. http://sdm.link/xeonphi
_______________________________________________
Samtools-help mailing list
Samtools-help@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/samtools-help
