基因数据处理50之cs-bwamem、bwa、snap、bwa-mem与art比较

(3) snap:

hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub/cs-bwamem$ cat G38L100c50Nhs20.snap.sam 
@HD VN:1.4  SO:unsorted
@RG ID:FASTQ    PL:Illumina PU:pu   LB:lb   SM:sm
@PG ID:SNAP PN:SNAP CL:single index G38L100c50Nhs20.fq -o G38L100c50Nhs20.snap.sam  VN:1.0beta.23
@SQ SN:chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38   LN:248956422
chr1-50 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  93465785    26  100M    *   0   0   TTCCACAATAGTTGAACTAATTTACAGTCCCACCAACAGTGTAAAAGTGTTCCTATTTCTCCACATCCTCTCCAGCACCTGTTGTTTCCTGACTTTTTAA    @@CDFDFDHFHGHIJH:IJJJ(JJE?JDIDEJIB@FGJIGBHJ()HG8(CIICGFFHEH=GI3@&@DD58FADDACHDDHFCD8D,DCC<CEFD<EDDCD    PG:Z:SNAP   NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-48 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  228133746   70  100M    *   0   0   ATCATTGTATGCCACAGAAATAATTAAATTTCCTTGTCAACTGACACATTATTATTAGGCACTCTCACCAGATCTTTACCCATGGCCATTTAAAGTGTGG    @>CFFFFFH<<GC1IIDCFJJHIGIHJ(IID7IJ,FJJJHJJJJ)GGBHIJFJFIFIHFE=HEIEE;CA)G0(D()HC@D(:EFDDC@;DDAC95(D?BD    PG:Z:SNAP   NM:i:1  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-47 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  13772988    70  100M    *   0   0   TTCAGTAATTCAGAATAACACATGAGGGAATGAATGAATGAATAAATAAAAAAAAACTGAATGAATAAATTACAAAAAATTGTGTTTCAGGGAAGAAAAA    CC@F(FFFDFH.HDHIGI(JIIIGGIEEJIIIHJJHHH3IJJIIJ3=EI>JDIGH((IBJCIEHGD>;J@HF+DC)CCCADBDBD+BDDDD5B5DDDE(C    PG:Z:SNAP   NM:i:2  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-46 16  chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  37474758    70  100M    *   0   0   GGGTCGGGGTCCTGTTCCCCGGTCCGCCGGGCCTCAGGACCCCTCCAACTTTGCCCAAGTTGGGAGAGCCGGGGAAGAGCACCAGGTTCCTGATCGGGAT    (5CBACDDD>FBDDDDDEC:CE(CBDFDDHEFH;FGEFHGHDGJJJDIGI:JEHJ=JJJJJH8CI?JJJG9JIII>IJIIGJ=EIJGAHHHHFFDFDCC?    PG:Z:SNAP   NM:i:1  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-45 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  29056657    70  100M    *   0   0   CTGGGATTACAGGTGCCCGCCACCATGCCCAGCTAATTTTTGTATTTTTGGTAGAGACAAGGTTTCACCATGTTGGCCGGGATTGTCTCGAACTCCTGAT    B@@FFFFFHHG)HIJJJJBJIJCJHGJIBFJJI3IIHDF@JIAJ9JJJIJJBIJJ?BJID8F:HFHA(+D>J>CG>7D=DDFF@EDC3D<CDDDC@BD@B    PG:Z:SNAP   NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-44 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  49993893    70  100M    *   0   0   CAATTTAGCCAAAACTGGCTAATCGTTTTACCAGAATCATTCCCATTGTTCAAGACCTATTTTAAGCTCCACTATCACCATAAAACTTTCCCGATCAGTT    C@CFFFFDHHHHH<IEJJ@JJJHI)IDIBIJA:HJHFJJJIJGGJJIIIIHGGJJGH<(IIIJI?ICDG;CDHFHCDCCB?FDDED:CD:>DD5C&DDCD    PG:Z:SNAP   NM:i:1  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-43 16  chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  194714506   70  100M    *   0   0   AATATGTTTTAATAATATCATATTTAAATTTGATGATACTTTAAAAATGGTTCCATGTGTGTTCTCTTGGGTTATTTCACAATCAATAAAAGGTCTGCAA    CCCCDDC@E>CDCDDC>D9CD=C)CGC>E@7.HF)DIBJBJJ.JEJEJ@JJIIIIGD?<IHH)FJJIIH*DJIBIIJHJIHHJFHIGHHHHDFDFDFC@B    PG:Z:SNAP   NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-42 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  35706203    70  100M    *   0   0   CAGGTTCAAGCGATTCTCCTGCCTCAGCCTCCTGAGTAGCTGGGATTACAGGCACGTGCCACCATGCCTGGCAATTTTTGTATTTTTAGTACAGATGGGG    CC@FDFFAHHGFHJIHJFJJII=@JEHIJIIJIJEJIJJHHGIJBBFJG6JJHJJG<F3JJHIFG(DCJDHDFDCHDDF7DDHBDFDDCDDCDCD;CBCB    PG:Z:SNAP   NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-41 16  chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  156482338   70  100M    *   0   0   GTGTGTGCATAGGCAGGTCTGCGTGTACATGCAACGTGGGCACGTGTCCATGTGGATGCAGGCGGGGGTATATCCTGGTGCCTGTGTGTATGGGCCCACC    D;CCDDCDCDDDCD:EDA@<C<E(GDDDGDJDDHJ@CJJI,=FHJJIGJ7GEC?IGJJIFBBICHJEIJJHHAIJIJI.IJGJJGJJHHGHHFFFFFB=B    PG:Z:SNAP   NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-40 16  chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  221779284   70  100M    *   0   0   CATGGCACATAGCACTTTGGTGATGGGGACTGCTTTGCTAATGTCAGGGTCAAGGGGTGCATGGACCATGGGCAGAGTGCTGGGCTCAGCCAAATGGTTC    DDBCDCDDDDDDD25F?DD@4I5HED?CAHGA?JJIIJB)IHFJJFCJII?@<HIIFGIJIJFG?JIJCIIJJ)IJJIJIJJIIEGJFHHHFFDFDFC@@    PG:Z:SNAP   NM:i:1  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-39 16  chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  3895605 70  100M    *   0   0   GTCCTCTCCGGATTGACAGGAGTCAAAACATGAGATCGGCTTAGCTTCAGTTTCGTCATGGATTAACCACCTCCAAGGTGTCAACTCCAAAATGTCAAGA    DD5CCAD&8DAD>D&FDDDCDBDD?6DD.FHDDIFE?@IDEGIBCGD?JFJ>JGBI,IJIF.JJIHJJJEIEGFJ=JJHJHHJFHIJHHHHHFFAFFC@B  PG:Z:SNAP NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-38 16  chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  33174926    70  100M    *   0   0   CACACATACATATATGTGTGTATATATATATATATATATATATACACACATATACATATATATGCACACACACATGTATGTATATGTATATGTATATGTG    CDDC(FBDC(AACBDDCBDDECEC5@H;HFDJFH>=FCHAHJFJ'H3JG9JFEHIJFDJJ9IJHEJIGJIJJJJJC;J?AJFJGEHFHDC<HDFFFDCC# PG:Z:SNAP NM:i:1 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-37 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 206124777 70 100M * 0 0 TCAGTCAATAGACATTTGGGTTGTTTCCACCTTGGGCAGGTTACAAATAATGCTGCTAGTGAACATTCATGTGCAAGTTTTTGTGTGGACATACGTTGTT CBCD8DDADC@DDDDB?CCHD;@AHECEEJIHAEII?E05GFJHHDJCJEJDBHJE7GJJJGJGGGJ=JIC(JIJHHIIIAGGJIIICHGHGFFFDF?CC PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-36 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 181673626 70 100M * 0 0 TCCACTGCCCAGAAAGAGGACATCCCTTATAGGACCAGCGGATGGAAGCCATGGGCTGGGCAGGACATTCCTGTCCCAACCCACATGGCAGCTAGAGTCC @@3DFFFFHHHFHJJIGGJHHFJJJIIJJJJDH*GJJJGJIJ6AIJIJDFDII=HFI2AH1AIEAAC?JEIEDJF.HDH@FAFDCDE2D:DDBD0DDBDD PG:Z:SNAP NM:i:1 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-35 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 152104780 70 100M * 0 0 ACCTCTATAAATACAATATCTTCAAATATGATTACATTCTGAGGTACTGAGGGTTAAGACTTCAACACGTGAACCTCTGTGGGGGTTGGGAGGTCACAAC +@?D>>D?DDBDFB)DDDDDC5(9>F;G)FB84/AJE3JJIJIGIGJBBIGCJCJGJGHJIDJ>IB7IGJGEGCIIGFJJJEFHIIJHHF=HFF8=F??= PG:Z:SNAP NM:i:1 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-34 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 13267477 3 100M * 0 0 ATACAGGGAAGGTTTTAATCTGTTTCAGACATTAGAAATACATATATTTATATATGGTATCTTTATTGGAGAACCTTTGGCCACATCAAAAGTATCAAAA @CCFFFFFAHHHHDEIJJJFI)<HBAJJIIGJFFGJI0IJJBIHFJ8JIIC=IDHJFBI8I5'3DDI.DDH@EDEH>EHED>F>CD(DDDECDED?BC@D    PG:Z:SNAP   NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-33 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  48968233    70  100M    *   0   0   TAATAGTAGGCAATAAACAAAGAGAGCAACTTAGGAGCCAGATCACATGTGGCCGCTCGAGCAATATGGTAAAAGTTCTGGACTTCATTCTAGGTGAATG    1CCB=FFFHHHHHEDHJIIAFG4JIFJIJB)JJI?(&JJIJJEE)HIJJBJ?HJ(=B(I@?I?8DC8C>JHJH>@EDFDD5DDDDDDDDCFD:=DCC(DD    PG:Z:SNAP   NM:i:2  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-32 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  88980623    70  100M    *   0   0   TAGTTCAGTAAACTATTTATCAAACAGGTGTCAGGTCATTTTAACATACTCCTTGCTTTGAACAATATTCATTCATACTTGGTACAAACTCTATATCCTA    B?CFDFFFHHFH3JIJJJIGJJJJJFJDEJGJ(EHFI>E=JIJ(GGJDFCH>>GJ=IHDJEHHDI>GEBJE@DD@HH'AA@ECC@BDEDDD@CDDADBDD PG:Z:SNAP NM:i:2 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-31 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 227005594 70 100M * 0 0 TCACCAGGCATCTTTACTGACTCACACCAATAGTAGTACTGGGATTAGAAATAAGACGCTGCAATACTCACAACCTAGGTGAAGTTAGTTAATTTGGGAA D@D5B=DACDDDDDBEFECBFDC5BCDDDCDFIDC8ICEIJ=DHIGHIJIJJB0HJJCDJHJGJIJI9GGHGGJ3@IJJAIJGGBGJ7HFHHDEBFF@CC PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-11 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 72085324 70 100M * 0 0 TACTAGCCTTGAAAATGTTTAAAATAATATTCCAGAGTTAATATTGTTGTCCCTGGTATGTTAAAGAGTATTTGTTATCATAGCCAATTCTTGAGTCTGC 8@DDCD4D>D?C3DF(DCCHDDDA;HDEIBFCHGHHHFFIFEG1JHIJIJCGEJIHJG)IH(IJ)BDJ??FHHJHCJJIFJHJJJGIGH)2HFFFA=<CC PG:Z:SNAP NM:i:2 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-10 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 214311330 70 100M * 0 0 GGTGATCCCTATTTGGTCCACTTTTGTTGGTAGTCTTCAAGCTTGATATCTGATTATCACTGTTGGAAGGTGTAAACTCACAGACTCAGAATTCTGGACT D&++DDEDDC3+CB(@D8DDEHD;?FDF?DBI9A7@JFHJ(I(AIJJ@DHJIIJHII>HFG1JJFII<IJJJ(<HIJGHIJJE=FH@FHHHHFDDDFC@C PG:Z:SNAP NM:i:1 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-9 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 152012629 70 100M * 0 0 AAGCAATTCTCTTGCTTTAGCCTCCCGAGAAGCTCGGATTACAGGCATGTCCACCACACCCAGCTAATTCTTTTGTATTTTTAGTAGACATGGGGTTTTG BC:FFDFFHHHFHJ9HIJHIJJJFIGJCI=H/IHH@IGJIJIJIGJJJJEEIB'JJDJJIJDIGICHEFD@D3:0A/(BECDDDDCBE>BD8DDDDDC8C    PG:Z:SNAP   NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-8  0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  79478960    70  100M    *   0   0   GGCAACACTTGAGAACACAAAGTGAGTTCTCACTTTGGGCGGTGGTTTCAGGCTTCAGGGTGGAGTTTTGTCAGGAACCCAACCTTTTCTGCCTAGAATT    @CCFDFDADHHAHFIH@ICJJHI5?JIJ)GCFEIJHG=II)HIGI9JJIJGHEJHFI8EIDG)GCI4FJF?I8HCDH;DD0&3CFDDDD@C4DCD6ADD>    PG:Z:SNAP   NM:i:1  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-7  0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  178190761   70  100M    *   0   0   GTAGCCGGAATAAACAGTCACTGTGAGTTGTCCATTTTAGAGCATAGGTTTTCAGGTGGTGAAGACCTGTCCTTAGTTGAATTTGTATGTGAATTAAACT    B?<;FDDFHHHHHJEGJFCJJIFJJJA=JHHGIIGJJIIGIGJ(D:DAFG7)&DJID9J)FCD/HHJEDFIJ<FJAF@D@DDADF?C@A@ADCDD@CDDD    PG:Z:SNAP   NM:i:1  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-6  0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  42572411    70  100M    *   0   0   AACCCTTTATCAGGTATGTATTATAAACATCGACTCTGTGGCTTGCATTTTCATTCTCCTTATATATCTTTTGATGAATCAAAGTTTTTAATTTGAATAT    BCCFFFFDAH)HHJJHIGG,HFH2JIJJ4IDI93IJJ<=JJ>IH7IJIJBIBG)CFH7DHHFAHFHEDIFBEFH;EBICA?3DD5D(DDBACDC(BADD:    PG:Z:SNAP   NM:i:1  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-5  0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  153186635   70  100M    *   0   0   GTCTTGACTCTTTATCCACTTTGCCAGTCTGTGTCTTGTAATTGGGGCATTTAGCCTATTTACATTTAAGGTTAATATTGTTATGTGTGAATTTGATCCT    C1@?FF=FHHHHHI?JJEJFIIIHG:.?>EEJEI(JG9J'IIHIJIHIJGJGJFJ9FJAG4EEC:DADE8DAEJFCCBBCDAEDDDD-DDDD@+DBC8D+ PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-4 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 145038405 3 100M * 0 0 AGTGGAAATAATACTCGTCAACATATGCCTTTCAAAAAAATTTTTTTTCATATTTTAAATTTACCTTTACTACCTATTTATTTGGTTCAAGGCTCCATTT C:CFFFDDHFHDHJIIJJJJ29CJ+JJJIIJIIFIG?JI08?CJJIFIFDEFDGBD>JAIDJDJ>JCBG(CG=DE5?(EDB3HDD>ED2:CCHDB<DDCC PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-3 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 84685021 70 100M * 0 0 TATTATTAAAACTATAAATGGACCAATTAAACAAACGTGTCATGAGCCAAGGAATATAAACTAATTCTTTACACCTGAAGTCCTTTAAAATGATTTAATT CCCDFFB=HHHFHJIIJ:E@J>JA2C<IEI2DGJHJGI8FJJIJJAHJCIJJJJ*JJ;F?HFDJCIJFDJ'CHECDEFD,DDFBCH<A7DCD-DDBDD(E    PG:Z:SNAP   NM:i:2  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-2  0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  62477842    70  100M    *   0   0   TAGGAAAATGGAGAAACTTTAATATGAAATCTTCCTGTTTTTCACATTATGTTTAGATTGTTACAGCATAAAATTTCAGAAACATTGCAAAAAGTTTTAA    @C=FFDFDHHHH>GJ@IEJGJIIJJJJF@JHHIIGGJ<IHJIG/J*?GD>ICIAJFJIH)H7E?GHEDI>HHFAHC@)(D>DDDED<DCDD=DBBC5DDE    PG:Z:SNAP   NM:i:2  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-1  0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  11355150    70  100M    *   0   0   ATTTATTGGCTGTCTTTCAGGCACATTTTAGCTGTCATCCAACATTCTCAACCTTAGTCCCCTTCTCTGGGCTAAGGGGAGAATGATGGTCCTACCCCAG    BC?DFFFFHH<FFJG(ICJIGJJIJJGIFJIJGJ(7FJJJJJFDHID)JCH=3DIJ5JGDI8@@I@A=>3<:IDCA9DDFFI(FADEBDCDCCDDB(DC>    PG:Z:SNAP   NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-30 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  9852129 70  100M    *   0   0   TGTGAAATGGAGTCAGCAGAGTGAGCCGGCCTCCACTCAGTGAGCCGGGTCTCCCCCACAGCCGGCATGTGCTGACCTCCTTCCAACTGCTCTACCAAGA    CBCDDFFDHGHHHIEJ+J<EFJIJIIJI><J(IIFJG)0JGIJ?8J5;J?D@9IJHDI=DI)DDHG@3FAI5FF?EDAHDC@DDGD3AA>D+?ECDDDDB  PG:Z:SNAP NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-29 16  chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  156397431   68  100M    *   0   0   TCAGCCTCCCGAGTAGCTGGGATTACAGGAACCTGCCACCACGCCCGGCTAATTTTTGTATTTTCAGTTGAGACGGGGTTTCACCATGTTGCCCAGGCTG    D1D(@9DDDC@D0C3=CDDJ;FDHDD@H2BDHIDAGDDDCDIFJ9GIFGIG@?)JJHJGFGJIB7JG>'IJIJJGJ+JJGIIHFIJIDHHHFFFFFFC@B PG:Z:SNAP NM:i:1 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-28 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 56986638 70 100M * 0 0 ACTCAGAACAGGTCTCCTTGTGGAACCATGGCCTTCCTTTTGGATCCTGGCCATGAGAGCCCATTCTTAGGAACCATGTTTCAATTCCAGTAGGTGATGT DD)DC@C<EDDD+DC0BDDBDDCECFDIJJ@)?HDACJGDFI?JGJH)JJJJJIEIJIIJJGIJIHIJHJCIGDHI>J@A)GIFJJJHHFFHFDDFFCCC PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-27 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 172015198 70 100M * 0 0 AGGTGTCAGTCCTCCAGCTTTGTTCTTCTTTTATATTGTGTTGGCTATCCTGGGCTCTTTGCTTCTCCATACAAAACTTAGAATCAGTTTGTTGATATCC B8BD>/D<BED@CCEBBEBCH,F?CCD.E;HGJBJ)IGD7HED5@6JJJCHIGHJIJFDJCIJJHGJIJJJIEF:FEJHBJ.JJJIGHBHCF2DDFFCC@ PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-26 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 233336763 70 100M * 0 0 AGATATACAGCAAAGTTTGAAAGCTACAGTTCTGAGGACCATATTTATGGATTCCTTCTTATATGTTATCTGGGTTGATATAGAAATTCTTCCATGGCTA CBCFDFF<H<G?AIBJJJGIEJIJIIJIIJJEGHEIJIIGI)GJHIJGF8JIIHED=DJH?IJFB;>;HHDDHHB?C9DE?DCE@D?B&5E>DDD7DD?D PG:Z:SNAP NM:i:1 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-25 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 196781488 0 100M * 0 0 GCCATTCTAACTGGTGTGAGATGGTATCTCATTGTGGTTTTGATTTGCATTTCTCTGATGGCCAGTGATGGTGAGCATTTTTTCATGTGTTTTTTGGCTG CCCDFFFFBGHHHHHCIGFJ:JAIGIJIJG)HCIIJGIHHJJJGEDHIHJHIII3J>JHJ?GDD?:;EFE(EDIJD?DDEAHCEDCDD?CDCF6D=>DDD PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-24 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 235841969 70 100M * 0 0 GTTGGCTACTAGCTTAGCAGAGGTGCAAAACCATGAATTTCTGGTGGTATGGATTTTTTCAGCTATTTCAGATTCACCAGCAGGATCCAGCTGCTTGGGT CCCFF?FFFHHHFGI,JEJIIG<JJ)I1GJG=ICJJEGJIJF<@IIBDJJIFDIEAIJB;JGADHJD,CBD@DEC;?DDHD<BEED&DD@DCDEDDDAD? PG:Z:SNAP NM:i:2 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-23 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 96545358 70 100M * 0 0 AGTGAAAAAGGCTGGCTGCCCTTCAATATCATCTTCAAATGTTAACAACACTGAATATTAATAAATTTCCTTTAGCGAATAATGAATCCAGCCTTCCTTA C@CF+FFFGGDHGJIBDJI2JGJIHHJJII?GJJJJGIJIJJGFJJG)IJ0HD0JIFJDJDFC;D7JGFFCEDFHADCDCCDEDDEAHDDD+9?<CA2:D PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-22 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 80270679 70 100M * 0 0 TTGTACACCCTATTTCTGACCAGAAGAAGGAGCATTTTGCTTTTTGCCAAATGAGAAGTGCATTCTGGAAACACTTGATGCCTGCACCACACCTCGAGTT ?@CFDDFFHFHHHJJJGC7J(GI8IJJJE?HHI>BJG*IJFJIDJHD0IEJIHDI>@H=EHGIAHJ33(EJCDEDA?FDG<ADDDCDDEF9DDDD@DDBD PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-21 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 35923261 70 100M * 0 0 CTAAGCAGCAGTGTTTTTGGATACTTTTTTTTTCTGTTTGTGAATAAGGCCAGCACTCAAGATGGGCAGCCAAGGGTGCACTGACTATTAGCTGGCCCAT =@@DFDFEHGHHH8JIJGJH1JJHHJIHJGH?IIFEJIIG87JI=IAJJJBJIJD(IIFI8JIHF=JDHEJHEHDDCEDCDEACDDCCAD<BDE+B8(DD PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-20 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 112489190 70 100M * 0 0 AGGGAATGAACTATGCACATCTATATAGTAACAGGGACAGATTTTTTTTTAACATGAGAGTGTAAAAAAAAGAAAAAGAAAAAAAAAGGCCAGGCACAGT DACDABD@DDDDDA7DDDC8GHI@EI(DC?FG'+8.FBDJIHIEGG=IIG=I@*DFIJJIBIIJIJIIHJCHBGFJJJI@F>HJIIIHHAHFAFDFFC@1    PG:Z:SNAP   NM:i:1  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-19 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  160371244   70  100M    *   0   0   AGGCCCTGGGCACAGGCAGAGAGCCCACCGGCTGGTCATGAGGGCCTCTTCCTTTCTCTGACCCAGGCACCTCGAGGGCTCTTCTCCTGGGTTCCTTCCG    @@:FDDFFCHHAHI:GEJFJGF@JJJFIC9JIIJJJ?IIEFHGJ'G?BFFBIIDIG,J)AJIHEGFBHCI&ECCD@EDD?)DED(D>3C?ABEEDDD4BD PG:Z:SNAP NM:i:1 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-18 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 179855835 70 100M * 0 0 AGCAATTAAAATAAATTAGGGTATCTTTAAAAGTTGTAAAATTATAGCAGTGAAGTACTGTTGACCAGGCACAGTGGCTCACACCTGTAATACCAGCACT DCEDBBDD/DD9DDD@DDFB(DDDHCHDF;C?;FJGC/IJ8DHEJ:DFGGIGHBIGIJDI(JDHGJJGJIHJII@HJJJ3JIJDIJBBHHFHFDFFF@C@ PG:Z:SNAP NM:i:2 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-17 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 207455995 70 100M * 0 0 GGTTCTTATGATTGGAAAGGTTAAAGAGTGACCTATAGGTCACTTTCCAATTATGAAAACAAAAAATTAAGAAATATATATATTTTCATTATTTCACTCC <DC>CBDDDCDD:&DDCFCFDDHDEJEDCFDJ;;EHGCD;CG?DIHGGCIJJJJ-GIJ7GIFHHHCGI)JJJJIJEGJIGJJJIH<GFGHHBFD@FDCC1 PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-16 0 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 114154603 70 100M * 0 0 TGTATCTTTCTGCTAAGCATAACAAGAAAGACAGAAAGCTCAACGGGAGGATTGAGGCTAGACTTAAAGTAGAGATCCCCTCAGAAACTGTGGAGTGAGG CCCF8FFDHHHH4JIJIGIJIIFJHJJ?JEDI9BG?I>GHJ7FJJJIF67EIIHD2C>?>DDHDE8E7@JEJ(IFDDC;EDCC:FD>@DBC>D5D>=<AB PG:Z:SNAP NM:i:0 RG:Z:FASTQ PL:Z:Illumina PU:Z:pu LB:Z:lb SM:Z:sm chr1-15 16 chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38 169767580 70 100M * 0 0 GGTGGGGGAGAGGAAAGGAAACGAGGGAGGAAAGGCCCTAATAGGGAGGATTTTGGAGTTTAGATTTTAAAATGATAAAGGTTGTTTGACACTCTAGGCA DEDD9DDD@DD4DDDAEDDDC@D7=D;DA)7;IIJFD(J?JJDGI(IDGD7D'3JIE;H?AC@EHJJE?JJHDFJIIIECG)GGJJECHFHHFDFDFC@C    PG:Z:SNAP   NM:i:2  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-14 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  117644126   70  100M    *   0   0   GCATTTCATTGTGGACTAATTTTCCCCCACTATTGAGGGAAGACCCTTTTGAGTACTCTATCTGATGCCCCATGAATGATAAAGTTTTATACTCTGGCTG    C?CFBDFFGHHH<JAGCIJJIIIGJJI8:JIJC(JGGJBJH-GIJIJJ;IIH;>JI5CJ=CD9DC-HGIJDCJHHDBEDDCC&DDDBD39DBCDDDDDCD    PG:Z:SNAP   NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-13 0   chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  104996994   70  100M    *   0   0   TTCTTGCTGGAACACATGTTTTCACCTTTACCTTCACCCACAGCCCAATGTGCATCAATATGGAGATAATGCAGTTCCATTTATACCTCTTTGTGGTTCA    =@?FFFDDHBHHHJHJIJ)JIJ<IJII++>HBHIJ*G:CJCJJJI?G)>GJI;JD3FJ8FJFGD;DDDDFBED7C<?E@C>7&A(ABC9CD+DCC&DDCA    PG:Z:SNAP   NM:i:3  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
chr1-12 16  chr1__AC:CM000663.2__gi:568336023__LN:248956422__rl:Chromosome__M5:6aef897c3d6ff0c78aff06ac189178dd__AS:GRCh38  108617705   70  100M    *   0   0   AGGTCGGGGAGATTGGGAAGAAGAATGAGCAAAGAAACCACCAGTGTGATCAGAGGAGGAAAGCAAAGCAGAGTCCTGTCCTGAAAACCAAATGAAGAAA    :=>+D(DCEC=@GHB(CDDDDHABDD+HBJJ9F?A35DDIE?JJHIHJJIEE?JFJ?7JBGJJI>JJGJBJIIBIJJJIIIIJGJGJHHDFHF3FFFCB@    PG:Z:SNAP   NM:i:0  RG:Z:FASTQ  PL:Z:Illumina   PU:Z:pu LB:Z:lb SM:Z:sm
hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub/cs-bwamem$ 

（编辑：晋中站长网）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!