Re: [Samtools-help] Zero Length Reads in SAM

2017-01-16 Thread Colin Hercus
Yes, I sometimes overlook the regular expression definitions. On 16 January 2017 at 19:12, James Bonfield wrote: > 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 work

Re: [Samtools-help] Zero Length Reads in SAM

2017-01-16 Thread James Bonfield
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 seque

Re: [Samtools-help] Zero Length Reads in SAM

2017-01-16 Thread Colin Hercus
Hi Rob, I tried that and it worked. Thanks. 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. Cheers, Colin On 16 January 2017 at 17:54, Robert Davies wrote: > On Mon, 16 Jan 2017, Colin Hercus wrote: > > Sometim

Re: [Samtools-help] Zero Length Reads in SAM

2017-01-16 Thread Robert Davies
On Mon, 16 Jan 2017, Colin Hercus wrote: > Sometimes we get reads to align that have been trimmed to zero length and > I'm wondering how these should be represented in SAM format. > > Here's a pair as reported by Novoalign that had been trimmed by cutadapt > and one read of the pair is zero length

[Samtools-help] Zero Length Reads in SAM

2017-01-15 Thread Colin Hercus
Hi, Sometimes we get reads to align that have been trimmed to zero length and I'm wondering how these should be represented in SAM format. Here's a pair as reported by Novoalign that had been trimmed by cutadapt and one read of the pair is zero length READID77*00**0