Changes

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]

Navigation menu