Fail to determine the sequence name

I am new to BAM file and samtools. I want to extract part of bam file by

samtools view _SRA_RUN_AUTO_ARCHIVE_TEST_vault_EGAR00001003032_5499_1.bam "Chr5:18000-45500" > output.bam

But I got

[bam_parse_region] fail to determine the sequence name.
[main_samview] region "Chr5:18000-45500" specifies an unknown reference name. Continue anyway.

What happened?

bam samtools

• 3,772 views • link

updated 4.1 years ago by Ram

• written 10.6 years ago by zwang10 •

Hello zwang10!

It appears that your post has been cross-posted to another site: http://seqanswers.com/forums/showpost.php?p=187067&postcount=1

This is typically not recommended as it runs the risk of annoying people in both communities.

• 0 views • link

written 10.6 years ago by Pierre Lindenbaum

1 answer

chech there is a really a Chr5 (case sensitive) in the bam by printing the header:

samtools view -H _SRA_RUN_AUTO_ARCHIVE_TEST_vault_EGAR00001003032_5499_1.bam

• 0 views • link

updated 6.6 years ago by Ram

• written 10.6 years ago by Pierre Lindenbaum

More posts like this

Error message when using samtools view
written 3.0 years ago by Shin Taguchi •

Hi, every one. I ran the following command. samtools view input.bam scaffold_list > subset_reads.sam I then received the following error message. [main_samview] region "scaffold_list" specifies …
I am trying to select some reference genome region of a bam file, but got an error
written 6.0 years ago by schlogl

Hi there, I am trying to use this command: samtools view -b sorted_SRR1972739.bwa.bam AF086833:3050-3199 > selected.bam And got this error: [main_samview] region "AF086833:3050" specifies an …
Failure in extracting reads from bam file due to unknown reference
written 6.8 years ago by marongiu.luigi

Hello, I need to extract reads from an alignment I have done of the whole human genome. I am looking for the gene [MACC1][1] which …
extracting spesific sequences from bam or bam.bai file using samtools
written 8.3 years ago by blooming.daisy333 •

im using samtool to extract specific regions from the sorted bam file/bam.bai file. e.g if i want to extract the sequence from 329729 to 329840 …
Extracting reads from given region
written 8.7 years ago by KVC_bioinfo

Hello all, I am extracting the reads from bam file which are falling in particular region using the following command. samtools view -h input.bam "Chr10:18000-45500" …
Samtools error in novoBreak " [main_samview] region "chrY:13119970-13120970" specifies an unknown r…
written 9.0 years ago by seocd74 •

Hi, I'm trying to run novoBreak with tumor and normal bam file. Tumor files are 'chr19.tumor.bam' and 'chr19.normal.bam' from ICGC-TCGA DREAM Mutation Calling challenge (URL: …
The novoBreak does not work well. Error : region "chr1:66670169-66671169" specifies an unknown refe…
written 9.0 years ago by seocd74 •

Hello. I'm trying to run novoBreak using human genome in AWS. NovoBreak is using samtools and bwa and so on for running and it include …
error in query a BAM file
written 9.3 years ago by shoujun.gu

I'm the beginner of using samtools. I downloaded: efetch -db=nuccore -format=gb -id=AF086833 > AF086833.gb and then try to query a bam file: samtools view -c …
Extracting Reads by position
written 10.4 years ago by michael.madans94 •

I have an alignment file for a protein encoded by about ~3000 bases (I have both bam and sam files), and I am looking to …
Region "Awk" Specifies An Unknown Reference Name
written 14.7 years ago by Zhshqzyc •

<p>Hello, can you look at what is wrong by using samtools command?</p> <pre><code>#!bin/bash samtools view -o raw_genome.txt my.bam awk '{print $3, $4, $10;}' \ raw_genome.txt …