Question: Does Last aligner only outputs the mapped reads to a bam file ?
0
gravatar for pinninti1991reddy
2.0 years ago by
pinninti1991reddy60 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 2.0 years ago by RamRS28k • written 2.0 years ago by pinninti1991reddy60
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 2.0 years ago • written 2.0 years ago by RamRS28k

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 2.0 years ago • written 2.0 years ago by h.mon30k

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

LAST ~= (BLAST - B)

ADD REPLYlink modified 2.0 years ago • written 2.0 years ago by Joe17k
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 2.0 years ago by RamRS28k
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: 1595 users visited in the last hour