This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Flawed Sam Header Regex?

According to the SAM format white paper (http://samtools.sourceforge.net/SAM1.pdf) header lines should be matched by /^@[A-Za-z][A-Za-z](\t[A-Za-z][A-Za-z0-9]:[ -~])+$/ or /^@CO\t.*/. However, this does not seem to be the case with the example in the format white paper - or with real world data.

ruby -e 'puts "@HD\tVN:1.3\tSO:coordinate" =~ /^@[A-Za-z][A-Za-z](\t[A-Za-z][A-Za-z0-9]:[ -~])+$/'

So is the SAM header or the regex flawed?

Cheers

Martin

sam format
• 2,824 views • link
updated by brentp brentp • written by Martin A Hansen Martin A Hansen

2 answers

I use rubular for testing regexes. You could add your test SAM string and then play with the regex to get the correct match?

• 1 views • link
written by Michael Barton • Michael Barton

Hm, looks like the regex is flawed. Whereas, it's currently:

/^@[A-Za-z][A-Za-z](\t[A-Za-z][A-Za-z0-9]:[ -~])+$/

It seems it should be:

/^@[A-Za-z][A-Za-z](\t[A-Za-z][A-Za-z0-9]:[ -~]+)+$/

Where the extra '+' allows more than 1 character following the :

• 1 views • link
written by brentp brentp

Funny, to my understanding [ -~]+ allows one or more of that group of chars space, dash and tilde. How it matches '1.3' and 'coordinate' baffles me.

• 0 views • link
written by Martin A Hansen Martin A Hansen

@masha it does look odd, but it's the same format as [A-Z] it means match any characters between (including endpoints) " " and "~". ord(" ") == 32 and ord("~") == 126

• 0 views • link
written by brentp brentp

Of cause. I submitted the question to the samtools mailing list. Waiting for an answer.

• 0 views • link
written by Martin A Hansen Martin A Hansen

Log in to answer this question.

More posts like this