Versions Compared

Key

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

...

Use the subcommands bedtools intersect and bedtools subtract we can find equal and different predictions between mappers. Try to figure out how to to do this on your own first. Hint: Remember that adding > output.vcf to the end of a command will pipe the output that is to the terminal into a file, so that you can save it.

Expand
I'm stuck? Show me the commands...
I'm stuck? Show me the commands...

Load Bedtools.

Code Block
module load bedtools

Finding common mutations.

Code Block
bedtools intersect -a bowtie.vcf -b bwa.vcf > common_bowtie_bwa.vcf

Finding mutations that are unique for each mapper.

Code Block
bedtools subtract -a bowtie.vcf -b common_bowtie_bwa.vcf > unique_bowtie.vcf
bedtools subtract -a bwa.vcf -b common_bowtie_bwa.vcf > unique_bwa.vcf

...