This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to align fastq files to a fasta file using biopython to produce VCF file

I have many fastq files that I want to align to a fast file using biopython package

fastq biopython fasta vcf
• 2,869 views • link
updated by Ram Ram • written by FadyNabil • FadyNabil

Using biopython for this task is unlikely to be the best approach, can you elaborate why you want to do that?

• 0 views • link
written by WouterDeCoster WouterDeCoster

I am running on a project that needs me to make alignment for fastq files

• 0 views • link
written by FadyNabil • FadyNabil

Yes, but BioPython is not what you want. Alignment of fastq files (NGS sequencing data, is it?) is done with specialized tools such as bwa or bowtie2. Please google for these tools, there are literally hundreds of posts, tutorials and manuals on how to do that. The entire process is called variant calling, there are pipelines for it, please search for these. Something like GATK best practices or bcftools mpileup as good keywords to start with.

• 0 views • link
written by ATpoint ATpoint

So what is the problem, you opened a new question? Feel free to add a comment if you need further help. Did you get stuck somewhere, if so where? There is no need to be ashamed asking for help, but please do not just open identical questions, that only annoys people :)

• 0 views • link
written by ATpoint ATpoint

Please do not delete posts that have received feedback.

• 0 views • link
written by Ram Ram

0 answers

No answers yet.

Log in to answer this question.

More posts like this