[BioC] pileup function in ShortRead package

Nora Rieber n.rieber at dkfz-heidelberg.de
Thu Oct 8 16:55:27 CEST 2009


Hi there,

I'm using the pileup function in the ShortRead package and it seems I
call it with the wrong arguments, or that there is some kind of bug in
the function. All my reads are piled up in forward direction only, even
if they were actually mapped to the backward strand.

Here's a small (nonfunctional since I assume I cannot append data to
this e-mail) example of how I call the pileup function:

#read in the MAQ alignment file:
aln_s_7_MCIP4<-readAligned("/home","sequenceALN.txt", type="MAQMapview")
aln_Pat_subset<-aln_Pat[1:10]
pos_Patient_subset<-position(aln_Pat_subset)
str_Patient_subset<-strand(aln_Pat_subset)
#use only a small subset (all positions are on the same chromosome 1)

subset_pileuptest<-pileup(pos_Patient_subset, 36,
max(pos_Patient,na.rm=TRUE)+37,dir=str_Patient_subset)

Is there anything wrong with my dir argument? I checked several places
all over the genome for the complete alignment and my reads all seem to
be considered on the + strand, even when they are not. Can anybody help?

This is my session info:

> sessionInfo()
R version 2.9.2 (2009-08-24)
x86_64-unknown-linux-gnu

locale:
LC_CTYPE=de_DE.UTF-8;LC_NUMERIC=C;LC_TIME=de_DE.UTF-8;LC_COLLATE=de_DE.UTF-8;LC_MONETARY=C;LC_MESSAGES=de_DE.UTF-8;LC_PAPER=de_DE.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=de_DE.UTF-8;LC_IDENTIFICATION=C

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] ShortRead_1.2.1   lattice_0.17-25   BSgenome_1.12.3   Biostrings_2.12.9
[5] IRanges_1.2.3

loaded via a namespace (and not attached):
[1] Biobase_2.4.1 grid_2.9.2    hwriter_1.1

Thanks,
Nora



More information about the Bioconductor mailing list