Changes

From Genome Analysis Wiki
Jump to navigationJump to search
157 bytes removed ,  11:10, 4 September 2011
Line 28: Line 28:  
* --reference
 
* --reference
 
** The reference genome is the same as karma reference genome. If the index files do not exist, qplot will create the index files using the input reference fasta file.
 
** The reference genome is the same as karma reference genome. If the index files do not exist, qplot will create the index files using the input reference fasta file.
** Or use the latest karma to generate reference index files using the command:
  −
karma create -i reference.fa
      
* --dbsnp
 
* --dbsnp
Line 37: Line 35:     
Although GC content can be calculated on fly each time, it is much more efficient to load a precomputed GC content from a file. To generate the file, use the following command
 
Although GC content can be calculated on fly each time, it is much more efficient to load a precomputed GC content from a file. To generate the file, use the following command
  /home/bingshan/code/calcGCContent/genGCContent -r referenceByKarma -w windowSize -o gccontent.out
+
  qplot -r reference.fa -w windowSize -o reference.gc
    
''Note'': Before running the qplot, it is critical to check how the chromosome numbers are coded. Some bam file use just numbers, others use chr + numbers. '''You need to make sure that the chromosome numbers from reference and dbsnp are consistent with the bam file.'''
 
''Note'': Before running the qplot, it is critical to check how the chromosome numbers are coded. Some bam file use just numbers, others use chr + numbers. '''You need to make sure that the chromosome numbers from reference and dbsnp are consistent with the bam file.'''
480

edits

Navigation menu