# HG changeset patch # User mvdbeek # Date 1522568296 14400 # Node ID 29734dbad6de1bcaac0b23a001d8563355333d7e # Parent bb09f4ee2ff830081c2ec1d2b39b0168f2857237 planemo upload for repository https://github.com/bardin-lab/readtagger/tree/master/galaxy commit 2e991a52c828b0234020e89748eafe1028e44a7d diff -r bb09f4ee2ff8 -r 29734dbad6de add_matesequence.xml --- a/add_matesequence.xml Sat Jan 06 07:21:07 2018 -0500 +++ b/add_matesequence.xml Sun Apr 01 03:38:16 2018 -0400 @@ -1,7 +1,7 @@ - + into tag field - readtagger + readtagger add_matesequence --version + modifies proper_pair flag in bam files - readtagger + readtagger + from multiple bam files macros.xml - readtagger + readtagger + based on softclips in control files + + readtagger + + confirm_insertions --version + + + + + + + + + + + + + + + + + + + + diff -r bb09f4ee2ff8 -r 29734dbad6de findcluster.xml --- a/findcluster.xml Sat Jan 06 07:21:07 2018 -0500 +++ b/findcluster.xml Sun Apr 01 03:38:16 2018 -0400 @@ -1,10 +1,10 @@ - + in bam files macros.xml - readtagger + readtagger findcluster --version - + + + + - - - + + make_bam + + + make_fasta + + + make_vcf + + + make_gff + + - + + + diff -r bb09f4ee2ff8 -r 29734dbad6de plot_coverage.xml --- a/plot_coverage.xml Sat Jan 06 07:21:07 2018 -0500 +++ b/plot_coverage.xml Sun Apr 01 03:38:16 2018 -0400 @@ -1,7 +1,7 @@ - + as area plot between BAM files - readtagger + readtagger + diff -r bb09f4ee2ff8 -r 29734dbad6de test-data/all_control_insertions.gff --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/all_control_insertions.gff Sun Apr 01 03:38:16 2018 -0400 @@ -0,0 +1,3 @@ +TEST find_softclip 5p_clip 16727068 16727111 6 + . ID=softclip_test2;consensus=TTTTTTT +X find_softclip 5p_clip 16727068 16727067 3 + . ID=softclip_0;consensus=ACTTAATA;max_mapq=25 +X find_softclip 3p_clip 16727112 16727111 4 + . ID=softclip_1;consensus=TTTAATAAAATTAAAGGACAACTAAACGGAAGGAAAATATCGTGAGCAG;max_mapq=28 diff -r bb09f4ee2ff8 -r 29734dbad6de test-data/all_treatment_insertions.gff --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/all_treatment_insertions.gff Sun Apr 01 03:38:16 2018 -0400 @@ -0,0 +1,13 @@ +TEST find_softclip 5p_clip 16727068 16727111 6 + . ID=softclip_test1;consensus=AAAAAA +X feature predicted_insertion 10766937 10766942 . + . Parent=R37_49;contig_support=2;fraction_full_length=0.03132530120481928;sbjct=RR48388_transposable_element_rover_Gypsy;sbjct_end=7308;sbjct_start=7074;type=predicted_insertion +X findcluster transposable_element_rover 10766937 10766942 24 + . ID=R37_49;genotype=heterozygous;genotype_likelihoods=6.86863481923e-26,0.999999994797,5.2033283093e-09;insert_reference_name=transposable_element_rover;left_inserts=0,GGGGTGCATATATAGCGCACATATGTACATATTTGTACATATATAGTAACATAATATGCTTCTCATATTACGTTTACATACTTACACTAATTGTACATACAATCTTGCACATGCA,1,CAACAATGTATTATCCTACCGACTGCGCCACTTACGTATTTTACTATTCGAGTCGAACAAATGATCCGTCGTTTGACTAAGATCAACGCCTTTAAAGAAGTTTCAGAATCAACACGAGGTATTTTTTTAGTTTGGA;left_mate_count=4;max_mapq=60;nref=12;right_inserts=0,GGCGTTGATCTTAGTCAAACGACGGATCATTTGTTCGACTCGAATAGTAAAATACGTAAGTATATACATATGCATATGCAGTCCGAATCAAAATCGAAAACAAGAACTGCACTTTTCGCTGCGATTGGGGCACCTTCAGC,1,GTTTACATACTTACACTAATTGTACATACAATCTTGCACATGCATAAACACATCAAACCAGCTTACATTTTTACTTACACTTAAGCGCATGATTTGTTGTGCATCCATACCGTTATTTTTCCGTTCTTTTTTTGTACACATATACTGATTAGACATTCCCGTTTAT,2,CACATATACTGATTAGACATTCCCGTTTATCGCGACTCACTTCGAGCCGATCAAAAACTCTGTACAGTCAGTCTTAAGCCGACAACGAAGAAATAAAGATCCAAACTAAAAAAATACCTCGTGTTGATTCTGAAACTTCTTTAAAGGCGTTGATCTTAGTCAAACGACGGATCATTTGTT;right_mate_count=14;softclip_clusters=;total_left_count=6;total_right_count=18;valid_TSD=True +X find_softclip 5p_clip 10766937 10766936 6 + . ID=softclip_0;consensus=GGCGTTGATCTTAGTCAAACGACGGATCATTTGTTCGACTCGAATAGTAAAATACGTAAG;max_mapq=60 +X find_softclip 3p_clip 10766943 10766942 3 + . ID=softclip_1;consensus=GTAACATAATATGCTTCTCATATTACGTTTACATACTTACACTAATTGTACATACAATCTTGCACATGCA;max_mapq=60 +X feature left_insert 16726713 16727261 . + . Parent=R37_26;fraction_full_length=0.003280519299225239;read_support=1;sbjct=FBti0059747_gtwin_Gypsy;sbjct_end=2536;sbjct_start=2489;type=left_insert +X findcluster gtwin 16726713 16727261 1 + . ID=R37_26;genotype=reference;genotype_likelihoods=1.0,1.93425522512e-14,6.88207144991e-72;insert_reference_name=gtwin;left_inserts=;left_mate_count=0;max_mapq=60;nref=47;right_inserts=;right_mate_count=1;softclip_clusters=;total_left_count=0;total_right_count=1;valid_TSD=False +X feature left_insert 16727068 16727111 . + . Parent=R37_27;fraction_full_length=0.6875;read_support=1;sbjct=FBti0062261_jockey2_Jockey;sbjct_end=63;sbjct_start=19;type=left_insert +X feature right_insert 16727068 16727111 . + . Parent=R37_27;fraction_full_length=0.0017704870646047116;read_support=1;sbjct=FBti0062666_HMS-Beagle_Gypsy;sbjct_end=10537;sbjct_start=10488;type=right_insert +X findcluster jockey2 16727068 16727111 6 + . ID=R37_27;genotype=heterozygous;genotype_likelihoods=8.55882087985e-05,0.999914411654,1.37666001367e-10;insert_reference_name=jockey2;left_inserts=;left_mate_count=0;max_mapq=25;nref=9;right_inserts=0,GATATGGACTTAATTGGAAATTAATTAATGCAATAATCCAATTAAGAATTTAATATTTAAAAAGGTGGACTGAATACTTAATACAAGAAGGATAATAGTACGCTAGTCAATTAAATAAAATAATTTCTTTAATAAAATTAAAGGAC;right_mate_count=2;softclip_clusters=softclip_2,softclip_3;total_left_count=3;total_right_count=4;valid_TSD=True +X find_softclip 3p_clip 16727068 16727111 6 + . ID=softclip_2;consensus=TTTAATAAAATTAAAGGAC;max_mapq=34 +X find_softclip 5p_clip 16727068 16727111 5 + . ID=softclip_3;consensus=GAATTTAATATTTAAAAAGGTGGACTGAATACTTAATA;max_mapq=34 +X find_softclip 3p_clip 16727110 16727109 1 + . ID=softclip_4;consensus=CCTTTAATAAAATTAAAGGACAACTAAACGGAAGGAAAATATCGTGAGCAGAAAATAATATTAAGTCTTTAAAAG;max_mapq=34 diff -r bb09f4ee2ff8 -r 29734dbad6de test-data/confirmed_insertions.gff --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/confirmed_insertions.gff Sun Apr 01 03:38:16 2018 -0400 @@ -0,0 +1,4 @@ +2R findcluster jockey2 16727068 16727111 6 + . ID=R37_27;softclip_clusters=softclip_test1 +TEST findcluster jockey2 16727068 16727111 6 + . ID=R37_27;softclip_clusters=softclip_test1 +X findcluster transposable_element_rover 10766937 10766942 24 + . ID=R37_49;genotype=heterozygous;genotype_likelihoods=6.86863481923e-26,0.999999994797,5.2033283093e-09;insert_reference_name=transposable_element_rover;left_inserts=0,GGGGTGCATATATAGCGCACATATGTACATATTTGTACATATATAGTAACATAATATGCTTCTCATATTACGTTTACATACTTACACTAATTGTACATACAATCTTGCACATGCA,1,CAACAATGTATTATCCTACCGACTGCGCCACTTACGTATTTTACTATTCGAGTCGAACAAATGATCCGTCGTTTGACTAAGATCAACGCCTTTAAAGAAGTTTCAGAATCAACACGAGGTATTTTTTTAGTTTGGA;left_mate_count=4;max_mapq=60;nref=12;right_inserts=0,GGCGTTGATCTTAGTCAAACGACGGATCATTTGTTCGACTCGAATAGTAAAATACGTAAGTATATACATATGCATATGCAGTCCGAATCAAAATCGAAAACAAGAACTGCACTTTTCGCTGCGATTGGGGCACCTTCAGC,1,GTTTACATACTTACACTAATTGTACATACAATCTTGCACATGCATAAACACATCAAACCAGCTTACATTTTTACTTACACTTAAGCGCATGATTTGTTGTGCATCCATACCGTTATTTTTCCGTTCTTTTTTTGTACACATATACTGATTAGACATTCCCGTTTAT,2,CACATATACTGATTAGACATTCCCGTTTATCGCGACTCACTTCGAGCCGATCAAAAACTCTGTACAGTCAGTCTTAAGCCGACAACGAAGAAATAAAGATCCAAACTAAAAAAATACCTCGTGTTGATTCTGAAACTTCTTTAAAGGCGTTGATCTTAGTCAAACGACGGATCATTTGTT;right_mate_count=14;total_left_count=6;total_right_count=18;valid_TSD=True +X findcluster jockey2 16727068 16727111 6 + . ID=R37_27;genotype=heterozygous;genotype_likelihoods=8.55882087985e-05,0.999914411654,1.37666001367e-10;insert_reference_name=jockey2;left_inserts=;left_mate_count=0;max_mapq=25;nref=9;right_inserts=0,GATATGGACTTAATTGGAAATTAATTAATGCAATAATCCAATTAAGAATTTAATATTTAAAAAGGTGGACTGAATACTTAATACAAGAAGGATAATAGTACGCTAGTCAATTAAATAAAATAATTTCTTTAATAAAATTAAAGGAC;right_mate_count=2;softclip_clusters=softclip_2,softclip_3;total_left_count=3;total_right_count=4;valid_TSD=True;FAIL=clip_seq_in_control diff -r bb09f4ee2ff8 -r 29734dbad6de test-data/dm6.bam Binary file test-data/dm6.bam has changed diff -r bb09f4ee2ff8 -r 29734dbad6de test-data/putative_insertions.gff --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/putative_insertions.gff Sun Apr 01 03:38:16 2018 -0400 @@ -0,0 +1,4 @@ +2R findcluster jockey2 16727068 16727111 6 + . ID=R37_27;softclip_clusters=softclip_test1 +TEST findcluster jockey2 16727068 16727111 6 + . ID=R37_27;softclip_clusters=softclip_test1 +X findcluster transposable_element_rover 10766937 10766942 24 + . ID=R37_49;genotype=heterozygous;genotype_likelihoods=6.86863481923e-26,0.999999994797,5.2033283093e-09;insert_reference_name=transposable_element_rover;left_inserts=0,GGGGTGCATATATAGCGCACATATGTACATATTTGTACATATATAGTAACATAATATGCTTCTCATATTACGTTTACATACTTACACTAATTGTACATACAATCTTGCACATGCA,1,CAACAATGTATTATCCTACCGACTGCGCCACTTACGTATTTTACTATTCGAGTCGAACAAATGATCCGTCGTTTGACTAAGATCAACGCCTTTAAAGAAGTTTCAGAATCAACACGAGGTATTTTTTTAGTTTGGA;left_mate_count=4;max_mapq=60;nref=12;right_inserts=0,GGCGTTGATCTTAGTCAAACGACGGATCATTTGTTCGACTCGAATAGTAAAATACGTAAGTATATACATATGCATATGCAGTCCGAATCAAAATCGAAAACAAGAACTGCACTTTTCGCTGCGATTGGGGCACCTTCAGC,1,GTTTACATACTTACACTAATTGTACATACAATCTTGCACATGCATAAACACATCAAACCAGCTTACATTTTTACTTACACTTAAGCGCATGATTTGTTGTGCATCCATACCGTTATTTTTCCGTTCTTTTTTTGTACACATATACTGATTAGACATTCCCGTTTAT,2,CACATATACTGATTAGACATTCCCGTTTATCGCGACTCACTTCGAGCCGATCAAAAACTCTGTACAGTCAGTCTTAAGCCGACAACGAAGAAATAAAGATCCAAACTAAAAAAATACCTCGTGTTGATTCTGAAACTTCTTTAAAGGCGTTGATCTTAGTCAAACGACGGATCATTTGTT;right_mate_count=14;total_left_count=6;total_right_count=18;valid_TSD=True +X findcluster jockey2 16727068 16727111 6 + . ID=R37_27;genotype=heterozygous;genotype_likelihoods=8.55882087985e-05,0.999914411654,1.37666001367e-10;insert_reference_name=jockey2;left_inserts=;left_mate_count=0;max_mapq=25;nref=9;right_inserts=0,GATATGGACTTAATTGGAAATTAATTAATGCAATAATCCAATTAAGAATTTAATATTTAAAAAGGTGGACTGAATACTTAATACAAGAAGGATAATAGTACGCTAGTCAATTAAATAAAATAATTTCTTTAATAAAATTAAAGGAC;right_mate_count=2;softclip_clusters=softclip_2,softclip_3;total_left_count=3;total_right_count=4;valid_TSD=True diff -r bb09f4ee2ff8 -r 29734dbad6de test-data/three_cluster_out.bam Binary file test-data/three_cluster_out.bam has changed diff -r bb09f4ee2ff8 -r 29734dbad6de test-data/three_cluster_out.gff --- a/test-data/three_cluster_out.gff Sat Jan 06 07:21:07 2018 -0500 +++ b/test-data/three_cluster_out.gff Sun Apr 01 03:38:16 2018 -0400 @@ -1,3 +1,8 @@ ##gff-version 3 -3R findcluster TE 13373515 13373524 22 + . ID=extended_and_annotated_roi.bam_0;genotype=homozygous;genotype_likelihoods=1.8828617809e-31,9.53673406912e-07,0.999999046327;left_insert=0,CTCGGAATGTATCTAACTAACAAACTCATATCAAATATAAGCAAGTGCGCCAATTCGTATGCATATGGACATATGGACATATACATATAGTAACATAATATGCTTCTCATATTACGTTTACATACTTACACTAATTGTACATACAATCTTGCACATGCATAAACACATCAAACCAGCTTACATTTTTACTTACACTTAAGCGCATGATTTGTTGTGCATCCATACCGTTATTTTTCC;left_mate_support=8;left_support=11;non_support=0;right_insert=0,GACAACGAAGAAATAAAGATCCAAACTAAAAAAATACCTCGTGTTGATTCTGAAACTTCTTTAAAGGCGTTGATCTTAGTCAAACGACGGATCATTTGTTCGACTCGAATAGTAAAATACGTAAGTATATAGATAGTCTATATTAATTTTAAAAGCTCAAAGGGGCGCAAGTCTCTCTCATCAGTTGTGTCTTTAGTTTTTTTGATTTGGTTTGGTG;right_mate_support=6;right_support=11;valid_TSD=False -3R findcluster TE 13374595 13374676 3 + . ID=extended_and_annotated_roi.bam_1;genotype=homozygous;genotype_likelihoods=0.00038864889331,0.199922270221,0.799689080885;left_insert=0,GTTCACCCGCGTCCGAGTTCCTGCTCCACTACTCCCTGGCTGCTGACTCACTGTTGTTATAGGGGTGGCTTCCCCTCTGTTCTTCCTGGGGGAATGCTGCATCTTCCCCAGCTCCAAAATGGCGG;left_mate_support=2;left_support=2;non_support=0;right_insert=;right_mate_support=1;right_support=1;valid_TSD=False +3R findcluster TE 13373206 13373206 2 + . ID=extended_and_annotated_roi.bam_0;genotype=homozygous;genotype_likelihoods=0.010684606252473286,0.32977179791584216,0.6595435958316843;insert_reference_name=None;left_inserts=;left_mate_count=2;max_mapq=60;nref=0;right_inserts=;right_mate_count=0;softclip_clusters=;total_left_count=2;total_right_count=0;valid_TSD=False +3R find_softclip 3p_clip 13373515 13373514 7 + . ID=SOFTCLIP_8392192716115976782;consensus=ACATATAGTAACATAATATGCTTCTCATATTACGTTTACATACTTACACTAATTGTAC;max_mapq=60 +3R findcluster TE 13373515 13373524 20 + . ID=extended_and_annotated_roi.bam_1;genotype=homozygous;genotype_likelihoods=6.2761999508986e-30,1.9073449948406318e-06,0.9999980926550052;insert_reference_name=None;left_inserts=GGAAAAATAACGGTATGGATGCACAACAAATCATGCGCTTAAGTGTAAGTAAAAATGTAAGCTGGTTTGATGTGTTTATGCATGTGCAAGATTGTATGTACAATTAGTGTAAGTATGTAAACGTAATATGAGAAGCATATTATGTTACTATATGTATATGTCCATATGTCCATATGCATACGAATTGGCGCACTTGCTTATATTTGATATGAGTTTGTTAGTTAGATACATTCCGAG;left_mate_count=6;max_mapq=60;nref=0;right_inserts=GACAACGAAGAAATAAAGATCCAAACTAAAAAAATACCTCGTGTTGATTCTGAAACTTCTTTAAAGGCGTTGATCTTAGTCAAACGACGGATCATTTGTTCGACTCGAATAGTAAAATACGTAAGTATATAGATAGTCTATATTAATTTTAAAAGCTCAAAGGGGCGCAAGTCTCTCTCATCAGTTGTGTCTTTAGTTTTTTTGATTTGGTTTGGTG;right_mate_count=6;softclip_clusters=SOFTCLIP_8392192716115976782,SOFTCLIP_3495167510329773485;total_left_count=9;total_right_count=11;valid_TSD=False +3R find_softclip 5p_clip 13373525 13373524 6 + . ID=SOFTCLIP_3495167510329773485;consensus=CTTAGTCAAACGACGGATCATTTGTTCGACTCGAATAGTAAAATACGTAAGT;max_mapq=60 +3R findcluster TE 13374539 13374676 1 + . ID=extended_and_annotated_roi.bam_2;genotype=homozygous;genotype_likelihoods=0.21259842519685054,0.3937007874015747;insert_reference_name=None;left_inserts=;left_mate_count=0;max_mapq=60;nref=0;right_inserts=;right_mate_count=1;softclip_clusters=;total_left_count=0;total_right_count=1;valid_TSD=False +3R findcluster TE 13374595 13374709 2 + . ID=extended_and_annotated_roi.bam_3;genotype=homozygous;genotype_likelihoods=0.010684606252473286,0.32977179791584216,0.6595435958316843;insert_reference_name=None;left_inserts=GTTCACCCGCGTCCGAGTTCCTGCTCCACTACTCCCTGGCTGCTGACTCACTGTTGTTATAGGGGTGGCTTCCCCTCTGTTCTTCCTGGGGGAATGCTGCATCTTCCCCAGCTCCAAAATGGCGG;left_mate_count=2;max_mapq=60;nref=0;right_inserts=;right_mate_count=0;softclip_clusters=;total_left_count=2;total_right_count=0;valid_TSD=False +3R find_softclip 5p_clip 13374677 13374676 1 + . ID=SOFTCLIP_2520711244791099460;consensus=TGGCAAAACATATGCAAAGTGCA;max_mapq=60 diff -r bb09f4ee2ff8 -r 29734dbad6de test-data/three_cluster_out.vcf --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/three_cluster_out.vcf Sun Apr 01 03:38:16 2018 -0400 @@ -0,0 +1,1916 @@ +##fileformat=VCFv4.1 +##FILTER= +##fileDate=2018-03-30 15:24:43.667621 +##source=readtagger-v1.2.3 +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##contig= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##ALT= +##ALT= +##ALT= +##ALT= +##ALT= +##ALT= +##ALT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +##FORMAT= +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT extended_and_annotated_roi.bam +3R 13373206 INS_2444456283915041025 N . . END=13373207;SVTYPE=INS:ME;SVLEN=.;MQ=60;EVENT=INS_2444456283915041025;MATEID GT:GL:AD:DP:SU:SU5:SU3:SR:SR5:SR3:MSP:MENAME:MESTART:MEEND:MEASSEMBLY5:MEASSEMBLY3 1/1:0.0106846,0.329772,0.659544:0,2:2:2:2:0:0:0:0:0::.:.:: +3R 13373515 INS_630831370416389201 N . . END=13373525;SVTYPE=INS:ME;SVLEN=.;MQ=60;EVENT=INS_630831370416389201;MATEID=SOFTCLIP_8392192716115976782,SOFTCLIP_3495167510329773485,SOFTCLIP_8392192716115976782,SOFTCLIP_3495167510329773485 GT:GL:AD:DP:SU:SU5:SU3:SR:SR5:SR3:MSP:MENAME:MESTART:MEEND:MEASSEMBLY5:MEASSEMBLY3 1/1:6.2762e-30,1.90734e-06,0.999998:0,20:20:20:9:11:8:3:5:0::.:.:GGAAAAATAACGGTATGGATGCACAACAAATCATGCGCTTAAGTGTAAGTAAAAATGTAAGCTGGTTTGATGTGTTTATGCATGTGCAAGATTGTATGTACAATTAGTGTAAGTATGTAAACGTAATATGAGAAGCATATTATGTTACTATATGTATATGTCCATATGTCCATATGCATACGAATTGGCGCACTTGCTTATATTTGATATGAGTTTGTTAGTTAGATACATTCCGAG:GACAACGAAGAAATAAAGATCCAAACTAAAAAAATACCTCGTGTTGATTCTGAAACTTCTTTAAAGGCGTTGATCTTAGTCAAACGACGGATCATTTGTTCGACTCGAATAGTAAAATACGTAAGTATATAGATAGTCTATATTAATTTTAAAAGCTCAAAGGGGCGCAAGTCTCTCTCATCAGTTGTGTCTTTAGTTTTTTTGATTTGGTTTGGTG +3R 13373515 SOFTCLIP_8392192716115976782 N . . END=13373514;SVTYPE=;MQ=60;EVENT=INS_630831370416389201;MATEID=INS_630831370416389201 CLIP_CONSENSUS:AD:SU ACATATAGTAACATAATATGCTTCTCATATTACGTTTACATACTTACACTAATTGTAC:0,7:7 +3R 13373525 SOFTCLIP_3495167510329773485 N . . END=13373524;SVTYPE=;MQ=60;EVENT=INS_630831370416389201;MATEID=INS_630831370416389201 CLIP_CONSENSUS:AD:SU CTTAGTCAAACGACGGATCATTTGTTCGACTCGAATAGTAAAATACGTAAGT:0,6:6 +3R 13374539 INS_2520711244791099460 N . . END=13374677;SVTYPE=INS:ME;SVLEN=.;MQ=60;EVENT=INS_2520711244791099460;MATEID GT:GL:AD:DP:SU:SU5:SU3:SR:SR5:SR3:MSP:MENAME:MESTART:MEEND:MEASSEMBLY5:MEASSEMBLY3 1/1:0.212598,0.393701,0.393701:0,1:1:1:0:1:0:0:0:0::.:.:: +3R 13374595 INS_7083183827793206730 N . . END=13374710;SVTYPE=INS:ME;SVLEN=.;MQ=60;EVENT=INS_7083183827793206730;MATEID GT:GL:AD:DP:SU:SU5:SU3:SR:SR5:SR3:MSP:MENAME:MESTART:MEEND:MEASSEMBLY5:MEASSEMBLY3 1/1:0.0106846,0.329772,0.659544:0,2:2:2:2:0:0:0:0:0::.:.:GTTCACCCGCGTCCGAGTTCCTGCTCCACTACTCCCTGGCTGCTGACTCACTGTTGTTATAGGGGTGGCTTCCCCTCTGTTCTTCCTGGGGGAATGCTGCATCTTCCCCAGCTCCAAAATGGCGG: +3R 13374677 SOFTCLIP_2520711244791099460 N . . END=13374676;SVTYPE=;MQ=60;EVENT;MATEID CLIP_CONSENSUS:AD:SU TGGCAAAACATATGCAAAGTGCA:0,1:1 diff -r bb09f4ee2ff8 -r 29734dbad6de update_mapq.xml --- a/update_mapq.xml Sat Jan 06 07:21:07 2018 -0500 +++ b/update_mapq.xml Sun Apr 01 03:38:16 2018 -0400 @@ -1,7 +1,7 @@ - + of supplementary alignments - readtagger + readtagger update_mapq --version + from SAM/BAM alignment files as FASTQ - readtagger + readtagger write_supplementary_fastq --version