This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error Import File With "Import" Of Rtracklayer

Hi, This is related to a the question found here (http://biostar.stackexchange.com/questions/3414?sort=newest#sort-top). I am trying to import a gff file using import("file.gff") which contains a set of chromosomal positions for all chrmosomes in mouse. When I use the rtracklayer command 

genes = import("file.gff")  
names(genes)= "chr1" "chr10" "chr11" "chr12" "chr13" "chr14" "chr15" "chr16" "chr17" "chr18" "chr19" "chr2" "chr3" "chr4" "chr5" "chr6" "chr7" "chr8" "chr9" "chrX" "chrY"

How do I make sure that the import doesn't cause the chromosomal order to change I want it to import it into teh dataRange genes in the order 

"chr1" "chr2" "chr3" "chr4" "chr5" "chr6" "chr7" "chr8" "chr9" "chr10" "chr11" "chr12" "chr13" "chr14" "chr15" "chr16" "chr17" "chr18" "chr19" "chrX" "chrY"

Thanks in Advance

chip-seq next-gen sequencing r
• 3,728 views • link
updated by Biostar Biostar • written by ska007 • ska007

The names are sorted in lexical order, however the order of chromosome names is irrelevant for further analysis e.g. overlap computation in GenomicRanges/IRanges. The other way around that also makes sense, if the order of names appearing makes any difference to your analysis, then there is possibly something wrong with your analysis.

• 0 views • link
written by Michael Dondrup Michael Dondrup

1 answer

I don't think the order in which the data are imported has any bearing on your problem.

First, I assume that rtracklayer imports entries in the order in which they appear in the GFF file.

Second, I think you are confused by the output of names(genes). In fact, the chromosomes are sorted in that output. However, the sorting is by the ASCII value of the characters in the strings, not by chromosome number - hence "chr10" comes after "chr1" but before "chr11". This is a quite normal way to sort strings and tells you nothing about how the data are stored.

• 1 views • link
written by Neilfws Neilfws

Hi, Well I thought so too at first. But then to check if this is what was happening. I extracted say genes$chr2 from the .gff file and the reads I had from the bam file for chr2 and I computed their overlaps. It gave me the correct overlap.

• 0 views • link
written by ska007 • ska007

I checked if the gff file was imported in that order (using the command spaces(genes)) and it was. The file.gff as such has the correct order .

• 0 views • link
written by ska007 • ska007

OK, so we're agreed that any issues you have with rtracklayer don't relate to chromosome order? I don't really understand your problem.

• 0 views • link
written by Neilfws Neilfws

Log in to answer this question.

More posts like this
  • How can I get this sorted by ascending order? (when the number and character is together)
    written by woojoy14 • woojoy14
    samtools view NA12878_dedup.bam | cut -f 3 | sort -u > chr.list8 > chr1 chr10 chr11 chr12 chr13 chr14 chr15 chr16 chr17 chr18 chr19 chr2 …
  • Removing a chromosome from a fasta file
    written by t86dan • t86dan
    Hello, I have a fasta file for human Genome GRCh37 Reference Assembly, but for some reason chromosome 20 is repeated and I would like to …
  • MAFTOOLS HELP - Mutational Signature
    written by David_emir David_emir
    Hi All, I am trying to run mutational signatures on my MAF file, using maftools (https://github.com/PoisonAlien/maftools). while trying to run this i am constantly My …
  • Error with correctGCBias in deepTools
    written by clmin • clmin
    Hello, I'm trying to correct the GC bias with deepTools. I successfully used computeGCBias > computeGCBias -b file.bam --effectiveGenomeSize 2150570000 -g mm9.2bit -l 200 --GCbiasFrequenciesFile …
  • modifying the header of VCF file
    written by zengtony743 zengtony743
    Thank you for whoever can help me!!! 1) When I ran $ `java -jar /hpf/tools/centos6/gatk/3.6.0/GenomeAnalysisTK.jar -R genome.fa -T SelectVariants --variant Newsnp.vcf -o testSNP.vcf -sn AKRJ …
  • coverageBed with -g option ?
    written by bioguy24 bioguy24
    I am running a coverageBed and getting the below error. ``` cmccabe@DTV-A5211QLM:~/Desktop/NGS$ coverageBed -d -sorted -a /home/cmccabe/Desktop/NGS/bed/bedtools/xgen_targets_sorted.bed -b /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam > /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_output.txt ERROR: Database file /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam …
  • problem with methylKit of the annotation
    written by hua.peng1314 hua.peng1314
    I am dealing with RRBS data follow [the manual][1]. I met some problem when doing annotation.The problem showed below. ``` > read.transcript.features("mm9.refseq.genes.bed",remove.unsual=TRUE,up.flank=1500,down.flank=500,unique.prom=TRUE) > diffAnn <- …
  • CONTRA problem. It does not stop.
    written by mafonso • mafonso
    Hi, I am trying to run an algorithm for CNV detection and was testing CONTRA. The test file works but using my own files the …
  • How To Retrieve Chromosome Lengths Using Bsgenome Package?
    written by James Ashmore • James Ashmore
    <p>Hello, I am following along with the <a href='http://www.ebi.ac.uk/training/online/sites/ebi.ac.uk.training.online/files/user/42/documents/embo_chipseq_course_1101026.pdf'>EBI's online course for ChIP-Seq analysis</a> and have encountered an error when trying to retrieve chromosome lengths:</p> …
  • What Is 3-P And 5-P Flanking Sequence? How To Find Them Using Pygr?
    written by Use • Use
    <p>For example, </p> <pre><code>&gt;&gt;&gt; hg17=worldbase.Bio.Seq.Genome.HUMAN.hg17() &gt;&gt;&gt; len(hg17) 46 &gt;&gt;&gt; hg17.keys() ['chr6_random', 'chr19_random', 'chr8_random', 'chrY', 'chrX', 'chr13', 'chr12', 'chr11', 'chr15_random', 'chr17', 'chr16', 'chr15', 'chr14', 'chr19', 'chr18', …