how to align many single end fastq files from the same study on each other without using any reference genome using python?
1
0
Entering edit mode
3.4 years ago
FadyNabil ▴ 20

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 • 1.5k views
ADD COMMENT
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode
3.4 years ago
JC 13k

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

ADD COMMENT

Login before adding your answer.

Traffic: 2516 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6