This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Alignment of Query-RNA Seq

I am trying to align query sequence with RNA-Seq data (100bp PE). I have tried Bowtie2 and BWA. However, it did not give me reads that are matching or mismatched with query sequence. (I am not interested in aligning RNA-Seq data with Reference genome).

rna-seq alignment
• 2,083 views • link
written by skevalkumar • skevalkumar

(I am not interested in aligning RNA-Seq data with Reference genome).

That is fine but why are you not using the query sequence as your reference to make an index? Then align your data against it.

Note: Doing this always runs into risk of having some reads align in locations that they did not originate from.

• 0 views • link
written by GenoMax GenoMax

I did indexing of my query sequence and tried (--end-to-end and --local parameters) to align with RNA-Seq data. However, overall alignment rate was 0.00%.

• 0 views • link
written by skevalkumar • skevalkumar

What is the length of the query sequence?

• 0 views • link
written by GenoMax GenoMax

I suggest that you try bbmap.sh from BBTools (https://sourceforge.net/projects/bbmap/ ). Something like this:

bbmap.sh -Xmx10g threads=4 in1=R1.fq.gz in2=R2.fq.gz out=file.bam maxindel=2000 intronlen=10 ambig=random ref=your_query.fa mappedonly=t

This will only write mapped reads to the bam file (this will require samtools to be in your path otherwise SAM format will be used). If you only want to see how many reads are aligning then omit out. All the stats will still be written to STDERR.

Was there no alignment even if you used command defaults for bwa and bowtie2?

• 0 views • link
written by GenoMax GenoMax

There was no alignment with bowtie2 (all default parameters with endtoend and local) and bwa mem.

• 0 views • link
written by skevalkumar • skevalkumar

Hi, I have tried bbmap.sh and it shows 0.0011% mapped reads. I am reading bbmap reference guide (it is new to me.) thank you for help. Please give me more suggestions.

• 0 views • link
written by skevalkumar • skevalkumar

Not something you want to hear but here goes :

  • Your data only has a tiny fraction of reads that map to the query you have
  • You could try mapping individual reads (R1 or R2) as single-end data to the query (remove in1= and in2=, just use in= with R1 and R2 files above command line) and see if the mapping improves (some other explanation can be explored there, if you see a good number of reads aligning)

Take 10 reads (from R1 and R2) and blast them at NCBI to confirm that you are looking at the correct sample set and the reads are aligning to right genome. Would eliminate the possibility that you have contamination of some kind.

• 0 views • link
written by GenoMax GenoMax

It is showing almost similar result (0.0008%). I have checked few reads in NCBI. It is matching with my plant species.

• 0 views • link
written by skevalkumar • skevalkumar

The the only plausible explanation is this:

Your data only has a tiny fraction of reads that map to the query you have
• 0 views • link
written by GenoMax GenoMax

0 answers

No answers yet.

Log in to answer this question.

More posts like this