...
Code Block |
---|
samtools mpileup -uf samtools_bowtie/REL606.5.fasta samtools_bowtie/sorted_REL606.5.bam \|bcftools view -vcg - \> samtools_bowtie/output.vcf
|
...
Results on Lonestar
Code Block |
---|
|
/corral-repl/utexas/BioITeam/ngs_course/mapping/samtools_bowtie/sorted_REL606.5.bam
/corral-repl/utexas/BioITeam/ngs_course/mapping/samtools_bowtie/output.vcf
|
Code Block |
---|
|
/corral-repl/utexas/BioITeam/ngs_course/mapping/samtools_bwa/sorted_REL606.5.bam
/corral-repl/utexas/BioITeam/ngs_course/mapping/samtools_bwa/output.vcf
|
Exercise.
VCF format has Allele Frequency tags denoted by AF1. Try the following command to see what value we have in our files.
...