1
$\begingroup$

I am new to RNA sequencing analysis. I have RNA samples for bacteria for different treatment. I did contaminant filtration, remove adapters and checked fastQC report: attached below for one sample (fwd and rev) fwd. reverse After mapping with the reference genome using bbsplit but the mapping is only 0.62%.

Any help/suggestions highly appreciated.

Thanks

$\endgroup$
1
  • 2
    $\begingroup$ Well, when you looked at the sequences of some of your reads, what species did they BLAST to? $\endgroup$
    – swbarnes2
    Mar 16, 2022 at 1:34

1 Answer 1

1
$\begingroup$

In addition to swbarnes2 comment, I'd suggest a more systematic approach.

Kraken2 is a tool for metagenomic analysis of reads, though it can be hard to set up.

mash screen is another method for metagenomic analysis, much easier to set up but noisier.

If you are dealing with contamination, treating your data as metagenomic might help.

Also, as always, ensure that your workflow and reference assembly are correct!!!

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.