Genewiz gave me bam files instead of fasta files- can I convert bam to fasta to have a genome assembly?
I had three bacterial genomes sequenced and Genewiz gave me fastq sequencing reads and bam files. I want the assembly fasta file- can I generate it myself using samtools bam to fasta? This sequencing project was done last year and it is unlikely that Genewiz has this data.
• 1,818 views
•
link
updated
by
Ram
•
written
by
kristina.mahan
•
written
by
kristina.mahan
1 answer
Use the FASTQ files with the program Flye to create a FASTA assembly. https://github.com/fenderglass/Flye
Alternatively, try SPADES.
Both are installable via miniconda I believe.
• 0 views
•
link
written
by
colindaven
Log in to answer this question.
More posts like this
-
assembly using CCS, CLR, CCS_CLR sequences together?
written by ycts •
-
Import a BAM file into SMRT Link
written by a.vincent.0312 •
-
Mauve alignment vs samtools for SNP calling
written by gabriel.jabud •
-
Conversion in fasta and sequence analysis
written by demaioflavio •
-
Issues with mapping reads onto Reference genome using bwa mem
written by kristina.mahan
-
How do we find the number of sequences with terminal mismatches?
written by MAPK
-
Variant calling and downstream analysis
written by arash.iranzadeh1980 •
-
Which is the best input file for denovo assembly with Illumina paired data: Fastq, Fasta or Fasta i…
written by minesrebollo •
-
kmergenie to estimate genome size and then sequencing coverage
written by cook284 •
-
Getting genbank annotation file from JGI
written by bioPiraten •
I'm not quite sure what you want. You want to make a de novo assembly based on your reads? You can do that with the fastq. You want to recreate the genome that the reads were aligned to to make the bam? I don't think you can do that easily.
Did you ask them to create an assembly? If not then you must assemble it yourself, for this please google for tutorials on genome assembly. The raw sequencing data are just the genome chopped into pieces and sequenced, this is how NGS works.