Missing phage genome in SRA metagenomic data
2
2
Entering edit mode
8 weeks ago

Hi all, I have 3 metagenomic assembled phage genomes from NCBI with their respective SRA reads downloaded. I just want to confirm the presence of these viruses in these SRA files. I was able to do that to 2 of them by just building the scaffolds with Spades (metaviral mode). I'm having trouble finding the third virus in its respective SRA data, I've tried: Spades, Bowtie (mapping the reads to the assembled genome) and VirSorter with no results. I do not know if I could do something else or there is some kind of error in the information provided by the NCBI. Any clues?

metagenome Spades SRA phage Bowtie • 325 views
ADD COMMENT
2
Entering edit mode
8 weeks ago

Sounds like a pretty thorough investigation to me. If there are no reads mapping to the assembled genome that was supposedly created from those reads, I would consider this highly suspicious and suggestive of an error in the database (e.g. the wrong file was submitted to SRA by mistake).

The only alternative explanations that come to my mind: You happen to be interested to find back Phi X 174 or some mapping issues that can be fixed with different settings (reads need to be quality-trimmed, RNA/DNA conversion required?). But at this point, I would also start considering a database error, in particular, since you could identify the other two phages.

ADD COMMENT
0
Entering edit mode

It seems like a database error indeed

ADD REPLY
1
Entering edit mode

Consider emailing SRA support and the submitter's with relevant information. This may get fixed eventually.

ADD REPLY
2
Entering edit mode
8 weeks ago
GenoMax 125k

You could do a blast search against SRA (you can only do this via NCBI web blast page) with that last phage and see what you get.

ADD COMMENT
0
Entering edit mode

Just tried that, had positive hits but with percentage identity below 100. No hits when using as query the phage gene that I'm investigating. Positive 100% percentage identity hits for the same gene from the other two phages in their SRA runs though

ADD REPLY

Login before adding your answer.

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