...
| Code Block |
|---|
cut -f 1 alignments/human_rnaseq_mem.sam | sort | uniq -c | less #This gives you a view where each read is listed next to the number of entries it has in the SAM file cut -f 1 alignments/human_rnaseq.sam | sort | uniq -c | awk '{print $1}' | sort | uniq -c | less #This gives essentially a histogram of the number of times each read aligned - a plurality of reads aligned twice, which seems reasonable since these are all reads crossing a junction, but plenty aligned more or less cut -f 1 alignments/human_rnaseq.sam | sort | uniq | wc -l #This gives a better idea of the alignment rate, which is how many reads aligned at least once to the genome. Divided by the number of reads in the original file, the real alignment rate is around 64.19 %. # PLEASE NOTE: some of these one-liners are only reasonably fast if the files are relatively small (around a million reads or less). For bigger files, there are better ways to get this information, mostly using samtools, which is a utility explicitly designed for manipulating SAM/BAM files. We'll cover samtools in the next section. |
...