Tip | ||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| ||||||||||||||
Use our summer school reservation (CoreNGS ) when submitting batch jobs to get higher priority on the ls6 normal queue idev -m 180 -N 1 -A OTH21164 -r CoreNGS-Fri .
Note that the reservation name (CoreNGS) is different from the TACC allocation/project for this class, which is OTH21164. |
Table of Contents |
---|
Overview
...
- Convert text-format SAM files into binary BAM files (samtools view -b) and vice versa (just samtools view)
- Sort BAM files by reference coordinates (samtools sort)
- Index BAM files that have been sorted (samtools index)
- Filter alignment records based on BAM flags, mapping quality or location (samtools view)
Since Ad we've seen, since BAM files are binary, they can't be viewed directly using standard Unix file viewers such as more, less and head. We have seen how, butsamtools view can be used to convert binary-format BAM files into text format for viewing. But samtools view also has options that let you do powerful filtering of the output. We focus on this filtering capability in this set of exercises.
...
Code Block | ||||
---|---|---|---|---|
| ||||
idev -m 180 -N 1 -A OTH21164 -r CoreNGS-Fri # or -A TRA23004 # or idev -m 90 -N 1 -A OTH21164 -p development # or -A TRA23004 |
Next set up a directory for these exercises, and copy an indexed BAM file there. This is the yeast_pe.sort.bam file from our The Basic Alignment Workflow exercises.
...