This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to align many single end fastq files from the same study on each other without using any reference genome using python?

I have many single end fastq files that I downloaded from SRA data base I want to align these files on each other to make a consensus genome that I will consider it as my reference genome how can I make this alignment by using python?

fastq fasta python vcf

I want to align these files on each other to make a consensus genome

That is not how you create a consensus for a genome. You would either need to align to a reference and/or assemble the data.

but is there is a way to make alignment on these files without using ref genome by using python?

Aligning millions of fastq reads to each other makes no logical sense. python is a programming language not a tool for analysis.

1 answer

That is called genome assembly, not alignment, also python is not used for that. Check assembly tools such as SPAdes.

Log in to answer this question.