From a vcf I can do already, but that will only tell me what was called.... basically I am trying to get the "true" # in the human reference as compared to say NS12911 or NA12878. Is this possible? Thank you :). Edit: what I trying to do is sequence NA12878 and output the ins or del by length from the vcf. I can get that data but is there a way to get the actual or close to actual insertion/deletion by length of that known sequence?
count indel by length in human reference
Is there a tool or script that can calculate the ins or del length by chromosome in the ucsc.hg19.fasta (human geneome)? Thanks.
• 1,924 views
•
link
updated
by
Sean Davis
•
written
by
bioguy24
•
written
by
bioguy24
1 answer
Step 1: Call variants using a variant caller such as platypus, freebayes, samtools, or GATK.
Step 2: Write a small custom script to calculate the indel lengths per chromosome.
• 1 views
•
link
written
by
Sean Davis
Well, the "true" number is not actually known, but for GIAB sample, you could get the consensus calls and use that as an approximation of the "truth".
• 2 views
•
link
written
by
Sean Davis
Log in to answer this question.
More posts like this
-
csv to custom manifest file
written by bioguy24
-
How to determine mutation type/subtype via pysam/samtools
written by mike.rightmire •
-
targets in a bed at or above a coverage
written by bioguy24
-
Extraction of read subsequences based on a region from a sam file
written by gordon.freeman •
-
gene functional analysis
written by blooming.daisy333 •
-
FreeBayes Variant Diversity capture
written by shubhra.bhattacharya
-
list variants that map to a genomic location
written by bioguy24
-
calculate nucleotide forward and reverse reads from vcf
written by bioguy24
-
convert gene name to entrez id
written by bioguy24
-
add strand info to a bed file
written by bioguy24
what 's your input ? VCF ? BAM ?
It would be a BAM comparing NS12911 or NA12878 to the hg19 reference. Is this possible? Thanks.