Question: How to convert fastq file with raw reads to draft genome fastq file ?
0
gravatar for gskbioinfo143
2.9 years ago by
India
gskbioinfo14360 wrote:

Hi sir/mam,

i have done mapping of Streptococcus pneumonia raw sequenced data with reference Streptococcus pneumonia ATCC 700669 using bowtie2. i have converted from .bam file to .fasta file. the fasta file it has the raw reads - that is the fasta file just stripped of the quality scores know how to get draft genome fasta file.

fasta read image https://ibb.co/cMDzRQ

Need help

thanks

fastq • 1.2k views
ADD COMMENTlink modified 2.8 years ago by Biostar ♦♦ 20 • written 2.9 years ago by gskbioinfo14360
2

It sounds like you are trying to assemble a library. In that case, I'd suggest preprocessing the data to clean it up, then assembling it with Spades to yield a draft genome. Mapping to a related genome is not usually part of this process unless you are trying to reduce contamination or you want to do a reference-guided assembly.

ADD REPLYlink modified 2.9 years ago • written 2.9 years ago by Brian Bushnell17k

Yes sir, im actually looking for reference based mapping of reads

ADD REPLYlink written 2.9 years ago by gskbioinfo14360
1

not sure what you mean by "how to get draft genome fasta file", please explain better.

ADD REPLYlink written 2.9 years ago by Michael Dondrup47k

You have to make an assembly "assembling", not mapping, if your reference is a very conserved specie you can make a guided assembly using spades, idba-hybrid or trinity, if it doesn´t, you have to make a "denovo" assembly in which case I would recommend you spades.

ADD REPLYlink written 2.9 years ago by Buffo1.8k

http://spades.bioinf.spbau.ru/release3.10.1/manual.html

we cant do reference based mapping (assembly) with spades as given in spades manuel

thanks need your suggestions

ADD REPLYlink written 2.9 years ago by gskbioinfo14360
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1130 users visited in the last hour