benplm
11/18/2018 - 8:52 PM
Bioinformatics
# quickly prints total number of reads from indexed bam file
samtools idxstats example.bam | cut -f3 | awk 'BEGIN {total=0} {total += $1} END {print total}'# quickly prints total number of reads from indexed bam file
samtools idxstats example.bam | cut -f3 | awk 'BEGIN {total=0} {total += $1} END {print total}'