Biostar Beta. Not for public use.
Targeted DNA seq analysis
0
Entering edit mode
17 months ago
kanwarjag • 960
United States

I have a DNA-seq data and have to count incidences of specific 10bp occurrence, But it is an invitro data and has no specific genome. I understand perhaps I can align two sequences, My question is since I dont have a genome how shoudl I handle these FAsTq files to take the fragments to aligningl blast.

Thanks. Any direction or pointer will be helpful

Targeted DNA-seq • 818 views
ADD COMMENTlink
0
Entering edit mode

If you only need to count the occurrence and doesn't allow for mismatch e.t.c (require exact match), then you can simply count the occurrence of the sub-string using your favorite programming language

ADD REPLYlink
0
Entering edit mode

It may have one or two bases of variations and I may have to use some kind of variation threshold.

ADD REPLYlink
0
Entering edit mode

I guess you can try converting the FASTQ into FASTA and use them as reference for BLAST, then align the FASTQ to the reference with BLAST allowing for multiple mapping and some degree of mutation. However, I am not particularly sure about that

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1