How to find average read length and number of reads in a bam file using samtools in linux?

How to find average read length and number of reads in a bam file using samtools in linux? I used first samtools stats file.bam | grep "average length" but it never stop running. What is the most efficient way of finding them?

samtools bam

• 3,177 views • link

updated 3.1 years ago by Ram

• written 4.1 years ago by MobiusT •

duplicate of your previous question : Find mean read length of a bam file dataset of reads?

I used first samtools stats file.bam | grep "average length" but it never stop running

what is the size of the bam, how long did you wait ?

• 0 views • link

written 4.1 years ago by Pierre Lindenbaum

The size is 24661059 and I waited aroun 10 mins

• 0 views • link

written 4.1 years ago by MobiusT •

0 answers

No answers yet.

More posts like this

Find mean read length of a bam file dataset of reads?
written 4.1 years ago by MobiusT •

I am working in a remote linux server and I have a dataset of reads in the file with extension bam extension. How to find …
how to identify where multi-mappers, or "non primary hits" are mapping to
written 5.8 years ago by bertb •

Hello, My RNAseq alignment QC is showing me that I have a fair number of "non primary hits", which, as I understand it, means the …
Why there are low assigned reads in featureCounts output and different outputs for quality checks?
written 5.9 years ago by newbie

I'm using rnaseq data (paired-end and strand-specific) which is rRNA depleted with Ribo-Zero Gold rRNA Removal Kit. I observed that most number of reads are …
Extract the alignments from a Bam file by name of the read
written 5.9 years ago by kojix2

Hello. I have a question about Bam files. How can I extract the alignment of a specific named read from a Bam file? I did …
Realign Aligned.sortedByCoord.out.bam to Aligned.toTranscriptome.out.bam for RSEM input
written 7.7 years ago by mkvalluru •

Hi, I was using STAR-RSEM to quantify isoforms (primary input fastq files). Recently I have received Aligned.sortedByCoord.out.bam from the collaborators. However, I can't use these …
calculate fragments depth with pysam
written 8.3 years ago by electronsource75 •

Hello, I am currently working with pysam to calculate sequencing “depth” of fragments. I use a bam file with Pair End reads and I managed …
samtools - Finding Average Fragment size for pair?
written 8.4 years ago by oars

I'm trying to find the average fragment size for reads in my .bam file and want to also exclude read pairs with a fragment size …
STAR alignment, RNA-seq
written 9.3 years ago by prasoonagar •

> I have a a query regarding STAR alignment. I used the following commands to convert the BAM files to fastq (as they were some …
Samtools idxstats fails to load index
written 9.6 years ago by lvande •

I recently used samtools in order to downsample some bam files I have from a sequencing run as such: Samtools view -h -s 0.5 file.bam …
Sam file header and truncated file
written 9.8 years ago by TriS

I posted this on the samtools groups with no luck, so hopefully here i get better luck I'm changing part of the header in a …