MALT extract reads and BWA aln
0
0
Entering edit mode
6 months ago

Hello, It is my first using of MALT.

I ran malt-build with all refseq bacteria genomes with mapping file --acc2taxonomy megan-nucl-Jan201.db

malt-run with ancient reads with such flags -id 85 -m BlastN -at SemiGlobal -top 1 -supp 0.01 -mq 100. open .rma6 file in Megan desktop, click at the species, then File --> Extract Reads

I got .fasta file with reads that assigned to the species.

Then I decided to validate this assignments. I’ve taken ID of the reference from .blast file (also download after clicking at the specie) and run BWA: bwa aln ref_index extracted_reads_from_megan.fasta > result.bwa

bwa samse ref_index result.bwa extracted_reads_from_megan.fasta > result.sam

And there are no hits in samtools flagstat. Then I tested with reference genome of this species from Refseq and also there are no hits. Should I have more strict parameters with build-run? When I tried to blastn some of these reads on the site, it showed the hits with [Eukaryotic synthetic construct chromosome 20]. The reads for malt-run was from unmapped of human. What could be the mistake?

A lot of thanks, Valery

ancient metagenomic megan malt • 169 views
ADD COMMENT

Login before adding your answer.

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