Bowtie-Build, Which File To Use Contigs.Fa Or Transcript.Fa
1
3
Entering edit mode
11.2 years ago
nbvasani ▴ 240

Hi Community users,

I have single-end reads fastq file. I have generated denovo assembly using velvet/oases. Now I need to check how good is my assembly. So I am using aligner, bowtie. But I am confused which fasta file to use for creating bowtie index i.e. shall I used contigs.fa from velveg or transcript.fa file from oases.

Do you guys prefer --best parameter option while running bowtie?

I would really appreciate your feedback and reponse.

Thanks,
Naresh

bowtie RNA-seq • 5.1k views
ADD COMMENT
3
Entering edit mode
11.2 years ago

Use both and compare the results.

Plus use bwa-mem instead of bowtie, that way your reads will align across junctions even for the contigs.fa file.

ADD COMMENT
0
Entering edit mode

Thanks Albert. I will try that.

ADD REPLY
0
Entering edit mode

Hi Albert, I tried bowtie with contigs.fa file, but I got only 6% alignment. Now I am trying to run with transcripts.fa. Any suggestion why alignment rate is so low. Fastq file used for alignment was trimmed to get better quality.

Thanks in advance. Naresh

ADD REPLY
0
Entering edit mode

that's why I said use bwa-mem

ADD REPLY
0
Entering edit mode

ohh ok. Thanks.

ADD REPLY
1
Entering edit mode

Hi Albert,

As per your advice I used bwa-mem for alignment it work out very well. It think I got 79% align. Can you please check my flagstat output. If I am interpreting result in correct manner.

If it is 79% align can you please help me how can I improve alignment above 90% by adding more parameter.

I used following cmd:

##bwa index contigs.fa
##bwa mem contigs.fa CombineIonXpressRNA_009_NareshPool_Chip1_2_WT1_fastx_trimmer_from_quality_trimmer_file.fastq > aln-se.sam
##samtools view -bS aln-se.sam > aligned_reads.bam
##samtools flagstat aligned_reads.bam
34569162 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 duplicates
27555361 + 0 mapped (79.71%:-nan%)
0 + 0 paired in sequencing
0 + 0 read1
0 + 0 read2
0 + 0 properly paired (-nan%:-nan%)
0 + 0 with itself and mate mapped
0 + 0 singletons (-nan%:-nan%)
0 + 0 with mate mapped to a different chr
0 + 0 with mate mapped to a different chr (mapQ>=5)

Appreciate your help.

Thanks,
Naresh

ADD REPLY
0
Entering edit mode

there is no reason to expect a 90% alignment rate, some or many of your reads may not have ended up in the assembly

ADD REPLY
0
Entering edit mode

Thanks Albert! That mean I understand correctly, 79% was aligned from flagstat output?

ADD REPLY

Login before adding your answer.

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