Count number of unique alignments to each chromosome from sam file

Count number of unique alignments to each chromosome from sam file

0

Hi,

I want to count the number of unique reads aligned to each chromosome after mapping to the genome. I don’t mind if something is multimapped, say if read A maps to both chr3 and chr6 that is fine, but if read A maps to chr6 twice I want to only count it once.

I’ve pieced together 2 commands which I think should both do this, but not sure whether they are right/ one is preferred.

option 1:
samtools view file.bam | cut -f 1,3 | sort | uniq -c > chr_stats.txt

option 2:
awk ‘ $1 !~ /@/ {print $3}’ samtools view file.bam | sort | uniq -c > chr_stats.txt

1) convert a bam to sam
2) pull out column 1 and 3 from the sam file (read ID and chr)
3) sort these
4) only keep unique lines
5) count the number of unique IDs per chromosome
6) output to result.txt file

I’m new to piping commands using bash so would appreciate if someone could clarify whther these are doing what I want them to do

Thanks 🙂


samtools


sort


unix


uniq


sam

• 36 views


updated 5 minutes ago by

107k

written 2 hours ago by

▴

30

Read more here: Source link