https://bioinformaticsonline.com/thewire/thread/40284? https://bioinformaticsonline.com/thewire/view/40284 Mon, 25 Nov 2019 09:23:51 -0600 https://bioinformaticsonline.com/thewire/view/40284 #Estimate the #insert #size with mapped reads $ head -10000 mappings.sam | awk '{ if ($9 > 0) { N+=1; S+=$9; S2+=$9*$9 }} END { M=S/N; print "n="N", mean="M", stdev="sqrt ((S2-M*M*N)/(N-1))}'

]]> LEGE