Hello, I need help in interpreting abricate (vfdb) result. Based on the result below can I assume that there is hilA gene in the analyzed strain? I mean if I can sum coverage results for these fragments and assume that I have 100% coverage for hilA gene?
How would you interpret these results?
Thank you in advance.
•
written
by
elina2410 •
1 answer
Yes, I think you are right. Looks like the two contigs have an edge “inside” of the hilA gene.
To check this guess, I probably would join the two contigs, then annotate them (e.g. using prokka), then check if prokka regards the junction as an hilA gene. Even if junction contains an unrecognized ORF, you can blast its sequence and check if it is similar to hilA.
To join the contigs correctly, you can investigate their “sticky ends” using the combinator-FQ tool.
Log in to answer this question.
-
Strand-specific quantification of RNA-seq data without strand-aware mapping
written by tdsone •
-
KaryoPlot: Gene coverage over centromeric regions
written by Mary •
-
How to interpret the Variant Allele Frequency??
written by hashim.rana11 •
-
Converting Abricate output (.tsv) to gff3 format
written by ghataksnd •
-
Top 3 gene clusters using MCODE in Cytoscape - Protein-protein interaction analysis
written by bioinformatics •
-
Combine gene annotations (from metagenomes) & gene coverage information into single file
written by lintonf •
-
Please help me with the MultiQC results
written by arr234 •
-
Interpreting plots from Cytoscape
written by xmsterl •
-
How to set different color in every gene and add legend using Chromomap?
written by amirandi1808 •
-
OMA gene vs splicing variants in the output
written by mateusz.konczal •
