This repository was archived by the owner on Jan 31, 2020. It is now read-only.
Description I'm interested in detecting internal tandem repeats in FLT3 with Pindel. Is there a way to determine (or even estimate) the variant allele frequency? I see an AD for the alt allele, but the ref portion of AD is always zero in my results, and the RD field is never used.
This is my VCF result for a 300bp TD from the Horizon HD829 reference.
chr13 28033909 . CTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCACTTTTCCAAAAGCACCTGATCCTAGTACCTTCCCTGCAAAGACAAATGGTGAGTACGTGCATTTTAAAGATTTTCCAATGGAAAAGAAATGCTGCAGAAACATTTGGCACATTCCATTCTTACCAAACTCTAAATTTTCTCTTGGAAACTCCCATTTGAGATCATATTCATATTCTCTGAAATCAACGTAGAAGTACTCATTATCTfGAGGAGCCGGTCACCTGTACCATCTGTAG CTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCACTTTTCCAAAAGCACCTGATCCTAGTACCTTCCCTGCAAAGACAAATGGTGAGTACGTGCATTTTAAAGATTTTCCAATGGAAAAGAAATGCTGCAGAAACATTTGGCACATTCCATTCTTACCAAACTCTAAATTTTCTCTTGGAAACTCCCATTTGAGATCATATTCATATTCTCTGAAATCAACGTAGAAGTACTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCACTTTTCCAAAAGCACCTGATCCTAGTACCTTCCCTGCAAAGACAAATGGTGAGTACGTGCATTTTAAAGATTTTCCAATGGAAAAGAAATGCTGCAGAAACATTTGGCACATTCCATTCTTACCAAACTCTAAATTTTCTCTTGGAAACTCCCATTTGAGATCATATTCATATTCTCTGAAATCAACGTAGAAGTACTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAG . PASS END=28034193;HOMLEN=0;SVLEN=284;SVTYPE=DUP:TANDEM;NTLEN=16 GT:AD 1/1:0,104
This is the relevant pindel output:
5 TD 284 NT 16 "TGCTTGAATCCGACTG" ChrID chr13 BP 28033909 28034194 BP_range 28033909 28034194 Supports 104 104 + 36 36 - 68 68 S1 2553 SUM_MS 3678 1 NumSupSamples 1 1 HD829 0 0 36 36 68 68
CCATTTGAGATCATATTCATATTCTCTGAAATCAACGTAGAAGTACTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAG tggattgagactcctgttttgctaattccataagctgttgcgttcatcacttttccaaaagcacctgatcctagtaccttccct
CCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCA
+ 28033951 0 HD829 @NB551050:207:HVWYGBGX9:2:12106:10347:17618:TTAAATTC
CCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCAC
- 28034252 70 HD829 @NB551050:207:HVWYGBGX9:2:12309:19545:19915:GGGCGGAG
CCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCAC
- 28034252 70 HD829 @NB551050:207:HVWYGBGX9:2:11303:1945:6705:GGGCGGAG
CCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCAC
- 28034251 70 HD829 @NB551050:207:HVWYGBGX9:1:13108:20254:4960:GGGCGGAG
ACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCA
+ 28033951 0 HD829 @NB551050:207:HVWYGBGX9:4:12611:14246:5494:TTAAATTC
ACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCA
+ 28033951 0 HD829 @NB551050:207:HVWYGBGX9:4:12502:13194:15006:TTAAATTC
ACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCA
+ 28033951 0 HD829 @NB551050:207:HVWYGBGX9:2:23104:6526:16936:TTAAATTC
ACCATCTGTAGGGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCA
+ 28033951 0 HD829 @NB551050:207:HVWYGBGX9:2:11305:26232:5971:TTAAATTC
ACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATCA
+ 28033951 0 HD829 @NB551050:207:HVWYGBGX9:1:13302:4351:14440:TTAAATTC
TACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATC
+ 28034043 0 HD829 @NB551050:207:HVWYGBGX9:3:12606:19024:6635:TTTATTAG
TACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATC
+ 28034043 0 HD829 @NB551050:207:HVWYGBGX9:2:21101:12036:3704:TTTATTAG
TACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATC
+ 28034043 0 HD829 @NB551050:207:HVWYGBGX9:1:21302:17014:13350:TTTATTAG
TACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCATC
+ 28034043 0 HD829 @NB551050:207:HVWYGBGX9:1:12108:14390:4730:TTTATTAG
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTC
- 28034439 70 HD829 @NB551050:207:HVWYGBGX9:4:11503:6235:16676:TCAATTTG
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCAT
- 28034088 70 HD829 @NB551050:207:HVWYGBGX9:3:22412:21934:14131:ATGCGCCC
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCAT
- 28034088 70 HD829 @NB551050:207:HVWYGBGX9:2:23210:10332:6756:ATGCGCCC
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCAT
- 28034087 70 HD829 @NB551050:207:HVWYGBGX9:2:22205:8591:13822:ATGCGCCC
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCAAAAGCTGTTGCGTTCAT
- 28034087 70 HD829 @NB551050:207:HVWYGBGX9:2:13308:7570:4608:ATGCGCCC
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCAT
- 28034441 70 HD829 @NB551050:207:HVWYGBGX9:1:23310:14372:13006:TCAATTTG
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCAT
- 28034087 70 HD829 @NB551050:207:HVWYGBGX9:1:23105:26150:10227:ATGCGCCC
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCAT
- 28034087 70 HD829 @NB551050:207:HVWYGBGX9:1:22212:10625:11351:ATGCGCCC
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCAT
- 28034087 70 HD829 @NB551050:207:HVWYGBGX9:1:22203:13586:19359:ATGCGCCC
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCAT
- 28034081 70 HD829 @NB551050:207:HVWYGBGX9:1:21307:3357:1342:TATTAATA
GTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCAT
- 28034087 70 HD829 @NB551050:207:HVWYGBGX9:1:21112:14535:7682:ATGCGCCC
TGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCA
+ 28034006 0 HD829 @NB551050:207:HVWYGBGX9:3:22408:11497:18158:GATGCTAA
TGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCA
+ 28034006 0 HD829 @NB551050:207:HVWYGBGX9:1:21312:5124:10200:GATGCTAA
TGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCA
- 28034118 70 HD829 @NB551050:207:HVWYGBGX9:3:22508:14456:12120:GGGTGAGG
TGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCA
- 28034118 70 HD829 @NB551050:207:HVWYGBGX9:3:11405:20064:5799:GGGTGAGG
TGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTCA
- 28034118 70 HD829 @NB551050:207:HVWYGBGX9:2:23203:2047:10968:GGGTGAGG
CTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTTC
- 28034178 70 HD829 @NB551050:207:HVWYGBGX9:3:23407:5988:11756:CTGTACCT
CCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTT
- 28034051 70 HD829 @NB551050:207:HVWYGBGX9:4:22512:18360:14598:CACACGCA
CCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTT
- 28034190 70 HD829 @NB551050:207:HVWYGBGX9:4:21608:17341:6781:TGATCTGC
CCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTT
- 28034051 70 HD829 @NB551050:207:HVWYGBGX9:4:13602:14377:6354:CACACGCA
CCTGTACCATCTGTAGTGCTTGAATCAGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGATGCGTT
- 28034247 70 HD829 @NB551050:207:HVWYGBGX9:4:12610:22280:10403:GTCTTCTA
CCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTT
- 28034051 70 HD829 @NB551050:207:HVWYGBGX9:3:13401:18674:3865:CACACGCA
CCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTG
- 28034247 24 HD829 @NB551050:207:HVWYGBGX9:2:22303:26516:10162:GAGGAGGT
CCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTG
- 28034244 14 HD829 @NB551050:207:HVWYGBGX9:2:22106:7501:14519:GAGGAGGT
CCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTT
- 28034395 70 HD829 @NB551050:207:HVWYGBGX9:2:21208:26053:12414:ACAAAGTA
CCTGTACCATCTGTAGTGCTTGAATCAGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGC
- 28034392 70 HD829 @NB551050:207:HVWYGBGX9:2:21201:10460:17099:ACAAAGTA
CCTGTACCATCTGTAGTGCTTGAATCAGACTGTGGATTGAGACTCCTGTTTTGATAATTCCATAAGCTGTTGCGTT
- 28034247 70 HD829 @NB551050:207:HVWYGBGX9:2:21105:20030:20197:GTCTTCTA
CCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTT
- 28034106 70 HD829 @NB551050:207:HVWYGBGX9:2:12307:9576:16463:ATTGTACC
CCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTT
- 28034139 70 HD829 @NB551050:207:HVWYGBGX9:2:12209:13951:20124:ACCATTTT
CCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCGTT
- 28034139 70 HD829 @NB551050:207:HVWYGBGX9:1:11307:1880:8527:ACCATTTT
CACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCG
+ 28034006 0 HD829 @NB551050:207:HVWYGBGX9:4:22511:7918:4406:TACGCAAT
CACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCG
+ 28034006 0 HD829 @NB551050:207:HVWYGBGX9:1:11111:2826:18615:TACGCAAT
CACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTGCG
+ 28034006 0 HD829 @NB551050:207:HVWYGBGX9:1:11107:10906:2941:TACGCAAT
GTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTG
+ 28033997 0 HD829 @NB551050:207:HVWYGBGX9:4:23401:14336:18003:TCGACCCC
GTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTG
+ 28033997 0 HD829 @NB551050:207:HVWYGBGX9:3:21403:13039:10860:TCGACCCC
GTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTG
+ 28033997 0 HD829 @NB551050:207:HVWYGBGX9:2:12308:9701:6442:TCGACCCC
GTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGCCTCCTGTTTTGCTAATTCCATAAGCTGTTG
- 28034168 70 HD829 @NB551050:207:HVWYGBGX9:3:22610:9239:6870:ACGTTTGG
GTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTG
- 28034168 70 HD829 @NB551050:207:HVWYGBGX9:1:22303:22003:12429:ACGTTTGG
GTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGTTG
- 28034168 70 HD829 @NB551050:207:HVWYGBGX9:1:22208:18967:12899:ACGTTTGG
CGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGT
- 28034002 70 HD829 @NB551050:207:HVWYGBGX9:4:13412:26560:16802:ATCTACAT
CGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGT
- 28034156 70 HD829 @NB551050:207:HVWYGBGX9:4:13409:5172:8987:CGTCCTCA
CGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGT
- 28034156 70 HD829 @NB551050:207:HVWYGBGX9:3:21502:14631:18475:CGTCCTCA
CGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGT
- 28034002 70 HD829 @NB551050:207:HVWYGBGX9:2:22102:6524:10394:ATCTACAT
CGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGT
- 28034002 70 HD829 @NB551050:207:HVWYGBGX9:2:13307:8316:2570:ATCTACAT
CGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGCTGT
- 28034002 70 HD829 @NB551050:207:HVWYGBGX9:1:22106:16579:18612:ATCTACAT
AGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGC
+ 28034031 0 HD829 @NB551050:207:HVWYGBGX9:2:23307:1615:7696:TGATTTTA
AGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGC
+ 28034031 0 HD829 @NB551050:207:HVWYGBGX9:2:22101:7039:2873:TGATTTTA
AGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGC
+ 28034031 0 HD829 @NB551050:207:HVWYGBGX9:2:12105:12815:19614:TGATTTTA
AGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAGC
- 28034110 70 HD829 @NB551050:207:HVWYGBGX9:4:23608:26285:17007:TCCTATAT
GAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAG
- 28034065 70 HD829 @NB551050:207:HVWYGBGX9:3:22605:19916:9995:TAACGCGC
GAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAG
- 28034065 70 HD829 @NB551050:207:HVWYGBGX9:2:13307:13714:13299:TAACGCGC
GAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCATAAG
- 28034065 70 HD829 @NB551050:207:HVWYGBGX9:1:11112:9549:17151:TAACGCGC
AGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCAT
+ 28033943 0 HD829 @NB551050:207:HVWYGBGX9:3:12501:15677:20386:CCCGCCGT
GAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGCCTCCTGTTTTGCTAATTCCAT
- 28034026 0 HD829 @NB551050:207:HVWYGBGX9:4:12612:1301:15164:GGAGGGGC
GAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCAT
- 28034026 0 HD829 @NB551050:207:HVWYGBGX9:3:22402:20735:3999:GGAGGGGC
GAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCAT
- 28034026 0 HD829 @NB551050:207:HVWYGBGX9:1:12301:6984:15578:GGAGGGGC
GAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCAT
+ 28033943 0 HD829 @NB551050:207:HVWYGBGX9:1:13107:4829:19376:CCCGCCGT
GAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCA
- 28034025 70 HD829 @NB551050:207:HVWYGBGX9:4:12611:26687:15177:GGAGGGGC
TGAGGAGCCGGTCACCTGTTCCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCA
+ 28034026 0 HD829 @NB551050:207:HVWYGBGX9:1:13201:22397:16718:GGACTAGC
TGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCA
- 28034034 0 HD829 @NB551050:207:HVWYGBGX9:3:23501:2694:6051:CTTATCAT
TGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCA
- 28034034 0 HD829 @NB551050:207:HVWYGBGX9:3:22503:14869:18783:CTTATCAT
TGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCA
- 28034034 0 HD829 @NB551050:207:HVWYGBGX9:3:22409:13338:3120:CTTATCAT
TGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCA
- 28034034 0 HD829 @NB551050:207:HVWYGBGX9:1:21302:23177:9270:CTTATCAT
TGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCA
+ 28033956 0 HD829 @NB551050:207:HVWYGBGX9:4:12612:1596:10321:TTTCACTT
TGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATTCCA
+ 28033956 0 HD829 @NB551050:207:HVWYGBGX9:3:21608:4123:7725:TTTCACTT
ATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
+ 28033943 70 HD829 @NB551050:207:HVWYGBGX9:3:21605:16607:1667:ATATTGTG
ATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATT
- 28034035 0 HD829 @NB551050:207:HVWYGBGX9:1:21206:13209:3018:TTCACTGG
ATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAATT
- 28034035 0 HD829 @NB551050:207:HVWYGBGX9:1:13212:22293:7979:TTCACTGG
TATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
+ 28033943 70 HD829 @NB551050:207:HVWYGBGX9:4:23406:14507:15113:ATATTGTG
TATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
+ 28033943 70 HD829 @NB551050:207:HVWYGBGX9:1:21312:7569:1792:ATATTGTG
TATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
- 28034055 0 HD829 @NB551050:207:HVWYGBGX9:4:23503:25434:14228:GGAACAGG
TATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
- 28034055 0 HD829 @NB551050:207:HVWYGBGX9:4:12402:23801:5936:GGAACAGG
TATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
- 28034055 0 HD829 @NB551050:207:HVWYGBGX9:1:22307:18490:2250:GGAACAGG
TATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
- 28034055 0 HD829 @NB551050:207:HVWYGBGX9:1:11204:5301:16962:GGAACAGG
TATCTGAGGAGCCGGTCACCTGGACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
+ 28033943 0 HD829 @NB551050:207:HVWYGBGX9:3:23502:18790:3055:ATATTGTG
TATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
- 28034011 0 HD829 @NB551050:207:HVWYGBGX9:4:11507:17931:9721:GCTTCTTT
TATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
- 28034011 0 HD829 @NB551050:207:HVWYGBGX9:3:11510:21491:13038:GCTTCTTT
TATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAAT
- 28034011 0 HD829 @NB551050:207:HVWYGBGX9:2:11107:10820:17669:GCTTCTTT
TTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCAGACTGTGGATTGAGACTCCTGTTTTGCTAA
- 28034151 0 HD829 @NB551050:207:HVWYGBGX9:4:22402:20785:10022:TCCTACTA
TTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAA
- 28034149 0 HD829 @NB551050:207:HVWYGBGX9:4:23412:21271:19233:TCCTACTA
TTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAA
- 28034148 0 HD829 @NB551050:207:HVWYGBGX9:3:13603:10611:8918:TCCTACTA
TTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAA
- 28034148 0 HD829 @NB551050:207:HVWYGBGX9:2:23209:17506:13728:TCCTACTA
TTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTGCTAA
- 28034148 0 HD829 @NB551050:207:HVWYGBGX9:1:23111:5615:15916:TCCTACTA
CTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTG
+ 28034050 70 HD829 @NB551050:207:HVWYGBGX9:4:23510:7714:16346:AGTAGTTT
CTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTG
+ 28034050 70 HD829 @NB551050:207:HVWYGBGX9:4:21411:18012:10549:AGTAGTTT
CTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTG
+ 28034037 70 HD829 @NB551050:207:HVWYGBGX9:3:23403:17619:11920:CCGACAAG
CTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTTTG
+ 28034050 70 HD829 @NB551050:207:HVWYGBGX9:2:21305:11394:16134:AGTAGTTT
TACTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCTGTTT
- 28034149 0 HD829 @NB551050:207:HVWYGBGX9:1:11305:22985:12447:AACGAATT
GAAGTACTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGACTCCT
+ 28033917 70 HD829 @NB551050:207:HVWYGBGX9:4:13610:11781:2495:GCCGAAAT
CGTAGAAGTACTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGAC
+ 28033748 70 HD829 @NB551050:207:HVWYGBGX9:2:22306:19889:11897:CGCGATCG
CGTAGAAGTACTCATTATCTGAGGAGCCGGTCACCTGTACCATCTGTAGTGCTTGAATCCGACTGTGGATTGAGAC
- 28034123 0 HD829 @NB551050:207:HVWYGBGX9:1:23110:25361:3815:ACTTTTCC
Reactions are currently unavailable
I'm interested in detecting internal tandem repeats in FLT3 with Pindel. Is there a way to determine (or even estimate) the variant allele frequency? I see an AD for the alt allele, but the ref portion of AD is always zero in my results, and the RD field is never used.
This is my VCF result for a 300bp TD from the Horizon HD829 reference.
This is the relevant pindel output: