No Sequencing Data at Low Positions
Opening a bam file for just chr22, I was surprised to see that there were no reads aligned until around position 16,050,000. In UCSC Genome Browser, looking at a window of chr22:16,049,420-16,050,420, I can see that there's "nothing," but then different tracks start. In chr1, I also think I remember seeing that alignments started only at 10,000.
Why is there nothing earlier in the chromosome? Do those positions not correspond to DNA? I've downloaded data on the length/end position of each chromosome; can I find a list of these "start" positions?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
mRNA positions from vcf
written by jindal.aashi21 •
-
FANTOM5 enhancers: intepretations
written by c.clarido
-
Assign a genomic region to a window
written by Ana
-
Problems uploading custom track to UCSC Genome Browser
written by rs •
-
Mutation position in COSMIC
written by ognjen011
-
how to find SNPs within a range of given windows?
written by Ana
-
BED alignments not positioned properly in IGV
written by bn •
-
Sliding window using gene position file
written by burnsro •
-
biomart gene coordinates do not correspond to genome browser
written by tonja.r
-
turning autoplot objects into a grob that can be arranged with grid.arrange
written by adeslatt •
Are you looking at the right genome build in UCSC? Also are you sure the data is aligned against UCSC genome build (which have a
chrprefix for chromosomes as opposed to other builds which may only have numbers.Beginnings of the chromosome sequence may only have N's since the ends of chromosomes are hard to sequence.
I think so. http://genome.ucsc.edu/cgi-bin/hgTracks?db=hg19&lastVirtModeType=default&lastVirtModeExtraState=&virtModeType=default&virtMode=0&nonVirtPosition=&position=chr22%3A16049522-16050522&hgsid=662132401_e7mUocboTcIMiTPiqCEGpqCVAIpq
22 is only 51304566 bp long. So the first ~16M bases is almost a third of it
Ah it does indeed like everything before that is "N" when I zoom into "base" on the browser.
So those N-nucleotides means that "we know there are nucleotides there, we are just not sure what they are"
Got it. Thank you! So is there somewhere I can find out how many of the first bases of each chromosome are N?
Hello,
you could use your language of choice to find the first position in each reference sequence which is not an N.
fin swimmer