Question: Chromosome numbering in sam files
0
gravatar for haleykgrant
7 months ago by
haleykgrant0 wrote:

I have a (probably stupid) question about how the chromosome positions are numbered in SAM files.

I used bowtie to align reads to hg19 and got SAM files as output. I pulled the columns indicating chromosome and position so that I have files that look like:

chr, position

chr1, {starting position for read 1}

chr1, {starting position for read 2}

chr1, {starting position for read 3}

... etc.

My question is: do the position indices reset at each chromosome (i.e position 1 on chromosome 2 would read "chr2, 1") or does the position indicate the position within the entire genome (i.e. if chromosome 1 has X positions, position 1 on chromosome 2 would read "chr2, X+1")?

alignment genome • 156 views
ADD COMMENTlink modified 7 months ago by ATpoint42k • written 7 months ago by haleykgrant0
0
gravatar for ATpoint
7 months ago by
ATpoint42k
Germany
ATpoint42k wrote:

It is a valid question, don't worry. Every chromosome (or contig or whatever you use to align against) is considered a separate unit so it starts at 1 and ends at the length of the chromosome. Start of 1 is for 1-based formats, for 0-based formats the start would be 0, see Cheat Sheet For One-Based Vs Zero-Based Coordinate Systems

ADD COMMENTlink modified 7 months ago • written 7 months ago by ATpoint42k
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: 1195 users visited in the last hour