Is it possible to do a variant call on a BLAT output?
1
0
Entering edit mode
3.8 years ago

Hey, i aligned reads (20 - 80 nt length) to the human genome. I won't get hits at the variable regions of the genome (eg. p-arm chr 13, 14, 15 ...), whereas the other parts are uniform covered. I am using a strict filter, to get exact matches. Is it possible to do a variant call or something like that to be more sensitive to this regions? Is it possible to map against the SNP-variants?

BLAT variant-call SNP snp • 761 views
ADD COMMENT
0
Entering edit mode

First, I assume your reads are Illumina? Is there a reason why you're BLATing such small reads and not use a read mapper (bwa/bowtie(2))? By limiting to exact matches, you will not get any variant position (since you're filtering them out). I'd use one of the read mappers I mentioned. Since these are typically used to find variants in reads, you should be able get something out of it. You could also leave out the exact match filter, but obviously then you need to take measures against false positives (ensure only trust alignments backed up by sufficient read coverage.)

ADD REPLY
0
Entering edit mode

Yes illumina reads. I can adjusted blat more sensitive than bwa (tile size, step size, ...). I also run a less strict filter and toke alignments containing 0 to 1 mismatch. Allowing one mismatch to a 20 nt read, means an identity of 95 % and i get more false positive results ... that's the reason why i am searching for a more sensitive way.

ADD REPLY
0
Entering edit mode
3.8 years ago
JC 13k

No directly, you will need to create some scripts to do that, but it's better to use standard tools and pipelines for variant calling.

ADD COMMENT

Login before adding your answer.

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