mkatari-bioinformatics-august-2013-gatknotes
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
mkatari-bioinformatics-august-2013-gatknotes [2014/06/11 13:46] – mkatari | mkatari-bioinformatics-august-2013-gatknotes [2014/07/02 15:19] – mkatari | ||
---|---|---|---|
Line 96: | Line 96: | ||
</ | </ | ||
- | Now we merge the bam files and then sort and index them | + | In some cases there may be a need to clean the sam/bam file(s) (soft-trimming the coordinates). To do this use CleanSam in Picard tools. You may want to just do it to all to avoid the error in a workflow, but it may not be necessary. |
< | < | ||
- | java -jar / | + | java -jar / |
| | ||
+ | | ||
+ | </ | ||
+ | |||
+ | Now we merge the bam files and then sort and index them. If you cleaned the bam file, remember to use the cleaned ones. | ||
+ | |||
+ | < | ||
+ | java -jar / | ||
+ | | ||
| | ||
| | ||
Line 110: | Line 118: | ||
- | Finall | + | Finally |
< | < | ||
Line 121: | Line 129: | ||
-glm SNP \ | -glm SNP \ | ||
-o PTC_human.gatk.vcf | -o PTC_human.gatk.vcf | ||
+ | |||
+ | </ | ||
+ | |||
+ | If you want to load the vcf file into IGV, remember to index it first. | ||
+ | < | ||
+ | module load igvtools | ||
+ | igvtools index PTC_human.gatk.vcf | ||
+ | </ | ||
+ | |||
+ | If you would like to generate a table of from the vcf file use the following command | ||
+ | < | ||
+ | java -jar / | ||
+ | -R PTC_Human.fasta | ||
+ | -T VariantsToTable \ | ||
+ | -V PTC_human.gatk.vcf \ | ||
+ | -F CHROM -F POS -F ID -F QUAL -F AC \ | ||
+ | -GF GT -GF GQ \ | ||
+ | -o PTC_human.gatk.vcf.table | ||
</ | </ |
mkatari-bioinformatics-august-2013-gatknotes.txt · Last modified: 2016/08/17 08:37 by mkatari