<p>There is nothing in the <a href="https://samtools.github.io/hts-specs/SAMtags.pdf" rel="nofollow noreferrer">SAMtags</a> specification that dictates the order of the entries in an SA tag.  It may be entirely dependent on the aligner that you used.</p>

<p>Looking at my own data that was aligned with <a href="https://github.com/lh3/bwa" rel="nofollow noreferrer">bwa</a>, I cannot see a definite pattern.  Longer matches tend to be first, but not always.</p>

<p>Your aligner documentation may have more information on how it uses SA tags.</p>

CC BY-SA 4.0
2020-04-17T16:12:18.713Z
2020-04-17T16:12:18.713Z
Andrew Whitwham
<p>I have several long-read sequences, which when aligned to a human reference genome, aligns to multiple chromosomal fragment as a result of chromothripsis shattering the linear DNA and randomly piecing them together. This creates for a single primary alignment multiple supplementary alignments.</p>

<p>With a single example of a read (only showing parts of all the bamfile information) being</p>

<p><code>cdaf1fcb-27bf-4890-8ac7-a097dfa12c38    16      chr11   15151487        60</code></p>

<p>With the primary alignment starting at chr11 at 15151487</p>

<p>Then the supplementary alignment can be seen as follows 
<code>SA:Z:chr13,49827969,+,1839S956M25D1584S,60,53;chr11,15152459,-,833S749M22D2797S,60,72;chrX,138324656,-,2540S404M2D1435S,60,20;chr11,15152066,+,278M11D4101S,60,20;</code></p>

<p>Where the alignments are to chr13, 11, X and 11. With both regions on chr11 being close to the initial region from the primary alignment, making the regions on chromosome 11 the same (I know I haven't calculated the full chromosomal interval using the Cigar string for this example).</p>

<p>I assume the order of that supplementary information above is in the actual order. But the problem is that my DNA should be circular which doesn't fit with the order of the alignments since the circle would then be chr 11 (primary, start of circle) -> chr 13 -> chr11 (same as primary) -> chrX -> chr11 (back to the start of circle) </p>

<p>My question is:
Q) So is the order of the supplementary information just pieces together randomly or is it true that I can assume the order is the true order of the sequencing read? </p>

<p>So thanks a lot for your help.</p>

<phylogenetics><sam><long-reads>
Chromosomal order of supplementary alignment in BAM file "SA" tag
# Chromosomal order of supplementary alignment in BAM file "SA" tag

(Asked by: username_0 on 2020-04-17)

I have several long-read sequences, which when aligned to a human reference genome, aligns to multiple chromosomal fragment as a result of chromothripsis shattering the linear DNA and randomly piecing them together. This creates for a single primary alignment multiple supplementary alignments.

With a single example of a read (only showing parts of all the bamfile information) being

`cdaf1fcb-27bf-4890-8ac7-a097dfa12c38 16 chr11 15151487 60`

With the primary alignment starting at chr11 at 15151487

Then the supplementary alignment can be seen as follows
`SA:Z:chr13,49827969,+,1839S956M25D1584S,60,53;chr11,15152459,-,833S749M22D2797S,60,72;chrX,138324656,-,2540S404M2D1435S,60,20;chr11,15152066,+,278M11D4101S,60,20;`

Where the alignments are to chr13, 11, X and 11. With both regions on chr11 being close to the initial region from the primary alignment, making the regions on chromosome 11 the same (I know I haven't calculated the full chromosomal interval using the Cigar string for this example).

I assume the order of that supplementary information above is in the actual order. But the problem is that my DNA should be circular which doesn't fit with the order of the alignments since the circle would then be chr 11 (primary, start of circle) -> chr 13 -> chr11 (same as primary) -> chrX -> chr11 (back to the start of circle)

My question is:
Q) So is the order of the supplementary information just pieces together randomly or is it true that I can assume the order is the true order of the sequencing read?

So thanks a lot for your help.

**Tags:** phylogenetics, sam, long-reads

---

## 1 answers

### Answer by username_1 on 2020-04-17. Score: 0 | Accepted answer

There is nothing in the [SAMtags](https://samtools.github.io/hts-specs/SAMtags.pdf) specification that dictates the order of the entries in an SA tag. It may be entirely dependent on the aligner that you used.

Looking at my own data that was aligned with [bwa](https://github.com/lh3/bwa), I cannot see a definite pattern. Longer matches tend to be first, but not always.

Your aligner documentation may have more information on how it uses SA tags.


---

