Refseq Number/Name In Bam File
I am using bamread in matlab to read a bam file that I downloaded from 1000G. In the website it says all the reads are mapped to reference sequence called "GRCh37" which make sense. But when I try to read it in matlab, it's syntax is
bamread(File,RefSeq,Range)
The problem is this RefSeq in bamread, doesn't recognize "GRCh37" whatsoever! then I tried to pass refseq code such as 1 or 2 or ... it works but I think I'm missing sth. How do I make sure that the reads are exactly mapped to reference GRCh37 while I don't see it in the bam file?
Thanks in advance!
• 3,269 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
Sami Salami •
•
written
by
Sami Salami •
1 answer
print the header of your bam with samtools:
samtools view -H file.bam | grep '^@SQ'
the chromosomes contained in your reference will be listed there.
• 0 views
•
link
updated
by
Ram
•
written
by
Pierre Lindenbaum
•
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
How To Sort Sam File
written by Chen Sun
-
The Reference Genome Of The Reading In Bam File
written by shl198
-
Is That Necessary To Remove The Chr M And Y Of Reference File And Bam File For The Downstream Work?
written by Tonyzeng •
-
Bam To Afg
written by bioLife •
-
Cufflinks Paired End Read Names
written by Frenkiboy
-
Length Of Contigs In Bam File Off By 1 Versus Reference Fasta
written by Tancata •
-
Merging Bam Files
written by rickyflintoff •
-
Incompatible Contigs Issue In Mitochondria Dna
written by ckr123 •
-
Samtools Idxstats
written by KCC
-
Samtools + Picard Markduplicates
written by Nery •
The point of confusion seems to be that GRCh37 is not a single sequence but multiple chromosomes which are usually named in bam files as 1, 2, ..., 22, X, Y