Thanks, Nils; I had already tried changing the prefix, but it still seems to fail -- I'll look into it again when I have time.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Hi nilshomer and thanks for your reply!
I've sorted and created the indexes but the samtools message remains the same.
Any other ideas? Do you know if there exists another way to extract a particular chromosomal region from a .bam file?
Thanks,
M.
Comment
-
Originally posted by Jeckow View PostHi nilshomer and thanks for your reply!
I've sorted and created the indexes but the samtools message remains the same.
Any other ideas? Do you know if there exists another way to extract a particular chromosomal region from a .bam file?
Thanks,
M.
Nils
Comment
-
Originally posted by tgolubch View PostThanks, Nils; I had already tried changing the prefix, but it still seems to fail -- I'll look into it again when I have time.
Nils
Comment
-
That sounds suspiciously similar to a wig file, although, the format is not quite what you're expecting. I assume what you're looking for would look like:
chr1 1 0
chr2 2 1
chr2 3 1
... etc
The problem is that such file would be incredibly inefficient, hence the wig file, which will provide a more compact version. (Presumably the "variable width" format is what you're looking for.)
Although not quite the same, FindPeaks 4.0 (Yes, sorry, I'm plugging my own software) will generate fixedstep format wig files from sam/bam files, which would give you a similar result, which would look like:
fixedStep chrom=chr22 start=16653550 step=1
6
7
8
9
10
11
......
where you get a single line that tells you the start point (chr and position), and then the coverage of each position after that.
There are plenty of other wig file generators out there as well - although I'm not sure how many of them read bam/sam yet.
CheersThe more you know, the more you know you don't know. —Aristotle
Comment
-
@Alberto Magi
What you are looking for is "pileup" (http://samtools.sourceforge.net/pileup.shtml). Please also see this page (http://samtools.sourceforge.net/pipe.shtml) about getting pileup in a specified region.
@apfejes
Pileup output is not part of Picard, you should be able to get similar things from GATK (http://www.broadinstitute.org/gsa/wi...alysis_Toolkit)
Comment
-
Thanks - pileup also fits that description. Variety is the spice of life. (=
Though, just to clarify, I was just looking to point Alberto in the right direction, and had forgotten about pileup. (I don't use pileup, but I do use Picard.)The more you know, the more you know you don't know. —Aristotle
Comment
-
I was curious if anyone has found a workaround to this error? It is coming up for me when viewing novoalign alignments. I am using samtools 0.1.6 (r453), and get the following error when loading a samtools converted bam file from novoalign's sam file
samtools: bam_lpileup.c:116: tview_func: Assertion `l == tv->n_pre' failed.
thanks
Originally posted by zee View PostI was wondering if I could get some help on the SAM specification.
I get the error:
When I page the viewer to my reads containing indels. A sample of my data is
And I'm using novoalign format where I've modified the novo2sam.pl script to produce a 'correct' CIGAR format. Most samtools operation work except for the tview
Comment
-
Pileup Base Quality
Hi, I was wondering if anyone could point me in the direction of any documentation that describes how to interpret the base quality information reported within a pileup format file.
I've checked here: http://samtools.sourceforge.net/pileup.shtml but can't find any information.
Cheers
Comment
Latest Articles
Collapse
-
by seqadmin
The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...-
Channel: Articles
05-06-2024, 07:48 AM -
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:35 AM
|
0 responses
14 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:35 AM
|
||
Started by seqadmin, 05-09-2024, 02:46 PM
|
0 responses
18 views
0 likes
|
Last Post
by seqadmin
05-09-2024, 02:46 PM
|
||
Started by seqadmin, 05-07-2024, 06:57 AM
|
0 responses
17 views
0 likes
|
Last Post
by seqadmin
05-07-2024, 06:57 AM
|
||
Started by seqadmin, 05-06-2024, 07:17 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
05-06-2024, 07:17 AM
|
Comment