This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can I generate VCF file with two different assembled genome fasta?

When I generate VCF file by pair of illumina sequence fastq files with reference assembled genome fasta file, I used BWA mem to assemble illumina fragment and PICARD to remove duplicates. And finally I used gatk4 haplotypecaller for VCF.

But I don't know how to generate VCF file by assembled genome fasta with reference assembled genome fasta. Is there any recommended pipeline to do this?

genomics snp vcf
• 2,245 views • link
written by kgwkk2 • kgwkk2

I used BWA mem to assemble illumina fragment

BWA is an aligner, not an assembler.

by assembled genome fasta with reference assembled genome fasta

You will have to explain this a little more. It is highly unclear what you have and what you want.

• 0 views • link
written by WouterDeCoster WouterDeCoster

I apologize for unclear question. I want to compare same species, different strain genomes by making each genome's VCF file. And type strain was used for standard genome for reference(https://www.ncbi.nlm.nih.gov/assembly/GCF_014117465.1). In case of illumina short read sequences (raw data, fastq), I could align short reads based on reference genome, and make BAM file for VCF analysis. But in case of some genome, there were only completely assembled genomic fasta available. I want to make similar form of VCF files from these assembled genomic fasta with reference genome.

Shortly, I want to compare VCF from illumina sequence and VCF from assembled genome equal position as possible.

• 0 views • link
written by kgwkk2 • kgwkk2

0 answers

No answers yet.

Log in to answer this question.

More posts like this