diff --git a/include/misc/statistics.hpp b/include/misc/statistics.hpp new file mode 100644 index 00000000..070ae97c --- /dev/null +++ b/include/misc/statistics.hpp @@ -0,0 +1,53 @@ +#pragma once + +#include + +#include + +/*! \brief Returns a list of pairs, where the first integer is a random chromosome and the second is a random position. + * The list is sorted by default. + * \param sample_size - The number of positions to generate. + * \param header - The header from the SAM file input. + * + * \return Returns a sorted forward_list of pairs of reference chromosome IDs and positions. + */ +std::forward_list> get_random_positions(uint64_t sample_size, seqan3::sam_file_header> & header) +{ + std::forward_list> output{}; + if (sample_size == 0) + { + output.push_front(std::pair{-1, -1}); + } + int32_t rand_chr{}; + int32_t rand_pos{}; + + for (int i = 0; i < sample_size; ++i) + { + rand_chr = std::rand() % header.ref_ids().size(); + rand_pos = std::rand() % std::get<0>(header.ref_id_info[rand_chr]); + + output.push_front(std::pair{rand_chr, rand_pos}); + } + + output.sort(); + return output; +} + +/*! \brief Get the read length of a read from its cigar string. + * \param cigar - The cigar according to the seqan3 representation (a vector of cigar alphabet). + * + * \return Returns the length of the cigar string. + */ + +uint32_t get_read_length(std::vector & cigar) +{ + using seqan3::get; + + uint32_t result{0}; + for (auto e : cigar) + { + result += get<0>(e); + } + + return result; +} diff --git a/include/variant_detection/variant_detection.hpp b/include/variant_detection/variant_detection.hpp index d94e2d74..3e04e0f9 100644 --- a/include/variant_detection/variant_detection.hpp +++ b/include/variant_detection/variant_detection.hpp @@ -23,6 +23,7 @@ * 2: sVirl_refinement_method) * \param min_var_length - minimum length of variants to detect (default 30 bp) * \param output_file_path output file - path for the VCF file + * \param sample_size - how many positions to sample for read depth estimation. * * \details Detects novel junctions from read alignment records using different detection methods. * The junctions are clustered using one of several clustering methods. @@ -36,4 +37,5 @@ void detect_variants_in_alignment_file(const std::filesystem::path & alignment_f const clustering_methods & clustering_method, const refinement_methods & refinement_method, const uint64_t & min_var_length, - const std::filesystem::path & output_file_path); + const std::filesystem::path & output_file_path, + const uint64_t & sample_size); diff --git a/src/iGenVar.cpp b/src/iGenVar.cpp index 71c85b2d..8710a5b7 100644 --- a/src/iGenVar.cpp +++ b/src/iGenVar.cpp @@ -64,6 +64,7 @@ struct cmd_arguments clustering_methods clustering_method{simple_clustering}; // default: simple clustering method refinement_methods refinement_method{no_refinement}; // default: no refinement uint64_t min_var_length = 30; + uint64_t sample_size{100000}; // default is to sample 100000 }; void initialize_argument_parser(seqan3::argument_parser & parser, cmd_arguments & args) @@ -115,6 +116,11 @@ void initialize_argument_parser(seqan3::argument_parser & parser, cmd_arguments parser.add_option(args.min_var_length, 'l', "min_var_length", "Specify what should be the minimum length of your SVs to be detected (default 30 bp).", seqan3::option_spec::advanced); + + // Options - Miscellaneous + parser.add_option(args.sample_size, 's', "sample_size", + "This value is used in sampling positions for average read depth, insert size, and read length. Setting to 0 will skip sampling.", + seqan3::option_spec::advanced); } int main(int argc, char ** argv) @@ -151,7 +157,8 @@ int main(int argc, char ** argv) args.clustering_method, args.refinement_method, args.min_var_length, - args.output_file_path); + args.output_file_path, + args.sample_size); return 0; } diff --git a/src/variant_detection/variant_detection.cpp b/src/variant_detection/variant_detection.cpp index 9aeef898..0d06edf3 100644 --- a/src/variant_detection/variant_detection.cpp +++ b/src/variant_detection/variant_detection.cpp @@ -4,6 +4,7 @@ #include // SAM/BAM support #include // FASTA support +#include "misc/statistics.hpp" // for calculating average read depth #include "modules/clustering/simple_clustering_method.hpp" // for the simple clustering method #include "modules/sv_detection_methods/analyze_cigar_method.hpp" // for the split read method #include "modules/sv_detection_methods/analyze_sa_tag_method.hpp" // for the cigar string method @@ -20,7 +21,8 @@ void detect_variants_in_alignment_file(const std::filesystem::path & alignment_f const clustering_methods & clustering_method, const refinement_methods & refinement_method, const uint64_t & min_var_length, - const std::filesystem::path & output_file_path) + const std::filesystem::path & output_file_path, + const uint64_t & sample_size) { // Open input alignment file using my_fields = seqan3::fields insertion_alleles{}; uint16_t num_good = 0; + // Obtain list of random chromosomes and positions to get read depth from. + std::forward_list> rand_sample = get_random_positions(sample_size, alignment_file.header()); + uint32_t avg_depth{0}; + uint32_t cur_depth{0}; + uint32_t num_pos{0}; + for (auto & rec : alignment_file) { const std::string query_name = seqan3::get(rec); // 1: QNAME @@ -60,6 +68,26 @@ void detect_variants_in_alignment_file(const std::filesystem::path & alignment_f continue; const std::string ref_name = ref_ids[ref_id]; + uint32_t read_length{get_read_length(cigar)}; + // If there are positions to be sampled, then sample them. + if (!rand_sample.empty() && sample_size > 0) + { + if (ref_id == std::get<0>(rand_sample.front()) && // Correct chromosome. + pos <= std::get<1>(rand_sample.front()) && // Read starts before or at the position. + std::get<1>(rand_sample.front()) <= pos + read_length) // Read ends after or at position. + { + ++cur_depth; + } + else if (ref_id > std::get<0>(rand_sample.front()) || // Record is at next chromosome + (ref_id == std::get<0>(rand_sample.front()) && // Record is at same chromosome BUT + pos > std::get<1>(rand_sample.front()))) // in a later position. + { + ++num_pos; + avg_depth += cur_depth; + cur_depth = 0; + rand_sample.pop_front(); + } + } for (uint8_t method : methods) { switch (method) { @@ -101,6 +129,14 @@ void detect_variants_in_alignment_file(const std::filesystem::path & alignment_f seqan3::debug_stream << num_good << " good alignments" << std::endl; } } + if (sample_size == 0 || num_pos == 0) // User requested no sampling, or all of the random positions generated were not covered. + { + avg_depth = 0; + } + else + { + avg_depth = avg_depth/num_pos; + } std::sort(junctions.begin(), junctions.end()); seqan3::debug_stream << "Start clustering...\n"; diff --git a/test/api/junction_detection_methods_test.cpp b/test/api/junction_detection_methods_test.cpp index eac47ba0..321749f0 100644 --- a/test/api/junction_detection_methods_test.cpp +++ b/test/api/junction_detection_methods_test.cpp @@ -8,6 +8,7 @@ std::filesystem::path tmp_dir = std::filesystem::temp_directory_path(); // g std::filesystem::path empty_path{}; const std::vector all_methods{cigar_string, split_read, read_pairs, read_depth}; const uint64_t sv_default_length = 30; +const uint64_t sample_size = 10; // Explanation for the strings: // Reference\tm2257/8161/CCS\t41972616\tForward\tRead\t0\t2294\tForward\tchr21 @@ -73,7 +74,8 @@ TEST(junction_detection, single_method_only) simple_clustering, no_refinement, sv_default_length, - empty_path); + empty_path, + sample_size); //TODO (eldariont): Currently, this (CLI-like) test compares the stdout of the method with an expected string in VCF format. //We need to replace this with real API tests that check inidividual parts of the pipeline. @@ -121,7 +123,8 @@ TEST(junction_detection, method_pairs) simple_clustering, no_refinement, sv_default_length, - empty_path); + empty_path, + sample_size); //TODO (eldariont): Currently, this (CLI-like) test compares the stdout of the method with an expected string in VCF format. //We need to replace this with real API tests that check inidividual parts of the pipeline. @@ -184,7 +187,8 @@ TEST(junction_detection, method_triples) simple_clustering, no_refinement, sv_default_length, - empty_path); + empty_path, + sample_size); //TODO (eldariont): Currently, this (CLI-like) test compares the stdout of the method with an expected string in VCF format. //We need to replace this with real API tests that check inidividual parts of the pipeline. diff --git a/test/api/junction_detection_test.cpp b/test/api/junction_detection_test.cpp index ba061a0e..60f42ca4 100644 --- a/test/api/junction_detection_test.cpp +++ b/test/api/junction_detection_test.cpp @@ -6,6 +6,7 @@ std::filesystem::path tmp_dir = std::filesystem::temp_directory_path(); // g std::filesystem::path empty_path{}; const std::vector default_methods{cigar_string, split_read, read_pairs, read_depth}; const uint64_t sv_default_length = 30; +const uint64_t sample_size = 10; // Explanation for the strings: // Reference\tchr21\t41972616\tForward\tRead\t0\t2294\tForward\t1 @@ -56,7 +57,8 @@ TEST(junction_detection, fasta_out_not_empty) simple_clustering, no_refinement, sv_default_length, - empty_path); + empty_path, + sample_size); // check_output_and_cleanup(expected_res); check_output_and_cleanup(empty_res); @@ -75,7 +77,8 @@ TEST(junction_detection, clustering_method_hierarchical) hierarchical_clustering, no_refinement, sv_default_length, - empty_path); + empty_path, + sample_size); // check_output_and_cleanup(""); check_output_and_cleanup(empty_res); @@ -92,7 +95,8 @@ TEST(junction_detection, clustering_method_self_balancing_binary_tree) self_balancing_binary_tree, no_refinement, sv_default_length, - empty_path); + empty_path, + sample_size); // check_output_and_cleanup(""); check_output_and_cleanup(empty_res); @@ -109,7 +113,8 @@ TEST(junction_detection, clustering_method_candidate_selection_based_on_voting) candidate_selection_based_on_voting, no_refinement, sv_default_length, - empty_path); + empty_path, + sample_size); // check_output_and_cleanup(""); check_output_and_cleanup(empty_res); @@ -128,7 +133,8 @@ TEST(junction_detection, refinement_method_sViper) simple_clustering, sViper_refinement_method, sv_default_length, - empty_path); + empty_path, + sample_size); // check_output_and_cleanup(expected_res); check_output_and_cleanup(empty_res); @@ -145,7 +151,8 @@ TEST(junction_detection, refinement_method_sVirl) simple_clustering, sVirl_refinement_method, sv_default_length, - empty_path); + empty_path, + sample_size); // check_output_and_cleanup(expected_res); check_output_and_cleanup(empty_res); diff --git a/test/cli/iGenVar_cli_test.cpp b/test/cli/iGenVar_cli_test.cpp index adecd31b..619840c4 100644 --- a/test/cli/iGenVar_cli_test.cpp +++ b/test/cli/iGenVar_cli_test.cpp @@ -82,6 +82,10 @@ const std::string help_page_advanced " -l, --min_var_length (unsigned 64 bit integer)\n" " Specify what should be the minimum length of your SVs to be detected\n" " (default 30 bp). Default: 30.\n" + " -s, --sample_size (unsigned 64 bit integer)\n" + " This value is used in sampling positions for average read depth,\n" + " insert size, and read length. Setting to 0 will skip sampling.\n" + " Default: 100000.\n" }; std::string expected_res diff --git a/test/data/datasources.cmake b/test/data/datasources.cmake index 7decf550..5aafa720 100644 --- a/test/data/datasources.cmake +++ b/test/data/datasources.cmake @@ -9,7 +9,7 @@ include (cmake/app_datasources.cmake) # head -20999 simulated.minimap2.hg19.coordsorted.sam | tail -4 > simulated.minimap2.hg19.coordsorted_cutoff.sam declare_datasource (FILE simulated.minimap2.hg19.coordsorted_cutoff.sam URL ${CMAKE_SOURCE_DIR}/test/data/simulated.minimap2.hg19.coordsorted_cutoff.sam - URL_HASH SHA256=af188e36dfab9be4d351985bca8d9b103d64cb859f422cb2c7671b220ccc9f7a) + URL_HASH SHA256=2423106a8d7d72508ab9009a07f2dc814fa3b01b3983d5975339004fd4b1f9df) # copies file to /data/detect_breakends_shorted.vcf declare_datasource (FILE detect_breakends_shorted.vcf diff --git a/test/data/simulated.minimap2.hg19.coordsorted_cutoff.sam b/test/data/simulated.minimap2.hg19.coordsorted_cutoff.sam index cda6face..42f292c3 100644 --- a/test/data/simulated.minimap2.hg19.coordsorted_cutoff.sam +++ b/test/data/simulated.minimap2.hg19.coordsorted_cutoff.sam @@ -1,3 +1,5 @@ +@HD VN:1.6 +@SQ SN:chr21 LN:46709983 m2257/8161/CCS 16 chr21 41970369 60 17M1I68M1I24M1D13M1I107M2I16M1I259M1I2M1I33M1I115M1I54M1I5M2D21M1I76M1I5M1I129M1I23M1I9M1I37M1I28M1I5M1I32M1D52M1I47M1I50M1I84M1I3M1D45M1I89M1I41M1I23M1I10M2I43M1I39M1I83M1I32M1I13M1I46M1I18M1I68M1I22M2I6M1D11M1I1M1I3M1D23M1I2M1I17M1I8M1I43M1I41M1I17M1I38M1I21M1I30M1D10M1I15M1I40M1I4M1D13M1I10M1681I7M1I1M1I3M2I14M1I20M1I12M1D40M1I12M1D4M1I2M1I10M1I31M1I5M1I3M1I22M1I13M1I23M1I7M1I5M1D47M2I2M1D9M1I5M1I13M1I6M1I4M1I1M1I49M1I12M1I8M1I5M1I9M1I48M1I6M1I17M1I4M1I6M1I8M1I6M2I18M1I14M1I16M1I9M1I4M1I30M1I20M1I11M1I29M1I3M1I16M1I11M1D12M1I59M1I2M1I6M1I51M1I13M1I10M1I3M1I21M1I28M1I18M1I1M1I23M1D10M1D3M1I20M1I2M1D31M1I9M1I8M1D20M2I49M1I8M1I2M1I41M1I14M1I28M1I47M1D29M1I2M1I11M1I14M1D10M1D34M1I11M1I16M1I11M1I45M1I5M1D13M1I19M1I21M1I5M1I2M1I59M1I21M1I20M1D4M1I3M1I4M1I24M1I10M1I4M1I26M1D8M1D7M1D13M1I18M1I14M1I18M1I24M1I21M1I22M2I26M1I11M1I30M1I3M1I5M1D37M1I8M1I14M1I49M1I27M1D9M1I17M1D7M1I6M1I3M2D39M1D16M1I4M1I40M1I5M1I54M1D41M1I25M1I5M1I16M1I135M1I2M1D80M1D34M1I24M2I6M1I78M1I4M1D44M2I1M1I10M1D6M1I29M1I2M1I6M1I10M1D14M1I9M1I44M1I8M1I3M1I19M1I26M1I18M1D15M2I16M1D5M1I3M1I11M1I4M1I21M1I12M1I26M1I15M1D19M1I15M1I15M1I44M1D10M1I11M1I2M1I13M1I16M1I4M1I11M1I24M1I6M1I11M1D15M1I3M2I1M1I43M1I4M1I17M1I15M1I4M1I11M1I7M1I4M1D3M1D6M1I6M1I42M1I12M1I17M1I51M1I20M1I1M1I24M1I14M1I5M1I8M1I12M1I39M1D17M1I1M1I10M1I5M1D54M1I32M1D15M2I3M1D5M1I29M1I9M1I20M1I13M1I4M1I45M1I44M1I5M1I7M1I23M1I16M1I2M2I10M2S * 0 0 TGATTGCTGGATCCTATAGGTAAGACTATGTTTAGTTTGTAAGAAACTGCCAAACTGTTTTTTTTTTTTTTTTTTTTTTTTTGAGATTGGAGTCTCACTCTGTCACCCAGGTACAGTGCAGTGGACACGATCTTGGCTCACTGAAACCTCCGCCTCCTGGGTTCAAGTGATTCTCCTGCCTCAGCCTCCCAAGTAGCTGGGATTACAGGTGCACGCCATGACACCCGGCTAAGTTTTTTGTATTTTTAGTGAGAGACAGGTTGTCATCATGTTGGCCAGGATGGACTTGATCTCTTGACCTTGTGATCTGCCCACCTCGGCCTCCCAAAGTGCTGGGATTACAGGCGTGAGCCACCGCGACTGGCCCAAACTGTCTTTTAAAGTGACTCTGCCGTCTCTGTTACTCCACATCTCATCAGCATTTGATATTGTAAGTTTTTTGAATTTTAGCCAGTTTAGTAGCTATAGAGTAGTATCTCGTTATTGTTTTAATTTAAAATTATTTAACGATCACTATAATGCTGATCATCTTTTCATATCTTTATCTCTCCACTTCTGTGTCTTCTTAGGTGAGGTGATCTTCCAGATCTTTTGTTCTTTTTTAACATTGGATCATTTGTTTTCACATTGTTGAGCTTTAAGGCTGCTTTATATATTTGGATAGCAAGTCTTTTTCCTATATGTGTTTCACAAAGATTTTCTCTCCATCTTTGGCTTGATCTTTACGTTCTCAATTGTTCATTCGATGTCTTTTGCTGAGCAGAAGGCTTTTTTTTATTATTATTATTTTACTTTAAGTTCTGGGATACAAGTGCAGAACATCGCAGGCTTTGTTACATAGGTGTACATGTACCATGGTGGTTTGCTGTACCTATCAACCTGTCATCTAGGTTTTAAGCCCCACATGCATTAGCTATTTATCCTAATGCTCTCCCTCCCTTGGCCCCCAACCCCCTGGCCTGGCCTCAGTATGTGTTGTTCTCCCTCCCTGTAGTCCATGTGTTATCATTGTTCTGGGACACAGCTAAAGTCAGTGTTAAGAGGGAAATTTATAGCACTAAAATGCCCCACATCAGAAAGCCAGAAAGATCTCAAATCACACCCTAACATCACAATTAAAAGAGCTAGAGAAGCAAGAGCAAACAAATCCTAAAAGCTAGCAGAAAACAATAAATAACTAAGATCAGAGTAGAACTGACAGGAGATAGAGACACAAAAAATCCTTCAAAAAAATCAATGAATCCAGGAGCCTGGTTTTTTGAAAATGTTAACAAAATAGATAGATTGCTAGATAGACTAGTAAAGAAGAAAAGAGAGAAGAATAGACACAATAATAACCGATAAAGGGGATATCACCACTGACCCCACAGAAATACAAACTACTCATCAGAGAATACTATAAACACCTCTATGCAAATAAACTAGAAAATCTAGAAGAAATGGATAAATTCCTGGACACATACACCCTCCCAAGGATGAAACCAGGAAGAAGTCAAATCTCTGGATAGACCAATACACAAGTTCTGAAATTGAGGCAGTACATTAATAGCCTCTACCAACCAAAAAAAGCCCAGAATCAGATGGATTCACAGCTGCAATTCTACAGAGTTACAAAGAGGAGCTGGTACCATCCTTGCTGAAACTATTTCAAATAATTGAAAAGGAGGGACTCCTCCCTAACTCATTTTATGAGGCCATCATCTCCTGATACCAAAACCTGGGCAGAGATACAACAATAAAAAGAGCAGAAGTCTTTTAAATTGAAGCGAAGTCAAACTATCTTAGCTAGCTTGGGTTGCTATAAGAAAATACCAATCAACTGGGTAATTTATGGAACGACTGACATTTCTCACAGTGCTGGAGACTGGATAGTCCAAGATCAGGGCACCAGCAGATTTGCTCGTCTTGTGTGGAGCTGCTTTCTGGGGTTCAAGATGGTGATTACGTTGTACGTCCTCACATGGTGGAAGAGCACAAAGGCAGCTCTCTGGGGATCTTTATTCATATGGGCACTAATCTCATTCACGAGGGCAGAGACCTCATGATCTTAATCACCTCTCAAAGGCCCCACCTCCTAAGACTATCACATACATATGAATGTTGGGGAAGACACATACATCCAGACCATAGCACCAACTTATTAATATTGTCCTTTCATTCGTCATACTTTTTGCGGTTGTATCTAAAAAGTCAAAACCAACCCAAGGTATACATAGACTTCCTCACTATATTATCTTTTAGTAGTTGTATTGTTTTGCACTTTATGATGTGGTCTCTAATCAATTTTTGAGTTTGAGTGGACCTCAGCAAACTCCCAGTAGAGCTGCAGCAGAGGGGGTCTGACTGTTAGGATGGAAAACTAACAAACAGAAGGCAATAGCATCAACAACAACAAAAAAAAACTGCCACACAAAAACCGCTATCCGAAGATCACCAACATCAAAGATCGAACAGGTAGACAAATGACGAAGAGAGGAAAAAACAGTGCAAAAAAGGCTGAAAAGCCCAGAACCACCTCTTCCCTCCAGAGGATCACAACTCCTCACCAGGAAAGGGAACAAAACTGCACAGAGAAGAGTTTGACCAATGACAGAAGTAGGCTTCAGCAGAATGGGTAATAACTCCTCTGAGCTAAAGGAGCATGTTCCTACCCTAATGCAAGGAAGCTAAAGATACTTGATAAAAGGTTACAGGAACTGCTAACTAAATAAACCAGTTCAGAGAAGAACATAAATGACCTAAATGGAGCTGAAAAACACAGCATGAGAACTTCATGAAGGAATACACAAGGTATCAACAGACAAGTCCATCAGGCAGAAGAAAGGGATATCAGAGATTGAAGATCAACTTAATGAAATAAAGCATGCAAGACGAGATTAGTGAGAAAAAAGAATTAAAAGAAATGAGCAAAGGCCTCAAGGAAATATGGGACTATGTGTAAAAGACCAAGCATACGGTTTGATTGGTGTATGTGAAAATGACAGGGAAAATGGAACCAAGTTGGAAAACACTCTTCAGGATATCATGCAGGAGAACCTCCCAACCTAGCAAGAGAAGCCAACATTCACATTCAGGAAATACAAGAGAACACCACCAAGATACTCCTTGAGAAGTAGCAAACCCCCAAGACACATAATTGTTCAGATTCAGGCAAGGGTGAAAATGAAGGAAAAAATGCTAAGAGAGCCAGAGAGAAAGGTATGGGTTATCCACAAAAGGGCCAGCCATCAGACTAAGAGCAATTCTCTGCAGAAACCCTACAACCAGAAGAGAGAAGGGGCCAATATTCAACATTCTTAAAGAAAAGAATTTTCAACCCAGAATTTCATATCCAGCCAAAACAAAGCTTCGTAAGTGAAGGAGAAATAAATTTCTTTACAGACAAGCAAATGACTGAAGAAGATTTTTGTCACCACCATGCCTGCCTTACAAGATCTCCTGAAGGAAGCACTAAGACATGGGAAGGAAAAAATCCAGTACCAAGCCACTGCTAAACCATACCAAAATGTAGAGACTCAATGCTTAGGATAGGAAACTGCATCAACTAGCAGTCAAAATAACCAGCTAGCATTCATAATGACAGGATCAAATTCAGACCACATACAATTATTAACCTTAAATGTAAATGGGCTAAATGCCGCAATTAAAAGACACATCACTGGCAAATTGGATAAAGAGTCAAGCCCAATCGGTGTGCTGTTATTCAAGGAGACACCACTCTCACGTGCAAGAGACACAGATAGGCTCGAAAATGAATAGGGATGAAGGAAGATTACCAAGCAAATGGAAAGCAAAAAAAAAAAAAGCAGGGGTTGCAAATCCTAGTCTCTGTAAAACACTACTTTAAACCAAGAAAGATCAAAAGAGACAAAGAGGGCTTTAATAATGGTAATAGGGGGATTAATTCAACAAGAAGAGTTAACTATCCTAAATATATATGCTGCCTAATACAGGCACACCCAGATTCATAAAGCATTTTGTTAGCTTGTTTTGTTAAAGCACGTGATTGTAAAATATTTGCTAAATTTTTGTGTCTAGTTCTATCACACATGGATTTTCAGGTGTTCCAGCATTATTTAGTTGAAAAGACTTCCTCTTCCCCCATTGAATTTGTCTTTGTTCCTTTGTCAATATTCAGTTGTAATATGATTGTGTGTGCATCTGTTTCTGAGCCACTCATTCTGTTCCAATTGACCTACTTGCCTGTTCTTCACACCAACAACCAGAGCTAGTGATTACCAAAAGTCTTGAAGTCAGGTAATGTCAGTTCTCTTCCAGGTTCTTCTCACAGTAATCTCCTTGCCTAGTTCTAGTGTCTGTGTTAGCTTTCCATGTAAGCTTTAGAATCAGTTTTTCAATGTGTAAAAATTGGGCTTGCTGAGATTTTTGATTGGGGATGGACTCTGAATTCTATAAATCAAGTTGAGAAAAACTGACATCGTAACAAAATTGAGTCCATGTCCATTGACAATAGAATCTCTCTTCCATTTTCTATTATATCTTCGTCCATCATTCATTAGAGTTTTGTCGTTTTTCTTCATGAAACATACTGTCCAGATGTATGTTAGATTATATATGTTAAGTGGGTTTTTTGGTGCTAATGGAAATCTATATTGGTTTTTCATTTCATAATTCCAATTGTTTTATGTCTGGTATATTGGAAAGCAATTGTACTGTTCATATATTAGTCTTAGTATCCTGCAACCTCTATGACTACCTTATTAATTCCAGAAATTTTTTTGTTAACTGTTTCAGATTTTCAAATAGACTACCATGTTCGATCTGTAGAACAAGGACCATTTTATTTATTTCTTCTCAATCTGTATATATTTTATTTCTCTTTTCTTGCCTTACCTGGATTAGCCTCTAGACTTCCAGCACAATGTTGAGAAAGGAGTGATGAAAAGGAAAATCTTTTCACTTGTTCCAATCTTAGATGCGGAAGAAACTTATTTTCTCATTATAAGTATGATTCTGGTTATAGGTTTTATGTAGATATCTTTTTCCAAGATGAGGAAGTTCCCCTCTATTCCCTAGTTTGTCTGAGACATTTTATCATAAATAAGTGCCCGGGTTTTACCAAATTCTATTTTGCATCTGTTGTTATGATCATATGATTTTCCATTCTTCAAGTCCTATTATAATGCGATAGGTTACATTTCCTCATTTTCTAATAGTTAAAGCATCCTTACAATACCTGGAATAAATCCAAATGCATCAATGGTACATTATTCTTTTTATACATTGTTGGATTTAATCTGCTAGTATTTATTGAGTAATTTTGTATCTAAATTCACTTGAGAGCTATTGATTCTGTAGTTTTCTTTGATGTCTTTTCTGTTTTGGTTTTAGGATAATACTGGCCTCATCGAAATGAGTTAGCCAAATGTTCCCTCTGCTTTCTACTTTTTGTTTTTCTTTTTTTTTTTTTTGAGATTGGGTCTCACTCTGTCACTTAAGTCTGAGTGCAGTGGCTAAGATCTCTGCTCACTACAACCCTCTGTCTCCCAGGCTGAAGTTAATTACTACCCACCTCAGCCTCCCAAGTAGCTGGGACTATAGCCATGTGCCTCCATGCTCAGATAATCTATTTGTATTTTTTGTAGAGATCAGGGTTTGTCACGTTGCCCGGCTGGGTGCTGGAAAATCCTGAGCTCAAGCAATCCGCACCGCTGCGGCTTTCTACAAAGTGCTGAGATCATAGTCATGAGCATAACGTAGCTCAGTTTTATTTTTAGGAAAGAGATTGCAAAGAATTCGGTATCATTTTTTTATCTTTAAGAGTTGGTAGATTTTCACTAGTGCAAACACTGGTGCGCTAGTGCTATCTTTTTTGGAAGGGTTATTAATTATTGAGTCAATCATTCTTTCATTGATATAGGCCCATCCACGAATATCTGTTATCTCCTTATACAAATCTTGATAGTCTGTGTACCTATCAAGAATTAGTTTATTTCATCCAAGTATTCAAATTTATGGCGTACAGGGCTTTTTCATAACCCTATTTATCCCTTTAATTATTATCCTTTTAACATCTATGAAAGCATTCCTGATGGAGTTCTCTTTTATTTCCAATATTAATATTCATGTCTTTCTTTTTGCTTGATAGCAGGTTAGGAGGTTATACATTTTATTGACATTTTTAAAGAACCAGATTTTAGTTTTGTGATTTTCTCTACTGGTTTTCGCTGGTCTCAATGTAGCTGATTTCTGGTCTAAGTTTTATTATTTTTCTTTTCTTCTGCTTGTTTTAGGCTGATATTGCACTTCTTACTCCAGTTTTCTAAGGGGAAGCTTCGACTATTGATTTCAAATCTTTTTTCTTTTCTATATCTATGCATTCAATGTTATAAGTTGTTCTGATGAAGCAGTGATTTCATTTGCATCCCACATTTTGATAGGTTATATTTTCATTTAGTTACAAATAATTTAAATTTCCCTTGAGATTTCTGCTTTGACTTATGTGTTATTTGGAAGTGTATTTTTATTCTCCAAATATTTAGAGATTTGCAGCTAGTTTTTTGTTATTAAATTTAATTCCATTGAAGTTTAATTTCATTATAATCAAGACCATATTTTATATGATGTCTATGATTTTATTTTCTTGAGATGTCTTTTATAGTCCAGAATATGAGTCTTGGTGAATAGACCATGTGAATATTTGTAAAAGAATGTGCATTCTGCTGTTGTTCATGAAGTATTCTATGAATGTAAATTAGATTCAGTTGGTTGATGGCATTGTTACCATAGTCAACTATAACCTTACTGATTTTCTGTCTGCTGGCTCTGTTAAACCATACTGACAGAGGGTGCATGAAGCCCCTAACTATAACAGCAAGTTTTATCTTAATTCGTGCTTATGCATTTATTTTACTGTCTATTTATTCCTTACGAAACTCCTCCTTTTGTTATGTAAATCTGAATTTCTGACCCAATGCTATTGTTCTTTCTTTCTGAATACATTCTCTTTAACCTTTCTTAAATGTTAAAGGCAAGGTTTACTGGTAACAAATCCCTCAATTTTTTTTCGTGTCTGACAAAAATCTTATTGTCTCCCTTCACTTTTGGAAGAAAACATTATGTAGACAGAGAATTTCTATCTTGGTGGGGTTTTTTGTTGCTGTTGTTACTGTATGTTGTTGTTCAATATGTACATTTCACTCCACTCGTCTTACTTGTCTAGCTTTGTGGCAAAAAGTCTCAGTATAATTCTTGCCTTTGTCTCTCTACAGGTAAGGTGTTTTTTCTGGCTTGTCTTCAATATTGTACTCCTTTTTCTTTAGACTTTCTGCAGTTTGGATAATGACCATGCCTAGGTGGTAGACTTTTTTGCATTTCTTTTTCTTTTTTGTTTTTTTTTCTAGCCTCTTGGATATGTTGGATTCTTAGTGTCTATCATTAATTTTGGAAAGCTGTCAGCTATTGTTATTTTCTAAATTTATCTTCTGCTCGCCTTTTCTCCTTTTTGCTGGCGATTCTCATTCATGCATGTAGTTAACCTTTGTATATTGTCGCCAGAGTACTTGAACAGTCTGTTTCTATTTTTCCATTCATTTCTTTTTCTCTTTGACATTTTGGTTTAAAAAGTTTTCTCCTGACGTATCTTCAAGCTCATTGATTCTTTCTTGTTTATATCTAGTTCTATTAATGAGCTCATCAAGAGGGCACTTTTTATTTCTGCTACATTAGCTTTTGAGTTCTAGGCACCTTCCTTTAATATTCTTCCTTAGAAGTTTCCATTTCTGTGCTGACATTATCCATCTGTTCTTGATGTCATCTACTTTTCCGTGATTAGAACCCTTAATTTGTTAATCATAGTTAAATTCCATATCCAATAATTCCAAAATCACTGCCATATATCGAGTCTGGTTCTGATGCTTATTCTGTCTCTTCTATGTATTTTTTCTTTGGCCTTCATCTTCTTGTCATTATTTGTTGAGAGTCAGACGACAGTATTTGGGGTAATGAAAGCTAATATACATTAGACTTTTAGTCATGAAGGTTTTACATGAATCTGGGCTGTATTTAATGTTTGATGTAGCTGCCAGGTGACACAGCCTTTAGTTTCTTTATCTCACTTTTCTCTCTTGGTTGTACTTTGTGCTTTCCCTAGGAACTCATTCTTGAGATAGAGGCTGCATCTTGGCAAAGCTCTCTTACT * NM:i:2019 ms:i:7104 AS:i:8765 nn:i:0 tp:A:P cm:i:311 s1:i:3786 s2:i:901 de:f:0.0496 MD:Z:109^C558A45^TA17A1T345^G236^A486A91^T15^C144T95^A13C55^G80^T4T47^A125^A49^T213T18A12T55A41T0G6C45^C89A157^T10^G25^T48^T36T172^T56^G10^A77T44^G152T7^A75^A8^C7^C89G0A62T51^A135^A26^C16^CA39^C119^T224^C3A76^A126T19^G55^G53^T141^T31^T97^C29G32A0G29^C108^G70A53^T3^T257^C0T32^C86^A18^T232 rl:i:945 m41327/11677/CCS 0 chr21 41972616 60 1574S9M1I4M1I2M1D16M2I5M1I5M1I4M1I6M1I2M1I3M1I6M1I2M1I9M1I18M1I3M1I1M1I6M1I8M1I16M1I11M1D6M1I17M1D9M1I4M1I25M1D6M1D4M1I15M2I4M1I23M1I4M2I15M1I4M1D3M1I11M2I9M1I4M1I7M1I17M2I2M1D11M1I32M1I5M1I3M2I15M1I13M1I2M1D20M1I6M1I9M1D5M1I9M1D10M1I5M1I3M1I1M1I9M1I19M1I5M1I63M2I5M2I4M1D5M2I20M1I2M1I19M1I10M1I5M1I14M1D3M1I6M1I17M1I9M1I34M1I13M2I1M2I3M1I9M1D6M2I16M2I5M1I2M1I2M1D2M1I1M1I10M1I6M1I2M1I8M1I9M1I12M1I7M2I45M1D11M1D6M1D4M3I20M2I4M1I3M1I2M1I26M2D21M1I13M1I8M1D49M1I20M1I8M1I15M1D8M1I2M2I5M1D8M1D5M1I22M1I8M2I3M1I8M1I3M1I9M4I15M1I6M1I30M2I7M1I2M1I24M2I18M2I2M2I16M1I5M1I1M1I2M1I7M1I31M1I6M2I2M1I24M1I3M1I6M1I6M1D5M1I9M1I12M1D3M1I11M1I7M1D15M1D5M1I4M1D11M1I3M1D5M1I6M1I8M1I5M1I4M1I7M1I2M1I14M2I5M1I8M1I8M1I16M1D2M1D17M1I3M1I2M1I12M1I18M1D3M1D17M1I10M1I5M1I28M1D16M1I14M1I5M1D35M1I15M5I12M2I11M1I4M1I4M1I7M1I6M1I10M1I15M1I5M1D16M1I9M1I11M1I8M1D17M1D10M1I8M1I5M1I1M1I8M1I5M1I2M1I6M2I12M1I18M1I3M1I3M1I7M1I7M1I11M1I4M1I5M1I12M1I9M1I3M1I20M1I3M1I6M1I13M1I2M1I5M1I4M1I15M1I9M2I12M1I5M1I3M1I2M1I14M1I7M1I9M1D6M1D18M1I7M1I5M2I6M1I2M1I2M1I5M3I1M1I39M1I26M1I5M1I5M1I27M1I13M1I8M1I6M1I1M1I6M1I16M1I7M1I17M1I2M4I17M1D9M2I4M1I1M1I36M1I12M1D6M1I27M1I3M1I5M1D9M1I4M1I3M3I14M1I14M1I13M1I15M1I14M1D5M1I3M1I4M1I2M1I9M1I13M1D15M2I29M1I7M1I7M2I22M1I15M1I16M1I4M1I1M1I20M1D5M1I1M2I7M1I2M1D7M2I15M1D7M2D6M1D8M1D23M1I7M1I5M1I8M1I2M1I15M4I4M1I5M2D5M1I30M1I6M1I16M1I5M1D15M1I26M1I6M1I7M1D11M1D5M1I23M1I15M1I24M1D6M1D2M1D6M1I6M1D15M1I2M1I11M1I3M1D7M1I20M1I21M1I5M1I10M1I12M1I4M1I6M1I17M1I5M1D35M1I2M1I17M1I4M1I13M1I11M3I10M2I3M1I3M2D4M1I10M1D11M1I10M1I2M1I10M2D7M1I6M2I1M1I3M1I7M1I17M1D5M1I6M1I8M1I2M2I4M1I4M1D7M1I2M1I3M1I1M1I19M1I2M1D32M2I8M2I29M1D29M1I19M1D3M1I19M1I15M1I8M1D4M1I11M2I25M1I10M1I15M1I5M1I2M1I4M1I7M2I9M1I2M1I4M1I13M1D7M2I19M1D9M1I11M2I4M1D20M1I14M1I6M1D10M1I3M1D28M2I6M2I1M1D14M1D6M2I7M1I20M1D11M2I2M1I2M1I3M1I1M1I2M1I7M1I5M1D3M1D16M1D13M1D3M1D8M1D6M1I3M1I7M1D12M1I5M1I11M1I6M1I2M1I6M1I4M1I7M1I3M1I8M2I3M2I3M1I11M1I1M1D3M1D3M1I15M1D16M1I3M1I5M1D45M1I4M1D6M1I9M1I8M1I3M1I14M1D20M1I22M1D4M1I13M1I20M1I31M1D5M1I4M1I3M2I2M1I22M1I3M2I3M2I19M1I22M1I14M2I27M1I4M2I7M1I26M2I58M1I1M1I13M1I3M1I54M1D4M1I1M1I2M2I17M2I6M1I2M1I12M1I3M1I5M1I16M2I5M1D16M2I8M1D3M1D10M1I77M1I15M1D5M1I18M1D6M1I3M3I21M1I32M1D5M1I37M1I2M1I3M2I35M1D5M1D12M1I16M1I24M1D8M1I2M1I6M1I5M1I7M1D47M1D3M1D11M3I12M1I4M1I15M1I24M1I3M1I7M1I9M1I4M1D21M1D3M1I14M1I26M1I20M1D29M1I1M1I31M1I8M1D13M1I30M1I9M1D16M1I13M2I10M1D2M1I17M1I5M1I8M1I4M1I13M1I36M1D28M1I7M2I24M2I32M1I35M2I13M2D8M1I46M1D8M1D3M1D12M1I3M1I10M1I34M1I4M1I18M1I3M1I7M1I15M1D7M1I4M1I12M1I12M1I27M1I9M2I40M1D5M1I26M1I32M1D4M1I6M1I5M1I17M1D8M1D46M1I1M1I19M1I10M1I13M1I15M1I10M1I10M1D2M2I3M1I5M1I1M1I17M1D12M1I4M1I12M1D6M1I24M1I8M1D15M1I18M1D31M1D2M1D30M1I13M1I5M1I6M1I31M1I15M1I5M1I14M1I25M1I49M1I15M1I3M1I3M1D37M1I6M1I2M1I17M2I22M1I42M1I45M1I8M2I17M1I16M1I39M1D33M1I16M1I6M1D14M1I5M1I6M1D11M1I3M1I2M1I2M1I5M1I16M1I10M1I8M1I11M1I25M1I4M1I21M1D30M1I2M1I10M1D24M1I18M1I12M1I10M1I11M1I20M1I36M1I4M1I14M1I68M1I11M1I55M1I6M1I15M1I5M2I17M1I3M1I19M1I2M1I6M1I61M1D18M1I21M1I7M1D69M1I7M2I15M1D10M1I6M1I5M1I16M1I19M1I8M1I28M1D17M1I7M1I3M1D1M1D5M1I15M3D4M1I3M1I5M1I12M1I10M1D33M1I13M1I8M1I6M1I29M1I47M1D9M1I7M1D34M1D11M1I17M1I11M1I8M1I48M1I9M1I9M2I7M1D20M1D5M1I10M1I37M1I24M1I10M1I24M1I5M2I7M1I20M1I4M2I2M1I13M1I12M1I43M1I34M1I1M1I17M1I11M1I36M3I9M1I40M1I41M1I40M1I4M1I3M1I9M1D5M2I23M1D23M3I16M1I3M1I36M1I8M1I15M1I18M1I7M1I5M1I28M2I14M1I9M2I16M1I11M1I36M1I70M1I5M1I35M2I3M1I30M1I3M1D7M1I42M1I26M1I13M1I4M1I19M1I7M1I3M1I17M1I7M1I7M3I6M1I69M1I6M1I42M1D5M1D18M1D7M1I34M * 0 0 GTCAAAAAGGGCTGAACAAGCCACGAAACACTTCTTCTCCTCCAGCAGGAGTTCACAACTCCCACACCAGAAAAGGTGAACAAAATCACAGAGGGAATGAGTTTGACCAACTGCGAAGTAGGCTATCAGCAGATGAGGTATAACTCCTCTTGAGCTAAAGAGAGATTGTTAACCCTAATGCAAGGAAGCTAAAGACCTTGATAAAAAGGTTACAGGAACTGCTATAGATATAACCAGTTCAGAGAAGAACATAAATGTACCTAATGGAGCTGAAAAACACAGCATGAGAACTTTAATGCAAAGAATACACAAGTATCATACAGATCAAGTGACCATCAGGAGCCAAGAAGGTATCAGGACTTGAAGATCAACTTCAATGAAAAATACATGAAGACAGAATTAGAGAAAGAAAGAGACATTAAAAGAAATGAGCAATAGCCTCAAGGAATAAGGGACTATGATGAAAAGACCAAGCTCGTTTGAATTGGGTGACTGATGAAAATGACAGGGGAAAATGAACCCAAGTTGTGCCAAAACACTACTTCAGGATATTCATGTCAGGAGAGACCTTCCCCAAGCTACAAGGAAGCCAACATTCAAATTCAGGAATACAGAGAACATCACACCGGAAGATACTCCTTGGAGAAGAAGCAAAAACCGCCACAAGCACCATAATTTGTCGGTAATTCCAGCAAGGGTGAAAATGCAGGAAAAAAAGTGCTAAAGCAGCCAGTAGACGAAAAGTCTTGGGTTATCCACAATAGGGCAGCCCATCAGACGTAAGAGGAGACCTCTCCTGCAGAAATCCCTTACAACTCAGAGAAGGAAAGGGGCCGCAATATTCAACATTCTTAAAGAAAGGAAATTTTCAACCCAGAATGTTCAATATCCAGCACGGAAACAAAGCTTCGTAAGTGAAGGAGAAAATAAATTCTTTAAGACAAGCAATGCTTTGAGAGATTTTTGTCACCAACGCGATGCCTGCCTGACAAGGATCTCCTGAAGGAAGGCACTAAACACTGGGAAGGAGAAAACAGTACAGCTCACTGCAGAAACCATAACCAAAATGTAGAGAACCTACTCAAATGCTAGGAAGAAACTGCCGATCAATAGCAGGTCAAAATAACCAGGTAGCATCATAAGTGACAGGATCAATATCAACACATAACAATTATTATACCTTATAACTTGTAAATGGGCTAAGATGCGGCAATTAAAAACAACACTACTGGCAAATTGGATAAAAGAGTCAAGACCCCGTCGGTGTTGCTGTATTCAGGAGACCCATCTCACAGTGCAAAGACACAGATACGGTCAAAATACATAGGATGAAGAAGATTTACCCAGAAAGGCAAAATGGATATAGCAAAATAAAAAAAAATGCCAGGGGTTGCAATCCGTGTACTCTTGATAAAATATCACTTTCAAACCAAAGAAAGGAACAATAGAGACAAAGCAGGGCCTTTACATAATGGGCCTTAAATGGGATTAATTCCATACCCAAGAAGAGGTTAACCATCCTAAAATATTATACTGTGTCTCCAATACAGGCACACCCTTGATTCATAAAGCTATTTTTGTTGCATGTCTGTAAAGAACGTATTGTATCAAATACTTTGCTTAAATTTTTGTCGTCCTATGTTTCTGTTGCACACATGGCATTTTCAGGTGTTCCAGCGATTGATTTTGTTAGAAAAGACTTATCCTTTCCCCATTGTAATTGTCTTTGTCCTTTCGTCAATATTCAGTTGAAATATTTGTGCTGCATTCTGTTTCTGAGACCCCATTCTGTTCATTGACTACTTTGCCTGTTCTTCCACCCCAAGCACCAGAAGCTAGTGATTACCAATAAGTGCCTTGAAGTCAGGTAAGTGTCGTTCCTCTCATGGTTCACTTCTACAGATATCATCCTTGCGCTATTCTAGGTCTTTTATCGCTTCCATGTAAGCCTTTAGAATCAGTTTTTCGATGTTAAAAAATTGGGCTTCGCTTGGAGTATTTGATTGGGTATGACTCTGAATCATAAAATCACGTTGAGAAAAACTAGACATCTGAAACAAAATGAGTCCCTGTCCATACAATAGAATGCTCTCATCCTATTTTCTATATAATCTTGTCTATTTCATTAGAAGTTTATGTAGTTTTCTTCATGTAAATACTGTCCAGATGTAGTTAGATTAATATTTAAGTGGGTTTTTTTGGGTGCAGTAATGAAATTATATATTGGTTATTCATTTCAGAATTTCGAATTGTTTATTGCAGAGTACATTGGACAAGCATATTGACTTTCATATTTATGTCTTGTTATCCTGCAACCCTCTATTGACTACTTAATTAATTCCCAAAAGTTTTTTGTTAACTATTTCAGGATTTTCAAATAGTTAGGCTACCCAAGTCATTGTGAAATCAAGGACAATTTTATTTATATTTACTATCCATAGTCTGTATATACTTTTATGTTCCCTTTTCTAGGCCTTTCTGGGATTAGGTCAGAACTTCCAAGGCACAATGTTGAAAAGGAGTGATGAAAAGAACAATCTTTTACTTGTCCAATCTTAGGGGAAGAACTTCATATTTTCTCATCGATAAGTAACTGATGGTCGTGATTATAGGTTTTATGTAGATTCTTTTTCAAGATGAGGAAGTTCCCCTCTTATTCCTATTCTGCCCGACACATTTTATCATAAAATAGTGCCGGTTTTACCAAATTCTTTTTTCCATCTGTGTGTTATGATCATATGATTTTCCATTCTTCTAGCCTATTATAATGCATAGGTTACCAGTTTTCCCATTTTCTATGTTCAAAGCATCCTTCAATACCTGGAGATAGATCCCAAAACTGCATCATGGGTCACATTATTCCGGTTTTTTATACATTGTTGGGATTTCAATCTGCTAGTATTTTATTCAGTAATTTTGGTTATCTAACATCTCATTAGAGCTATTGATCTGTAGTGCTTTCTGGTGATGTCTTTAGCTCTGTGTTTTGGTTTTAGGATTAATATCATGTGCCTCATCCAGATGAGTTAGCAAATGTTCCCTCGTCTTCGTACTTTTATTATTTTTTTTTTTTTTTTTTTTGAGATTTGTGGTCTCTACTCTGCACTTAAGTCTGGAGTTGCAGTGGCTAGTCTGCTGCTCACTACAAACCTCTTCTCCCAGGCTGAAGAATTCGTCCCCCTCAGCCTCCGCAATAGCTCGGGACTCATAACCATGGTCCCTTCCACTGCTCAGAATCAATTATTTGTATTTTATTTGTGAGAGACAGAGGTTTGTCAATGTTGCCCAGGCTGGCTGAACTCCTGAGCCAAAGCCAATTCGCGCCCGCTGCGGACTTCTCAAAGTGCTGAGACATGTCATGAGACATAGCGTGCAGCTCAGCTCTTAATATTTTAGGAAGAGGTTGCAACGAACTGGTTCATTTTTTTTCTTTAGAGAGTTCGTAGATTCTCACTGTGCAAACACTGGTGCCTGACGCTATCTTCTTTGGAAAGGTAATTAATTTTGAACGTGAGTCAATTTCAGTTTCATTGATATTAGGCCCCATTTCAGAATTATCTGTTTTCTCCTTATTACAAATCTTGATAGATCTGTTCCTTCAAGAAATTAGTTTTATTTCATTCCAAGTATTCAAAATTTATGCTACAGGGCTTTTCATACCCTATTATCCCCTTTAACTTATTGATTCCTTTTAGACATCTTAGTGAAAGGCCATTCCTGATGATGTTCTCTTTTTATTCCAAGTATGTAAGTAATTCACTGTCTTCATTTTTGCTTGATTAGCTCAGGTATAGAGGTTTACCGAATTTTATTCGACGATTTTTAAAGAACCAGATTTGTAGTTTTTGCTTGATTTTTCCTACATGTGTTTCAATGGTTCTCAATGGTGCTGAGTTTCTGGTCCATAGGTTTTATTATTTTTCATTTCCTCTCTGCTTGTTTCAGCGCTGATAGTTGCACTTCTACTCCGTTTTCTAAGGTGGAAGCTTTCGACTAATTGATCTTTCAAGATCCTCTTTTTGACCGTTTTCTAATCTATGCATTCAATGTTATAAGTTTTCTGTGTAAGCAGTGATTTCATTGCATCCCACAGTTTTGAATAGGATTATATTCTCATGTAGTTACAAATAATGTTAAATGTCCCTTGGAGATTTCATGCTTTAGCACTTATGGTGTTATTTGGAAGTGTTATTTTTGATTCTCCAAATATTTAGAAGAAATATCTGCAGCTGTCTTTAGTTATTAAAACTTTACACTTCATTTGAAGTTCAATTTCATTGTAATCAAGACCACTATTTTATATGAGTCTATGGATTTAATTTTTCTTGAGATGTCTTTTCATACGTCCAAATATAGTCGCTGGCTGAACTATAGACCATGTGAATTTTGTAAAGAATGTGGCAATCTGCTGTTGGTTCATGAAGTATTCTTATGAATGTAAATTGATTCTTGTGTGTTCTGTATGGCATTGCTTACATAGGTCAATATAACCTTACTGATAATTTCTGTCTGCTGGCTCTGTTAACTACTGTACAGAGGCGGTGCTGATAAGCCCCTAATCATAACAGCAATGTTTTTCTAATTCTGCCTTATGCATTTTATTTGTACTCGTTCATTTATTCCTAAGCAACTCTCCTATCGTTGTTATAGTAATCTGATTATTTCTGACCCATGCATTTTCTTTTCTGATACATTCTTTAACCTTTCTTAAATGTTAAAGGGCAGGTATTACTCGGTAACAATATGTCCCTCAATTGTTTTGTTAGTGTCCTGACAAATCTTTTATTTCTCCTTCACTTTTGAAGAAACATCTATGTATGACAGATAATTCTTTCGTTGGTGGTTTTTTGTTGCTGTTTGTTGCTGTTGTTGTTGTTCAATACGTGTACAGTTTCACTCACTCTCTTACTGTGTAAGATTGTGGCAAAAAGTCCAGTAGTAATTCTTAGTTTTGCTCTCTCTACAGGTAAGGTGTTTTTTCCTGCTGTTCAATAATTGTCCTTTTTATTTAGATTATCTTGCAGTTTGGACTATACATGCCGTAGGTGTAGACTTTTTTGCATTTTCTTTTTTTTTTTTTTTTTATTTCTAGAGCCTCTTGAGATATGTGGTTTAAGTGTTCTATCTAATAATTTTGAAGAGCTGGTCAGTATTGTGATTTCAAATATATCTTCTGCTCCCTTTTACTACCTTTTTCTGGGATTCTCCATTCCTGCATGTGTTATGACCTTTTGTAACATTTGTCCCAGACGGTATCTTACAGACCTGTTTCTATTTTCCATTCATTTTTTTTTCTCCTATTGCATTTTGTTAAAAACGTTTCTGACGCTGTACGAATCTTTCAAGCTCATTGATTCTTCTTGGTTTATCATCTAGTCGTAATTTAATTGAGTCATCAAGAGTGCATCATTTTTATTTCTGCTACATTCGCTTTGAGTTCTAGCACCTCCTTTAAGATCTTCCTATTAGAGTTCGTCCTTTCCTGTGCTGACATTATCCATCTGTCTTGCTTGCCATCTACTTTTCCTATTAGCAACCTAAATTCTGTGAATCTAGTTTAAATTCCATATCCAATAGGTTCCAAAATTACTGTCTATATATAGTCTTGGTTCTGATGCCCTTATTCTGTCTCTGCATATGTATTGTTTTCTTGCCATTTCATTTCTTAGCATTTATTCTGCTTGACGAGTCAGCAACACAGTATATTGGGGTAAATAAGAGCTAATTAATTAGATACTTTTAGCATGAGGTTTTCATGAATCTCGGGCTGTATTTTCAATGTTGATGTAGCTGCAGGTGACTACAGCCTTTAGTTTACTTTATTCAATTTTCTACTCTGTTGTCTTTGTGTTTCCCTAGGAACTCAAATTCTTCCGATAGAGGCTGCATCTGCAGCTTTGTCTTACGCTGTAATTCACGGGTATTATTTGGAGCCTTCGGTGTGAATGATAGATTGCTGAGAGTGTTGGGGAGGAAGCTCTCCATATCTTAAGATAAAACTCGCTTGCAGATACCCGAAGCTCCCTGGTTGTGACTTTCACGATGTAGCTTTCTGGTTAGTGTTTAATCTTGTTTAGTTTGGTTTTTTGAGCTTTTTCTTCCGCTCTTCTCTCTATGTGAGGCAAGCAGATCAGAAGGGGGCACCAGGGCTAATTTTCCTTTCTCTCAAAGTTGACAATGTTCTGGTAAGGTACCGTTTCTTGAAGGACATGGTTTTCTTTATGAGAATCGCAAGGGGCTTTGCTTCAGAATCACTTACCTTCTCACTCCAATGCTCAGAGATGGAGCGGGATTTTTCTCTGACCTTCATTGAGCAATCTGATGAGCTGTCCTGGAGGAAAAACTCAAGCAAAGCTGTAGGGGCTCCTAGGAAATATTAACCTCAAGGCCTCCACCAATGACTTAGTCCCGAGAAGTTCATTAGAATTCTTTAACCATTTACGTGTTCCTACTCAATCACTGGCTCCGCAAGCTTCCTGCCCCCTGTAAGGACTGACCTTGGCTCTCATTCTCTGTGTTCGACCCTTGTCTCTACCAGATTGCACACCAGTGGTTTGCCCCTTGTGACCTCAATCCTCTGATGAATCCAGAAAAGTCATAGATTTTCATTCTGTTCAGCTCTATTTTTCATGTTATGGAGCCCCAGGAGTAATCACTTCCAAGCTCTCTGGGTGCCCTCGCTGAAAATCAAAGTCGCAAGACCTGATTTTCTTACGAAGGGCCGAAAATGAGTACAATCATTTTGGCGTTTCGAGAGAATACAGGGTCTCTATTAGGCAGCCCTGAAGTGTCCTATTGCGTAATATGAAGAGCCATAAACAAAAATGCGAATAAGTGAACATGGCTGCGTCCCAATACAGCTTCATTTACAAAAAAAGGGGTTGGGCCACATTTGGCCATGTGAGCCCCATTTTCTGATCCCCTGCTCTATACAGGCTCAGAATGGAATAGAGCCCATCGGGGCCCAAGACACTTAGTAAACATGAGCAAAGGTGAATGGGTTTTAACAAGAAAATATGTTCTTCACTTATATATGTGCTTAATGCACTATATCTTGTTTTTTTTTTTTTTTTTTTTTTTTTTTTAATACAGCTCACGTATTACCCAGGTCTGGAATGCAGAGGTGACAGTCTCAGCTCACTGCAACCTCCCCTCCCAGAGTGTCAAGCTTATTCATTCTGCCCAGCCTCCCGAGTAGCTGGGATTACAGGCATGTGCCACCACGCCCAGTAATTTTGTATTTTTTATAGTAGAGACAGTGGTTCTCACCATGCTGGCCAGGGCTGGAGTCAAACTCCTGAACTCTAGGTCGAACCACCCCACCCTGCGCCTCCAAAGTGCTGGGATTACAGAGTGCAGCCACTGTGCCCAGGCCAATGCACATTCTTGTTTCCATAACAAATGTTATATTTCATTTTTTTTAAACAAATACTCTACTTTAGATAGGATATTATGAATATCATAAATCATTTAATGCTTCTGTAATTATAAAAAGAGACAGCCACCACAATTGGCAGCTGGTACTCTGTCTCAGTCGTTGGTGTTGGGCCACGAGGGGGTGAATGGGGACATTGAATTGAGCTGGAGGACAGAGCACCAAAGCCACTGGCTGGAAGCTGCCAGGGTGCGGAATGAGAGCCTCGCATTGCCAGATTTATAATTCTTCCAGAGGAGCTGGAATCCAGATTTTATGTTACGTTTAAAATGGTTTAGTGCGTTGGCAACTATTCAAAATAAATTTATTTAAACATTGCCCCGGTCAACAGTCAAACACGAAATATTTGTCAAGTACAAATCCTCTAGCTGGCTATTACATTGGCTGGCTTGGTTGAATTTTGGCCTCCAACCCAGGCCAGCCCCAACGTTTACCAAACGGGTGCTGGGAGTTAAGAGGTGCCCAAGGCTCATCACGCCCAGGAAGACAGCAGAGTGTGACCTGGGCTGTAGTAGCTGAAGCATGCCCACTTCTGTATGCGATGCTCCACAGCTCCAGAGCAAATGTACCTTTAGCTCTCCACTGGAATGAACATCCTCTCTGCTCTGTTTTTCAGTCTGATATGGTTTGGCTGTGTTCATCACCCCTAAAACTTATCTTGAATTCCCACATATTTTGGGAGAGACCTCTGGGTAGGTAATTGAATCATGGGGGCAGGTCTTTTCCTGTGCTGTTCTTGTGACAATGAATAAGCTCACCAAGATACTTATGGGTTTTATAAAGGGCGATTTCCCTGACAAGCTCTGTCTTTGCCTGCTGCCATCCTTGTAAAACATGACTTGCGTTCCTCCTTGCCTTCCACCACTGTTTGTGAGAGCCTTGCCAGCCACCGTGGAACTCTGAGTGCCATTAAACCATTTTTCCTGTTAACAATATACCCGACGTCTTGGGTATGTCTTTTTAGCAGTGTGAGAAACAACAGTCCCACCATGAGTCGATTTTCTTCCACTGATATTAAGAATTCCTCGTGTTACATTCTTACATACCCTGTTCAGAAATCAACTACTTCACACTCTTTTTTCAAAAAGCACAATTTGTTCAGTTAATCAGCAAAAAGAACTAATAAACACAATCACAATTAGTGCCTAGTGCTGTTGAAATCACAGCAAGATAAATGGTTCCAATCAAGTCAAAGGAATTCTAACAAAGGAACATCATATATAGAGAAATGTGCTTCAAGTTATTTGGATTCCAAATTTCATATGCCTCATTCTCCCTTTAATTCTGATTTTTTTTTTTATTTTTTTTTAGGAAAGTAATATTCATCTTTTCAGGGATTTCATCTTTTTGTAAAGCGATATTTCTTTGCTATTAGGGGAAAAATGAAAATAGAAGCATGATATGGTTTGGCTGGTGTCCCCACCCAAATCTCATCGTGAATTGTACTGCCATAATCTCCCACGTGTCATGGGAGGGGCCCGGTGGGAGGTAATGGAAACATGGGGCGGGTGAGGGGGTTTCCCTCACACTTGTTCTCATGCTACTGAAATAAGTCTCATGAGATCTGATGTTTTTATAAGGGGTTTCCCATTCACTTGGCTCTCATTCTCTCTTGCCTGCCTGCCATGTAAGACGTAGCCTTTACCTTCTGCCATGAGGTGTGAAGGCCTCCCAGCCACGTCGGATACATGTTGAGTCCCATTAAACCTCTTTTATCTTTATAAATTTACCCTGTTCTCGGGTATGCTCTTTATCAGCAGCAAGGTAATGAAGCTAAGTACAAAGCAATTGGTTAGAATTTCCCACATTGATAGGCGCCATTCCATAGGTTGACCACTTTGGAGGAGAAGGAAGTTGGTACAGGGCCAAAGGCTCCTTCTCAGAGGACTGCAGTGAAGAGGTGAGTGAGGTAAGGGGGGAGTATTCAGGAGGGCCCTGGCTTGCCTCTTGGTGTTGCCATCACTCTTGTTCCTGTGATCCAAGAGTATCTTCAAAAGACACGAATTATCACCTCTAAAGAGCACCATAAATGGTAGGAGGAGACTGATTATAAAATATACATAATCCTCAATGAATTCTGCGTGAAGTATTTTTCTAGTGGCAGTGAGAAATTATCCTAAAAACGTGATGACATCACAATCTCACTTCTATTTCATGGCTCACAATTATTTGAAATGCCATAGCAGGTGCTATTTTTCTGGCATTAGGGCATAAAAAATATAATATTTTGAAAATATATTTGAAACTGCAAAGAACAGCAAAATATCCTGTGATTATATGCAAATTAAAATTAGGGAGAGAATGGTTAGGTGTGTTGAGCAGGAGTAACAAAAAGAAATCGTGCCATTGTATACCAGGTAACATTGTATATTTTCTCATTTTCCAAACGGTGTTCTAATGGAAAAAATTAAGAGATACAGCTCAAATAAACTGCTTTTGTCGCCCACGATAAGCTCTGGAAAACGGTTAGCTGAAGTGAAACAGGGTGTGTTGGACAACAGAATATCTCAGCACTTTTAGTTGCTTGTATGCGTTGTAGTTCTTAGCTAGTTCTACATATGAGGTATTTTTCTCTATTTGAACCATACAAAGCCCATATCTCAAATGCAGACTTTCAAACTCTCCCTGAGAAGTATTTAATGGCAACACAATTTGGGTAAACATTGAATTCATATTTAATTGTTGGGAAATTTAAATAAATTCTACATATTTAGAAAGTTAACTTTGCAATACTCTTAAAGGATACAACAATTTGGTATGCAAGTGTACAGGGACCTCTTCTCTTCTTTAAATTTAAGGTTATAATTCCAAGAGGCTACCACTCGTTTTCACCCCTGGGTAATCTTGCAAGTATCTTTCTGGGTAAAGTGTTTTGAATAAGTTGAATTGCTTCAGATTGTTAATTAATTAATTTAATTAATTATTTGGGAAGGGTCTCACTCTGTCACCCGGCTGTGAGTGCAGTGCGCACAATCATGGCTCACTGCAGCCTCGAGCTCCTGGGCCTCAAGCAGTCCTCCATCCTCAGGCCCCCTGCATATGCTTGGACTGCAGGTATGTGACACTCATGAGTCCAGCTATTTTTTTGTATTTTTAGAAAACGATGGCGGCTTTTTGCCATGTTCGCCCAGGCTGGACCTGGAACTCCTGGACTCAAGTCATCTACCTGCCTTGGCCTCACGAAAGTGCTGGTACTACAGGCATGAGCCACCGCGCGCACAGCCTGCTTCAGATTTGTTAAATCACGATTTTTGCCGTTGTTGTTTCCCCAGATGGGTTTAGCAATACGAGTGAAACGTATGCACATCATAAGAAATTTTCAATGCCAGAAGAAGCAATTTCAACAAAATTGTTGCCATAAAATTTGTACTGGAGCAAGATTAGAACACCTGTTTAGTTATCAATTAATACCTTTTCATTACTGTTGTTGCCCTGGCTTTTGTCTTCCCTTTGCCTGGCTTTCTGTGCAAGATACCTGGGCTAGGGTTAATTATACTGATTTGTAGCAAGGAGCGGTCAATTAAGATTTATGCTTGTGGGGAAATGCTTGGCACTCTCTGCTGCTGTTAGAATAACACAAAGGTAGTGATCTTTCTAGATCCTGCAAGTGTGGGAAGGGATGAACTTCTGTGAGAAAAGGCTAAAGATGACTCAGGAAAGGGGAGAGCGCATTGTTAGGACTGATATTGTTGGGGAGAGATGGTAAAAAAGATTTGGGGAAATATCAACTGGAGGGAGAGGTAGGTGGATAAAATAGAAACCCCTTTTAGAATTAGGAGCTAATGGGAGGCCAGCATTTGCTTTGGACAGCAGGACAGGCAGGCCTTGCCTGGGAAAGTCTAGCTCGACAGAGATGGGCAGATGGTGGTGAGACTATCACGGCCTGCTCAAGCTTCTGGTTAGATTCTTATTACTCAGTGCTCTTATATCAAGATGAATTCTGGTGACAGTCCAACAAAGAGGCTCACAGTTCACGTAGGGAGTTAATTTAAAAGAGGCCACAAATCCTGCCCTGCATAGGGAGACGGACTATTTTCTTACGGTTCACCTAAGGATACCCTTTCTTGACAACCAAACACCTTTAAAACAAGCACGGGGCTTGGGACATAGTCCACACGAACATCTCGAAGGATTGCTGGGCTGAGAAACAGCATATTTAAATCTCTCATCAAAAATGAATTCGAAAGCCCAGCGAGACAAATCTACTGCTTTACTGTAGAAAGAAAAAAATGCCCACATTTC * NM:i:1057 ms:i:12862 AS:i:12862 nn:i:0 tp:A:P cm:i:235 s1:i:3069 s2:i:51 de:f:0.0929 SA:Z:chr22,17456973,+,2S1444M126I10105S,60,229; MD:Z:9T5^T5C115^T23^T25C2T9^C6^C69^A53^T29A4G0T18A0T25^T6A20T7^T14^G124^A15T14C13T30^A44A0G3T13G25T5^C31^T29T5C9C40G1A13^T11^A6^A15T0A0T41^TC29G1T3T2G3^T12T0A27G50^G15^T8^A30A67G48T47A22T0G51^T26^A21^G15^T9^A14^G14G6G66^T2^G3A8T0C38^T3^A13A15T14A6A3T4^A22G12^A20T8T10T7A75^G44^G17^A67A0T106C0T11C11T0A16A30T18^T6^A128T4T20T71T14^T17C0A8T9A24^T11T1A27^G9T37T38^A5A30^C68C0T59C1A4^C15^A22^T7^TC6^A8^T55T13^AA47G6A7^G20A33^C11^T7C30C28^C2T0G2^T2^T12^T6C24^G86T8G1A9^C6T91^GA4T9^T33^GT20T20^T29^C34^T24A0T17A2T22^T9T23C0T8T4^A22A9C5C6^G29T31G0T38G1A7^A26^A24^T40^C3C9^T35^A14^T7C17T7^T33^A3^A16^C13^T3^A8^T16^A22G59^G3^G18^T24^G20G1C2C23^A10T29^T3C9A0C27^G10G57^T22T1C2C1A115G144^A35A0T21A14^T24^A3^T102^G23^A11G5A3G40^A77G4^T5^T5G17T28^A28^T47^C3^T89^C21^T33A29^C69^T52^G39^T85^A19C53A65^TG36G0T16^A8^A3^A18A0G86^T111^G63^T32^G8^A124^A28^A28^C38^T18T14^A31^A2^T214^C26A0G223^C55^C25^C118^A12G29^T28T0C387^A46^A9C42C1T36^A92^G27^A1^C20^AAG34^A136^T6T9^C34^T8G111^C20^A148A163A24T123^C28^T19T70G134G1T26A19T1A114^C275^A5^A18^T3C37 rl:i:474 m21263/13017/CCS 0 chr21 41972617 60 1487S8M1D13M1I4M1I19M1I7M1D8M1I11M2D6M1I12M1D4M1I19M1I1M1I18M1I3M1I4M1I6M1I6M1D22M1I4M1I6M3I12M1I17M1I7M1I16M1I6M2I24M1I6M1I2M1I8M1I9M1D12M1I8M1I1M1I7M2I7M1I12M2D19M1I4M2I8M2D6M1D10M1I14M1I5M1D11M1I13M2I12M1I27M1I29M1I8M1I2M1I4M2I5M1D3M2D11M1I26M1I2M1I10M1I1M1I19M1D3M2I12M2I6M1I43M4I3M1D2M2D2M1D6M1D37M1I3M1I20M1I4M1D2M2D4M2I3M1D41M1I1M1I2M1I9M1I1M1I16M2I6M1I10M2I9M1I2M1I16M1I1M1I11M1D15M1I1M1I13M1I3M1D15M1I4M1D11M1I9M1I2M1D21M1D10M1I3M1I2M1I3M2I7M1I16M1D25M1I12M1I4M1I18M1I32M1D3M1D6M1I1M2D11M1I7M1D12M1I5M3I5M1I8M1I2M1I7M1D9M2I27M1I5M1I1M1I1M1I6M1I15M1I2M1I17M1I5M1I8M1D23M2I10M2I2M1I7M1I3M1I3M1D22M1D4M2I6M3I21M1D13M1D4M1I3M1I35M1I21M1D29M2I1M1I1M1I6M1I3M1I11M1D6M1I6M1I12M1D3M1I15M3I5M2I12M1I8M1D6M3I8M1I26M1I6M3I8M1D3M1I5M1I54M1I13M1I15M1I9M2I29M1I16M1I11M1D7M1I3M1I5M1I5M1I13M1I8M1I5M2I6M1I3M1I19M1D3M2I17M1I17M1I14M1I29M1I6M1I16M2I10M1I14M1I3M1I8M1I5M3I14M2I14M1I15M1I10M1I3M2I10M2I12M1I21M1D27M1I1M1I60M1D11M1I24M1D3M1D8M1I2M1I2M1I11M1I11M1I4M1I2M1I6M2I4M1D3M2I7M1D10M1I26M1I4M1I4M1I26M1D14M1I6M3I17M1D4M1I20M1I33M1I19M1I4M1I31M1I6M2I8M1D17M1I7M1I14M1I2M1I16M1I32M1I4M1D12M1I7M1I3M1D3M1I3M2I32M1I24M1I2M1I16M1I4M1I22M1I4M1D2I4M1D3M2I28M2I11M1I2M1I18M2I4M1I4M1I3M1I11M2I8M1D9M4I23M1D22M2D9M2I1M1I6M1I23M1I11M1I10M1D13M1D4M1I5M2I3M1I9M1I3M1I4M1I6M1I12M1I5M1I12M2I9M2I2M1I1M1I20M1I5M1I21M1I16M1I6M1D10M1I7M1I20M1I1M1I1M1I10M1I26M1D2M1D10M1I8M4I15M1I4M1I4M1D30M1I7M2I7M1I5M1D9M1D13M1I13M1I6M1I12M1I9M1I6M1I23M1I3M1I4M1D12M1I10M1I6M2I3M1I8M1D22M2I12M1I5M3I8M2I3M1I4M1I13M1I3M1I5M1I34M1I25M1I2M1I19M2I4M1I3M1D2M1D2M1I5M1I14M1I1M1I19M1I81M1I20M2I4M1D5M1I7M1I9M1I7M1I17M1D5M1I6M1I14M1I3M1D22M1I1M1I11M1I14M1I5M1I11M1I2M1I5M1D22M1D9M1I13M1I13M1D5M1D10M1I7M1I6M1I13M1D5M3I10M1I13M1I7M1I22M1I12M1I4M1I3M1I10M2I26M2I2M1I8M1I5M1I27M1I1M1I22M4I9M1I14M1I18M1I7M1D6M1I5M1I4M1D19M1D17M1I8M1D10M1I25M2D4I12M1I6M1D16M1I11M2I2M1D2M1I22M1I9M1I23M1I3M1I4M1I11M1I2M2I7M1I13M1D5M1D10M1D3M1I25M1I7M1I38M2I1M1I21M1I3M1I5M1I10M1D4M1I14M1I13M1I5M1I1M1I6M1I2M1I16M1D2M1I15M1I2M1I9M1I1M1I4M1D16M1I2M1I12M1D9M1D22M3I1M1I10M2I5M1I3M1I52M1I2M1I14M2I2M1I3M1I5M1D28M1I5M1I4M1I2M1I12M3I11M1I3M1I17M1D23M2I13M1I4M1I18M2D6M1I4M1I28M1D13M2I8M1I12M1I8M1I24M1I1M1I10M1I4M2I2M1I3M1I7M1I27M1I1M1I16M1I19M1I15M1I9M1I2M1I21M1I9M1I20M1I13M2I8M1I7M2I29M1I5M1I28M1I2M1I7M1D82M1I34M1D45M1I14M1I7M1I4M1D6M1I68M1I3M1I7M1I33M1D5M1I24M1I1M1I7M1D5M1D4M1I14M1I17M1I24M2I13M1D11M1I29M1D4M1I22M1I7M1I11M1I7M1I18M1D26M1I3M1I24M1I3M1I12M1I1M1I9M2I1M1I23M1I23M2I10M1I19M1I11M1I31M1D10M1D4M1I32M1I18M1I16M2I4M1I9M1I7M1D6M1I14M1I51M1D9M1D11M1I5M2D5M1I7M1I20M1I14M2I59M1I8M1I5M1I16M1I5M1I33M1I19M1I5M1I18M1I35M1I21M1I6M1I1M1I31M1I56M1I18M1I37M1I28M1I21M1I46M1D8M1I4M1D36M1D31M1I23M1D8M1I76M1I8M1I3M1I4M2I67M1I4M1I29M1I6M1I66M1I5M1I12M1I5M1I10M1I20M1D25M1D12M1I9M1I13M1I57M1I3M1I44M1I9M1I22M1I17M1I3M1I5M1I17M1D4M1I16M1I18M1I3M1I6M2I6M1D2M1D23M1D37M1I3M1I4M1I2M1I3M1D14M1I13M1I62M2I18M2I34M1I10M1D1M1D22M1I25M1I30M2D45M1I3M1I4M1I58M1I5M1D5M1D42M1D5M1I3M1I1M1I48M2I9M1D11M1I5M2D16M1I13M1I8M1D21M1I6M1I19M1I3M1I3M1I36M1D2I14M1I49M2I5M1I16M1I8M1I8M1I11M1I4M1I34M1I4M1I32M3I21M1I31M1I9M1I17M1I12M1I23M1I19M1I7M1I10M1I33M1D78M1I10M1D22M1I5M1I14M1I10M1I23M1I41M1I9M1I11M1I24M1D4M1I2M2I12M1I43M1I14M1I26M1I6M1I14M1I18M1I2M1I33M1I42M1I35M1I39M1D6M1D29M1I61M1I8M1D5M1I20M1I16M1I17M1I5M1I17M2I3M1D19M1D18M1I11M2I19M1D9M1D30M1I6M1I38M1I18M3I25M1I18M1I57M1I11M1D18M1D15M1I12M1D12M2I12M1D12M1I31M2I21M2I1M1I9M1I3M1I7M1I6M1D4M1I4M1D1M1D8M1I7M1I19M1I35M1D20M1I13M1D8M1I14M1I27M1I11M2I27M1I33M1I10M1I2M1I5M1I3M1I5M1I10M1I2M1I11M1I1M1I36M4I8M1I25M1I2M1D4M1I11M1I6M1I18M1I4M1I5M1I41M1I6M1I38M2I9M1I4M1I20M2I32M1I1M1I8M1I2M1I45M1I30M2I16M1I9M1I22M1I4M2I12M1I9M1I34M1D60M1D6M1I1M1I1M1I34M1I15M1I9M1I16M2I25M1I6M1I9M1I14M1I63M1I14M1I36M1D14M1I12M1I9M1I10M1I11M1I6M * 0 0 GAACAAACTGCACAAGACGAATGAGTTTGACCCAAGCTGACAGAAGTAGGCTCAGCGAAGATGGGTAATAACTCCTCTGAGCTAAGGACCATGTTTAACCTAATGCCAAGGAAGCTAAAGACCGTGATAAAAGGTTAACATGGAACTGCTAACTAGAATAATCCAGTATACAAGAGAAAGAACATAAATGAACTCTATAATGGGGCTGAAAAACACAGCATGAGCACTTCATGAAAGAATACACAAGTATCAACAGACTATAGTCCATCAAGGGTCTAGAAGAAAAAGGATATCAGAGATTGAAGTCACATTATGAAATAAAAGCATGAAACGATACTAGAGTATAAAAAGAATTAAAAGAAATGACGCAAAGCCTCAAGGAAATAGGGCACTATTGAAAAAGGACAAGCTAACCTTTGATTGGTCGTATGTGAAAATGACAGGCGAAATGGCGAAACCAGTTGAGAAAATCACTCTTCAGTGATATCATGCAGGCAGAAGCCTCCCCAACCTGCAAGGGAAGCCAAATTCAAAAGATTCAGGATAATACAAGAACATACTCCAGATCTCCTTGAGACGAAGAGCAAACGCCCAGACACGATAATTAGATCGAGTATCCGCTCAAGGCGTGAAATGAAGGGGAAAAAATGCTAACAGATGGCAGGGCCAGAGAGAAAGGTTTGGGTTATCCCACAAAGGGCCAGACGCCATCATGCACTTGGAGCAGATCTCTTGCAGAAATCCTACACACTCAGAAGAGAGAACGGGGCCAATTTCAACATTCTTAAAGTAAAGAATCTTTCAACAAGAATTTCATATCCAGACGCCAAACAACAGCTTCGTAAAGTGAAGGAGAAATAAATTTTTTACAGACAAGCAAATATGCGAGATGCAGTTTTTCTCACCACCATGCCCTGCCCTAGCCAAGAATTTTCTGACAGGGGAAGCACTAAACATGGGAAGGATAAACCAGTACCAGCCACTGCAAAACCATACCAAAATGTAGAGACCATCAATGCTAAGGAGAAACTGCACTCAACGTCAGCAGTCAAAATAAACCAGCTAGCATCATAATGACAGGATCAAATTCTACACATAACAATATTGTAACTCTTAAATGTAAATGGGCTAAATGCCCAATTAAAACGACACATACTGGGCAACATTGGATAAAGACGTCAAACCCATCGGTGTTAGCTGTATTAGGAGACCCATCTTCACGTGCACAAGACAACACGATAGGCTCAAAATACATAGGATGGAAGGAAGATTTACCAAGCAAATGGAAAAGCAAATAAAAAAAAAAGCAGGGGTTGCCAATTCCTAGTCTCTGATCAAAACACACTTTAAAACCAAGAAAGATCAAAAGAGACAGTAAGAGGGCTTTACAAATGGCTATAAGGGAATTAATACCAAACAAGAAGAGTTACACTACCTATAATATATATGTGCTAAGTTACAGAGCCACACGCTAGTTCAATAAAGCATTTGGTTGTGTTGTTAAAGCAGCGTATTTGTAAAATATTTGCTAAATTTTTGTGCTAGTTTCTTTTCACACATGTTTTCAGGGTGTTCCAGCATATTCGGTCGAAAAGACTATCCTTGTCCCCCATTGAATTGTCTTTCGTTACCTTATGTCAATTATTCAATGAATATATTTGTGTGCATCTGGTTTGCTGAGCACACCTCATTCTGTTCCCACTGACGTAGCTGCCTTGTTCTTACCACAACCACACGAAGACTAGTGGAATTACCAAAAGTCTTGATGTCAGGTTAATGTGCACGTTCTCTCTATGGTTCTTTACAGTACCTCCATTGCCTATCTACTAGGTCTGTTTTTGCTTTTCCATGTAAGTTATTATCAGGTTTTCAATAGTGTACAAAAATTGTTGCTGGATTTTGATTCGGGATGACTGTGAACTCTATAATCAAGTTGATGAAAAACTGACATGCCGTAACAAAATTCGAGTCCTGTCCATGAACATAGAATCTCTTCTCCATTTCTATATATCTTGTCTATTTCGATTAGAGTCTTATGTACGGTTTTATCGTAAATACTGTCCCAGATGTAGTTAGATTAATATTTAATGTGGGGTTTTTTGCGCTGCTAATAGAAATTATATTGTTGCTTTCGATTCAAAATTTCCAACTTGTTTATAGCTGGTACATTTGGAAGCATTTGACTTACATATAGATTTTATCGTTCCTGCACCCTCTATGACTACTTATTAATTCCAGATATTTTTTATGTCTAACTGTTTCAGATTTTCAATATAGCTGATGACTCACTGTGAACAAGGACAAATTTATTTATTTCTTCTCAATCTGTCAGTATTATTTTTTTGTCCCTTTTCTTGCCTTCCTGTGGATTAAGCTCATACTGATCCAGCACAGATCGTTGAAAAAGGGTGATCGAAAAAGGAAAATTTTTACTTGTTCCAACTCCTTAGAGGGAAGATAACTATTTTCTCATTATTAAAGTTGATGTCTGTTAATAGGTTTTCATTAGATTCTTTTTTGCAAGATGGGAAGTTCCCTCTCATACTTCCTCTAGTTTGGCTGAGAAATTTTTATATAAATAAGTGCCGGTTTTACCAAACTTCTTTTTTGCATTCTGTTTGTTATGATCATATGATATTTCATTATTCAGCCTATTATAATCCGATAGGTACTTTCCTGCTTTCTAATGTTAAAAGCATCTTCAATACCTGTGAATATCGAATCCGAAATGCATACAGTGGTACATATTCTTTTTCTATACATTGTTGGATTTAATCTGCTAGTTATTTCTGATTTGAGTCAATTATGTATCTAAAGTTACATTAGAGCTATTGATCAAGTAGTTTTTCTTGGATGTCTTTATCTGTTTTGGTTTTATAGGATAATAACCTGGGCCTGACTCAGAATAGTTAGCAAATGTTCCCTCTGCTCTAGCCTTTTTCTGTTTCTGTTTTTTTTTTTTTTGGATTGGGTCTCACCTGTACACTTTAGTCTGGAGTGCAGTGACATGATCTCTGCTCACATACAACCTCTGTCTCCCAGGCGAAGTAATTCTCCCACCTCAGCCTCCCAAAGGCTCAGCTGGATACGTAACGCCATGTCCTCCATTGCTCACGATAATTATTTGATTCTTTGGTAGAGACAGGAAAGTTTGGCTCATGTTGCCCACGGCTGGTCGGAAATTCACTTGAGCTAGAAGCAATCCGCCCGCTGCGGCTTCTACAAAGTAGCGCTGAGATATAGGTCATGGAGACATAACGTCAGCTCAGCTTTTATTTTTGAGAAGAGATTGCAAAGAATTGGTTATCATTTTTTTTGCTTTAAGAGTTGGTAGGATTTCACTTGAGTGCAAACATCGGTGCCTGTAGGTATCTGTTTTTGGAAGGTTATTGAATTAATGAGTAATTTCTGTTCAATTGATTATAGAGCCCATTCAGAATCATCTGTTTCCTCCTCGAAGACAGAATGCTTGATAGTCTGTGTCCTTAAGAAAAATTAGTTTATTTCATACCAAGTATTCAAATTTACTGGCTACAGGGCTTCTTCATAACCCTAATATCCCTTTAATTATTCATCCTTCGTAACATCTATGAAAGACCATTCCTGATAGAGTTCTGTTTTATATTCTCAATATTAGATATTACATCATGTCTTCTTTTTCTGCTTGATAGCCAGAGTTAGGGCTTTACCACATTTTATTGATCATTATTTTAAAGAAACCCAGATTTTAGTGTTTGCTGATTATTTCTACTGGTTCCTGGTCTCAATGTAGCTGATTTCTCGCGTCAAAGTTTTATTATTTTCTTTCTTCTGCTTGTTATAGGCTGATATTGCACTTCTTACTCAGTTTTCTAATGGTGGAAGCTTCGACTATTGATTTAAACTTTTTTCGTTATTACTAATCTATGCGATTCAATGTTACTAAGATTATTCTGTGGGAAGAGTCAGATTTCATGCATCCCACAATTTTGATAGGTTATATTTTCATTTAGGTTAGCAAATTAATTTAAATTTCCCTTGAGATTTCTCTTTGACTTATGTGCTTATTTCGAGCAAGTGTATTTTTACTTCCAGAATATTTAGAGATTTGCAGCGTGTCTTTATGTTATTAAATTTAATTCCATTGAATGTTTAATTTCATTATAATCAAAGATCCATAATTTATATGATGTCTAGGATTTTAATATTTATTGCGAGATGTCTTTATAGTCAAGAATATAAGTCTTGTGTGAATAGACCATGTTGTAATTTGTAAAGAATGTCGCATTCTGTTGTTGTTCATGAAGTATTCTATGAAATGCAATTAGTTTCAGGTTGTTTGGATGCATTTGCTTTTCATAGGTCATCTATAACCTTACTGATTTTCCTGTCTGCTGGCTCTGTTAACTACTTGCACAGAGGGGTGCTGAAGGCCCACTAACTATAACAGCAAGTTTTTCCTAATCTCTGTTAATTGAATTTTATTTTACTGTCAATTATTCCCATACGAAACTCCGTCACTTTTGTTATGTAAATCTTAGAATGTTCTAGACGCCATGCTATTTGCTCTTCTTTTGAATACATGTGCTATTTTAACCTTTCTTAAATGTTAAGGCAGGTTTACTGGTAACAATCCCTCAATCCTCTTTTTGTTGTCTGACAAAAATCTTTAATTCGTCCTTCACTTTCTGAAGAAACATATGTAGACAGAGATTCTTATCTAGTGGATGGGTTTTTGTGTATGCTCGTTGTTAACTGTTGTTGTTGGTTCACATACTGTACATTTCTCACTCCTCCTTCGTCCTTACTTGTGTAGCTTGTGGCCAAAACAGTCCAGTATAATTCTTAGCTGTTGTCTCTCTACAGGTCAAGCTGTTTTCTCTGGCCTTGTTTCCAATATTGTCTCTTATTCTTGTCATGATTTCTGCAGGTTTGGATATGACATGCCTAGGTGTAACTTTTTGCATTATCTTTTTTGTTGTTTATTTTTTTTTTCGTGAGACCTCTGGATATGTGGTTTAGTGTCTATCATTAATGTTTGGAAATAGCATTCTAGCTATGTTATTTCAATATATCTTCTGCCTCCCTTTTCTCCTTTTTTCGTGGGATTCTCATATCTGCATGTAGTTATAGCCTTTTGTAATTGTCCCAGAGTACCTTGGAACGTCAGTTTCTTTATTCTCCATTCAATTTTTGCTTTCCTCTTTGCTTTTGGCTTAAAAAGTTTCTCCATTGACGTATCTTCTAAGCTTTGCATTGATTTCCTTCTCTTAGTTTATATCTAGTACTAGTTAATAGAGCTCATCAAAGGCACTTTTTATTTCTGCTACATTTGCTTTTGAGTTCTAGCACCTCCTGTTCAAATTCTTCCTTACAGTTTACCCATCTTCGTCTCGACATCTTTCCATCTGTTCTATCGCTTGTCATCTACTTTTCCTTATTAGAACCTTAATTCTGTTAATCATAGTTAAATTCCATATCCAATAATTCCCAAATCACTGCCATATATGAGTCTGGTTCCTGATGCTTATTCTGTCTCTTGTCATTGTATATTTTTTTGTGCCTTTCACTTTCTTGTTCATTATTTGCTGAGAGCAGACGACAGTATTTTGGGTAATGAAAGGCTATATAATTAGACTTTTAGCATGTAGGGTTTTACATGCAATCGTGGCTGTATCTTAATGGATTGATGTAGACTGGCAGGGACACAGCCTTTAGTTTCTTTACTCACTTTTGCTCTCTGGTTGTCTTTTGTGTTTCCCTGGAACCATTTTTCAAGTAGAGGCTTGCATCTTTGCAGCTCTCTTGCTGTTACATTTCACTGGATATTTATTTGGAGCCCTTGTTGGATGTGATGCTGAGTGGTTAGGCGTGAGGAAGCTCATCCAGTACTTCTTCAGATACTAAATCTCAGCTTGCAGTATACCCAAGGTTCTCCTGGATTGGTGACTTTTCAGATGTGCTTTGTGGTTGTGTTTCATTTTGTTTGTTTGTTTTTTAGCTGTGCTTCCTTCCTGCCCCCTATGTGAGGCCAGGCAGGATAGAAGGGGAGCCACAGGGCTAAGATTTCGCATTCTCACGTTGGACAATGTTCGTTAAGGTAGCCTTCCTATGAAGGACTGGTTTTCTTGATAGAAGATGCAATGGGCTTGCTTCATTAATACTTACCTTCGTCACTTCCATGCTCAGAACTGGAAGGGGATTTTTACCTTGTACCTTCATGTGAGAATCTGGTGAAGCTTCCTGCGAGGAAAAACTCAAGAAAGCTGTCAGGCGGCTCCCTAGGAAATAGTTAGAACTCTCGAAGCCTCACATACTAGTCCCAGCAATTCTTGCAACTTACCATTTACGTGTTCATACCAAATCACTCGGCTCCGCAAGCTTCTGGCCCCTGTAAGCTGACCTTGGGCCATCTCATTCTCTGTGTTGACCTCGTCATCTCCAAGATTGCAGACAGTGGGTTTGCCCTGTGACCCTAATTCCTCTGTATGAAGTCCCAGAAGAAGGTGATAGATTTTCATTTGATTCAGCTTTTTTTCATTGCTTATGAGCCACTAGGATACTCACTCCCAAGCTCCTACTGTGTGCCCTCCTGGAAATCAAGTCAGCAAACTTTTTCTTACTCGGTAAGGTCCAAACTTGATATAATGATTTTTGGCTTTGAGAGATACAGAGTCTCTATTGCAGCTCCTGAAGTGTCCTAATATGTAATATGAAAGTTAAGACCATTAAACAAATTCGAATAAGTGAACATGGCTGCGTGCCCAATTACAGGCGTTCATTTACAAATAGAAAGGGGGTTGAGGCTCACATTTGGCCTGTGAGCCCATTTTGCTGATCCCTGCTCTACATACAGGCTACAGAAATGAAAAGGCCATCAGGGGCCAACATTACGTAATACATGCGCAAAGGTGAATGGGTTTTAAAAAAAAATATGTTCCTTTCACTTATTATATGTGCTTATATGCACATGTCTTTTTTTTTTTTTTTTTTTTTTCTCTTTTTTAAGAGCAGTACCTGCACATGTATTGCCCCAGGGTGGAATGCAGTGGTGCAGTCGTCCAGCTCACTGCAACCTACCACCTCCCAGGTTCAAGCATATTCTTCTGCCTCAGGCCTCCCGACGTAAGCTGGGATTACAGGCATGTGCCCACCACGCTCCAGCTAATTTTTGTATTTTGTAGTAGAGACAGGGAGTTTCACCAATGCTGGGTCCAGTCTGGAGTCAAACTCCTGAACTCAGAGCGAACCCACCCACCCTGGCCTCCCAAAGTGCTGAGGTATTACAGTGTGAGCCACTGTGCCCAGCCAATGCAGCTTCTTATTTCCATAAAAATGTTATATTTCATTTTTCTTTAAACAAATACTCTAGCTTTAGATAGGAATATGAATATCATAAATCATTTATGCTTCTGAATTATAATAAAGAGACAGCCACACAATTGGCAGCTGGGTACTCTGTCTCAGGTGTTGGCTGTTGGGCCACCGAGGGGGTGATGGGGACATTGATGAGCTGGAGTGAAGAGCACCAAAGCCACTGCTGGAGCTGCCAGGCTGGAGAATGAGTAGCCTGCATTGCCAGATTTATAATTCTTCCAGAGAGCTGGGAAATCCAGATTTTATGTTACCTATATAAAATGTTAGTTTGGCCAACTATTCAAAATAAAATTTTTAAACATTGCTCCAGGTCAACAGTCAAACACAAATAGATTTGTCAAGTACAATCCTCTAGCGTGGCTAACATTGGCTGGCTTGTGGTTGAATTGGTCCTCCAACCCAGGCCAGCCCCACAGTTTTACCCAAACGGGTGCCATGGGAGGTTAAAGAAGGTGCCCAAGTCACACAGCCAGGAGACAGCAGAGTAGTGAACCTGGGCTGTAGTAGCTGAGCAGACCCGACTTCTGTATGCTGCAGCTCACAGACCTTCAGAGCAAATGTGACTTTTAGCTTTCCATGGAATGAACATCTCTCTGGGCTCTGTTTTCCAGTCTGATATGGTTTGGCTTGTGTCATCACTCCAAAACTTATCTTGAATTCCCACATATTTTGGAGAGACCTCTGGAGAGGTAATTGAATCATGGGGGCAGGTCTTTCCATGTGCTGTTCTTGTGACAGATGAATAAGTCTCACACGAGATCCTTATGGTTATTATATAGGCGAGTTTTCCCTGAACAAGTCTCTGTCTTTGCCTGCTGCCATCCTTGTAAAACATGACTTGGTCCTCCTTGCTTCCACCAGTTTGTGAGGCTCTTGCGCCACAGTGGAACATCTGAGTCCATTAAACCTTTATTCCAGTATAAATTGGACCCAGTCTTGGGTATGTCTTTATTAGCAGTGTGAAAACACAGTCCCACCACTGAGTCACTTTTCTTCACACTGAATATTAAGAATCCTCGTTGTTTAACATTCTTACATACCTTTTCAGAAATCAACTAACCTTCACACTCTTTTTTTAACAAAGCAACAATTATGTTTCAGTTAGATCAGCAAAAAGAACTAATAAAACAATCACCATTATTGCCTGTGCTGTGAAATCACAGGCAAGAGTGAAATGGTTCCAATCAGTCAAAGGAATTCTACAAAAGAACATCATATATAAGAAATGTGCTTCAAGTTATTTGGATCCAAATTTCATATCGCCTCATTCTCCCTTTAATTTCTGATTTTTTTTTTTTTTTTTTTTAGGAAATAAATGTCCATCTTTTCAGGGATTTCATCTTTTTGAGAAAGCGATATTCTTTCTTTCAGGGAAATAATGAAAAGAAGCATGATATGGTTTGGCTGTGTCCCCACCAAATCTCCATCTAATTTGTACTGCCATAATTCCCACGTGTCATGGGAGGGCCCGGTGGGAGGTAATGGAAACATGGGCAGGGTGGGGGTTTCCCTCACACTGTCTCATGCTTACTGAATAAGTCTCATGAGATCTGATGTTTTTATAAGGGGTTTCCCCATTCACTTGGCTCTCATTCTCTCTTGCCTTGCTGCCATTGTGAAGATTGGTGCCTTTCACCTTCTGCCATGAGTGTGAGGGCTCCCCAGCCACGTGGAACTGTGAGTCCATTAAACCCTCTTTTTTCTTTATAAATTACCCTGTCTCGGGTTATGTCATTGATCAGCAGCAAGGTAATGAACTAATACAAAGCAATTGGTTAGAATATTCCCACATTGAGAGGCTGCCATCTCCATAGGTGCCAACTTTGGGATGGAGAATGGAAGTTGGTACAGGGCCAAGTCCCTTCTCAGAGGCTGCAGTGAAAGTGAGTGAGGAGAGGGGGGAGATTTCAGGAGGGCCTCTGGCTTGCTCTTGGTTGTGCCATCACTCTTGTTCCTGTGATCCAGAGATCTTCAAACAGACCAGGAATATCACCTCTAAAGAGCACCATAAATGGTAGGAGCGGATCTGATTAAAGAGATATGCATAACCCTCATGAATTTCTGCTGAAGTATTTTGTCTGAGTGGCCAGTGAGAAATTATCCTAAAACCGTGATGACATCACATCCTCATTCTATTTCATGGCTTACAAATTTTGAGAAATGCAAGCAGTGTATTTTTCTGGCATTAGGATAAAATGTAATATTTTGAAAATATATTTGAAACTCGCAGAAGACACTAGCAAATATCCTGTGATATATATGCAAAATTGAAAATTAGGGAGAAATGGTTAGGTGTGTTGAGCAGAGTAACAAAAAAGAAACCGTGCCATTGGCTATACCAGGTAACATTGTTAATATTTTCTCATTTCCTAAACGGTGTTCTAAGGAGAAAATAAGAGTCAGCTCAAAATAAACTGTTTTGCTGACCAGATAAGCTCTGGAAAGGGTCTAGCTGAAGTGAAACAGGTGTGTTGACAACAATATCTCAGCACTTTTAGTGTGCTTGTATGCGTTGTATTCTTAGATAGAATCTCCTACTATGAGGTATTTTTCAAGTCTATTCACCATCAAAGCCCATATTCAAATGCAGAATCTTTCAACTCCCCTGAGAAGTATTTAATGGAACACAATTTGGGAAACATTGATCATTCTTAGATTTGTTGGGAAATTTAAATAAATTCACATATTTAGAAAGTTAACTTTGCATAATACTCTTAAGGATACAACTAATTTGTTTGAAGTGTACCAGAGGACCTCTTCTCTATCTTTAAATTAAGGTTATTAATTCCAGGAGGCTACCTACTCGTTTTCACCCTGGGTAAATATTGCCAAGTTCTTTCTGGGTAAGGTGTTTTGAATAAGTTGATATTGCTTCAGATTGTTTATTAATTAATTAATTAATTTTGGGACAGGGTCTCACTCTGTCACCCAAGGGCTGAGAGTGCAGTGGCACAATTCATGGCTGCACTGCAGCCCTCGAGCTCCATGGGACTCAAGCAGTCCTCCATCCTCAGGCCCCTGCATATGCTTGGGACTGCAGGTATGTGACACCATGACCAGCTATTTTTTTTTGTATTTTTAGAAAAGAATGGGATTTTGCCATGTTGCCCAGGCTGGACTTGGAACTCCTTGGACTCAAGTCATCTACCCTGCCTTGGCCGTCACAAAGTGCTGGTACTACAGGCCATGAGCCACCGCGCCCAGCCCTGCTTCCAGATTTTTAGAATCACGATTTTGCCGTTGTTGTTTCCCCAGATGGTTTAGAAAGAGTGAAACTATGCACATCATATGAAATTTTCAATGCCAGAAGAAGCTAATCAACAAAATTGTTGCCATTAAAATTTGTCTGGGGCAAGATAGAACACCTGTTTTAGTTTATCAATTAATACTCTTTTCATTATTGTTTTGCCTGGCTTTCTGTCTCACTTTGCCTGGCTTTCTGTGCAATGATACCTGGGCTAGGGTTCATTTATGATCTTGTAGCAAGGAACGGCAATTAAGATTTATGCTTGTGGGAGAATTTGCTTGGCATCTCCTGCTCTGTGAGAATAACACAAGGTAGTGATCTTTCTAATCCTAGCAGTGTGGAAGGGAATGAACTTCTGTGAGAAAAGGAAAGAGTGACTCAAGGAAACGGAGAGCTATTGTTAGGACTGAATTGCTTGGGGGAAGATGGTAAAAAAGATGTTGGGAAATATGCAACTGGAGGAGAGGTAAGTGGATAAAATAGAAACCCTTATTAAGAATTAGGAGCTAATGGGAGCGCAGCATTTGCTTCTGACAGAGGACAGGCAGGCCTTGCCTGGGAAAGTCTAGCCCAGGAGGGCAGATGGTGGTGAGACTATCACGGCTGGCCTCAAGCTCTGGTTAGATTCTTATTACTCAGTGCTCTTATATCAAGATGATTCTGGTGAGCAGTCCAAAAAGACGGTCACAGTTCACGTGGGATATAATTTAAAAGAGGCCTACAATCCTGCCTCATAGCGGAGAGCGGACTATTTCTTAGGTGCTCAAGGATACCTTTCTTGACAAAAAACACCTTTAAAACAAGGCACGGGGCTTGGGGGACATAGTCCACACGAAATCTCGAGGTTCTGGGCTGAGAAACAGCATATTTAAATCTTCTCATACAAAAATGACATTCAGAAAGCCCAGCGAGACAATATCCGACTCTTTACTGTAGAAAGAATAAAAAAATGCCCACATTTCTGGGAGTGAAATTGCAGACAGTCTGACTGCCCTGACGGCCACTCGGTGTACATTTCTCCCCTCCATCTCGGACACTGCACCCAGGCCGTTGTGGACTGAACTTGTGTCCTACCTTGGTTTGCCCCTTCGAGCAAACTGAGCAACTAGAGAAATTGAAGGCAAACCAAGAACAGCCCTAGCATGTCTGTAAAGAGCCCAGCACGGGCTGCCCCTATGACCCTGGGGCTGGCAAGGACCGCCCATCCCCTTCCAATCAGAGCTAGAAGGGCTATCCACTCCGGTGGTGAACATAGTGTGCTTATACGACTTTATAGTGTGGGGCTGTGTGATTCTGCCAAAGTGCTTGAATTGGCACTTCTCCTGTGTCTAGATTCCATCCCCTCCCTCCCCGCAATAAAAATGCACATTTCATAGGCTCTGCCTCTAATCTGCATTAATAAGATTCGAAGTTTTACATGCACCATGTCCACGACCTTCAGGGCCAGGCATCAGGGTGGCCGTCACCTAATTCATGTTTCCACTAAAGTGATGTATTCATGGGGAGGTGAATTATTCCAAGGTCTTGAGATAAGATGAGACCTATGCAAGGTCCAGTCAGCTTTGGTGCAAATGCACTTGGCTTAGAGTTAATAGCAGTTTTGACTCAGGCACTTCTCAGGGCCTTTTGAATATTTTCATTGCAAACCTCTCTAAAAGACGTGAGAGTGCTTGCAGAATTTCTCCAAATACAGTTCACCCTCAGAACCGTTCTCTTCAGAGCAGCTTTGTCCCACAAAATATAACATTAGAAACATAACATGTGGAACTGCATATGTGGTTTCACGGATTTTCAAAAAGTCATATAAAAATGTAAAAATGCAAACAGGTGGATAATTAGTTTTCAATTTTTTTTTTTTTTTTTTTGGGATGGAGACTCAACCGTGTCGTCCAGGCTGGAGTACAGTGGCAGGCGATCTCGGCTCACTGGCAACCTCTCGCCTCCCGGGTTCAAGCGATTCTTCCTCGGCCTCAGCCTCCGCGGGTAGCTGGGGACTACAGGTGCACGCCTCCACCCACAGCTAATTTTGTAATTTTAGTAGAGACTGGGTTTTACCATGCTGGCCAGGATGGGCTCAATCTCTTACTTCACTCGCATTTGCCTGCCTTGGCCTCACAAGGTGTTGGGATCTACAGGTGCGAGCCATCTGTGCCCGAGCCTAGTTTTCAATTTTATAATAATACTGTTAGGTAACACAGCAGTATCATAAACTATTAGTCATTTCCACATGTCAATTCATATGAAATATCAATGAGACAATTGATACCCTTTGATGTGACACCATGTCTTTGAAAAGAGGGTATGTATTTTAATGGCACCCTCAATTCAGACCAGCCATACTTTGACTCCCAACACTGTTCGATGCCTAGGGCTACCCCTACTGGTATGGTGCAGACTCCAGCATGTCTTCAGGA * NM:i:1158 ms:i:15084 AS:i:15084 nn:i:0 tp:A:P cm:i:299 s1:i:4031 s2:i:61 de:f:0.0882 SA:Z:chr22,17457038,+,1380M107I11530S,60,197; MD:Z:3T4^T43^T19^GA18^T4T2T53^G0T46T4C2C0T15C1A3C0A27A31^C7T24A14^CT3G0A5T20^GC6^A19C9^A51C0A58^C0T2^AT69^G7A0T20T11G1A5A7T9^G2^TT2^A6^A29A34^A2^AC0C6^T16T33A31G22G1A17^C32^T19^A22^G13C7^A32C8^C66C16G7^T3^A7^AT18^C39^T53T29T12^T40C7^G22^T13T1T15^A13^T25G1T0A34^T37G4T0A7^G24^T6T36^T7C6C39^C39A0G68C0T8A0T1C26T5^C46T1T25^C63T22T32C20A34A1G52T1C8^T31T31T24^C35^C3^T50^C10^T70^G21G13T1^C85T15T12C10^T9C54C27^T0A6A14^G8A9A91^T4^C5C17T68^C10C21^A22^AT35T24^T13^A70A69G2^T31T43^G2^T21T19^T40T0G7^T9^A89^A3T6A3T24^A6T142G12^T2^G8A42C42A51^A10C27T6^T28^A38T0G14T16^T22^T15T19^A5^T4C2G28^A6A12A0T28G0T7A21A98T8T31A0C3^T6T5T2^C4A14^T1G23^A15G0A18^AA18^C29^C96^T5^T10^A23C28C60^C21C1A23C13^C33^G2A5T12G8^G3A5^A27G91^A4G49A27^C58^AG15A22^C98C164G66^A27C0A87^A70^G117^G37^T5^G72^A40^T69^G92A0C102^G10^G88A1^G71^C9^T16^CA36T151C55A227G2A38^C12^G1G1A32^G54^T99C31C75T115^G25^G50G0T45A1T34G0A9T2A4A5T51^A53^C2^T23^G8A0T39^A151^A1^A24C18A33^AG80G34^A5^T42^T66^A16^TG37^T70A17^A367^G57A1T28^A4A154^G112G0G123G0C51^T6^T98^C83^C19^C0C47^C9^A45A83C0A72^C18^T27^A24^G90^C8^G1^T69^T33^C24C46A98A56C12^G25G49G4T54G254^T60^G31A217^G11G0T3G1C43 rl:i:661