Question: How can I know that the blast hits of my NGS data mapping against my query sequence?
0
gravatar for apl00028
14 months ago by
apl0002870
apl0002870 wrote:

I did a blastn of my NGS sequences against a NCBI database, following the next steps:

firstly make a database: makeblastdb -in hk_db.fasta -title "HK" -dbtype nucl -out hk_db -parse_seqids

then, I did the blast analysis: blastn -query PV001_1_trim.fa -db hk_db -evalue 1e-4 -outfmt "6 qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore qlen qcovs" -out PV001_1_blast.out

I got the results and I analyzed by R. I got the number of hits that match among my query sequence and my NGS data, I would like to know if there is some way to know if these hits mapping along the query sequence or not.

Thanks in advance.

alignment ngs • 386 views
ADD COMMENTlink written 14 months ago by apl0002870

May I ask why you use blast instead of a dedicated NGS aligner? Maybe there is a better way of achieving your analysis goal.

ADD REPLYlink written 14 months ago by ATpoint46k

Yeah, I though that It could be useful make an aligment with Bowtie2. But I do not what paramether shows that my hits mapping in an equitative way with my query sequence using Bowtie2

ADD REPLYlink written 14 months ago by apl0002870

Can you give a representative example what exactly you want to see. Blast is going to take a very long time for an entire NGS run.

ADD REPLYlink written 14 months ago by ATpoint46k

I would like to know that the positive hits of my NGS data mapping in more than one position, I mean that it covers a long position of my query sequence

ADD REPLYlink written 14 months ago by apl0002870
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1973 users visited in the last hour
_