diff --git a/include/stellar/app/stellar.main.tpp b/include/stellar/app/stellar.main.tpp index 5e75a2d1d..ecb30d3e3 100644 --- a/include/stellar/app/stellar.main.tpp +++ b/include/stellar/app/stellar.main.tpp @@ -155,14 +155,11 @@ _stellarMain( _postproccessQueryMatches(databaseStrand, refLen, options, forwardMatches, disabledQueryIDs); }); // measure_time - if (_shouldWriteOutputFile(databaseStrand, forwardMatches)) + stellar_runtime.forward_strand_stellar_time.output_eps_matches_time.measure_time([&]() { - stellar_runtime.forward_strand_stellar_time.output_eps_matches_time.measure_time([&]() - { - // output forwardMatches on positive database strand - _writeAllQueryMatchesToFile(forwardMatches, queryIDs, databaseStrand, options.outputFormat, outputFile); - }); // measure_time - } + // output forwardMatches on positive database strand + _writeAllQueryMatchesToFile(forwardMatches, queryIDs, databaseStrand, options.outputFormat, outputFile); + }); // measure_time outputStatistics = _computeOutputStatistics(forwardMatches); }); // measure_time @@ -219,14 +216,11 @@ _stellarMain( _postproccessQueryMatches(databaseStrand, refLen, options, reverseMatches, disabledQueryIDs); }); // measure_time - if (_shouldWriteOutputFile(databaseStrand, reverseMatches)) + stellar_runtime.reverse_strand_stellar_time.output_eps_matches_time.measure_time([&]() { - stellar_runtime.reverse_strand_stellar_time.output_eps_matches_time.measure_time([&]() - { - // output reverseMatches on negative database strand - _writeAllQueryMatchesToFile(reverseMatches, queryIDs, databaseStrand, options.outputFormat, outputFile); - }); // measure_time - } + // output reverseMatches on negative database strand + _writeAllQueryMatchesToFile(reverseMatches, queryIDs, databaseStrand, options.outputFormat, outputFile); + }); // measure_time outputStatistics.mergeIn(_computeOutputStatistics(reverseMatches)); }); // measure_time diff --git a/include/stellar/stellar_launcher.hpp b/include/stellar/stellar_launcher.hpp index 3780acd4a..198f83141 100644 --- a/include/stellar/stellar_launcher.hpp +++ b/include/stellar/stellar_launcher.hpp @@ -6,14 +6,6 @@ namespace stellar { -template -bool _shouldWriteOutputFile(bool const databaseStrand, StringSet const, TId> > > const & matches) -{ - // if databaseStrand == true always outputs - // if databaseStrand == false only outputs if TAlphabet == Dna5 or TAlphabet == Rna5 - return databaseStrand || IsSameType::VALUE || IsSameType::VALUE; -} - template void _postproccessQueryMatches(bool const databaseStrand, uint64_t const & refLen, StellarOptions const & options, @@ -37,7 +29,7 @@ void _postproccessQueryMatches(bool const databaseStrand, uint64_t const & refLe // adjust length for each matches of a single query (only for dna5 and rna5) // TODO: WHY? This seems like an arbitrary restriction - if (_shouldWriteOutputFile(databaseStrand, matches)) + if (databaseStrand || IsSameType::VALUE || IsSameType::VALUE) _postproccessLengthAdjustment(refLen, matches); } diff --git a/test/cli/gold_standard/dna_both/minLen20.gff b/test/cli/gold_standard/dna_both/minLen20.gff index 2e8498193..bf6449c0f 100644 --- a/test/cli/gold_standard/dna_both/minLen20.gff +++ b/test/cli/gold_standard/dna_both/minLen20.gff @@ -602,3 +602,100 @@ seq1 Stellar eps-matches 739905 739923 95 + . seq2;seq2Range=623901,623920;cigar seq1 Stellar eps-matches 841474 841492 95 + . seq2;seq2Range=619616,619635;cigar=9M1I10M;mutations=10A seq1 Stellar eps-matches 912127 912145 95 + . seq2;seq2Range=281234,281253;cigar=3M1I16M;mutations=4A seq1 Stellar eps-matches 997497 997515 95 + . seq2;seq2Range=570509,570528;cigar=15M1I4M;mutations=16G +seq1 Stellar eps-matches 765718 765740 95.6521 - . seq2;seq2Range=480465,480487;cigar=23M;mutations=15T +seq1 Stellar eps-matches 259870 259892 95.6521 - . seq2;seq2Range=242055,242077;cigar=23M;mutations=14T +seq1 Stellar eps-matches 986804 986825 95.4545 - . seq2;seq2Range=255704,255725;cigar=22M;mutations=9G +seq1 Stellar eps-matches 795116 795137 95.4545 - . seq2;seq2Range=100168,100188;cigar=15M1D6M;mutations= +seq1 Stellar eps-matches 590053 590074 95.4545 - . seq2;seq2Range=245228,245249;cigar=22M;mutations=17A +seq1 Stellar eps-matches 563210 563231 95.4545 - . seq2;seq2Range=355253,355274;cigar=22M;mutations=4A +seq1 Stellar eps-matches 353237 353258 95.4545 - . seq2;seq2Range=905108,905129;cigar=22M;mutations=3G +seq1 Stellar eps-matches 971517 971537 95.4545 - . seq2;seq2Range=55261,55282;cigar=9M1I12M;mutations=10T +seq1 Stellar eps-matches 960203 960223 95.238 - . seq2;seq2Range=617579,617598;cigar=19M1D1M;mutations= +seq1 Stellar eps-matches 887924 887944 95.238 - . seq2;seq2Range=527360,527380;cigar=21M;mutations=7T +seq1 Stellar eps-matches 814505 814525 95.238 - . seq2;seq2Range=546149,546168;cigar=8M1D12M;mutations= +seq1 Stellar eps-matches 721801 721821 95.4545 - . seq2;seq2Range=739398,739419;cigar=10M1I11M;mutations=11C +seq1 Stellar eps-matches 576858 576878 95.238 - . seq2;seq2Range=491562,491582;cigar=21M;mutations=3C +seq1 Stellar eps-matches 558950 558970 95.238 - . seq2;seq2Range=683633,683652;cigar=12M1D8M;mutations= +seq1 Stellar eps-matches 479518 479538 95.238 - . seq2;seq2Range=83080,83100;cigar=21M;mutations=19C +seq1 Stellar eps-matches 433635 433655 95.238 - . seq2;seq2Range=879014,879033;cigar=2M1D18M;mutations= +seq1 Stellar eps-matches 215664 215684 95.238 - . seq2;seq2Range=212347,212367;cigar=21M;mutations=10G +seq1 Stellar eps-matches 182569 182589 95.238 - . seq2;seq2Range=29771,29791;cigar=21M;mutations=13A +seq1 Stellar eps-matches 995415 995434 95 - . seq2;seq2Range=228306,228325;cigar=20M;mutations=8G +seq1 Stellar eps-matches 930403 930422 95 - . seq2;seq2Range=832382,832400;cigar=15M1D4M;mutations= +seq1 Stellar eps-matches 924969 924988 95 - . seq2;seq2Range=932201,932219;cigar=18M1D1M;mutations= +seq1 Stellar eps-matches 915539 915558 95.238 - . seq2;seq2Range=799352,799372;cigar=8M1I12M;mutations=9C +seq1 Stellar eps-matches 886471 886490 95 - . seq2;seq2Range=337574,337592;cigar=16M1D3M;mutations= +seq1 Stellar eps-matches 875106 875125 95 - . seq2;seq2Range=127516,127535;cigar=20M;mutations=9G +seq1 Stellar eps-matches 860711 860730 95 - . seq2;seq2Range=137910,137928;cigar=6M1D13M;mutations= +seq1 Stellar eps-matches 857248 857267 95.238 - . seq2;seq2Range=911396,911416;cigar=16M1I4M;mutations=17A +seq1 Stellar eps-matches 805362 805381 95 - . seq2;seq2Range=46863,46882;cigar=20M;mutations=4G +seq1 Stellar eps-matches 782079 782098 95 - . seq2;seq2Range=744396,744414;cigar=8M1D11M;mutations= +seq1 Stellar eps-matches 763601 763620 95 - . seq2;seq2Range=245161,245179;cigar=10M1D9M;mutations= +seq1 Stellar eps-matches 761138 761157 95 - . seq2;seq2Range=849566,849585;cigar=20M;mutations=7C +seq1 Stellar eps-matches 715611 715630 95 - . seq2;seq2Range=769274,769293;cigar=20M;mutations=11A +seq1 Stellar eps-matches 698067 698086 95 - . seq2;seq2Range=872432,872451;cigar=20M;mutations=4G +seq1 Stellar eps-matches 687547 687566 95 - . seq2;seq2Range=870365,870383;cigar=10M1D9M;mutations= +seq1 Stellar eps-matches 677977 677996 95 - . seq2;seq2Range=439433,439452;cigar=20M;mutations=11C +seq1 Stellar eps-matches 673645 673664 95 - . seq2;seq2Range=874398,874416;cigar=14M1D5M;mutations= +seq1 Stellar eps-matches 645038 645057 95 - . seq2;seq2Range=234191,234209;cigar=8M1D11M;mutations= +seq1 Stellar eps-matches 644858 644877 95 - . seq2;seq2Range=676277,676295;cigar=7M1D12M;mutations= +seq1 Stellar eps-matches 644813 644832 95 - . seq2;seq2Range=865801,865819;cigar=17M1D2M;mutations= +seq1 Stellar eps-matches 641510 641529 95 - . seq2;seq2Range=607675,607693;cigar=13M1D6M;mutations= +seq1 Stellar eps-matches 620345 620364 95 - . seq2;seq2Range=670708,670727;cigar=20M;mutations=16G +seq1 Stellar eps-matches 589984 590003 95 - . seq2;seq2Range=562730,562749;cigar=20M;mutations=3T +seq1 Stellar eps-matches 584612 584631 95.238 - . seq2;seq2Range=427382,427402;cigar=19M1I1M;mutations=20G +seq1 Stellar eps-matches 579568 579587 95 - . seq2;seq2Range=894330,894348;cigar=19M1D;mutations= +seq1 Stellar eps-matches 563144 563163 95 - . seq2;seq2Range=603709,603727;cigar=13M1D6M;mutations= +seq1 Stellar eps-matches 559003 559022 95.238 - . seq2;seq2Range=545195,545215;cigar=15M1I5M;mutations=16A +seq1 Stellar eps-matches 555342 555361 95 - . seq2;seq2Range=244036,244054;cigar=4M1D15M;mutations= +seq1 Stellar eps-matches 554010 554029 95 - . seq2;seq2Range=500180,500199;cigar=20M;mutations=19A +seq1 Stellar eps-matches 529470 529489 95.238 - . seq2;seq2Range=261759,261779;cigar=15M1I5M;mutations=16G +seq1 Stellar eps-matches 494432 494451 95 - . seq2;seq2Range=109070,109088;cigar=8M1D11M;mutations= +seq1 Stellar eps-matches 486491 486510 95 - . seq2;seq2Range=60990,61009;cigar=20M;mutations=14G +seq1 Stellar eps-matches 479793 479812 95 - . seq2;seq2Range=853178,853197;cigar=20M;mutations=12C +seq1 Stellar eps-matches 478520 478539 95 - . seq2;seq2Range=100372,100390;cigar=12M1D7M;mutations= +seq1 Stellar eps-matches 472552 472571 95 - . seq2;seq2Range=303200,303218;cigar=7M1D12M;mutations= +seq1 Stellar eps-matches 471158 471177 95.238 - . seq2;seq2Range=32620,32640;cigar=13M1I7M;mutations=14G +seq1 Stellar eps-matches 455364 455383 95.238 - . seq2;seq2Range=669817,669837;cigar=5M1I15M;mutations=6C +seq1 Stellar eps-matches 442608 442627 95 - . seq2;seq2Range=641357,641375;cigar=7M1D12M;mutations= +seq1 Stellar eps-matches 428970 428989 95 - . seq2;seq2Range=791789,791807;cigar=6M1D13M;mutations= +seq1 Stellar eps-matches 408714 408733 95 - . seq2;seq2Range=574104,574123;cigar=20M;mutations=6A +seq1 Stellar eps-matches 404450 404469 95 - . seq2;seq2Range=236052,236070;cigar=4M1D15M;mutations= +seq1 Stellar eps-matches 402092 402111 95 - . seq2;seq2Range=928807,928825;cigar=6M1D13M;mutations= +seq1 Stellar eps-matches 377549 377568 95 - . seq2;seq2Range=920264,920282;cigar=16M1D3M;mutations= +seq1 Stellar eps-matches 362963 362982 95 - . seq2;seq2Range=618522,618540;cigar=12M1D7M;mutations= +seq1 Stellar eps-matches 316447 316466 95.238 - . seq2;seq2Range=962072,962092;cigar=9M1I11M;mutations=10A +seq1 Stellar eps-matches 277755 277774 95 - . seq2;seq2Range=477270,477288;cigar=13M1D6M;mutations= +seq1 Stellar eps-matches 252314 252333 95 - . seq2;seq2Range=518921,518939;cigar=9M1D10M;mutations= +seq1 Stellar eps-matches 238268 238287 95.238 - . seq2;seq2Range=83582,83602;cigar=16M1I4M;mutations=17A +seq1 Stellar eps-matches 123831 123850 95 - . seq2;seq2Range=791401,791420;cigar=20M;mutations=15A +seq1 Stellar eps-matches 90438 90457 95 - . seq2;seq2Range=877934,877953;cigar=20M;mutations=6G +seq1 Stellar eps-matches 77396 77415 95 - . seq2;seq2Range=785298,785316;cigar=7M1D12M;mutations= +seq1 Stellar eps-matches 66376 66395 95 - . seq2;seq2Range=842457,842475;cigar=9M1D10M;mutations= +seq1 Stellar eps-matches 60971 60990 95 - . seq2;seq2Range=760087,760106;cigar=20M;mutations=13T +seq1 Stellar eps-matches 53898 53917 95 - . seq2;seq2Range=949225,949244;cigar=20M;mutations=10G +seq1 Stellar eps-matches 31264 31283 95 - . seq2;seq2Range=993852,993871;cigar=20M;mutations=4C +seq1 Stellar eps-matches 26503 26522 95 - . seq2;seq2Range=842546,842565;cigar=20M;mutations=9C +seq1 Stellar eps-matches 19458 19477 95 - . seq2;seq2Range=30538,30556;cigar=7M1D12M;mutations= +seq1 Stellar eps-matches 954644 954662 95 - . seq2;seq2Range=566182,566201;cigar=3M1I16M;mutations=4A +seq1 Stellar eps-matches 905577 905595 95 - . seq2;seq2Range=528544,528563;cigar=16M1I3M;mutations=17G +seq1 Stellar eps-matches 805470 805488 95 - . seq2;seq2Range=41851,41870;cigar=12M1I7M;mutations=13G +seq1 Stellar eps-matches 795190 795208 95 - . seq2;seq2Range=361614,361633;cigar=8M1I11M;mutations=9T +seq1 Stellar eps-matches 763547 763565 95 - . seq2;seq2Range=469456,469475;cigar=4M1I15M;mutations=5C +seq1 Stellar eps-matches 732502 732520 95 - . seq2;seq2Range=446091,446110;cigar=2M1I17M;mutations=3C +seq1 Stellar eps-matches 643822 643840 95 - . seq2;seq2Range=244584,244603;cigar=16M1I3M;mutations=17T +seq1 Stellar eps-matches 603126 603144 95 - . seq2;seq2Range=192740,192759;cigar=15M1I4M;mutations=16G +seq1 Stellar eps-matches 508105 508123 95 - . seq2;seq2Range=400207,400226;cigar=8M1I11M;mutations=9T +seq1 Stellar eps-matches 447396 447414 95 - . seq2;seq2Range=828387,828406;cigar=17M1I2M;mutations=18A +seq1 Stellar eps-matches 424852 424870 95 - . seq2;seq2Range=882879,882898;cigar=13M1I6M;mutations=14T +seq1 Stellar eps-matches 418601 418619 95 - . seq2;seq2Range=159496,159515;cigar=17M1I2M;mutations=18G +seq1 Stellar eps-matches 399065 399083 95 - . seq2;seq2Range=932045,932064;cigar=14M1I5M;mutations=15C +seq1 Stellar eps-matches 346392 346410 95 - . seq2;seq2Range=246662,246681;cigar=4M1I15M;mutations=5G +seq1 Stellar eps-matches 313455 313473 95 - . seq2;seq2Range=993116,993135;cigar=6M1I13M;mutations=7G +seq1 Stellar eps-matches 303067 303085 95 - . seq2;seq2Range=935085,935104;cigar=8M1I11M;mutations=9C +seq1 Stellar eps-matches 293473 293491 95 - . seq2;seq2Range=421802,421821;cigar=18M1I1M;mutations=19C +seq1 Stellar eps-matches 285955 285973 95 - . seq2;seq2Range=517350,517369;cigar=5M1I14M;mutations=6G +seq1 Stellar eps-matches 273607 273625 95 - . seq2;seq2Range=9351,9370;cigar=6M1I13M;mutations=7G +seq1 Stellar eps-matches 216605 216623 95 - . seq2;seq2Range=669253,669272;cigar=9M1I10M;mutations=10T +seq1 Stellar eps-matches 71324 71342 95 - . seq2;seq2Range=218249,218268;cigar=6M1I13M;mutations=7C +seq1 Stellar eps-matches 3331 3349 95 - . seq2;seq2Range=631237,631256;cigar=9M1I10M;mutations=10T diff --git a/test/cli/gold_standard/dna_both/minLen20.txt b/test/cli/gold_standard/dna_both/minLen20.txt index 5af85a133..b200706f1 100644 --- a/test/cli/gold_standard/dna_both/minLen20.txt +++ b/test/cli/gold_standard/dna_both/minLen20.txt @@ -13075,3 +13075,1264 @@ Query positions: 570508..570528 ---------------------------------------------------------------------- +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 765740..765717 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 480464..480487 + + 0 . : . : + CCATTCAACAGGGGGCAGCAGTC + |||||||||||||| |||||||| + CCATTCAACAGGGGTCAGCAGTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 259892..259869 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 242054..242077 + + 0 . : . : + CTTTAACCCTGGTATATCAGCCT + ||||||||||||| ||||||||| + CTTTAACCCTGGTTTATCAGCCT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 986825..986803 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 255703..255725 + + 0 . : . : + CGTACTCGTCCGGGGTTGCTTA + |||||||| ||||||||||||| + CGTACTCGGCCGGGGTTGCTTA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 795137..795115 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 100167..100188 + + 0 . : . : + TAAGGTCCGTTGATATCGCGCA + ||||||||||||||| |||||| + TAAGGTCCGTTGATA-CGCGCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 590074..590052 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 245227..245249 + + 0 . : . : + GCGGCGCGTAGCTCAAGAGCCC + |||||||||||||||| ||||| + GCGGCGCGTAGCTCAAAAGCCC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 563231..563209 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 355252..355274 + + 0 . : . : + GGCCCCACACTCTTTCAATGCG + ||| |||||||||||||||||| + GGCACCACACTCTTTCAATGCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 353258..353236 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 905107..905129 + + 0 . : . : + AACATTTTAGGATCTGCTTAAC + || ||||||||||||||||||| + AAGATTTTAGGATCTGCTTAAC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 971537..971516 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 55260..55282 + + 0 . : . : + GGTTCAATA-TCCTTAGTAACC + ||||||||| |||||||||||| + GGTTCAATATTCCTTAGTAACC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 960223..960202 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 617578..617598 + + 0 . : . : + CTACGTTCAGCCGCCGTACAC + ||||||||||||||||||| | + CTACGTTCAGCCGCCGTAC-C + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 887944..887923 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 527359..527380 + + 0 . : . : + CATATACAAAGGTCTCCGACA + |||||| |||||||||||||| + CATATATAAAGGTCTCCGACA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 814525..814504 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 546148..546168 + + 0 . : . : + AACTTAAGCAGGACGAAGCGT + |||||||| |||||||||||| + AACTTAAG-AGGACGAAGCGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 721821..721800 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 739397..739419 + + 0 . : . : + TTGGACGATA-TTTCAGATACG + |||||||||| ||||||||||| + TTGGACGATACTTTCAGATACG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 576878..576857 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 491561..491582 + + 0 . : . : + TGAAAGTAGTTAACCAACCGA + || |||||||||||||||||| + TGCAAGTAGTTAACCAACCGA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 558970..558949 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 683632..683652 + + 0 . : . : + ATTATGCTCCTAGAGTGAGGG + |||||||||||| |||||||| + ATTATGCTCCTA-AGTGAGGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 479538..479517 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 83079..83100 + + 0 . : . : + TGGAGTGCGGAGGTCCTGAAC + |||||||||||||||||| || + TGGAGTGCGGAGGTCCTGCAC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 433655..433634 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 879013..879033 + + 0 . : . : + GCAGGTCCGCCGAAGGAGTTC + || |||||||||||||||||| + GC-GGTCCGCCGAAGGAGTTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 215684..215663 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 212346..212367 + + 0 . : . : + GGCGCCCATTCAAGCGTTATA + ||||||||| ||||||||||| + GGCGCCCATGCAAGCGTTATA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 182589..182568 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 29770..29791 + + 0 . : . : + GACGTTGAGATGGCCGACCCG + |||||||||||| |||||||| + GACGTTGAGATGACCGACCCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 995434..995414 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 228305..228325 + + 0 . : . : + CCGAGCAAAAAACGGTCCCA + ||||||| |||||||||||| + CCGAGCAGAAAACGGTCCCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 930422..930402 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 832381..832400 + + 0 . : . : + GCGATGCACCTACGTTGTCA + ||||||||||||||| |||| + GCGATGCACCTACGT-GTCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 924988..924968 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 932200..932219 + + 0 . : . : + CTCTACCTCTTCCACCTAAG + |||||||||||||||||| | + CTCTACCTCTTCCACCTA-G + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 915558..915538 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 799351..799372 + + 0 . : . : + TCATCAAT-CCAGAGTAGCGT + |||||||| |||||||||||| + TCATCAATCCCAGAGTAGCGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 886490..886470 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 337573..337592 + + 0 . : . : + TGTACTGTAAAAAACCGTGC + |||||||||||||||| ||| + TGTACTGTAAAAAACC-TGC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 875125..875105 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 127515..127535 + + 0 . : . : + CGAAGCATTAACAGCATTGT + |||||||| ||||||||||| + CGAAGCATGAACAGCATTGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 860730..860710 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 137909..137928 + + 0 . : . : + CCCGCAGAGTAGAGGGAATA + |||||| ||||||||||||| + CCCGCA-AGTAGAGGGAATA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 857267..857247 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 911395..911416 + + 0 . : . : + GCCATGTCGGCATAGG-CTAC + |||||||||||||||| |||| + GCCATGTCGGCATAGGACTAC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 805381..805361 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 46862..46882 + + 0 . : . : + GTTAAGAGCTGTGGCCTAGA + ||| |||||||||||||||| + GTTGAGAGCTGTGGCCTAGA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 782098..782078 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 744395..744414 + + 0 . : . : + GTCAGTATCCAACAAGTTTG + |||||||| ||||||||||| + GTCAGTAT-CAACAAGTTTG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 763620..763600 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 245160..245179 + + 0 . : . : + CGTAACTCATCTCTTCCTTC + |||||||||| ||||||||| + CGTAACTCAT-TCTTCCTTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 761157..761137 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 849565..849585 + + 0 . : . : + TGTATTATTGGTTCCAAGAG + |||||| ||||||||||||| + TGTATTCTTGGTTCCAAGAG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 715630..715610 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 769273..769293 + + 0 . : . : + CGCGCATTTATCCTGGGGCA + |||||||||| ||||||||| + CGCGCATTTAACCTGGGGCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 698086..698066 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 872431..872451 + + 0 . : . : + GAACAAGCGACGCGGCCGGG + ||| |||||||||||||||| + GAAGAAGCGACGCGGCCGGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 687566..687546 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 870364..870383 + + 0 . : . : + GGTCAAAATGGAATCTGATT + |||||||||| ||||||||| + GGTCAAAATG-AATCTGATT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 677996..677976 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 439432..439452 + + 0 . : . : + CCGGTCCTGTATGTAGCGAA + |||||||||| ||||||||| + CCGGTCCTGTCTGTAGCGAA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 673664..673644 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 874397..874416 + + 0 . : . : + ACCTCAAGTTCGGTATGTGT + |||||||||||||| ||||| + ACCTCAAGTTCGGT-TGTGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 645057..645037 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 234190..234209 + + 0 . : . : + TGCGTGGAGCGTTTCAAAAG + |||||||| ||||||||||| + TGCGTGGA-CGTTTCAAAAG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 644877..644857 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 676276..676295 + + 0 . : . : + GTCAACAGCATGGATGGCGG + ||||||| |||||||||||| + GTCAACA-CATGGATGGCGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 644832..644812 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 865800..865819 + + 0 . : . : + AGGCAATGCGGTGATCGCAT + ||||||||||||||||| || + AGGCAATGCGGTGATCG-AT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 641529..641509 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 607674..607693 + + 0 . : . : + ATCCTGGGAGCCCCGATACT + ||||||||||||| |||||| + ATCCTGGGAGCCC-GATACT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 620364..620344 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 670707..670727 + + 0 . : . : + GCTTGTCTGCATTGTTGCCG + ||||||||||||||| |||| + GCTTGTCTGCATTGTGGCCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 590003..589983 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 562729..562749 + + 0 . : . : + GCGTATCCTCCGGCGTTTAG + || ||||||||||||||||| + GCTTATCCTCCGGCGTTTAG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 584631..584611 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 427381..427402 + + 0 . : . : + ATACGTCCGCTCGTCAGTG-C + ||||||||||||||||||| | + ATACGTCCGCTCGTCAGTGGC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 579587..579567 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 894329..894348 + + 0 . : . : + TATGTTATGTGCGGTGTCAA + ||||||||||||||||||| + TATGTTATGTGCGGTGTCA- + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 563163..563143 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 603708..603727 + + 0 . : . : + CACGCGTCCGGCGCGATGAA + ||||||||||||| |||||| + CACGCGTCCGGCG-GATGAA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 559022..559002 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 545194..545215 + + 0 . : . : + ACACAGTCGGGACAC-TAACT + ||||||||||||||| ||||| + ACACAGTCGGGACACATAACT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 555361..555341 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 244035..244054 + + 0 . : . : + TCTTCGGTGGCGTCAATAGT + |||| ||||||||||||||| + TCTT-GGTGGCGTCAATAGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 554029..554009 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 500179..500199 + + 0 . : . : + GAAGCGGAATAGGACGGTTG + |||||||||||||||||| | + GAAGCGGAATAGGACGGTAG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 529489..529469 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 261758..261779 + + 0 . : . : + GGAAGAGCCGTCGAA-TCATG + ||||||||||||||| ||||| + GGAAGAGCCGTCGAAGTCATG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 494451..494431 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 109069..109088 + + 0 . : . : + CTGTTATATCGCGGACCTCA + |||||||| ||||||||||| + CTGTTATA-CGCGGACCTCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 486510..486490 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 60989..61009 + + 0 . : . : + CTCGGTAATTTGGAAGTAGA + ||||||||||||| |||||| + CTCGGTAATTTGGGAGTAGA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 479812..479792 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 853177..853197 + + 0 . : . : + AGCCAAAGGACAGTGGCTAG + ||||||||||| |||||||| + AGCCAAAGGACCGTGGCTAG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 478539..478519 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 100371..100390 + + 0 . : . : + ATGACTGCGTCCGTGGTTTC + |||||||||||| ||||||| + ATGACTGCGTCC-TGGTTTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 472571..472551 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 303199..303218 + + 0 . : . : + ACCTCCACCCTTCGAGCTCA + ||||||| |||||||||||| + ACCTCCA-CCTTCGAGCTCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 471177..471157 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 32619..32640 + + 0 . : . : + ATAAACCAGCCTG-AACTAGA + ||||||||||||| ||||||| + ATAAACCAGCCTGGAACTAGA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 455383..455363 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 669816..669837 + + 0 . : . : + TACTA-AAGTCTGTTAAGGCG + ||||| ||||||||||||||| + TACTACAAGTCTGTTAAGGCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 442627..442607 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 641356..641375 + + 0 . : . : + CTTGGGGTAGGGGGAGCTCT + ||||||| |||||||||||| + CTTGGGG-AGGGGGAGCTCT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 428989..428969 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 791788..791807 + + 0 . : . : + CCTCCGCGAAGGACTCCATC + |||||| ||||||||||||| + CCTCCG-GAAGGACTCCATC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 408733..408713 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 574103..574123 + + 0 . : . : + TTGTGGTTCAATTTCGTGGC + ||||| |||||||||||||| + TTGTGATTCAATTTCGTGGC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 404469..404449 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 236051..236070 + + 0 . : . : + AGCGTGAGCTCCGTAGAGGA + |||| ||||||||||||||| + AGCG-GAGCTCCGTAGAGGA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 402111..402091 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 928806..928825 + + 0 . : . : + CTTCACTTATGTGGTAAAAT + |||||| ||||||||||||| + CTTCAC-TATGTGGTAAAAT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 377568..377548 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 920263..920282 + + 0 . : . : + GCGGCTCGACTTTGAGATCA + |||||||||||||||| ||| + GCGGCTCGACTTTGAG-TCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 362982..362962 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 618521..618540 + + 0 . : . : + ACAGCTCGTTTCTGTCCCAC + |||||||||||| ||||||| + ACAGCTCGTTTC-GTCCCAC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 316466..316446 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 962071..962092 + + 0 . : . : + GGTAAGTTC-AGCTTGAAACG + ||||||||| ||||||||||| + GGTAAGTTCAAGCTTGAAACG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 277774..277754 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 477269..477288 + + 0 . : . : + TATGATTTCGGCGGATTATA + ||||||||||||| |||||| + TATGATTTCGGCG-ATTATA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 252333..252313 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 518920..518939 + + 0 . : . : + CTTCGCACACACATGGTCTT + ||||||||| |||||||||| + CTTCGCACA-ACATGGTCTT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 238287..238267 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 83581..83602 + + 0 . : . : + GTCAATCTTGAATATA-CTTG + |||||||||||||||| |||| + GTCAATCTTGAATATAACTTG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 123850..123830 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 791400..791420 + + 0 . : . : + CCCATATATGCAGCCCAGTT + |||||||||||||| ||||| + CCCATATATGCAGCACAGTT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 90457..90437 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 877933..877953 + + 0 . : . : + CAAGAAAACAGAGAGGGCTC + ||||| |||||||||||||| + CAAGAGAACAGAGAGGGCTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 77415..77395 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 785297..785316 + + 0 . : . : + TCACCTTAATGTGTTCTGCA + ||||||| |||||||||||| + TCACCTT-ATGTGTTCTGCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 66395..66375 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 842456..842475 + + 0 . : . : + CACGGTTATCAACCTTTGCC + ||||||||| |||||||||| + CACGGTTAT-AACCTTTGCC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 60990..60970 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 760086..760106 + + 0 . : . : + AGTTGCCTACGGATACAACT + |||||||||||| ||||||| + AGTTGCCTACGGTTACAACT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 53917..53897 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 949224..949244 + + 0 . : . : + CCTGCCTAGTTTCACGCTCA + ||||||||| |||||||||| + CCTGCCTAGGTTCACGCTCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 31283..31263 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 993851..993871 + + 0 . : . : + TAATGTGCCAGATAATAGTC + ||| |||||||||||||||| + TAACGTGCCAGATAATAGTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 26522..26502 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 842545..842565 + + 0 . : . : + ATATTCTAGCCGTCTCCAGG + |||||||| ||||||||||| + ATATTCTACCCGTCTCCAGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 19477..19457 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 30537..30556 + + 0 . : . : + ATAGTTCAGCATAGGGGTCT + ||||||| |||||||||||| + ATAGTTC-GCATAGGGGTCT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 954662..954643 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 566181..566201 + + 0 . : . : + TAG-GACGACTGTTAGGCCA + ||| |||||||||||||||| + TAGAGACGACTGTTAGGCCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 905595..905576 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 528543..528563 + + 0 . : . : + ATTCTTTCTCTCTCTT-CGG + |||||||||||||||| ||| + ATTCTTTCTCTCTCTTGCGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 805488..805469 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 41850..41870 + + 0 . : . : + ACATTCCTACTC-ACCGTGC + |||||||||||| ||||||| + ACATTCCTACTCGACCGTGC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 795208..795189 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 361613..361633 + + 0 . : . : + CCGTATGA-AGTTAAGTGTC + |||||||| ||||||||||| + CCGTATGATAGTTAAGTGTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 763565..763546 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 469455..469475 + + 0 . : . : + CTCA-TTAAATGGCGCCATC + |||| ||||||||||||||| + CTCACTTAAATGGCGCCATC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 732520..732501 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 446090..446110 + + 0 . : . : + GA-TTCGGCCACTTCTATGG + || ||||||||||||||||| + GACTTCGGCCACTTCTATGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 643840..643821 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 244583..244603 + + 0 . : . : + GACTTCATGGGAGCTT-CAT + |||||||||||||||| ||| + GACTTCATGGGAGCTTTCAT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 603144..603125 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 192739..192759 + + 0 . : . : + TGCGCGAAATGGGTC-CACG + ||||||||||||||| |||| + TGCGCGAAATGGGTCGCACG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 508123..508104 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 400206..400226 + + 0 . : . : + TCCCACCA-TATCCTCAGTC + |||||||| ||||||||||| + TCCCACCATTATCCTCAGTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 447414..447395 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 828386..828406 + + 0 . : . : + TGCCACCGTGCTTCCCT-TG + ||||||||||||||||| || + TGCCACCGTGCTTCCCTATG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 424870..424851 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 882878..882898 + + 0 . : . : + CACCGTACAAGAC-AATATT + ||||||||||||| |||||| + CACCGTACAAGACTAATATT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 418619..418600 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 159495..159515 + + 0 . : . : + AGGGTCGGCCGTTCTGG-CG + ||||||||||||||||| || + AGGGTCGGCCGTTCTGGGCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 399083..399064 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 932044..932064 + + 0 . : . : + GCTGTAAAGACCCC-TACGT + |||||||||||||| ||||| + GCTGTAAAGACCCCCTACGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 346410..346391 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 246661..246681 + + 0 . : . : + TCTT-GGTCTAGGGCGTGGT + |||| ||||||||||||||| + TCTTGGGTCTAGGGCGTGGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 313473..313454 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 993115..993135 + + 0 . : . : + TATGGT-TCACCAATCTTGT + |||||| ||||||||||||| + TATGGTGTCACCAATCTTGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 303085..303066 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 935084..935104 + + 0 . : . : + CACGTAAA-CTCTAAATTCC + |||||||| ||||||||||| + CACGTAAACCTCTAAATTCC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 293491..293472 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 421801..421821 + + 0 . : . : + CGACGTCCGCGGCATGAT-T + |||||||||||||||||| | + CGACGTCCGCGGCATGATCT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 285973..285954 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 517349..517369 + + 0 . : . : + GTAAA-AGCTGCTGTTGGTT + ||||| |||||||||||||| + GTAAAGAGCTGCTGTTGGTT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 273625..273606 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 9350..9370 + + 0 . : . : + CAAATA-GGGGTATGACAAT + |||||| ||||||||||||| + CAAATAGGGGGTATGACAAT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 216623..216604 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 669252..669272 + + 0 . : . : + TGCAGTAAA-AATATACACC + ||||||||| |||||||||| + TGCAGTAAATAATATACACC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 71342..71323 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 218248..218268 + + 0 . : . : + TGCCCT-CCTTTAGCACCCG + |||||| ||||||||||||| + TGCCCTCCCTTTAGCACCCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 3349..3330 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 631236..631256 + + 0 . : . : + CGGACCTGG-AGTGTCTGAC + ||||||||| |||||||||| + CGGACCTGGTAGTGTCTGAC + + +---------------------------------------------------------------------- + diff --git a/test/cli/gold_standard/dna_reverse/minLen20.gff b/test/cli/gold_standard/dna_reverse/minLen20.gff index e69de29bb..55cf3fae4 100644 --- a/test/cli/gold_standard/dna_reverse/minLen20.gff +++ b/test/cli/gold_standard/dna_reverse/minLen20.gff @@ -0,0 +1,97 @@ +seq1 Stellar eps-matches 765718 765740 95.6521 - . seq2;seq2Range=480465,480487;cigar=23M;mutations=15T +seq1 Stellar eps-matches 259870 259892 95.6521 - . seq2;seq2Range=242055,242077;cigar=23M;mutations=14T +seq1 Stellar eps-matches 986804 986825 95.4545 - . seq2;seq2Range=255704,255725;cigar=22M;mutations=9G +seq1 Stellar eps-matches 795116 795137 95.4545 - . seq2;seq2Range=100168,100188;cigar=15M1D6M;mutations= +seq1 Stellar eps-matches 590053 590074 95.4545 - . seq2;seq2Range=245228,245249;cigar=22M;mutations=17A +seq1 Stellar eps-matches 563210 563231 95.4545 - . seq2;seq2Range=355253,355274;cigar=22M;mutations=4A +seq1 Stellar eps-matches 353237 353258 95.4545 - . seq2;seq2Range=905108,905129;cigar=22M;mutations=3G +seq1 Stellar eps-matches 971517 971537 95.4545 - . seq2;seq2Range=55261,55282;cigar=9M1I12M;mutations=10T +seq1 Stellar eps-matches 960203 960223 95.238 - . seq2;seq2Range=617579,617598;cigar=19M1D1M;mutations= +seq1 Stellar eps-matches 887924 887944 95.238 - . seq2;seq2Range=527360,527380;cigar=21M;mutations=7T +seq1 Stellar eps-matches 814505 814525 95.238 - . seq2;seq2Range=546149,546168;cigar=8M1D12M;mutations= +seq1 Stellar eps-matches 721801 721821 95.4545 - . seq2;seq2Range=739398,739419;cigar=10M1I11M;mutations=11C +seq1 Stellar eps-matches 576858 576878 95.238 - . seq2;seq2Range=491562,491582;cigar=21M;mutations=3C +seq1 Stellar eps-matches 558950 558970 95.238 - . seq2;seq2Range=683633,683652;cigar=12M1D8M;mutations= +seq1 Stellar eps-matches 479518 479538 95.238 - . seq2;seq2Range=83080,83100;cigar=21M;mutations=19C +seq1 Stellar eps-matches 433635 433655 95.238 - . seq2;seq2Range=879014,879033;cigar=2M1D18M;mutations= +seq1 Stellar eps-matches 215664 215684 95.238 - . seq2;seq2Range=212347,212367;cigar=21M;mutations=10G +seq1 Stellar eps-matches 182569 182589 95.238 - . seq2;seq2Range=29771,29791;cigar=21M;mutations=13A +seq1 Stellar eps-matches 995415 995434 95 - . seq2;seq2Range=228306,228325;cigar=20M;mutations=8G +seq1 Stellar eps-matches 930403 930422 95 - . seq2;seq2Range=832382,832400;cigar=15M1D4M;mutations= +seq1 Stellar eps-matches 924969 924988 95 - . seq2;seq2Range=932201,932219;cigar=18M1D1M;mutations= +seq1 Stellar eps-matches 915539 915558 95.238 - . seq2;seq2Range=799352,799372;cigar=8M1I12M;mutations=9C +seq1 Stellar eps-matches 886471 886490 95 - . seq2;seq2Range=337574,337592;cigar=16M1D3M;mutations= +seq1 Stellar eps-matches 875106 875125 95 - . seq2;seq2Range=127516,127535;cigar=20M;mutations=9G +seq1 Stellar eps-matches 860711 860730 95 - . seq2;seq2Range=137910,137928;cigar=6M1D13M;mutations= +seq1 Stellar eps-matches 857248 857267 95.238 - . seq2;seq2Range=911396,911416;cigar=16M1I4M;mutations=17A +seq1 Stellar eps-matches 805362 805381 95 - . seq2;seq2Range=46863,46882;cigar=20M;mutations=4G +seq1 Stellar eps-matches 782079 782098 95 - . seq2;seq2Range=744396,744414;cigar=8M1D11M;mutations= +seq1 Stellar eps-matches 763601 763620 95 - . seq2;seq2Range=245161,245179;cigar=10M1D9M;mutations= +seq1 Stellar eps-matches 761138 761157 95 - . seq2;seq2Range=849566,849585;cigar=20M;mutations=7C +seq1 Stellar eps-matches 715611 715630 95 - . seq2;seq2Range=769274,769293;cigar=20M;mutations=11A +seq1 Stellar eps-matches 698067 698086 95 - . seq2;seq2Range=872432,872451;cigar=20M;mutations=4G +seq1 Stellar eps-matches 687547 687566 95 - . seq2;seq2Range=870365,870383;cigar=10M1D9M;mutations= +seq1 Stellar eps-matches 677977 677996 95 - . seq2;seq2Range=439433,439452;cigar=20M;mutations=11C +seq1 Stellar eps-matches 673645 673664 95 - . seq2;seq2Range=874398,874416;cigar=14M1D5M;mutations= +seq1 Stellar eps-matches 645038 645057 95 - . seq2;seq2Range=234191,234209;cigar=8M1D11M;mutations= +seq1 Stellar eps-matches 644858 644877 95 - . seq2;seq2Range=676277,676295;cigar=7M1D12M;mutations= +seq1 Stellar eps-matches 644813 644832 95 - . seq2;seq2Range=865801,865819;cigar=17M1D2M;mutations= +seq1 Stellar eps-matches 641510 641529 95 - . seq2;seq2Range=607675,607693;cigar=13M1D6M;mutations= +seq1 Stellar eps-matches 620345 620364 95 - . seq2;seq2Range=670708,670727;cigar=20M;mutations=16G +seq1 Stellar eps-matches 589984 590003 95 - . seq2;seq2Range=562730,562749;cigar=20M;mutations=3T +seq1 Stellar eps-matches 584612 584631 95.238 - . seq2;seq2Range=427382,427402;cigar=19M1I1M;mutations=20G +seq1 Stellar eps-matches 579568 579587 95 - . seq2;seq2Range=894330,894348;cigar=19M1D;mutations= +seq1 Stellar eps-matches 563144 563163 95 - . seq2;seq2Range=603709,603727;cigar=13M1D6M;mutations= +seq1 Stellar eps-matches 559003 559022 95.238 - . seq2;seq2Range=545195,545215;cigar=15M1I5M;mutations=16A +seq1 Stellar eps-matches 555342 555361 95 - . seq2;seq2Range=244036,244054;cigar=4M1D15M;mutations= +seq1 Stellar eps-matches 554010 554029 95 - . seq2;seq2Range=500180,500199;cigar=20M;mutations=19A +seq1 Stellar eps-matches 529470 529489 95.238 - . seq2;seq2Range=261759,261779;cigar=15M1I5M;mutations=16G +seq1 Stellar eps-matches 494432 494451 95 - . seq2;seq2Range=109070,109088;cigar=8M1D11M;mutations= +seq1 Stellar eps-matches 486491 486510 95 - . seq2;seq2Range=60990,61009;cigar=20M;mutations=14G +seq1 Stellar eps-matches 479793 479812 95 - . seq2;seq2Range=853178,853197;cigar=20M;mutations=12C +seq1 Stellar eps-matches 478520 478539 95 - . seq2;seq2Range=100372,100390;cigar=12M1D7M;mutations= +seq1 Stellar eps-matches 472552 472571 95 - . seq2;seq2Range=303200,303218;cigar=7M1D12M;mutations= +seq1 Stellar eps-matches 471158 471177 95.238 - . seq2;seq2Range=32620,32640;cigar=13M1I7M;mutations=14G +seq1 Stellar eps-matches 455364 455383 95.238 - . seq2;seq2Range=669817,669837;cigar=5M1I15M;mutations=6C +seq1 Stellar eps-matches 442608 442627 95 - . seq2;seq2Range=641357,641375;cigar=7M1D12M;mutations= +seq1 Stellar eps-matches 428970 428989 95 - . seq2;seq2Range=791789,791807;cigar=6M1D13M;mutations= +seq1 Stellar eps-matches 408714 408733 95 - . seq2;seq2Range=574104,574123;cigar=20M;mutations=6A +seq1 Stellar eps-matches 404450 404469 95 - . seq2;seq2Range=236052,236070;cigar=4M1D15M;mutations= +seq1 Stellar eps-matches 402092 402111 95 - . seq2;seq2Range=928807,928825;cigar=6M1D13M;mutations= +seq1 Stellar eps-matches 377549 377568 95 - . seq2;seq2Range=920264,920282;cigar=16M1D3M;mutations= +seq1 Stellar eps-matches 362963 362982 95 - . seq2;seq2Range=618522,618540;cigar=12M1D7M;mutations= +seq1 Stellar eps-matches 316447 316466 95.238 - . seq2;seq2Range=962072,962092;cigar=9M1I11M;mutations=10A +seq1 Stellar eps-matches 277755 277774 95 - . seq2;seq2Range=477270,477288;cigar=13M1D6M;mutations= +seq1 Stellar eps-matches 252314 252333 95 - . seq2;seq2Range=518921,518939;cigar=9M1D10M;mutations= +seq1 Stellar eps-matches 238268 238287 95.238 - . seq2;seq2Range=83582,83602;cigar=16M1I4M;mutations=17A +seq1 Stellar eps-matches 123831 123850 95 - . seq2;seq2Range=791401,791420;cigar=20M;mutations=15A +seq1 Stellar eps-matches 90438 90457 95 - . seq2;seq2Range=877934,877953;cigar=20M;mutations=6G +seq1 Stellar eps-matches 77396 77415 95 - . seq2;seq2Range=785298,785316;cigar=7M1D12M;mutations= +seq1 Stellar eps-matches 66376 66395 95 - . seq2;seq2Range=842457,842475;cigar=9M1D10M;mutations= +seq1 Stellar eps-matches 60971 60990 95 - . seq2;seq2Range=760087,760106;cigar=20M;mutations=13T +seq1 Stellar eps-matches 53898 53917 95 - . seq2;seq2Range=949225,949244;cigar=20M;mutations=10G +seq1 Stellar eps-matches 31264 31283 95 - . seq2;seq2Range=993852,993871;cigar=20M;mutations=4C +seq1 Stellar eps-matches 26503 26522 95 - . seq2;seq2Range=842546,842565;cigar=20M;mutations=9C +seq1 Stellar eps-matches 19458 19477 95 - . seq2;seq2Range=30538,30556;cigar=7M1D12M;mutations= +seq1 Stellar eps-matches 954644 954662 95 - . seq2;seq2Range=566182,566201;cigar=3M1I16M;mutations=4A +seq1 Stellar eps-matches 905577 905595 95 - . seq2;seq2Range=528544,528563;cigar=16M1I3M;mutations=17G +seq1 Stellar eps-matches 805470 805488 95 - . seq2;seq2Range=41851,41870;cigar=12M1I7M;mutations=13G +seq1 Stellar eps-matches 795190 795208 95 - . seq2;seq2Range=361614,361633;cigar=8M1I11M;mutations=9T +seq1 Stellar eps-matches 763547 763565 95 - . seq2;seq2Range=469456,469475;cigar=4M1I15M;mutations=5C +seq1 Stellar eps-matches 732502 732520 95 - . seq2;seq2Range=446091,446110;cigar=2M1I17M;mutations=3C +seq1 Stellar eps-matches 643822 643840 95 - . seq2;seq2Range=244584,244603;cigar=16M1I3M;mutations=17T +seq1 Stellar eps-matches 603126 603144 95 - . seq2;seq2Range=192740,192759;cigar=15M1I4M;mutations=16G +seq1 Stellar eps-matches 508105 508123 95 - . seq2;seq2Range=400207,400226;cigar=8M1I11M;mutations=9T +seq1 Stellar eps-matches 447396 447414 95 - . seq2;seq2Range=828387,828406;cigar=17M1I2M;mutations=18A +seq1 Stellar eps-matches 424852 424870 95 - . seq2;seq2Range=882879,882898;cigar=13M1I6M;mutations=14T +seq1 Stellar eps-matches 418601 418619 95 - . seq2;seq2Range=159496,159515;cigar=17M1I2M;mutations=18G +seq1 Stellar eps-matches 399065 399083 95 - . seq2;seq2Range=932045,932064;cigar=14M1I5M;mutations=15C +seq1 Stellar eps-matches 346392 346410 95 - . seq2;seq2Range=246662,246681;cigar=4M1I15M;mutations=5G +seq1 Stellar eps-matches 313455 313473 95 - . seq2;seq2Range=993116,993135;cigar=6M1I13M;mutations=7G +seq1 Stellar eps-matches 303067 303085 95 - . seq2;seq2Range=935085,935104;cigar=8M1I11M;mutations=9C +seq1 Stellar eps-matches 293473 293491 95 - . seq2;seq2Range=421802,421821;cigar=18M1I1M;mutations=19C +seq1 Stellar eps-matches 285955 285973 95 - . seq2;seq2Range=517350,517369;cigar=5M1I14M;mutations=6G +seq1 Stellar eps-matches 273607 273625 95 - . seq2;seq2Range=9351,9370;cigar=6M1I13M;mutations=7G +seq1 Stellar eps-matches 216605 216623 95 - . seq2;seq2Range=669253,669272;cigar=9M1I10M;mutations=10T +seq1 Stellar eps-matches 71324 71342 95 - . seq2;seq2Range=218249,218268;cigar=6M1I13M;mutations=7C +seq1 Stellar eps-matches 3331 3349 95 - . seq2;seq2Range=631237,631256;cigar=9M1I10M;mutations=10T diff --git a/test/cli/gold_standard/dna_reverse/minLen20.txt b/test/cli/gold_standard/dna_reverse/minLen20.txt index e69de29bb..323c9a590 100644 --- a/test/cli/gold_standard/dna_reverse/minLen20.txt +++ b/test/cli/gold_standard/dna_reverse/minLen20.txt @@ -0,0 +1,1261 @@ +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 765740..765717 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 480464..480487 + + 0 . : . : + CCATTCAACAGGGGGCAGCAGTC + |||||||||||||| |||||||| + CCATTCAACAGGGGTCAGCAGTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 259892..259869 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 242054..242077 + + 0 . : . : + CTTTAACCCTGGTATATCAGCCT + ||||||||||||| ||||||||| + CTTTAACCCTGGTTTATCAGCCT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 986825..986803 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 255703..255725 + + 0 . : . : + CGTACTCGTCCGGGGTTGCTTA + |||||||| ||||||||||||| + CGTACTCGGCCGGGGTTGCTTA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 795137..795115 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 100167..100188 + + 0 . : . : + TAAGGTCCGTTGATATCGCGCA + ||||||||||||||| |||||| + TAAGGTCCGTTGATA-CGCGCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 590074..590052 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 245227..245249 + + 0 . : . : + GCGGCGCGTAGCTCAAGAGCCC + |||||||||||||||| ||||| + GCGGCGCGTAGCTCAAAAGCCC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 563231..563209 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 355252..355274 + + 0 . : . : + GGCCCCACACTCTTTCAATGCG + ||| |||||||||||||||||| + GGCACCACACTCTTTCAATGCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 353258..353236 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 905107..905129 + + 0 . : . : + AACATTTTAGGATCTGCTTAAC + || ||||||||||||||||||| + AAGATTTTAGGATCTGCTTAAC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 971537..971516 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 55260..55282 + + 0 . : . : + GGTTCAATA-TCCTTAGTAACC + ||||||||| |||||||||||| + GGTTCAATATTCCTTAGTAACC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 960223..960202 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 617578..617598 + + 0 . : . : + CTACGTTCAGCCGCCGTACAC + ||||||||||||||||||| | + CTACGTTCAGCCGCCGTAC-C + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 887944..887923 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 527359..527380 + + 0 . : . : + CATATACAAAGGTCTCCGACA + |||||| |||||||||||||| + CATATATAAAGGTCTCCGACA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 814525..814504 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 546148..546168 + + 0 . : . : + AACTTAAGCAGGACGAAGCGT + |||||||| |||||||||||| + AACTTAAG-AGGACGAAGCGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 721821..721800 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 739397..739419 + + 0 . : . : + TTGGACGATA-TTTCAGATACG + |||||||||| ||||||||||| + TTGGACGATACTTTCAGATACG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 576878..576857 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 491561..491582 + + 0 . : . : + TGAAAGTAGTTAACCAACCGA + || |||||||||||||||||| + TGCAAGTAGTTAACCAACCGA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 558970..558949 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 683632..683652 + + 0 . : . : + ATTATGCTCCTAGAGTGAGGG + |||||||||||| |||||||| + ATTATGCTCCTA-AGTGAGGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 479538..479517 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 83079..83100 + + 0 . : . : + TGGAGTGCGGAGGTCCTGAAC + |||||||||||||||||| || + TGGAGTGCGGAGGTCCTGCAC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 433655..433634 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 879013..879033 + + 0 . : . : + GCAGGTCCGCCGAAGGAGTTC + || |||||||||||||||||| + GC-GGTCCGCCGAAGGAGTTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 215684..215663 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 212346..212367 + + 0 . : . : + GGCGCCCATTCAAGCGTTATA + ||||||||| ||||||||||| + GGCGCCCATGCAAGCGTTATA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 182589..182568 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 29770..29791 + + 0 . : . : + GACGTTGAGATGGCCGACCCG + |||||||||||| |||||||| + GACGTTGAGATGACCGACCCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 995434..995414 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 228305..228325 + + 0 . : . : + CCGAGCAAAAAACGGTCCCA + ||||||| |||||||||||| + CCGAGCAGAAAACGGTCCCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 930422..930402 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 832381..832400 + + 0 . : . : + GCGATGCACCTACGTTGTCA + ||||||||||||||| |||| + GCGATGCACCTACGT-GTCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 924988..924968 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 932200..932219 + + 0 . : . : + CTCTACCTCTTCCACCTAAG + |||||||||||||||||| | + CTCTACCTCTTCCACCTA-G + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 915558..915538 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 799351..799372 + + 0 . : . : + TCATCAAT-CCAGAGTAGCGT + |||||||| |||||||||||| + TCATCAATCCCAGAGTAGCGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 886490..886470 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 337573..337592 + + 0 . : . : + TGTACTGTAAAAAACCGTGC + |||||||||||||||| ||| + TGTACTGTAAAAAACC-TGC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 875125..875105 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 127515..127535 + + 0 . : . : + CGAAGCATTAACAGCATTGT + |||||||| ||||||||||| + CGAAGCATGAACAGCATTGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 860730..860710 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 137909..137928 + + 0 . : . : + CCCGCAGAGTAGAGGGAATA + |||||| ||||||||||||| + CCCGCA-AGTAGAGGGAATA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 857267..857247 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 911395..911416 + + 0 . : . : + GCCATGTCGGCATAGG-CTAC + |||||||||||||||| |||| + GCCATGTCGGCATAGGACTAC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 805381..805361 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 46862..46882 + + 0 . : . : + GTTAAGAGCTGTGGCCTAGA + ||| |||||||||||||||| + GTTGAGAGCTGTGGCCTAGA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 782098..782078 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 744395..744414 + + 0 . : . : + GTCAGTATCCAACAAGTTTG + |||||||| ||||||||||| + GTCAGTAT-CAACAAGTTTG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 763620..763600 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 245160..245179 + + 0 . : . : + CGTAACTCATCTCTTCCTTC + |||||||||| ||||||||| + CGTAACTCAT-TCTTCCTTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 761157..761137 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 849565..849585 + + 0 . : . : + TGTATTATTGGTTCCAAGAG + |||||| ||||||||||||| + TGTATTCTTGGTTCCAAGAG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 715630..715610 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 769273..769293 + + 0 . : . : + CGCGCATTTATCCTGGGGCA + |||||||||| ||||||||| + CGCGCATTTAACCTGGGGCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 698086..698066 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 872431..872451 + + 0 . : . : + GAACAAGCGACGCGGCCGGG + ||| |||||||||||||||| + GAAGAAGCGACGCGGCCGGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 687566..687546 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 870364..870383 + + 0 . : . : + GGTCAAAATGGAATCTGATT + |||||||||| ||||||||| + GGTCAAAATG-AATCTGATT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 677996..677976 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 439432..439452 + + 0 . : . : + CCGGTCCTGTATGTAGCGAA + |||||||||| ||||||||| + CCGGTCCTGTCTGTAGCGAA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 673664..673644 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 874397..874416 + + 0 . : . : + ACCTCAAGTTCGGTATGTGT + |||||||||||||| ||||| + ACCTCAAGTTCGGT-TGTGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 645057..645037 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 234190..234209 + + 0 . : . : + TGCGTGGAGCGTTTCAAAAG + |||||||| ||||||||||| + TGCGTGGA-CGTTTCAAAAG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 644877..644857 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 676276..676295 + + 0 . : . : + GTCAACAGCATGGATGGCGG + ||||||| |||||||||||| + GTCAACA-CATGGATGGCGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 644832..644812 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 865800..865819 + + 0 . : . : + AGGCAATGCGGTGATCGCAT + ||||||||||||||||| || + AGGCAATGCGGTGATCG-AT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 641529..641509 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 607674..607693 + + 0 . : . : + ATCCTGGGAGCCCCGATACT + ||||||||||||| |||||| + ATCCTGGGAGCCC-GATACT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 620364..620344 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 670707..670727 + + 0 . : . : + GCTTGTCTGCATTGTTGCCG + ||||||||||||||| |||| + GCTTGTCTGCATTGTGGCCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 590003..589983 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 562729..562749 + + 0 . : . : + GCGTATCCTCCGGCGTTTAG + || ||||||||||||||||| + GCTTATCCTCCGGCGTTTAG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 584631..584611 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 427381..427402 + + 0 . : . : + ATACGTCCGCTCGTCAGTG-C + ||||||||||||||||||| | + ATACGTCCGCTCGTCAGTGGC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 579587..579567 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 894329..894348 + + 0 . : . : + TATGTTATGTGCGGTGTCAA + ||||||||||||||||||| + TATGTTATGTGCGGTGTCA- + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 563163..563143 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 603708..603727 + + 0 . : . : + CACGCGTCCGGCGCGATGAA + ||||||||||||| |||||| + CACGCGTCCGGCG-GATGAA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 559022..559002 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 545194..545215 + + 0 . : . : + ACACAGTCGGGACAC-TAACT + ||||||||||||||| ||||| + ACACAGTCGGGACACATAACT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 555361..555341 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 244035..244054 + + 0 . : . : + TCTTCGGTGGCGTCAATAGT + |||| ||||||||||||||| + TCTT-GGTGGCGTCAATAGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 554029..554009 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 500179..500199 + + 0 . : . : + GAAGCGGAATAGGACGGTTG + |||||||||||||||||| | + GAAGCGGAATAGGACGGTAG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 529489..529469 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 261758..261779 + + 0 . : . : + GGAAGAGCCGTCGAA-TCATG + ||||||||||||||| ||||| + GGAAGAGCCGTCGAAGTCATG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 494451..494431 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 109069..109088 + + 0 . : . : + CTGTTATATCGCGGACCTCA + |||||||| ||||||||||| + CTGTTATA-CGCGGACCTCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 486510..486490 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 60989..61009 + + 0 . : . : + CTCGGTAATTTGGAAGTAGA + ||||||||||||| |||||| + CTCGGTAATTTGGGAGTAGA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 479812..479792 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 853177..853197 + + 0 . : . : + AGCCAAAGGACAGTGGCTAG + ||||||||||| |||||||| + AGCCAAAGGACCGTGGCTAG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 478539..478519 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 100371..100390 + + 0 . : . : + ATGACTGCGTCCGTGGTTTC + |||||||||||| ||||||| + ATGACTGCGTCC-TGGTTTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 472571..472551 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 303199..303218 + + 0 . : . : + ACCTCCACCCTTCGAGCTCA + ||||||| |||||||||||| + ACCTCCA-CCTTCGAGCTCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 471177..471157 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 32619..32640 + + 0 . : . : + ATAAACCAGCCTG-AACTAGA + ||||||||||||| ||||||| + ATAAACCAGCCTGGAACTAGA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 455383..455363 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 669816..669837 + + 0 . : . : + TACTA-AAGTCTGTTAAGGCG + ||||| ||||||||||||||| + TACTACAAGTCTGTTAAGGCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 442627..442607 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 641356..641375 + + 0 . : . : + CTTGGGGTAGGGGGAGCTCT + ||||||| |||||||||||| + CTTGGGG-AGGGGGAGCTCT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 428989..428969 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 791788..791807 + + 0 . : . : + CCTCCGCGAAGGACTCCATC + |||||| ||||||||||||| + CCTCCG-GAAGGACTCCATC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 408733..408713 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 574103..574123 + + 0 . : . : + TTGTGGTTCAATTTCGTGGC + ||||| |||||||||||||| + TTGTGATTCAATTTCGTGGC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 404469..404449 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 236051..236070 + + 0 . : . : + AGCGTGAGCTCCGTAGAGGA + |||| ||||||||||||||| + AGCG-GAGCTCCGTAGAGGA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 402111..402091 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 928806..928825 + + 0 . : . : + CTTCACTTATGTGGTAAAAT + |||||| ||||||||||||| + CTTCAC-TATGTGGTAAAAT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 377568..377548 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 920263..920282 + + 0 . : . : + GCGGCTCGACTTTGAGATCA + |||||||||||||||| ||| + GCGGCTCGACTTTGAG-TCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 362982..362962 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 618521..618540 + + 0 . : . : + ACAGCTCGTTTCTGTCCCAC + |||||||||||| ||||||| + ACAGCTCGTTTC-GTCCCAC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 316466..316446 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 962071..962092 + + 0 . : . : + GGTAAGTTC-AGCTTGAAACG + ||||||||| ||||||||||| + GGTAAGTTCAAGCTTGAAACG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 277774..277754 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 477269..477288 + + 0 . : . : + TATGATTTCGGCGGATTATA + ||||||||||||| |||||| + TATGATTTCGGCG-ATTATA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 252333..252313 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 518920..518939 + + 0 . : . : + CTTCGCACACACATGGTCTT + ||||||||| |||||||||| + CTTCGCACA-ACATGGTCTT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 238287..238267 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 83581..83602 + + 0 . : . : + GTCAATCTTGAATATA-CTTG + |||||||||||||||| |||| + GTCAATCTTGAATATAACTTG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 123850..123830 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 791400..791420 + + 0 . : . : + CCCATATATGCAGCCCAGTT + |||||||||||||| ||||| + CCCATATATGCAGCACAGTT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 90457..90437 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 877933..877953 + + 0 . : . : + CAAGAAAACAGAGAGGGCTC + ||||| |||||||||||||| + CAAGAGAACAGAGAGGGCTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 77415..77395 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 785297..785316 + + 0 . : . : + TCACCTTAATGTGTTCTGCA + ||||||| |||||||||||| + TCACCTT-ATGTGTTCTGCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 66395..66375 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 842456..842475 + + 0 . : . : + CACGGTTATCAACCTTTGCC + ||||||||| |||||||||| + CACGGTTAT-AACCTTTGCC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 60990..60970 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 760086..760106 + + 0 . : . : + AGTTGCCTACGGATACAACT + |||||||||||| ||||||| + AGTTGCCTACGGTTACAACT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 53917..53897 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 949224..949244 + + 0 . : . : + CCTGCCTAGTTTCACGCTCA + ||||||||| |||||||||| + CCTGCCTAGGTTCACGCTCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 31283..31263 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 993851..993871 + + 0 . : . : + TAATGTGCCAGATAATAGTC + ||| |||||||||||||||| + TAACGTGCCAGATAATAGTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 26522..26502 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 842545..842565 + + 0 . : . : + ATATTCTAGCCGTCTCCAGG + |||||||| ||||||||||| + ATATTCTACCCGTCTCCAGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 19477..19457 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 30537..30556 + + 0 . : . : + ATAGTTCAGCATAGGGGTCT + ||||||| |||||||||||| + ATAGTTC-GCATAGGGGTCT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 954662..954643 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 566181..566201 + + 0 . : . : + TAG-GACGACTGTTAGGCCA + ||| |||||||||||||||| + TAGAGACGACTGTTAGGCCA + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 905595..905576 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 528543..528563 + + 0 . : . : + ATTCTTTCTCTCTCTT-CGG + |||||||||||||||| ||| + ATTCTTTCTCTCTCTTGCGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 805488..805469 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 41850..41870 + + 0 . : . : + ACATTCCTACTC-ACCGTGC + |||||||||||| ||||||| + ACATTCCTACTCGACCGTGC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 795208..795189 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 361613..361633 + + 0 . : . : + CCGTATGA-AGTTAAGTGTC + |||||||| ||||||||||| + CCGTATGATAGTTAAGTGTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 763565..763546 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 469455..469475 + + 0 . : . : + CTCA-TTAAATGGCGCCATC + |||| ||||||||||||||| + CTCACTTAAATGGCGCCATC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 732520..732501 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 446090..446110 + + 0 . : . : + GA-TTCGGCCACTTCTATGG + || ||||||||||||||||| + GACTTCGGCCACTTCTATGG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 643840..643821 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 244583..244603 + + 0 . : . : + GACTTCATGGGAGCTT-CAT + |||||||||||||||| ||| + GACTTCATGGGAGCTTTCAT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 603144..603125 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 192739..192759 + + 0 . : . : + TGCGCGAAATGGGTC-CACG + ||||||||||||||| |||| + TGCGCGAAATGGGTCGCACG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 508123..508104 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 400206..400226 + + 0 . : . : + TCCCACCA-TATCCTCAGTC + |||||||| ||||||||||| + TCCCACCATTATCCTCAGTC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 447414..447395 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 828386..828406 + + 0 . : . : + TGCCACCGTGCTTCCCT-TG + ||||||||||||||||| || + TGCCACCGTGCTTCCCTATG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 424870..424851 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 882878..882898 + + 0 . : . : + CACCGTACAAGAC-AATATT + ||||||||||||| |||||| + CACCGTACAAGACTAATATT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 418619..418600 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 159495..159515 + + 0 . : . : + AGGGTCGGCCGTTCTGG-CG + ||||||||||||||||| || + AGGGTCGGCCGTTCTGGGCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 399083..399064 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 932044..932064 + + 0 . : . : + GCTGTAAAGACCCC-TACGT + |||||||||||||| ||||| + GCTGTAAAGACCCCCTACGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 346410..346391 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 246661..246681 + + 0 . : . : + TCTT-GGTCTAGGGCGTGGT + |||| ||||||||||||||| + TCTTGGGTCTAGGGCGTGGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 313473..313454 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 993115..993135 + + 0 . : . : + TATGGT-TCACCAATCTTGT + |||||| ||||||||||||| + TATGGTGTCACCAATCTTGT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 303085..303066 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 935084..935104 + + 0 . : . : + CACGTAAA-CTCTAAATTCC + |||||||| ||||||||||| + CACGTAAACCTCTAAATTCC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 293491..293472 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 421801..421821 + + 0 . : . : + CGACGTCCGCGGCATGAT-T + |||||||||||||||||| | + CGACGTCCGCGGCATGATCT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 285973..285954 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 517349..517369 + + 0 . : . : + GTAAA-AGCTGCTGTTGGTT + ||||| |||||||||||||| + GTAAAGAGCTGCTGTTGGTT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 273625..273606 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 9350..9370 + + 0 . : . : + CAAATA-GGGGTATGACAAT + |||||| ||||||||||||| + CAAATAGGGGGTATGACAAT + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 216623..216604 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 669252..669272 + + 0 . : . : + TGCAGTAAA-AATATACACC + ||||||||| |||||||||| + TGCAGTAAATAATATACACC + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 71342..71323 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 218248..218268 + + 0 . : . : + TGCCCT-CCTTTAGCACCCG + |||||| ||||||||||||| + TGCCCTCCCTTTAGCACCCG + + +---------------------------------------------------------------------- + +Database sequence: seq1 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 (complement) +Database positions: 3349..3330 +Query sequence: seq2 length=1000000 numMatches=500 errorRate=0.05 matchMinLength=50 matchMaxLength=200 +Query positions: 631236..631256 + + 0 . : . : + CGGACCTGG-AGTGTCTGAC + ||||||||| |||||||||| + CGGACCTGGTAGTGTCTGAC + + +---------------------------------------------------------------------- +