This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to find snp changes between pairs of samples that have genomes sequenced

Hi everyone I'm new to bioinformatics and I want to identify significant snp changes between pairs of samples (patient isolates that are resistant and non-resistant to fluconazole) of cryptococcus neoformans which have had their genomes sequenced and [put on ebi database[1] I've tried looking for a way to do this on fungidb but am unable to find the right tool. Would anybody be able to give me some advice on how to proceed? Thanks.

variant genome
• 1,187 views • link
written by bryanchang2016 • bryanchang2016

Here is an outline of what you can do:

  1. Get original fastq sequence for samples from EBI
  2. Get reference genome sequence for Cryptococcus neoformans.
  3. Index the genome with an aligner of your choice (bwa, bowtie2, bbmap etc)
  4. Align your samples to the reference.
  5. Sort/index alignment files.
  6. Call variants on the two alignments
  7. Use bcftools isec to find changes that are sample specific.
• 0 views • link
written by GenoMax GenoMax

I have compared each sample to reference on veupathdb and have obtained vcf files for each sample. But is there a way I can do bcftool on veupathdb/fungidb to find changes between two vcf files (resistant/nonresistant)?

• 0 views • link
written by bryanchang2016 • bryanchang2016

0 answers

No answers yet.

Log in to answer this question.

More posts like this