From Genome Analysis Wiki
Jump to navigationJump to search
92 bytes removed
, 05:23, 8 February 2012
Line 179: |
Line 179: |
| </source> | | </source> |
| | | |
− | | + | To complete our example analysis, we could run: |
− | Here is how that might work:
| |
| | | |
− | bin/TrioCaller --shotgun result/chr20.mpileup.vcf --pedfile ped/triocaller.ped --rounds 20 --prefix result/chr20.triocaller | + | TrioCaller --shotgun result/chr20.mpileup.vcf --pedfile ped/triocaller.ped --rounds 20 --prefix result/chr20.triocaller |
| | | |
− | The format of output file is same as the input file. Here, you might replace chr20.mpileup.vcf with other vcf file (e.g. chr20.polymutter.vcf). <br /> | + | The format of output file is same as the input file. Again, you can review the contents of the updated VCF file using the more command: |
− | Again, you can review the contents of the updated VCF file using the more command: | |
| | | |
| more result/chr20.triocaller.vcf | | more result/chr20.triocaller.vcf |
| | | |
| + | All right. Congratulations! You have come to the end and learned basic skills for the variant calling in trios. |
| | | |
− | All right. Congratulations! You have come to the end and learned basic skills for the analysis of sequence data. <br />
| |
| If you have any question, feel free to contact Wei Chen at [mailto:weichen.mich@gmail.com weichen.mich@gmail.com] or Goncalo Abecasis at [mailto:goncalo@umich.edu goncalo@umich.edu] | | If you have any question, feel free to contact Wei Chen at [mailto:weichen.mich@gmail.com weichen.mich@gmail.com] or Goncalo Abecasis at [mailto:goncalo@umich.edu goncalo@umich.edu] |