Question: Does Last aligner only outputs the mapped reads to a bam file ?
0
gravatar for pinninti1991reddy
16 months ago by
pinninti1991reddy30 wrote:

Hi, I used Last on Illumina genome dataset, after generating bam file its showing only mapped reads. I'm not able to figure it out. Why it is not showing any other statistics ?

Aligner - Last
Data-set - Human genome
Read-length - Illumina 250 bp (PE reads)

##CMD##
./lastal -Q1 -P10 hg38 lastl/SRR1q26.fastq lastl/SRR2q26.fastq > SRR.maf
After converting SRR.maf ----> SRR.last.bam

##CMD##
`samtools flagstat SRR.last.bam`
1067609571 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 secondary
0 + 0 supplementary
0 + 0 duplicates
**1067609571 + 0 mapped (100.00% : N/A)**
0 + 0 paired in sequencing
0 + 0 read1
0 + 0 read2
0 + 0 properly paired (N/A : N/A)
0 + 0 with itself and mate mapped
0 + 0 singletons (N/A : N/A)
0 + 0 with mate mapped to a different chr
0 + 0 with mate mapped to a different chr (mapQ>=5)

Is it problem with the alignment ? How to sort it out. Thanks!

ADD COMMENTlink modified 16 months ago by RamRS24k • written 16 months ago by pinninti1991reddy30
1

I'm curious - is this a valid use for lastal? The website says it does local alignments from FASTA files to genomes, why are you using this on FASTQ reads that come with a lot of possibilities for sequencing error and whatnot?

ADD REPLYlink modified 16 months ago • written 16 months ago by RamRS24k

To put it differently, why not use one of the tried and extensively tested programs, such as bwa, bowtie2, hisat2, star, bbmap?

edit: never mind my comment, you indeed have been testing several different programs:

Comparing alignment stats BWAmem vs. Soap2 aligner

gSNAP aligner generating truncated sam file ?

ADD REPLYlink modified 16 months ago • written 16 months ago by h.mon28k

Definitely not. Last is an aligner like Blast and Clustal are aligners. It's not a mapper.

LAST ~= (BLAST - B)

ADD REPLYlink modified 16 months ago • written 16 months ago by Joe15k
1

Which is why I asked OP. Technically, local alignment for reads makes a little sense, as reads are sequences after all. But one should be aligning reads semi-globally, so they are aligned globally, but the alignment is local to the database sequence. I don't understand OP's logic behind using local alignments on reads.

ADD REPLYlink written 16 months ago by RamRS24k
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: 2114 users visited in the last hour