Targeted DNA seq analysis
0
0
Entering edit mode
8.8 years ago
kanwarjag ★ 1.2k

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 • 1.6k views
ADD COMMENT
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 REPLY
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 REPLY
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 REPLY

Login before adding your answer.

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