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
Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.
Somatic Genomics
“We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...-
Channel: Articles
05-24-2024, 01:16 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 06-03-2024, 06:55 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
06-03-2024, 06:55 AM
|
||
Started by seqadmin, 05-30-2024, 03:16 PM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
05-30-2024, 03:16 PM
|
||
Comprehensive Sequencing of Great Ape Sex Chromosomes Yields Insights into Evolution and Genetic Variability
by seqadmin
Started by seqadmin, 05-29-2024, 01:32 PM
|
0 responses
29 views
0 likes
|
Last Post
by seqadmin
05-29-2024, 01:32 PM
|
||
Started by seqadmin, 05-24-2024, 07:15 AM
|
0 responses
216 views
0 likes
|
Last Post
by seqadmin
05-24-2024, 07:15 AM
|
Comment