BLAST and GFF script to find locations of SNPs
1
0
Entering edit mode
17 months ago
Kiki_01 • 0

Hi! I found SNPs in transcript using GATK and I need to find their locations in genes. I did the BLAST trascripts to the genome. Does anyone know how should python script looks like if I'm trying to connect together my BLAST results with GFF file to find in which genes SNPs are located? I'd appreciate any help!

BLAST gff script SNP • 319 views
ADD COMMENT
0
Entering edit mode
17 months ago
colindaven ★ 3.2k

The only way to go is to start again - realign your reads to the genome and call the SNVs using GATK there.

Transcriptome != genome (and this is really important when mapping short reads).

BLAST, GFFs etc is all a wild goosechase with no hope of success in my opinion.

ADD COMMENT
0
Entering edit mode

Dear colindaven, thank you for your respond, but that's what i have to do according to my thesis. I have found SNPs in transcript and i need to find their locations in genome using script with BLAST results and GFF file together.

ADD REPLY

Login before adding your answer.

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