genome / pindel

Pindel can detect breakpoints of large deletions, medium sized insertions, inversions, tandem duplications and other structural variants at single-based resolution from next-gen sequence data. It uses a pattern growth approach to identify the breakpoints of these variants from paired-end short reads.
GNU General Public License v3.0
162 stars 89 forks source link

how to get allele frequencies with Pindel? #117

Open myourshaw opened 4 years ago

myourshaw commented 4 years ago

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