bedtools questions ! can't sort/index by samtools after bedToBam command

Hi all ! After I use STAR to alignment , I got the BAM file , this BAM file could index successfully by samtools , but after converting this BAM file to BED with command

 bedtools bamtobed -i sample.bam > sample.bed

and then,I use bedtobam convert it to bam

 bedtools bedtobam -i /sample.bed -g genome > new_sample.bam

I want to index this new BAM ,so I run

samtools index new_sample.bam

But I got this error finally

[E::hts_idx_push] Unsorted positions on sequence #1: 248926319 followed by 127352

samtools index: failed to create index for "new_sample.bam"

Why? Is there any wrong that makes this happen? How I can fix this problem ?

And I try to use samtools view to see what happen,but I got this,here is the command and error

samtools view -h new_sample.bam | head -10

@HD VN:1.0 SO:unsorted

@PG ID:BEDTools_bedToBam

VN:Vv2.29.2 [main_samview] truncated file.

And samtools sort failed either

samtools sort: truncated file. Aborting

BTW:samtools version=1.9 bedtools version=2.29.2

I new to NGS ,any reply is greatly appreciated!

rna-seq software error

• 1,505 views • link

written 5.8 years ago by zhaizhaoyu7 •

the command you posted does not look right, that slash in front of sample.bed indicates a file in the root of you computer, most likely it is not what you ran:

bedtools bedtobam -i /sample.bed -g genome > new_sample.bam

overall something is not quite right with your sam file. Investigate the SAM file you've created.

samtools view new_sample.bam

the follow the errors if you get any

• 0 views • link

written 5.8 years ago by Istvan Albert

Thanks for your reply!

I cut my absolute path because it's too long ,sorry that I left a '/'.

Apart from this,I didn't got any SAM file , I only have the BAM file in coordinate after I ran STAR.And I can index this BAM successfully. That's what I'm confused,there are only two steps:bedtools bamtobed and bedtobam,then I got the wrong BAM .

Here is samtools view I got without -h

[main_samview] truncated file.

Any reply will be appreciated.

• 0 views • link

written 5.8 years ago by zhaizhaoyu7 •

0 answers

No answers yet.

More posts like this

samtools index error
written 5.5 years ago by pt.taklifi

Hello everyone I am trying to extract uniquely aligned reads and sort them with `samtools` with the following command samtools view -h -f 2 -q …
cant make index file for my bam file that sorted by -n
written 5.5 years ago by meisam.radfar •

hi i have a bam file (from tophat2) i use : samtools view -b -h -T Triticum_aestivum.fa 19358.bam -o 119358.bam then: samtools sort -n 119358.bam …
Samtools error -1
written 5.8 years ago by rin •

Hi all, I am analyzing some data from [GSE156255][1] and I keep getting an error when ordering the bam files. I tried to look it …
getting error when indexing bam files
written 7.1 years ago by Sara

I have aligned my fastq files to hg19 using this command: bwa mem hg19.fa /DATA/myfile.fastq.gz then I made bam files from sam files using: samtools …
samtools sort: truncated file. Aborting
written 7.3 years ago by m.t.lorenc •

Hi, I got an error by using this the below command: samtools sort -O SAM -n -@ 2 -o OZBenth2_.fastp.fq.gz.name-sorted.sam OZBenth2_.fastp.fq.gz.sam [W::sam_read1] Parse error at …
After sorting bam file not able to index ?
written 7.6 years ago by sunnykevin97

Hi after sorting not able to index the bam file showing error, how to sort it out ? samtools index -b -@ 4 A-DS1546_1__R1_dedup_2_ReadGroupssorted.bam >[E::hts_idx_push] …
samtools arror: [E::hts_idx_push] NO_COOR reads not in a single block at the end 0 -1
written 8.1 years ago by marongiu.luigi

Dear all, I have aligned a file with the command: bwa mem -t 10 <refFile> <input_1.fq.gz> <input_2.fq.gz> -o <output>.sam then I converted the SAM to …
Error while generating index files using samtools
written 8.5 years ago by bisht20diksha •

Hello. I am trying to identify differential alternative splicing between samples using splAdder. It requires sorted bam files and index files. I have sorted bam …
bgzf_blocked error when using samtools
written 8.7 years ago by milanlove •

when i used samtools sort for sorting .bam file, following error was appeared. any body know that how i can solve it? [E::bgzf_read] Read block …
How to specify the sort based on name in samtools sort?
written 8.7 years ago by anc.informatics •

I would like to know how to issue the command when sorting a BAM file based on name. The problem is that when I used …