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.''' |