I am willing to make a bunch of blast queries in ncbi hosted SRA datasets.
The blast query is fast, I want to search for two short sequences X and Y in SARS-CoV-2 datasets, Y is a mutated version of X, and stop when 100 reads are found, in order to find if sometimes Y appears as a minority variant of X, and conversely.
It seems this is not allowed (anymore?), neither by the blast Web API nor by the blastn command line executable.
So this fails blastn -db SRX10000025 -query XY.fasta -remote -out blast_output.out
as well as github/srablast and blast/Web API.
On the other hand the blast interface allows to do it for one SRA.
Is there any way to get over this?
I made a demo there showing that ncbi is blocking SRA requests outside of the Web interface:
http://cov2.infinityfreeapp.com/blast.html nt request, works
http://cov2.infinityfreeapp.com/blast_sra.html SRA request, is blocked.
Sometimes the SRA request works too, but more frequently it fails with