Question: The Meaning Of "Seq Being Reverse Complemented" In Sam Format
3
gravatar for Pascal
7.0 years ago by
Pascal1.4k
Barcelona
Pascal1.4k wrote:

Hi

After reading the SAM specs it is still not very clear what is the meaning of these two fields in the FLAG:

  • SEQ being reverse complemented
  • SEQ of the next segment in the template being reversed

Could you please explain it to me?

UPDATE: if I understood it well when "SEQ being reverse complemented" is NOT SET it means that the read has been aligned to the forward strand, when SET it means that it has been aligned on the reverse strand. Is that correct?

And can we use these fields to detect an incorrect orientation of a paired-end read ?

Regards, Pascal

sam structural • 2.8k views
ADD COMMENTlink modified 4.6 years ago by Biostar ♦♦ 20 • written 7.0 years ago by Pascal1.4k
3
gravatar for Istvan Albert
7.0 years ago by
Istvan Albert ♦♦ 78k
University Park, USA
Istvan Albert ♦♦ 78k wrote:

As you mention it means that the reverse complement of the read aligns to the position that is indicated; and that is equivalent to saying that it aligns to the complementary strand in reverse direction.

The reason for the wording is probably to remind you that the position field is still the leftmost position of the alignment and not the actual 5' start of the read.

You may use this information to detect properly paired reads, although there is more to these than just being in the correct orientation, for example the distance between them and other considerations may need to be taken into account. Flag 2, read mapped in proper pair should be filled in by the aligner.

ADD COMMENTlink written 7.0 years ago by Istvan Albert ♦♦ 78k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1736 users visited in the last hour