Filtered reads from assembly and raw reads?
1
0
Entering edit mode
5 weeks ago
el97004 ▴ 80

Hi all!

Anyone know of a software tool that can be used to get reads from a de-novo assembly if I have the raw reads? the raw reads went through several steps of preprocessing and filtering and I no longer have the final reads that were used to make the assembly.

Thank you!

sequencing denovo assembly • 297 views
ADD COMMENT
4
Entering edit mode
5 weeks ago
GenoMax 119k

can be used to get reads from a de-novo assembly if I have the raw reads?

If you have an assembly then the individual reads are no longer identifiable since they are now part of a long sequence.

If you are looking to see which part of an assembly a read may belong to then you can simply align the reads to the assembly.

ADD COMMENT
0
Entering edit mode

Thanks GenoMax ! Once aligned and obtained bam file could I not then just convert that to fastq?

ADD REPLY
0
Entering edit mode

Not sure then what your question is. What format is your query data in? That must already be in fastq format? Are you simply looking to extract reads that align to the assembly then (sometimes all reads may not align)?

ADD REPLY
0
Entering edit mode

Yeah I was just trying to find the set of reads that were used to make up the final assembly. Some of those raw reads will definitely not align since they were filtered out (e.g. contaminants).

Format information: assembly: fasta; raw reads: fastq

ADD REPLY
1
Entering edit mode

In that case use an aligner of your choice and extract reads that map: How To Filter Mapped Reads With Samtools

Be sure to get primary alignments (otherwise) there may be multiple copies of reads in extract.

ADD REPLY

Login before adding your answer.

Traffic: 1797 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