Handling of Reverse Strand in GAM file
In the gam file, for a reverse strand how is the CIGAR string represented? Is the read sequence already in reverse compliment form in gam file and the CIGAR string follows this read seq and the Flag is raised for reverse strand?
• 1,579 views
•
link
updated
by
glenn.hickey
•
written
by
ved_vyas •
•
written
by
ved_vyas •
1 answer
There is some documentation on the GAM format here: https://github.com/vgteam/libvgio/blob/master/deps/vg.proto
Note that in vg we are phasing out GAM in favour of GAF: https://github.com/lh3/gfatools/blob/master/doc/rGFA.md#the-graph-alignment-format-gaf with "cs" cigar tags as defined here: https://lh3.github.io/minimap2/minimap2.html
Most (but not yet all) vg tools have options to input or output GAF instead of GAM
• 0 views
•
link
written
by
glenn.hickey
Log in to answer this question.
More posts like this
-
Clarification regarding SAM flags "mate reverse strand" (flag 16/0x10) and "read reverse strand" (f…
written by kalavattam •
-
gene strand determined by flag in read1 and read2 should be intepreted differently?
written by tomas4482 •
-
Meaning of the column of sam file
written by Jimpix •
-
Force BlastN to align only to minus strand of subject
written by mwittep •
-
Flags in bam file and gene strand orientation
written by skylinesky •
-
Extracting CIGAR string from a specific region of a read in a bam file using pysam?
written by kiki4133 •
-
How to use scanbamparam and CIGAR from Rsamtools to count specific reads?
written by MAPK
-
Translating DNA to protein with perl
written by ksi216 •
-
Bowtie2 local alignment SAM file SEQ output - reverse complement or not ?
written by Roger •
-
Reverse Read Direction In A Bam File
written by lin.barnum •
