I just ran an assembly on yeast genomes using Flye and I want to polish those assemblies with Pilon but it requires a sorted BAM file.
How do I make a BAM file of the resulting assembled.fasta?
1
-
1$\begingroup$ What reads do you have? Only long read or short read or whatever? You need to align the starting reads to your assembly and then supply that bam, as implied by the user172818 answer. $\endgroup$– Maximilian PressSep 19, 2022 at 22:54
Add a comment
|
1 Answer
$\begingroup$
$\endgroup$
1
Run a short-read mapper. For example:
bwa index assembled.fasta
bwa mem -pt16 assembled.fasta read1.fq.gz read2.fq.gz \
| samtools sort -m4G -@4 -o align.bam -
answered Sep 19, 2022 at 15:51
-
$\begingroup$ This is good if the user has Illumina, that isn't clear (though it is implied by using pilon). $\endgroup$ Sep 19, 2022 at 22:53