First genomic coordinate that Nanopore read aligns to
I aligned Nanopore reads to a genome of interest and outputted a sorted BAM file. I want to determine the genomic coordinate that corresponds to the first base of each aligned read aka where in the genome the start of each read aligns. Any suggestions?
1 answer
From the SAM file specification (section 1.4):
- POS: 1-based leftmost mapping POSition of the first CIGAR operation that “consumes” a reference base (see table below). The first base in a reference sequence has coordinate 1. POS is set as 0 for an unmapped read without coordinate. If POS is 0, no assumptions can be made about RNAME and CIGAR.
In other words, the first position of the alignment should be the fourth field in the SAM/BAM file.
• 0 views
•
link
written
by
Dave Carlson
Log in to answer this question.
More posts like this
-
extracting data from single end aligned bam file - chromosome, strand, start and stop from cigar
written by Jalil Sharif •
-
How to extract paired (PE) reads where the single (SE) reads align to different genes from a BAM fi…
written by chrisk •
-
Read depths on whole genome (WGS)
written by Jean-Charles •
-
whole exome sequence fragment size analysis
written by adR
-
Illumina Overlapping R1/R2 reads, error-correction in python
written by Mick •
-
Replacing reference with aligned reads
written by elena.syp •
-
What is the definition of "read depth" vs "coverage"? (again...)
written by ariel.balter
-
counting the reads mapped to specific coordinate in the genome
written by alirezamomeni707 •
-
Determine relative position (start and end) of chimeric secondary alignment within original long re…
written by mattia.bosio85 •
-
tool simulate paired-end reads and report where each read come from
written by Chen Sun
