```r
  require(GenomicRanges)
  require(GenomicAlignments)
  ```
  
  ```
  ## Loading required package: GenomicAlignments
  ## Loading required package: Rsamtools
  ```
  
  ```r
  require(Rsamtools)
  
  ## our 'region of interest'
  roi <- GRanges("chr14", IRanges(19653773, width=1)) 
  ## sample data
  require('RNAseqData.HNRNPC.bam.chr14')
  ```
  
  ```
  ## Loading required package: RNAseqData.HNRNPC.bam.chr14
  ```
  
  ```r
  bf <- BamFile(RNAseqData.HNRNPC.bam.chr14_BAMFILES[[1]], asMates=TRUE)
  ## alignments, junctions, overlapping our roi
  paln <- readGAlignmentsList(bf)
  j <- summarizeJunctions(paln, with.revmap=TRUE)
  j_overlap <- j[j %over% roi]
  
  ## supporting reads
  paln[j_overlap$revmap[[1]]]
  ```
  
  ```
  ## GAlignmentsList of length 8:
  ## [[1]] 
  ## GAlignments with 2 alignments and 0 metadata columns:
  ##       seqnames strand      cigar qwidth    start      end width njunc
  ##   [1]    chr14      -  66M120N6M     72 19653707 19653898   192     1
  ##   [2]    chr14      + 7M1270N65M     72 19652348 19653689  1342     1
  ## 
  ## [[2]] 
  ## GAlignments with 2 alignments and 0 metadata columns:
  ##       seqnames strand     cigar qwidth    start      end width njunc
  ##   [1]    chr14      - 66M120N6M     72 19653707 19653898   192     1
  ##   [2]    chr14      +       72M     72 19653686 19653757    72     0
  ## 
  ## [[3]] 
  ## GAlignments with 2 alignments and 0 metadata columns:
  ##       seqnames strand     cigar qwidth    start      end width njunc
  ##   [1]    chr14      +       72M     72 19653675 19653746    72     0
  ##   [2]    chr14      - 65M120N7M     72 19653708 19653899   192     1
  ## 
  ## ...
  ## <5 more elements>
  ## ---
  ## seqlengths:
  ##                   chr1                 chr10 ...                  chrY
  ##              249250621             135534747 ...              59373566
  ```