Skip to content

Commit

Permalink
Document SAM ambiguity of 1bp long seq with QUAL "*".
Browse files Browse the repository at this point in the history
"*" is either QUAL 9, or QUAL unavailable.  Made a recommendation in
a footnote, mainly as an indication that the ambiguity exists.  In
practice it's vanishingly unlikely to matter.

Fixes samtools#715
  • Loading branch information
jkbonfield committed Jan 29, 2024
1 parent 2543789 commit a8550ef
Showing 1 changed file with 4 additions and 0 deletions.
4 changes: 4 additions & 0 deletions SAMv1.tex
Original file line number Diff line number Diff line change
Expand Up @@ -651,6 +651,10 @@ \subsection{The alignment section: mandatory fields}\label{sec:alnrecord}
wrong}\}$. This field can be a `*' when quality is not stored. If
not a `*', {\sf SEQ} must not be a `*' and the length of the quality string
ought to equal the length of {\sf SEQ}.
\footnote{There is a small ambiguity with a sequence exactly one
base-pair long having quality 9. This is ASCII `*' so it could be
interpreted as either QUAL 9 or as quality unavailable. Given
this ambiguity, we recommend interpreting it as unavailable.}
\end{enumerate}
\subsection{The alignment section: optional fields}\label{sec:alnaux}
Expand Down

0 comments on commit a8550ef

Please sign in to comment.