Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Output BCF file. This is a binary form of the text Variant Call Format (VCF).

Code Block
titleThis is *one* command. Put it all on one line.
samtools mpileup -u -f samtools_bowtie/NC_012967.1.fasta samtools_bowtie/SRR030257.sorted.bam > samtools_bowtie/SRR030257.bcf

...

Convert BCF to human-readable VCF:

Code Block
titleThis is *one* command. Put it all on one line.
bcftools view -v -c -g samtools_bowtie/SRR030257.bcf > samtools_bowtie/SRR030257.vcf

...