Skip to content

Instantly share code, notes, and snippets.

@sirusb

sirusb/GAlignments_example.r

Created February 19, 2015 14:51
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save sirusb/33826333fae542855e35 to your computer and use it in GitHub Desktop.
Save sirusb/33826333fae542855e35 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
GAlignments_example.r
require(pasillaBamSubset)
require(GenomicAlignments)
# BAM نقراء الملف
bamFile <- untreated1_chr4()
gReads <- readGAlignmentsFromBam(bamFile)
#GAlignments object with 204355 alignments and 0 metadata columns:
# seqnames strand cigar qwidth start end width njunc
# <Rle> <Rle> <character> <integer> <integer> <integer> <integer> <integer>
# [1] chr4 - 75M 75 892 966 75 0
# [2] chr4 - 75M 75 919 993 75 0
# [3] chr4 + 75M 75 924 998 75 0
# ... ... ... ... ... ... ... ... ...
# [204353] chr4 + 75M 75 1348268 1348342 75 0
# [204354] chr4 - 75M 75 1348449 1348523 75 0
# [204355] chr4 - 75M 75 1350124 1350198 75 0
# -------
# seqinfo: 8 sequences from an unspecified genome
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment