Originally posted by bioinfo816
View Post
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
-
Originally posted by bioinfo816 View PostI need suggestion to extract the chromosome or chromosome region (e.g. chr1: 100,000,000-200,000,000) without sorting bam file. Any idea related with this would be appreciated.
Having said that, a small python/perl/whatever script could do this (you could even put together a shell script, if you're so inclined). Just filter the output of samtools view through it.Last edited by dpryan; 09-24-2013, 05:27 AM.
Leave a comment:
-
chromosome extraction without sorting
I need suggestion to extract the chromosome or chromosome region (e.g. chr1: 100,000,000-200,000,000) without sorting bam file. Any idea related with this would be appreciated.
Leave a comment:
-
In a BAM file, the actual chromosome name is only stored in the header (each read has an associated index number that is used by samtools or other programs to get the actual name of the chromosome/contig).
If you have a SAM file then this won't work, you'll have to write a little script in python/perl/whatever.
Leave a comment:
-
Thanks,
I thought that samtools reheader only changes the actual header info, but what about in all of the individual alignments? I do know that the same genome build was used for everything so that shouldn't be too much of a problem.
Leave a comment:
-
You can use samtools reheader. Just samtools view -H file.bam > header.sam, then edit that to include the "chr" (don't change the ordering, though!). You can then use that for the reheader command. This assumes, of course, that this is the only difference between the genomes.
Leave a comment:
-
the "chr" in the chromosome name
Hi,
I have been given a few exome datasets (already mapped and I don't have the fastq files). The problem is that when these were mapped, they used a reference annotation that did not have the "chr" in with the chromosome ID. All of my other exomes that I have do have this annotation. As you can image, this is being a little problamatic when I am trying to do the analysis with all of the samples at the same time. Is there a simple way to add the "chr" to these samples without having to do all of the remapping?
Thanks
Latest Articles
Collapse
-
by seqadmin
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
160 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|
Leave a comment: