Comments (3)
In the SAM/BAM format spec, it's confusing to describe "TLEN" as a "signed length", because lengths are non-negative quantities. TLEN could be described more clearly as an "offset". It is the observed offset from one end of the template to the other. For the leftmost segment, the offset is positive. For the rightmost, it is negative. In other cases it is undefined or 0.
from hts-specs.
You are correct about the flags in the example, and I am in the process of fixing it.
IMHO it's barely reasonable to describe something as a signed length, to mean a length (which is of course unsigned) together with a sign attached to it denoting something related. However in any case, we've been reworking the description of TLEN to hopefully make it clearer -- please have a look at #23 and comment there on whether you think it's an improvement.
from hts-specs.
The same change to the flags needs to be made for the padded SAM version of the same example that is provided in section 3.2. That is, the flags 163 and 83 for the first and last reads of qname=r001 should be 99 and 147, respectively. Sorry for not mentioning this previously-- I hadn't finished working all the way through the spec.
from hts-specs.
Related Issues (20)
- VCF: "Genotype fields" vs "FORMAT" and per-sample HOT 1
- primary, secondary, and supplementary alignments with optional MM tags HOT 5
- Modified base single letter codes update HOT 7
- test/sam: Duplicate aux field tags
- test/vcf: Duplicate contig header record ID
- FAIRsharing Record Query - BED format
- CRAM: Need to improve feature positions description HOT 1
- is `*` better than `\*`? HOT 1
- cram: interpretation of "unmapped" flag in a pseudocode seems incorrect HOT 1
- SVCLAIM: VCF4.4 and backward compatibility with VCF4.3 HOT 1
- How to retrieve the primary alignment for secondary and supplementary reads HOT 8
- Is there a semantic difference between GT=./. and GT=0/0 + GQ=0 ? HOT 15
- cram: Inconsistent descriptions of auxiliary tag types HOT 1
- SA tag CIGAR format
- vcf: Handling structured header records with missing IDs in VCF 4.1/4.2 HOT 1
- bcf: First phasing indicators not set in genotype (GT) value examples
- CSI file is BGZF compressed but this is not mentioned in the CSV1 spec HOT 2
- Questions about third-party use of test data HOT 6
- VCF Draft 4.5 and Modified Bases HOT 27
- VCF4.4 SVLEN requirement across different variant representations HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from hts-specs.