Usage:
  last-pair-probs --help
  last-pair-probs [options] interleaved-alignments
  last-pair-probs [options] alignments1 alignments2

Read alignments of paired DNA reads to a genome, and: (1) estimate the
distribution of distances between paired reads, (2) estimate the probability
that each alignment represents the genomic source of the read.

Options:
  -h, --help            show this help message and exit
  -r, --rna             assume the reads are from potentially-spliced RNA
  -e, --estdist         just estimate the distribution of distances
  -m M, --mismap=M      don't write alignments with mismap probability > M
                        (default: 0.01)
  -f BP, --fraglen=BP   mean distance in bp
  -s BP, --sdev=BP      standard deviation of distance
  -d PROB, --disjoint=PROB
                        prior probability of disjoint mapping (default: 0.02
                        if -r, else 0.01)
  -c CHROM, --circular=CHROM
                        specifies that chromosome CHROM is circular (default:
                        chrM)
  -V, --version         show program's version number and exit
# distance sample size: 35
# distance quartiles: 222 248 277
# estimated mean distance: 248
# estimated standard deviation of distance: 40.7715
# fraglen=248 sdev=40.7715 disjoint=0.01 genome=16571
a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.000222
s chrM 11 50 + 16571 TATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTC
s 4/2   0 50 +    50 TATCACCCTATTAACCACTCACGAAAACTCTCCATACATTTAATATTTTC
q 4/2                ?BCC@?CCBCC@CBCACCCCCCCC@ABCCAABBCCBC@A=B@CBCACCC@

a score=180 EG2=0.59 E=3.6e-13 mismap=8.89e-07
s chrM 9112 50 + 16571 CTAATTCTACTGACTATCCTAGAAATCGCTGTCGCCTTAATCCAAGCCTA
s 6/1     0 50 -    50 CTAATTCTACTAACTATCCTAAAAATCGCTATCACCTTAATCCAAACCTA
q 6/1                  FFFDFECECEBCDCA:FAEFDCBD:EEBEAAECCEGGGGCFCDCFGFGGG

a score=220 EG2=6.3e-05 E=2.3e-17 mismap=9.34e-11
s chrM 8903 47 + 16571 CCCACTTCTTACCACAAGGCACACCTACACCCCTTATCCCCATACTA
s 6/2     3 47 +    50 CCCACTTCTTACCACAAAACACACCTACACCCCTTATCCCCCTACTA
q 6/2                  ?C<@BA;5AB:AB@B?B<AB>A<AB:@AB9;=@AA43A>@;,5=@@A

a score=166 EG2=15 E=9.9e-12 mismap=2.31e-05
s chrM 12726 38 + 16571 CTTAGTTACCGCTAACAACCTATTCCAACTGTTCATCG
s 7/2      0 38 +    50 CTTAATTACCACTAACAACCTATTCCAACTATTCATCG
q 7/2                   BCC@>CCBCBBCBCBB;B?6,:>@?BB<?:+4>89ABB

a score=180 EG2=0.59 E=3.6e-13 mismap=9.02e-07
s chrM 14235 50 + 16571 CATACAAAGCCCCCGCACCAATAGGATCCTCCCGAATCAACCCTGACCCC
s 9/1      0 50 -    50 CATACAAAACCCCCACACCAATAAAATCCTCCCGAATCAACCCTAACCCC
q 9/1                   FFFFFFFFFFFFFFFDFFFFEEEEEDEFFFFBFFFEEECDFFFFFFFBEF

a score=253 EG2=3.4e-08 E=6e-21 mismap=4.26e-14
s chrM 14024 50 + 16571 CTAAAACAATTTCACAGCACCAAATCTCCACCTCCATCATCACCTCAACC
s 9/2      0 50 +    50 CTAAAACAATTTCACAACACCAAAGCTCCACCTCCATCATCACCTCAACC
q 9/2                   B=9=B+BBBB>BBB=:BB@;<@@B:<A4@?A:?=BABAAAA>>@<A>@AA

a score=252 EG2=4.2e-08 E=7.7e-21 mismap=5.68e-14
s chrM 14182 50 + 16571 TACACCAACAAACAATGTTCAACCAGTAACCACTACTAATCAACGCCCAT
s 13/1     0 50 -    50 TACACCAACAAACAATATTCAACCAATAACCACTACTAATCAACGCCCAT
q 13/1                  D=D?BCDBBDCDCAEEE?DEFFFFFCCCC?CCCC6EEBEDDDDFFFDDFD

a score=202 EG2=0.0039 E=1.8e-15 mismap=5.85e-09
s chrM 13997 49 + 16571 ACCTAACCTGACTAGAAAAGCTATTACCTAAAACAATTTCACAGCACCA
s 13/2     1 49 +    50 ACCTAACCTAACTAAAAAAACTATTACCTAAAACAATTTCACAACACCA
q 13/2                  @BBAB?ACBA@@BAA>ABAB:>AA@9?@AA9=@BB;@>@A8><2?99A@

a score=210 EG2=0.00062 E=2.7e-16 mismap=2.86e-05
s chrM 4992 39 + 16571 CTACGCAAAATCTTAGCATACTCCTCAATTACCCACATA
s 14/2    2 39 +    50 CTACGCAAAATCTTAACATACTCCTCAATTACCCACATA
q 14/2                 CACCC<BBA5=CCCCCB><ACBC@C+BCBBBCC<@B9BB

a score=194 EG2=0.024 E=1.3e-14 mismap=0.00111
s chrM 2765 47 + 16571 AACAAACCCACAGGTCCTAAACTACCAAACCTGCATTAAAAATTTCG
s 15/2    0 47 +    50 AACAGACCCACAAATCCTAAACTACCAAACCTACATTAAAAATTTCG
q 15/2                 2@1>(?@)@C>(>CCC>)@1:15=>9@=ACCC@>A?*=B?8,>4@BC

a score=228 EG2=1e-05 E=3.3e-18 mismap=2.13e-11
s chrM 349 50 + 16571 CCAAACCCCAAAAACAAAGAACCCTAACACCAGCCTAACCAGATTTCAAA
s 16/1   0 50 -    50 CCAaaccccaaaaacaaaaaacccTAACACCAACCTAACCAAATTTCAAA
q 16/1                GGGGFGGGGGEBGGGFGGGGCFFDFGGFGGEGGFEGGGGEGGFGGGDFFB

a score=204 EG2=0.0025 E=1.1e-15 mismap=4.24e-09
s chrM 5387 50 + 16571 CATATCTAACAACGTAAAAATAAAATGACAGTTTGAACATACAAAACCCA
s 18/1    0 50 -    50 CATATCCAACAACGTAAAAATAAAATAacaatTtaaACATACAAAACCCA
q 18/1                 DDCDDFDEEEGEFFFGGFGGGGGGGFFFEEEEGGFGGGBGDGGGGGGGGG

a score=226 EG2=1.6e-05 E=5.4e-18 mismap=2.75e-11
s chrM 5164 50 + 16571 GCACCTGAAACAAGCTAACATGACTAACACCCTTAATTCCATCCACCCTC
s 18/2    0 50 +    50 GCACCTAAAACAAACTAACATAACTAACACCCTTAATTCCATCCACCCTC
q 18/2                 (>CCCCCCBCC=ABCBB=ACCBC3*@B@C=9CBC,CA=-C@CB=CA>9?/

a score=174 EG2=2.3 E=1.5e-12 mismap=3.82e-06
s chrM 8013 41 + 16571 TACTCCCGATTGAAGCCCCCATTCGTATAATAATTACATCA
s 19/1    9 41 -    50 TACTCCCAATTAAAACCCCCATTCGTATAATAATTACATCA
q 19/1                 FGGFGFGFGGGFCGFGFGGGGDGGGGFGGGGDGGGAGFGFF

a score=154 EG2=2.3e+02 E=1.7e-10 mismap=0.000372
s chrM 7826 50 + 16571 CTACGCATCCTTTACATAACAGACGAGGTCAACGATCCCTCCCTTACCAT
s 19/2    0 50 +    50 CTACACATACTTTACATAACAAACGAAATCAAGAATCCCTCCCTTACCAT
q 19/2                 =@9@B>C<';>)*8?'>67@1;(?<;A<@:@@(9A@@)(38>9+@<8?6@

a score=210 EG2=0.00062 E=2.7e-16 mismap=9.12e-10
s chrM 4644 47 + 16571 TATTTCCTCACGCAAGCAACCGCATCCATAATCCTTCTAATAGCTAT
s 20/1    3 47 -    50 TATTTCCTCACGCAAACAACCACATCCATAATCCTTCTAATAACTAT
q 20/1                 FFF?FFF?EEFFEFFFFEEDEEEEDEEEEEDEEEEE5FFFFFGDGGD

a score=207 EG2=0.0012 E=5.5e-16 mismap=1.81e-09
s chrM 4452 50 + 16571 TATACCCTTCCCGTACTAATTAATCCCCTGGCCCAACCCGTCATCTACTC
s 20/2    0 50 +    50 TATACCCTTCCCGTACTAATTAGTCCCCTAACCCAACCCACCATCTACTC
q 20/2                 )@,27@CBC:1?>C?-2AA?59'<)?B7A46)@>:*?8,5)<24?)5>?;

a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.00026
s chrM 607 46 + 16571 ACACTGAAAATGTTTAGACGGGCTCACATCACCCCATAAACAAATA
s 24/2   0 46 +    50 ACACTAAAAATATTTAAACGAACTCACATCACCCCATAAACAAATA
q 24/2                BACBCCCACCBCCCCCCCCCACBCCABABACBCCCCCCBBCCCBCB

a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.000269
s chrM 9781 50 + 16571 TCTACGGCTCAACATTTTTTGTAGCCACAGGCTTCCACGGACTTCACGTC
s 25/1    0 50 -    50 TCTACGACTCAACATTTTTTATAACCACAAACTTCCACGAACTTCACGTC
q 25/1                 @5CCCBEBEEEAEEBEBEE?GGGGGGGGGGFGGGGGGGGGDGFGGA=FEF

a score=199 EG2=0.0077 E=3.8e-15 mismap=1.42e-08
s chrM 9554 49 + 16571 GCCCCCAACAGGCATCACCCCGCTAAATCCCCTAGAAGTCCCACTCCTA
s 25/2    1 49 +    50 GCCCCCAACAGACATCACCCCACTAAATCCCCTACAAATCCCACTACTA
q 25/2                 B=BBCCC@25)=CB9B7;BC7>?A4A;5@BA+;6)()1BCA748@)=1?

a score=240 EG2=6.6e-07 E=1.7e-19 mismap=2.97e-08
s chrM 3458 48 + 16571 ACGCCATAAAACTCTTCACCAAAGAGCCCCTAAAACCCGCCACATCTA
s 26/2    2 48 +    50 ACTCCATAAAACTCTTCACCAaaaaacccctaaaaccCGCCACATCTA
q 26/2                 1@'(=CBBABA,AC664*>0?B@.9,BA<:=?>-B6.)*=>0?59;'3

a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.000595
s chrM 16102 50 + 16571 TACTGCCAGCCACCATGAATATTGTACGGTACCATAAATACTTGACCACC
s 31/1     0 50 -    50 TACTACCAACCACCATAAATATTATACGATAccataaataCTTAACCACC
q 31/1                  GGGGGGGGGGGGAGGGGGGEFFFFFEFFFFEEEEEEEEEDEEEEEECEEE

a score=206 EG2=0.0016 E=7e-16 mismap=6.34e-09
s chrM 15963 50 + 16571 AAAAGTCTTTAACTCCACCATTAGCACCCAAAGCTAAGATTCTAATTTAA
s 31/2     0 50 +    50 AAAAATCTTTAACTCCACCATTAACACCCAAAACTAAAATTCTAATTTAA
q 31/2                  BBBB@,AB>B=BB:@>AA8?B<>7@A@?A@A<@@BAA>?.9,9A?A<@?A

a score=300 EG2=7.3e-13 E=9.6e-27 mismap=0
s chrM 459 50 + 16571 CCTCCCACTCCCATACTACTAATCTCATCAATACAACCCCCGCCCATCCT
s 34/2   0 50 +    50 CCTCCCACTCCCATACTACTAATCTCATCAATACAACCCCCGCCCATCCT
q 34/2                BCABBBBCBBBCAACB;BCAB@BBBBBCABBA?AB>B@5@B@BA5BCB7?

a score=204 EG2=0.0025 E=1.1e-15 mismap=4.31e-09
s chrM 15490 50 + 16571 CCTCCTAGGCGACCCAGACAATTATACCCTAGCCAACCCCTTAAACACCC
s 35/1     0 50 -    50 CCTCCTAAACGACCCAAACAATTATACCCTAACCAACCCCTTAAACACCC
q 35/1                  AC?C<CEBEE@@?CC:A>C>EEEEEE?=DECC=CADECEEEEEEEDEEEE

a score=150 EG2=5.6e+02 E=4.4e-10 mismap=0.00101
s chrM 15318 49 + 16571 CCCTAGCAGCACTCCACCTCCTATTCTTGCACGAAACGGGATCAAACAA
s 35/2     1 49 +    50 CCCTAACAACACTCCACCTCCTATTCTTACACAAAACGAAATCAAACAA
q 35/2                  BC?CAAB@B=CCCCBBB@B@>@@BBACB?@B@CB@CBACC>A4BBACAA

a score=186 EG2=0.15 E=8.5e-14 mismap=2.62e-07
s chrM 9210 43 + 16571 ACCCACCAATCACATGCCTATCATATAGTAAAACCCAGCCCAT
s 36/1    4 43 -    50 ACCCACCAATCACATACCTATCAtataATAAAACCCAACCCAT
q 36/1                 C@DEEFBD7EEE?FFFEEE:EFFFEE??D?B9EB?D5CEEE?C

a score=176 EG2=1.5 E=9.2e-13 mismap=2.59e-06
s chrM 9030 50 + 16571 CCTACTCATGCACCTAATTGGAAGCGCCACCCTAGCAATATCAACCATTA
s 36/2    0 50 +    50 CCTACTCATACACCTAATTAAAAACGCCACCCTAACAATATCAACCATTA
q 36/2                 BB:@@BCB4BBCB?*@A>BB>99BCBBCC>53B<<(>)(<9(87:3B??C

a score=276 EG2=1.8e-10 E=1.2e-23 mismap=0
s chrM 10929 50 + 16571 CCTCCGACCCCCTAACAACCCCCCTCCTAATACTAACTACCTGACTCCTA
s 39/1     0 50 -    50 CCTCCGACCCCCTAACAACCCCCCTCCTAATACTAACTACCTAACTCCTA
q 39/1                  ?CBDD3;@4@B>>,>C7@>>.?=BB>7?4@C8>C3EFFFFFFFFFFDFFF

a score=205 EG2=0.002 E=8.9e-16 mismap=3.73e-09
s chrM 10700 50 + 16571 ACTAGTCTCAATCTCCAACACATATGGCCTAGACTACGTACATAACCTAA
s 39/2     0 50 +    50 ACTAATCTCAATCTCCAACACATATGACCTAAACTACATACATAACCTAA
q 39/2                  B?=BBBBB?B@BABBB=@BABBAAA@BBAA>;0===8>=;7.5:<;@@@A

a score=168 EG2=9.2 E=6.2e-12 mismap=0.000191
s chrM 15071 48 + 16571 TACTCAGAAACCTGAAACATCGGCATTATCCTCCTGCTTGCAACTATA
s 43/1     0 48 -    50 TACTCAAAAACCTAAAACATCGACATTATCCTCCTACTTACAACTATA
q 43/1                  FBGGGFGGGGGEEGGGFGEEGGGEGBDGEFGGFFGGGFGFGGGGGGGG

a score=281 EG2=5.6e-11 E=2.9e-24 mismap=0
s chrM 14777 47 + 16571 ATAAAATTAATTAACCACTCATTCATCGACCTCCCCACCCCATCCAA
s 43/2     0 47 +    50 ATAAAATTAATTAACCACTCATTCATCGACCTCCCCACCCCATCCAA
q 43/2                  BC@BCCCC=B=BBB@B@A@B?B@B9AAB>>BBBA):B<>??05:?7@

a score=264 EG2=2.7e-09 E=3.3e-22 mismap=1.42e-14
s chrM 3543 48 + 16571 CACCATCGCTCTTCTACTATGAACCCCCCTCCCCATACCCAACCCCCT
s 44/1    0 48 -    50 CACCATCGCTCTTCTACTATAAACCCCCCTCCCCATACCCAACCCCCT
q 44/1                 FFFFFFDFFDFFFDFDDDD:FFFFEEEBEEEFFFFGGGEGGGFGGGGG

a score=154 EG2=2.3e+02 E=1.7e-10 mismap=0.00082
s chrM 3288 27 + 16571 AATTCCTCTTCTTAACAACATACCCAT
s 44/2    0 27 +    50 AATTCCTCTTCTTAACAACATACCCAT
q 44/2                 >*>C+3(>:'<CC?B@-7'':0%'@-B

a score=204 EG2=0.0025 E=1.1e-15 mismap=4.61e-09
s chrM 358 50 + 16571 AAAAACAAAGAACCCTAACACCAGCCTAACCAGATTTCAAATTTTATCTT
s 47/1   0 50 -    50 AAAAACAAAAAAGCCTAACACCAACCTAACCAAATTTCAAATTTTATCTT
q 47/1                EE?EEBEEEEBCEEDBEE1EDDEEEEEB6A?=@CC;AAAADDCD@CB;;@

a score=215 EG2=0.0002 E=7.9e-17 mismap=3.71e-10
s chrM 130 50 + 16571 TCTTTGATTCCTGCCTCATTCTATTATTTATCGCACCTACGTTCAATATT
s 47/2   0 50 +    50 TCTTTAATTCCTACCTCATTCTATTATTTATCGCACCAACATTCAATATT
q 47/2                AC>=ABBAA@CB6BCCB9?BAB=CCB@<B;C<CA@C@,9CCBAB>?=BBA

a score=222 EG2=4e-05 E=1.4e-17 mismap=5.74e-11
s chrM 419 41 + 16571 CACTTTTAACAGTCACCCCCCAACTAACACATTATTTTCCC
s 48/1   9 41 -    50 CACTTTTAACAATCACCCCCCAACTAACACATTATTTTCCC
q 48/1                GGGGGFGGGGGGGFDGFFGGEGEGGGEEFFFDGDFGGGGEG

a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.00021
s chrM 207 50 + 16571 TTAATtaattaaTGCTTGTAGGACATAATAATAACAATTGAATGTCTGCA
s 48/2   0 50 +    50 TTAATTAATTAATACTTATAAAACATAAtaataaCAATTAAATATCTGCA
q 48/2                BCBBB@BBBB@BBB?BB>B?BBBAB@>A?B?<AAAB>B9BAB<@A@>ABA

a score=213 EG2=0.00031 E=1.3e-16 mismap=5.32e-10
s chrM 13205 46 + 16571 TCTGCGCCCTTACACAAAATGACATCAAAAAAATCGTAGCCTTCTC
s 55/1     4 46 -    50 TCTACGCCCTTACACAAAATAACATCAAAAAAATCATAGCCTTCTC
q 55/1                  0*3(7%:??47@@@B?BCC9EECCB6;EAEEGFGFGGEFFDGGFDG

a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.000248
s chrM 13027 50 + 16571 CCCTGACTCCCCTCAGCCATAGAAGGCCCCACCCCAGTCTCAGCCCTACT
s 55/2     0 50 +    50 CCCTAACTCCCCTCAACCATAAAAAACCCCACCCCAATCTCAGCCCTACT
q 55/2                  BBCBBBBC@C>B@B?@BBBBABAAB?A@B@BB>BBB?AACA@@BABBBBA

a score=252 EG2=4.2e-08 E=7.7e-21 mismap=5.68e-14
s chrM 10892 50 + 16571 ATCAACAACAACCTATTTAGCTGTTCCCCAACCTTTTCCTCCGACCCCCT
s 59/1     0 50 -    50 ATCAACAACAACCTATTTAACTATTCCCCAACCTTTTCCTCCGACCCCCT
q 59/1                  GEGFGDFFFFEBDBEAGGGFFFEFFGFGGGGGGGGDGGGGEGGGFGGGGG

a score=157 EG2=1.1e+02 E=8.4e-11 mismap=0.000167
s chrM 10691 34 + 16571 CCTAGCCCTACTAGTCTCAATCTCCAACACATAT
s 59/2    16 34 +    50 CCTAACCCTACTAATCTCAATCTCCAACACATAT
q 59/2                  ?A0*:BABB?B?=6><?<BA@5A99>024)?B>A

a score=168 EG2=9.2 E=6.2e-12 mismap=1.37e-05
s chrM 6611 40 + 16571 ATTTTTCGGTCACCCTGAAGTTTATATTCTTATCCTACCA
s 60/1    0 40 -    50 ATTTTTCGATCACCCTAAAATTTATATTCTTATCCTACCA
q 60/1                 FFFFDFFGFFGGGFGGGDGGGFGGGEEEEEEGGGGGGEGG

a score=224 EG2=2.5e-05 E=8.8e-18 mismap=3.69e-11
s chrM 6412 48 + 16571 ATATAAAACCCCCTGCCATAACCCAATACCAAACGCCCCTCTTCGTCT
s 60/2    0 48 +    50 ATATAAAACCCCCTACCGTAACCCAATACCAAACACCCCTCTTCGTCT
q 60/2                 B>B>=)=(@3:87/;7;(;:B8((:+@<)@7<B@:27(<ABAABAAB@

a score=198 EG2=0.0097 E=4.8e-15 mismap=4.95e-08
s chrM 11907 41 + 16571 TAACCACGTTCTCCTGATCAAATATCACTCTCCTACTTACA
s 63/1     1 41 -    50 TAACCACATTCTCCTAATCAAATATCACTCTCCTACTTACA
q 63/1                  C>?;@CAC:>@@>6?DDCDBD@>@CC@CAB1=;5DD:DDDD

a score=200 EG2=0.0061 E=3e-15 mismap=3.13e-08
s chrM 11765 50 + 16571 CGCACTCACAGTCGCATCATAATCCTCTCTCAAGGACTTCAAACTCTACT
s 63/2     0 50 +    50 CGCTCTCCCAATCGCATCATAATCCTCTCTCAAAAAGTTCAAACACTACT
q 63/2                  3)>(<@))@A*<*>(0(94<:0)<+BA:>:7>B8'?%66-?>B5).A6@0

a score=204 EG2=0.0025 E=1.1e-15 mismap=4.69e-09
s chrM 13054 42 + 16571 CCCACCCCAGTCTCAGCCCTACTCCACTCAAGCACTATAGTT
s 64/1     0 42 -    50 CCCACCCCAATCTCAGCCCTACTCCACTCAAACACTATAGTT
q 64/1                  GCDDGDF?GDDGDGG@FDGGAFFDGECCE@GDGDFEGDGBFG

a score=198 EG2=0.0097 E=4.8e-15 mismap=1.85e-08
s chrM 12879 44 + 16571 TTTCATCCTCGCCTTAGCATGATTTATCCTACACTCCAACTCAT
s 64/2     3 44 +    50 TTTCATCCTCACCTTAACATAATTTATCCTACACTCCAACTCAT
q 64/2                  CAC@CCACCCC=CCCC0CCCCCCACCCBCCCCCCC@CCBC>C@C

a score=167 EG2=12 E=7.8e-12 mismap=3.79e-05
s chrM 1767 46 + 16571 TGGCGCAATAGATATAGTACCGCAAGGGAAAGATGAAAAATTATAA
s 65/2    2 46 -    50 TGGCGTAACAGATATAGTATTGTAAGGGAAAGATGAAAAATTATAA
q 65/2                 ??@B@<2;+A:B=>BB/CB?BBBB?@A??>BBB<60:>;>B;B;AA

a score=228 EG2=1e-05 E=3.3e-18 mismap=2.23e-11
s chrM 16146 46 + 16571 ACCACCTGTAGTACATAAAAACCCAACCCACATCAaaccccccccc
s 67/1     4 46 -    50 ACCACCTATAATACATAAAAACCCAACCCACATCAAaccccccccC
q 67/1                  GGGGGFGGGEGGGGGBGGGFGGGGFGGGGGGGGGFGGGGGGFFFFE

a score=258 EG2=1.1e-08 E=1.6e-21 mismap=4.81e-10
s chrM 4154 47 + 16571 ACCAACTCATACACCTCCTATGAAAAAACTTCCTACCACTCACCCTA
s 73/1    2 47 -    50 ACCAACTCATACACCTCCTATAAAAAAACTTCCTACCACTCACCCTA
q 73/1                 BCDBBBBDBBDBCBCCCCDDBDBBCCCC0CCCAFFDGBDDDCDDDD?

a score=210 EG2=0.00062 E=2.7e-16 mismap=9.03e-10
s chrM 12100 47 + 16571 ATCCCTCAACCCCGACATCATTACCGGGTTTTCCTCTTGTAAATATA
s 74/1     0 47 -    50 ATCCCTCAACCCCGACATCATTACCGAATTTTCCTCTTATAAATATA
q 74/1                  FGGGGGGGFGGGGEGFGGGGGGF?GGEGFEGGGGGGGGGGGGGGGGE

a score=153 EG2=2.8e+02 E=2.2e-10 mismap=0.000421
s chrM 11892 48 + 16571 AACTCTCTGTGCTAGTAACCACGTTCTCCTGATCAAATATCACTCTCC
s 74/2     2 48 +    50 AACTCTCTATACTAATAACCACATTCTCCTAATCAAATATCACTGTCC
q 74/2                  C>*3B6)<CBA;=@3AC:B3BBB;A>)@A?CCBBCB:3>)*@B=+??C

a score=228 EG2=1e-05 E=3.3e-18 mismap=5.56e-11
s chrM 4854 46 + 16571 CTTCTTCTCACATGACAAAAACTAGCCCCCATCTCAATCATATACC
s 75/1    4 46 -    50 CTTCTTCTCACATAACAAAAACTAACCCCCATCTCAATCATATACC
q 75/1                 GEGGGGCCEBEEDEEAFGBFEFGGDGFGGGEBFGGFGGDDFFFEFF

a score=199 EG2=0.0077 E=3.8e-15 mismap=4.26e-08
s chrM 4585 44 + 16571 CTTTTATTCCAGTTCTAACCAAAAAAATAAACCCTCGTTCCACA
s 75/2    3 44 +    50 CTTTTATTCCAATTCTAACCAAAAAAATAAAACCTCATTCCACA
q 75/2                 ?CA==A=BCBCABAAA9BBBB;@B?B@B@7@.CA9AACBB==.B

a score=276 EG2=1.8e-10 E=1.2e-23 mismap=7.79e-12
s chrM 5469 50 + 16571 ACGCTACTCCTACCTATCTCCCCTTTTATACTAATAATCTTATAGAAATT
s 78/1    0 50 -    50 ACGCTACTCCTACCTATCTCCCCTTTTATACTAATAATCTTATAAAAATT
q 78/1                 EEEEEE?BAEEBEEEBE:BEEEDDEEE:DBEEEDEEEBEEBEEEEEBDEE

a score=193 EG2=0.03 E=1.6e-14 mismap=0.0014
s chrM 1833 46 + 16571 ACTAACCCCTATACCTTCTGCATAATGAATTAACTAGAAATAACTT
s 82/2    4 46 +    50 ACTAACCCCTATACCTTCTACAAAATAAATTAACTAAAAATAACTT
q 82/2                 B5BCBCBC?CBCBBAA<BAB@4'AC@ABCCBCCBACBA@B@BCB7<

a score=228 EG2=1e-05 E=3.3e-18 mismap=4.63e-07
s chrM 10231 50 + 16571 TAGCTATTACCTTCTTATTATTTGATCTAGAAATTGCCCTCCTTTTACCC
s 84/1     0 50 -    50 TAGCTATTACCTTCTTATTATTTAATCTAAAAATTACCCTCCTTTTACCC
q 84/1                  EEEEEEEEEEAEEEEEEEEEDBDB@EEEEEEEEEEEEEE=EEEBEEAE:D

a score=203 EG2=0.0031 E=1.4e-15 mismap=4.42e-09
s chrM 11684 50 + 16571 TTCACCGGCGCAGTCATTCTCATAATCGCCCACGGACTCACATCCTCATT
s 86/2     0 50 +    50 TTCACCGACACAATCATTCTCATAATCGCCCACGAACTCACATCCTCATT
q 86/2                  ,A;B<1@?BC>48C@B?C9BB.<?BC=:CB=BBC=BAB*):C>=='<CBB

a score=204 EG2=0.0025 E=1.1e-15 mismap=3.54e-09
s chrM 8530 42 + 16571 AACGAAAATCTGTTCGCTTCATTCATTGCCCCCACAATCCTA
s 87/1    0 42 -    50 AACGAAAATCTATTCGCTTCATTCATTACCCCCACAATCCTA
q 87/1                 EGDFEGFGGGGGDGGGFGGGFGGCGGGEGGEFBGGGGFGGGG

a score=175 EG2=1.9 E=1.2e-12 mismap=2.71e-06
s chrM 8329 49 + 16571 TTAAAGATTAAGAGAACCAACACCTCTTTACAGTGAAATGCCCCAACTA
s 87/2    1 49 +    50 TTAAAAATTAAAAGAACCAACACCTCTTTACAATAAAATACCCCAACTA
q 87/2                 ACC@=BBBC<C=C@BBA@C=AB@BB=?@BB?BBBBABACAA4A?BBB?A

a score=210 EG2=0.00062 E=2.7e-16 mismap=2.86e-05
s chrM 1543 47 + 16571 TTATATAGAGGAGACAAGTCGTAACATGGTAAGTGTACTGGAaagtg
s 89/1    0 47 +    50 TTATATAGAGGAGATAAGTCGTAATATGGTAAGTGTATTGGAAAGTG
q 89/1                 EEEBEEEEEAFFFFFEGGGGGGGAEGGADGFGFGAGGGGGGGGEGGA

a score=180 EG2=0.59 E=3.6e-13 mismap=1.1e-06
s chrM 4414 50 + 16571 TAAATAAGCTATCGGGCCCATACCCCGAAAATGTTGGTTATACCCTTCCC
s 91/1    0 50 -    50 TAAATAAACTATCGAACCCATACCCCGAAAATATTGATTATACCCTTCCC
q 91/1                 BEEEEEEEBEEBECACGGFGFAAFFGDDGEGGDGGEGEGGFGGGFGGGGG

a score=183 EG2=0.3 E=1.7e-13 mismap=5.55e-07
s chrM 4245 46 + 16571 CATTCCCCCTCAAACCTAAGAAATATGTCTGATAAAAGAGTTACTT
s 91/2    4 46 +    50 CATTCCCCCTCAAACCTAAGAAATATATCTAATAAAAAAATTACTT
q 91/2                 @6@>7@4:BABBB==B@BBB@(@59A@BAA@B9>B??32B<A<B=;

a score=157 EG2=1.1e+02 E=8.4e-11 mismap=0.000375
s chrM 7286 50 + 16571 TCTAACAGCAGTAATattaataattttcatgatTTGAGAAGCCTTCGCTT
s 95/1    0 50 -    50 TCTAACAACAAtaatattaataattttcataatttaAAAAACCTTCGCTT
q 95/1                 A06+0@@;>@C:@@@;FFFFGFGGFFFFFFFEFECGGGGGFGFGGGEGGD

a score=251 EG2=5.3e-08 E=1e-20 mismap=1.71e-13
s chrM 7140 50 + 16571 CATTTCACTATCATATTCATCGGCGTAAATCTAACTTTCTTCCCACAACA
s 95/2    0 50 +    50 CATTTCACTATCATATTCATCGACATAAATCTAACTTTCTTCCCACAACA
q 95/2                 B@BABA>BB=.=?*@@?2<A?6BA=*?@>>8:?0;=6>@A8?76@7<?@?

a score=180 EG2=0.59 E=3.6e-13 mismap=8.7e-07
s chrM 9555 50 + 16571 CCCCCAACAGGCATCACCCCGCTAAATCCCCTAGAAGTCCCACTCCTAAA
s 98/1    0 50 -    50 CCCCCAACAAACATCACCCCACTAAATCCCCTAAAAATCCCACTCCTAAA
q 98/1                 EDEEBEBCEBFDFEFDGGGGEGGGGFGEGGGBGGGGGGFGGGGGG?EEEE

a score=154 EG2=2.3e+02 E=1.7e-10 mismap=0.000335
s chrM 9350 26 + 16571 AACCAACACACTAACCATATACCAAT
s 98/2    0 26 +    50 AACCAACACACTAACCATATACCAAT
q 98/2                 9BBBA6=?8*<(;BAB@7B+AA<<A<

# distance sample size: 35
# distance quartiles: 222 248 277
# estimated mean distance: 248
# estimated standard deviation of distance: 40.7715
# fraglen=248 sdev=40.7715 disjoint=0.01 genome=16571
156	chrM	11	50	+	16571	4/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000222
180	chrM	9112	50	+	16571	6/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.89e-07
220	chrM	8903	47	+	16571	6/2	3	47	+	50	47	EG2=6.3e-05	E=2.3e-17	mismap=9.34e-11
166	chrM	12726	38	+	16571	7/2	0	38	+	50	38	EG2=15	E=9.9e-12	mismap=2.31e-05
180	chrM	14235	50	+	16571	9/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=9.02e-07
253	chrM	14024	50	+	16571	9/2	0	50	+	50	50	EG2=3.4e-08	E=6e-21	mismap=4.26e-14
252	chrM	14182	50	+	16571	13/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
202	chrM	13997	49	+	16571	13/2	1	49	+	50	49	EG2=0.0039	E=1.8e-15	mismap=5.85e-09
210	chrM	4992	39	+	16571	14/2	2	39	+	50	39	EG2=0.00062	E=2.7e-16	mismap=2.86e-05
228	chrM	349	50	+	16571	16/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=2.13e-11
204	chrM	5387	50	+	16571	18/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.24e-09
226	chrM	5164	50	+	16571	18/2	0	50	+	50	50	EG2=1.6e-05	E=5.4e-18	mismap=2.75e-11
174	chrM	8013	41	+	16571	19/1	9	41	-	50	41	EG2=2.3	E=1.5e-12	mismap=3.82e-06
154	chrM	7826	50	+	16571	19/2	0	50	+	50	50	EG2=2.3e+02	E=1.7e-10	mismap=0.000372
210	chrM	4644	47	+	16571	20/1	3	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.12e-10
207	chrM	4452	50	+	16571	20/2	0	50	+	50	50	EG2=0.0012	E=5.5e-16	mismap=1.81e-09
156	chrM	607	46	+	16571	24/2	0	46	+	50	46	EG2=1.4e+02	E=1.1e-10	mismap=0.00026
156	chrM	9781	50	+	16571	25/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000269
199	chrM	9554	49	+	16571	25/2	1	49	+	50	49	EG2=0.0077	E=3.8e-15	mismap=1.42e-08
240	chrM	3458	48	+	16571	26/2	2	48	+	50	48	EG2=6.6e-07	E=1.7e-19	mismap=2.97e-08
156	chrM	16102	50	+	16571	31/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000595
206	chrM	15963	50	+	16571	31/2	0	50	+	50	50	EG2=0.0016	E=7e-16	mismap=6.34e-09
300	chrM	459	50	+	16571	34/2	0	50	+	50	50	EG2=7.3e-13	E=9.6e-27	mismap=0
204	chrM	15490	50	+	16571	35/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.31e-09
186	chrM	9210	43	+	16571	36/1	4	43	-	50	43	EG2=0.15	E=8.5e-14	mismap=2.62e-07
176	chrM	9030	50	+	16571	36/2	0	50	+	50	50	EG2=1.5	E=9.2e-13	mismap=2.59e-06
276	chrM	10929	50	+	16571	39/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=0
205	chrM	10700	50	+	16571	39/2	0	50	+	50	50	EG2=0.002	E=8.9e-16	mismap=3.73e-09
168	chrM	15071	48	+	16571	43/1	0	48	-	50	48	EG2=9.2	E=6.2e-12	mismap=0.000191
281	chrM	14777	47	+	16571	43/2	0	47	+	50	47	EG2=5.6e-11	E=2.9e-24	mismap=0
264	chrM	3543	48	+	16571	44/1	0	48	-	50	48	EG2=2.7e-09	E=3.3e-22	mismap=1.42e-14
154	chrM	3288	27	+	16571	44/2	0	27	+	50	27	EG2=2.3e+02	E=1.7e-10	mismap=0.00082
204	chrM	358	50	+	16571	47/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.61e-09
215	chrM	130	50	+	16571	47/2	0	50	+	50	50	EG2=0.0002	E=7.9e-17	mismap=3.71e-10
222	chrM	419	41	+	16571	48/1	9	41	-	50	41	EG2=4e-05	E=1.4e-17	mismap=5.74e-11
156	chrM	207	50	+	16571	48/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.00021
213	chrM	13205	46	+	16571	55/1	4	46	-	50	46	EG2=0.00031	E=1.3e-16	mismap=5.32e-10
156	chrM	13027	50	+	16571	55/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000248
252	chrM	10892	50	+	16571	59/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
157	chrM	10691	34	+	16571	59/2	16	34	+	50	34	EG2=1.1e+02	E=8.4e-11	mismap=0.000167
168	chrM	6611	40	+	16571	60/1	0	40	-	50	40	EG2=9.2	E=6.2e-12	mismap=1.37e-05
224	chrM	6412	48	+	16571	60/2	0	48	+	50	48	EG2=2.5e-05	E=8.8e-18	mismap=3.69e-11
198	chrM	11907	41	+	16571	63/1	1	41	-	50	41	EG2=0.0097	E=4.8e-15	mismap=4.95e-08
200	chrM	11765	50	+	16571	63/2	0	50	+	50	50	EG2=0.0061	E=3e-15	mismap=3.13e-08
204	chrM	13054	42	+	16571	64/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=4.69e-09
198	chrM	12879	44	+	16571	64/2	3	44	+	50	44	EG2=0.0097	E=4.8e-15	mismap=1.85e-08
167	chrM	1767	46	+	16571	65/2	2	46	-	50	46	EG2=12	E=7.8e-12	mismap=3.79e-05
228	chrM	16146	46	+	16571	67/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=2.23e-11
258	chrM	4154	47	+	16571	73/1	2	47	-	50	47	EG2=1.1e-08	E=1.6e-21	mismap=4.81e-10
210	chrM	12100	47	+	16571	74/1	0	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.03e-10
153	chrM	11892	48	+	16571	74/2	2	48	+	50	48	EG2=2.8e+02	E=2.2e-10	mismap=0.000421
228	chrM	4854	46	+	16571	75/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=5.56e-11
199	chrM	4585	44	+	16571	75/2	3	44	+	50	44	EG2=0.0077	E=3.8e-15	mismap=4.26e-08
276	chrM	5469	50	+	16571	78/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=7.79e-12
228	chrM	10231	50	+	16571	84/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=4.63e-07
203	chrM	11684	50	+	16571	86/2	0	50	+	50	50	EG2=0.0031	E=1.4e-15	mismap=4.42e-09
204	chrM	8530	42	+	16571	87/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=3.54e-09
175	chrM	8329	49	+	16571	87/2	1	49	+	50	49	EG2=1.9	E=1.2e-12	mismap=2.71e-06
210	chrM	1543	47	+	16571	89/1	0	47	+	50	47	EG2=0.00062	E=2.7e-16	mismap=2.86e-05
180	chrM	4414	50	+	16571	91/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=1.1e-06
183	chrM	4245	46	+	16571	91/2	4	46	+	50	46	EG2=0.3	E=1.7e-13	mismap=5.55e-07
157	chrM	7286	50	+	16571	95/1	0	50	-	50	50	EG2=1.1e+02	E=8.4e-11	mismap=0.000375
251	chrM	7140	50	+	16571	95/2	0	50	+	50	50	EG2=5.3e-08	E=1e-20	mismap=1.71e-13
180	chrM	9555	50	+	16571	98/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.7e-07
154	chrM	9350	26	+	16571	98/2	0	26	+	50	26	EG2=2.3e+02	E=1.7e-10	mismap=0.000335
# distance sample size: 35
# distance quartiles: 222 248 277
# estimated mean ln[distance]: 5.51343
# estimated standard deviation of ln[distance]: 0.16408
# fraglen=5.51343 sdev=0.16408 disjoint=0.02 genome=16571
156	chrM	11	50	+	16571	4/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000213
180	chrM	9112	50	+	16571	6/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=9.4e-07
220	chrM	8903	47	+	16571	6/2	3	47	+	50	47	EG2=6.3e-05	E=2.3e-17	mismap=9.88e-11
166	chrM	12726	38	+	16571	7/2	0	38	+	50	38	EG2=15	E=9.9e-12	mismap=2.49e-05
180	chrM	14235	50	+	16571	9/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=9.6e-07
253	chrM	14024	50	+	16571	9/2	0	50	+	50	50	EG2=3.4e-08	E=6e-21	mismap=5.68e-14
252	chrM	14182	50	+	16571	13/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
202	chrM	13997	49	+	16571	13/2	1	49	+	50	49	EG2=0.0039	E=1.8e-15	mismap=5.63e-09
210	chrM	4992	39	+	16571	14/2	2	39	+	50	39	EG2=0.00062	E=2.7e-16	mismap=1.44e-05
194	chrM	2765	47	+	16571	15/2	0	47	+	50	47	EG2=0.024	E=1.3e-14	mismap=0.00056
228	chrM	349	50	+	16571	16/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=1.97e-11
204	chrM	5387	50	+	16571	18/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.66e-09
226	chrM	5164	50	+	16571	18/2	0	50	+	50	50	EG2=1.6e-05	E=5.4e-18	mismap=3.02e-11
174	chrM	8013	41	+	16571	19/1	9	41	-	50	41	EG2=2.3	E=1.5e-12	mismap=3.6e-06
154	chrM	7826	50	+	16571	19/2	0	50	+	50	50	EG2=2.3e+02	E=1.7e-10	mismap=0.000351
210	chrM	4644	47	+	16571	20/1	3	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=8.91e-10
207	chrM	4452	50	+	16571	20/2	0	50	+	50	50	EG2=0.0012	E=5.5e-16	mismap=1.77e-09
156	chrM	607	46	+	16571	24/2	0	46	+	50	46	EG2=1.4e+02	E=1.1e-10	mismap=0.000241
156	chrM	9781	50	+	16571	25/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000297
199	chrM	9554	49	+	16571	25/2	1	49	+	50	49	EG2=0.0077	E=3.8e-15	mismap=1.57e-08
240	chrM	3458	48	+	16571	26/2	2	48	+	50	48	EG2=6.6e-07	E=1.7e-19	mismap=1.49e-08
156	chrM	16102	50	+	16571	31/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000635
206	chrM	15963	50	+	16571	31/2	0	50	+	50	50	EG2=0.0016	E=7e-16	mismap=6.77e-09
300	chrM	459	50	+	16571	34/2	0	50	+	50	50	EG2=7.3e-13	E=9.6e-27	mismap=0
204	chrM	15490	50	+	16571	35/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.01e-09
150	chrM	15318	49	+	16571	35/2	1	49	+	50	49	EG2=5.6e+02	E=4.4e-10	mismap=0.00094
186	chrM	9210	43	+	16571	36/1	4	43	-	50	43	EG2=0.15	E=8.5e-14	mismap=2.44e-07
176	chrM	9030	50	+	16571	36/2	0	50	+	50	50	EG2=1.5	E=9.2e-13	mismap=2.41e-06
276	chrM	10929	50	+	16571	39/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=0
205	chrM	10700	50	+	16571	39/2	0	50	+	50	50	EG2=0.002	E=8.9e-16	mismap=4.13e-09
168	chrM	15071	48	+	16571	43/1	0	48	-	50	48	EG2=9.2	E=6.2e-12	mismap=0.000127
281	chrM	14777	47	+	16571	43/2	0	47	+	50	47	EG2=5.6e-11	E=2.9e-24	mismap=0
264	chrM	3543	48	+	16571	44/1	0	48	-	50	48	EG2=2.7e-09	E=3.3e-22	mismap=1.42e-14
154	chrM	3288	27	+	16571	44/2	0	27	+	50	27	EG2=2.3e+02	E=1.7e-10	mismap=0.000861
204	chrM	358	50	+	16571	47/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=5.09e-09
215	chrM	130	50	+	16571	47/2	0	50	+	50	50	EG2=0.0002	E=7.9e-17	mismap=4.1e-10
222	chrM	419	41	+	16571	48/1	9	41	-	50	41	EG2=4e-05	E=1.4e-17	mismap=5.94e-11
156	chrM	207	50	+	16571	48/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000218
213	chrM	13205	46	+	16571	55/1	4	46	-	50	46	EG2=0.00031	E=1.3e-16	mismap=4.97e-10
156	chrM	13027	50	+	16571	55/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000232
252	chrM	10892	50	+	16571	59/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
157	chrM	10691	34	+	16571	59/2	16	34	+	50	34	EG2=1.1e+02	E=8.4e-11	mismap=0.000171
168	chrM	6611	40	+	16571	60/1	0	40	-	50	40	EG2=9.2	E=6.2e-12	mismap=1.34e-05
224	chrM	6412	48	+	16571	60/2	0	48	+	50	48	EG2=2.5e-05	E=8.8e-18	mismap=3.61e-11
198	chrM	11907	41	+	16571	63/1	1	41	-	50	41	EG2=0.0097	E=4.8e-15	mismap=5.78e-08
200	chrM	11765	50	+	16571	63/2	0	50	+	50	50	EG2=0.0061	E=3e-15	mismap=3.65e-08
204	chrM	13054	42	+	16571	64/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=4.34e-09
198	chrM	12879	44	+	16571	64/2	3	44	+	50	44	EG2=0.0097	E=4.8e-15	mismap=1.72e-08
167	chrM	1767	46	+	16571	65/2	2	46	-	50	46	EG2=12	E=7.8e-12	mismap=4.04e-05
228	chrM	16146	46	+	16571	67/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=2.46e-11
258	chrM	4154	47	+	16571	73/1	2	47	-	50	47	EG2=1.1e-08	E=1.6e-21	mismap=2.42e-10
210	chrM	12100	47	+	16571	74/1	0	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.41e-10
153	chrM	11892	48	+	16571	74/2	2	48	+	50	48	EG2=2.8e+02	E=2.2e-10	mismap=0.000439
228	chrM	4854	46	+	16571	75/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=5.37e-11
199	chrM	4585	44	+	16571	75/2	3	44	+	50	44	EG2=0.0077	E=3.8e-15	mismap=4.11e-08
276	chrM	5469	50	+	16571	78/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=3.92e-12
193	chrM	1833	46	+	16571	82/2	4	46	+	50	46	EG2=0.03	E=1.6e-14	mismap=0.000704
228	chrM	10231	50	+	16571	84/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=2.33e-07
203	chrM	11684	50	+	16571	86/2	0	50	+	50	50	EG2=0.0031	E=1.4e-15	mismap=4.48e-09
204	chrM	8530	42	+	16571	87/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=3.52e-09
175	chrM	8329	49	+	16571	87/2	1	49	+	50	49	EG2=1.9	E=1.2e-12	mismap=2.7e-06
210	chrM	1543	47	+	16571	89/1	0	47	+	50	47	EG2=0.00062	E=2.7e-16	mismap=1.44e-05
180	chrM	4414	50	+	16571	91/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=1.02e-06
183	chrM	4245	46	+	16571	91/2	4	46	+	50	46	EG2=0.3	E=1.7e-13	mismap=5.14e-07
157	chrM	7286	50	+	16571	95/1	0	50	-	50	50	EG2=1.1e+02	E=8.4e-11	mismap=0.000372
251	chrM	7140	50	+	16571	95/2	0	50	+	50	50	EG2=5.3e-08	E=1e-20	mismap=1.71e-13
180	chrM	9555	50	+	16571	98/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=9.06e-07
154	chrM	9350	26	+	16571	98/2	0	26	+	50	26	EG2=2.3e+02	E=1.7e-10	mismap=0.000349
# fraglen=248 sdev=40.7715 disjoint=0.01 genome=16571
156	chrM	11	50	+	16571	4/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000222
180	chrM	9112	50	+	16571	6/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.89e-07
220	chrM	8903	47	+	16571	6/2	3	47	+	50	47	EG2=6.3e-05	E=2.3e-17	mismap=9.34e-11
166	chrM	12726	38	+	16571	7/2	0	38	+	50	38	EG2=15	E=9.9e-12	mismap=2.31e-05
180	chrM	14235	50	+	16571	9/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=9.02e-07
253	chrM	14024	50	+	16571	9/2	0	50	+	50	50	EG2=3.4e-08	E=6e-21	mismap=4.26e-14
252	chrM	14182	50	+	16571	13/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
202	chrM	13997	49	+	16571	13/2	1	49	+	50	49	EG2=0.0039	E=1.8e-15	mismap=5.85e-09
210	chrM	4992	39	+	16571	14/2	2	39	+	50	39	EG2=0.00062	E=2.7e-16	mismap=2.86e-05
194	chrM	2765	47	+	16571	15/2	0	47	+	50	47	EG2=0.024	E=1.3e-14	mismap=0.00111
228	chrM	349	50	+	16571	16/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=2.13e-11
204	chrM	5387	50	+	16571	18/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.24e-09
226	chrM	5164	50	+	16571	18/2	0	50	+	50	50	EG2=1.6e-05	E=5.4e-18	mismap=2.75e-11
174	chrM	8013	41	+	16571	19/1	9	41	-	50	41	EG2=2.3	E=1.5e-12	mismap=3.82e-06
154	chrM	7826	50	+	16571	19/2	0	50	+	50	50	EG2=2.3e+02	E=1.7e-10	mismap=0.000372
210	chrM	4644	47	+	16571	20/1	3	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.12e-10
207	chrM	4452	50	+	16571	20/2	0	50	+	50	50	EG2=0.0012	E=5.5e-16	mismap=1.81e-09
156	chrM	607	46	+	16571	24/2	0	46	+	50	46	EG2=1.4e+02	E=1.1e-10	mismap=0.00026
156	chrM	9781	50	+	16571	25/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000269
199	chrM	9554	49	+	16571	25/2	1	49	+	50	49	EG2=0.0077	E=3.8e-15	mismap=1.42e-08
240	chrM	3458	48	+	16571	26/2	2	48	+	50	48	EG2=6.6e-07	E=1.7e-19	mismap=2.97e-08
156	chrM	16102	50	+	16571	31/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000595
206	chrM	15963	50	+	16571	31/2	0	50	+	50	50	EG2=0.0016	E=7e-16	mismap=6.34e-09
300	chrM	459	50	+	16571	34/2	0	50	+	50	50	EG2=7.3e-13	E=9.6e-27	mismap=0
204	chrM	15490	50	+	16571	35/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.31e-09
150	chrM	15318	49	+	16571	35/2	1	49	+	50	49	EG2=5.6e+02	E=4.4e-10	mismap=0.00101
186	chrM	9210	43	+	16571	36/1	4	43	-	50	43	EG2=0.15	E=8.5e-14	mismap=2.62e-07
176	chrM	9030	50	+	16571	36/2	0	50	+	50	50	EG2=1.5	E=9.2e-13	mismap=2.59e-06
276	chrM	10929	50	+	16571	39/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=0
205	chrM	10700	50	+	16571	39/2	0	50	+	50	50	EG2=0.002	E=8.9e-16	mismap=3.73e-09
168	chrM	15071	48	+	16571	43/1	0	48	-	50	48	EG2=9.2	E=6.2e-12	mismap=0.000191
281	chrM	14777	47	+	16571	43/2	0	47	+	50	47	EG2=5.6e-11	E=2.9e-24	mismap=0
264	chrM	3543	48	+	16571	44/1	0	48	-	50	48	EG2=2.7e-09	E=3.3e-22	mismap=1.42e-14
154	chrM	3288	27	+	16571	44/2	0	27	+	50	27	EG2=2.3e+02	E=1.7e-10	mismap=0.00082
204	chrM	358	50	+	16571	47/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.61e-09
215	chrM	130	50	+	16571	47/2	0	50	+	50	50	EG2=0.0002	E=7.9e-17	mismap=3.71e-10
222	chrM	419	41	+	16571	48/1	9	41	-	50	41	EG2=4e-05	E=1.4e-17	mismap=5.74e-11
156	chrM	207	50	+	16571	48/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.00021
213	chrM	13205	46	+	16571	55/1	4	46	-	50	46	EG2=0.00031	E=1.3e-16	mismap=5.32e-10
156	chrM	13027	50	+	16571	55/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000248
252	chrM	10892	50	+	16571	59/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
157	chrM	10691	34	+	16571	59/2	16	34	+	50	34	EG2=1.1e+02	E=8.4e-11	mismap=0.000167
168	chrM	6611	40	+	16571	60/1	0	40	-	50	40	EG2=9.2	E=6.2e-12	mismap=1.37e-05
224	chrM	6412	48	+	16571	60/2	0	48	+	50	48	EG2=2.5e-05	E=8.8e-18	mismap=3.69e-11
198	chrM	11907	41	+	16571	63/1	1	41	-	50	41	EG2=0.0097	E=4.8e-15	mismap=4.95e-08
200	chrM	11765	50	+	16571	63/2	0	50	+	50	50	EG2=0.0061	E=3e-15	mismap=3.13e-08
204	chrM	13054	42	+	16571	64/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=4.69e-09
198	chrM	12879	44	+	16571	64/2	3	44	+	50	44	EG2=0.0097	E=4.8e-15	mismap=1.85e-08
167	chrM	1767	46	+	16571	65/2	2	46	-	50	46	EG2=12	E=7.8e-12	mismap=3.79e-05
228	chrM	16146	46	+	16571	67/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=2.23e-11
258	chrM	4154	47	+	16571	73/1	2	47	-	50	47	EG2=1.1e-08	E=1.6e-21	mismap=4.81e-10
210	chrM	12100	47	+	16571	74/1	0	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.03e-10
153	chrM	11892	48	+	16571	74/2	2	48	+	50	48	EG2=2.8e+02	E=2.2e-10	mismap=0.000421
228	chrM	4854	46	+	16571	75/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=5.56e-11
199	chrM	4585	44	+	16571	75/2	3	44	+	50	44	EG2=0.0077	E=3.8e-15	mismap=4.26e-08
276	chrM	5469	50	+	16571	78/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=7.79e-12
193	chrM	1833	46	+	16571	82/2	4	46	+	50	46	EG2=0.03	E=1.6e-14	mismap=0.0014
228	chrM	10231	50	+	16571	84/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=4.63e-07
203	chrM	11684	50	+	16571	86/2	0	50	+	50	50	EG2=0.0031	E=1.4e-15	mismap=4.42e-09
204	chrM	8530	42	+	16571	87/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=3.54e-09
175	chrM	8329	49	+	16571	87/2	1	49	+	50	49	EG2=1.9	E=1.2e-12	mismap=2.71e-06
210	chrM	1543	47	+	16571	89/1	0	47	+	50	47	EG2=0.00062	E=2.7e-16	mismap=2.86e-05
180	chrM	4414	50	+	16571	91/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=1.1e-06
183	chrM	4245	46	+	16571	91/2	4	46	+	50	46	EG2=0.3	E=1.7e-13	mismap=5.55e-07
157	chrM	7286	50	+	16571	95/1	0	50	-	50	50	EG2=1.1e+02	E=8.4e-11	mismap=0.000375
251	chrM	7140	50	+	16571	95/2	0	50	+	50	50	EG2=5.3e-08	E=1e-20	mismap=1.71e-13
180	chrM	9555	50	+	16571	98/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.7e-07
154	chrM	9350	26	+	16571	98/2	0	26	+	50	26	EG2=2.3e+02	E=1.7e-10	mismap=0.000335
# distance sample size: 35
# distance quartiles: 222 248 277
# estimated mean distance: 248
# estimated standard deviation of distance: 40.7715
# fraglen=248 sdev=40.7715 disjoint=0.001 genome=16571
156	chrM	11	50	+	16571	4/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000222
180	chrM	9112	50	+	16571	6/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.89e-07
220	chrM	8903	47	+	16571	6/2	3	47	+	50	47	EG2=6.3e-05	E=2.3e-17	mismap=9.34e-11
166	chrM	12726	38	+	16571	7/2	0	38	+	50	38	EG2=15	E=9.9e-12	mismap=2.31e-05
180	chrM	14235	50	+	16571	9/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=9.02e-07
253	chrM	14024	50	+	16571	9/2	0	50	+	50	50	EG2=3.4e-08	E=6e-21	mismap=4.26e-14
252	chrM	14182	50	+	16571	13/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
202	chrM	13997	49	+	16571	13/2	1	49	+	50	49	EG2=0.0039	E=1.8e-15	mismap=5.85e-09
210	chrM	4992	39	+	16571	14/2	2	39	+	50	39	EG2=0.00062	E=2.7e-16	mismap=0.000288
228	chrM	349	50	+	16571	16/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=2.13e-11
204	chrM	5387	50	+	16571	18/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.24e-09
226	chrM	5164	50	+	16571	18/2	0	50	+	50	50	EG2=1.6e-05	E=5.4e-18	mismap=2.75e-11
174	chrM	8013	41	+	16571	19/1	9	41	-	50	41	EG2=2.3	E=1.5e-12	mismap=3.82e-06
154	chrM	7826	50	+	16571	19/2	0	50	+	50	50	EG2=2.3e+02	E=1.7e-10	mismap=0.000372
210	chrM	4644	47	+	16571	20/1	3	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.12e-10
207	chrM	4452	50	+	16571	20/2	0	50	+	50	50	EG2=0.0012	E=5.5e-16	mismap=1.81e-09
156	chrM	607	46	+	16571	24/2	0	46	+	50	46	EG2=1.4e+02	E=1.1e-10	mismap=0.00026
156	chrM	9781	50	+	16571	25/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000269
199	chrM	9554	49	+	16571	25/2	1	49	+	50	49	EG2=0.0077	E=3.8e-15	mismap=1.42e-08
240	chrM	3458	48	+	16571	26/2	2	48	+	50	48	EG2=6.6e-07	E=1.7e-19	mismap=2.99e-07
156	chrM	16102	50	+	16571	31/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000595
206	chrM	15963	50	+	16571	31/2	0	50	+	50	50	EG2=0.0016	E=7e-16	mismap=6.34e-09
300	chrM	459	50	+	16571	34/2	0	50	+	50	50	EG2=7.3e-13	E=9.6e-27	mismap=0
204	chrM	15490	50	+	16571	35/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.31e-09
150	chrM	15318	49	+	16571	35/2	1	49	+	50	49	EG2=5.6e+02	E=4.4e-10	mismap=0.00101
186	chrM	9210	43	+	16571	36/1	4	43	-	50	43	EG2=0.15	E=8.5e-14	mismap=2.62e-07
176	chrM	9030	50	+	16571	36/2	0	50	+	50	50	EG2=1.5	E=9.2e-13	mismap=2.59e-06
276	chrM	10929	50	+	16571	39/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=0
205	chrM	10700	50	+	16571	39/2	0	50	+	50	50	EG2=0.002	E=8.9e-16	mismap=3.73e-09
168	chrM	15071	48	+	16571	43/1	0	48	-	50	48	EG2=9.2	E=6.2e-12	mismap=0.000191
281	chrM	14777	47	+	16571	43/2	0	47	+	50	47	EG2=5.6e-11	E=2.9e-24	mismap=0
264	chrM	3543	48	+	16571	44/1	0	48	-	50	48	EG2=2.7e-09	E=3.3e-22	mismap=1.42e-14
154	chrM	3288	27	+	16571	44/2	0	27	+	50	27	EG2=2.3e+02	E=1.7e-10	mismap=0.00082
204	chrM	358	50	+	16571	47/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.61e-09
215	chrM	130	50	+	16571	47/2	0	50	+	50	50	EG2=0.0002	E=7.9e-17	mismap=3.71e-10
222	chrM	419	41	+	16571	48/1	9	41	-	50	41	EG2=4e-05	E=1.4e-17	mismap=5.74e-11
156	chrM	207	50	+	16571	48/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.00021
213	chrM	13205	46	+	16571	55/1	4	46	-	50	46	EG2=0.00031	E=1.3e-16	mismap=5.32e-10
156	chrM	13027	50	+	16571	55/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000248
252	chrM	10892	50	+	16571	59/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
157	chrM	10691	34	+	16571	59/2	16	34	+	50	34	EG2=1.1e+02	E=8.4e-11	mismap=0.000167
168	chrM	6611	40	+	16571	60/1	0	40	-	50	40	EG2=9.2	E=6.2e-12	mismap=1.37e-05
224	chrM	6412	48	+	16571	60/2	0	48	+	50	48	EG2=2.5e-05	E=8.8e-18	mismap=3.69e-11
198	chrM	11907	41	+	16571	63/1	1	41	-	50	41	EG2=0.0097	E=4.8e-15	mismap=4.95e-08
200	chrM	11765	50	+	16571	63/2	0	50	+	50	50	EG2=0.0061	E=3e-15	mismap=3.13e-08
204	chrM	13054	42	+	16571	64/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=4.69e-09
198	chrM	12879	44	+	16571	64/2	3	44	+	50	44	EG2=0.0097	E=4.8e-15	mismap=1.86e-08
167	chrM	1767	46	+	16571	65/2	2	46	-	50	46	EG2=12	E=7.8e-12	mismap=3.8e-05
228	chrM	16146	46	+	16571	67/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=2.23e-11
258	chrM	4154	47	+	16571	73/1	2	47	-	50	47	EG2=1.1e-08	E=1.6e-21	mismap=4.85e-09
210	chrM	12100	47	+	16571	74/1	0	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.03e-10
153	chrM	11892	48	+	16571	74/2	2	48	+	50	48	EG2=2.8e+02	E=2.2e-10	mismap=0.000421
228	chrM	4854	46	+	16571	75/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=5.56e-11
199	chrM	4585	44	+	16571	75/2	3	44	+	50	44	EG2=0.0077	E=3.8e-15	mismap=4.26e-08
276	chrM	5469	50	+	16571	78/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=7.86e-11
228	chrM	10231	50	+	16571	84/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=4.67e-06
203	chrM	11684	50	+	16571	86/2	0	50	+	50	50	EG2=0.0031	E=1.4e-15	mismap=4.42e-09
204	chrM	8530	42	+	16571	87/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=3.54e-09
175	chrM	8329	49	+	16571	87/2	1	49	+	50	49	EG2=1.9	E=1.2e-12	mismap=2.71e-06
210	chrM	1543	47	+	16571	89/1	0	47	+	50	47	EG2=0.00062	E=2.7e-16	mismap=0.000288
180	chrM	4414	50	+	16571	91/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=1.1e-06
183	chrM	4245	46	+	16571	91/2	4	46	+	50	46	EG2=0.3	E=1.7e-13	mismap=5.55e-07
157	chrM	7286	50	+	16571	95/1	0	50	-	50	50	EG2=1.1e+02	E=8.4e-11	mismap=0.000375
251	chrM	7140	50	+	16571	95/2	0	50	+	50	50	EG2=5.3e-08	E=1e-20	mismap=1.71e-13
180	chrM	9555	50	+	16571	98/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.7e-07
154	chrM	9350	26	+	16571	98/2	0	26	+	50	26	EG2=2.3e+02	E=1.7e-10	mismap=0.000335
# distance sample size: 35
# distance quartiles: 222 248 277
# estimated mean distance: 248
# estimated standard deviation of distance: 40.7715
# fraglen=248 sdev=40.7715 disjoint=0.01 genome=16571
156	chrM	11	50	+	16571	4/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000222
180	chrM	9112	50	+	16571	6/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.89e-07
220	chrM	8903	47	+	16571	6/2	3	47	+	50	47	EG2=6.3e-05	E=2.3e-17	mismap=9.34e-11
166	chrM	12726	38	+	16571	7/2	0	38	+	50	38	EG2=15	E=9.9e-12	mismap=2.31e-05
180	chrM	14235	50	+	16571	9/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=9.02e-07
253	chrM	14024	50	+	16571	9/2	0	50	+	50	50	EG2=3.4e-08	E=6e-21	mismap=4.26e-14
252	chrM	14182	50	+	16571	13/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
202	chrM	13997	49	+	16571	13/2	1	49	+	50	49	EG2=0.0039	E=1.8e-15	mismap=5.85e-09
210	chrM	4992	39	+	16571	14/2	2	39	+	50	39	EG2=0.00062	E=2.7e-16	mismap=2.86e-05
194	chrM	2765	47	+	16571	15/2	0	47	+	50	47	EG2=0.024	E=1.3e-14	mismap=0.00111
228	chrM	349	50	+	16571	16/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=2.13e-11
204	chrM	5387	50	+	16571	18/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.24e-09
226	chrM	5164	50	+	16571	18/2	0	50	+	50	50	EG2=1.6e-05	E=5.4e-18	mismap=2.75e-11
174	chrM	8013	41	+	16571	19/1	9	41	-	50	41	EG2=2.3	E=1.5e-12	mismap=3.82e-06
154	chrM	7826	50	+	16571	19/2	0	50	+	50	50	EG2=2.3e+02	E=1.7e-10	mismap=0.000372
210	chrM	4644	47	+	16571	20/1	3	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.12e-10
207	chrM	4452	50	+	16571	20/2	0	50	+	50	50	EG2=0.0012	E=5.5e-16	mismap=1.81e-09
156	chrM	607	46	+	16571	24/2	0	46	+	50	46	EG2=1.4e+02	E=1.1e-10	mismap=0.00026
156	chrM	9781	50	+	16571	25/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000269
199	chrM	9554	49	+	16571	25/2	1	49	+	50	49	EG2=0.0077	E=3.8e-15	mismap=1.42e-08
240	chrM	3458	48	+	16571	26/2	2	48	+	50	48	EG2=6.6e-07	E=1.7e-19	mismap=2.97e-08
156	chrM	16102	50	+	16571	31/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000595
206	chrM	15963	50	+	16571	31/2	0	50	+	50	50	EG2=0.0016	E=7e-16	mismap=6.34e-09
300	chrM	459	50	+	16571	34/2	0	50	+	50	50	EG2=7.3e-13	E=9.6e-27	mismap=0
204	chrM	15490	50	+	16571	35/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.31e-09
150	chrM	15318	49	+	16571	35/2	1	49	+	50	49	EG2=5.6e+02	E=4.4e-10	mismap=0.00101
186	chrM	9210	43	+	16571	36/1	4	43	-	50	43	EG2=0.15	E=8.5e-14	mismap=2.62e-07
176	chrM	9030	50	+	16571	36/2	0	50	+	50	50	EG2=1.5	E=9.2e-13	mismap=2.59e-06
276	chrM	10929	50	+	16571	39/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=0
205	chrM	10700	50	+	16571	39/2	0	50	+	50	50	EG2=0.002	E=8.9e-16	mismap=3.73e-09
168	chrM	15071	48	+	16571	43/1	0	48	-	50	48	EG2=9.2	E=6.2e-12	mismap=0.000191
281	chrM	14777	47	+	16571	43/2	0	47	+	50	47	EG2=5.6e-11	E=2.9e-24	mismap=0
264	chrM	3543	48	+	16571	44/1	0	48	-	50	48	EG2=2.7e-09	E=3.3e-22	mismap=1.42e-14
154	chrM	3288	27	+	16571	44/2	0	27	+	50	27	EG2=2.3e+02	E=1.7e-10	mismap=0.00082
204	chrM	358	50	+	16571	47/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.61e-09
215	chrM	130	50	+	16571	47/2	0	50	+	50	50	EG2=0.0002	E=7.9e-17	mismap=3.71e-10
222	chrM	419	41	+	16571	48/1	9	41	-	50	41	EG2=4e-05	E=1.4e-17	mismap=5.74e-11
156	chrM	207	50	+	16571	48/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.00021
213	chrM	13205	46	+	16571	55/1	4	46	-	50	46	EG2=0.00031	E=1.3e-16	mismap=5.32e-10
156	chrM	13027	50	+	16571	55/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000248
252	chrM	10892	50	+	16571	59/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
157	chrM	10691	34	+	16571	59/2	16	34	+	50	34	EG2=1.1e+02	E=8.4e-11	mismap=0.000167
168	chrM	6611	40	+	16571	60/1	0	40	-	50	40	EG2=9.2	E=6.2e-12	mismap=1.37e-05
224	chrM	6412	48	+	16571	60/2	0	48	+	50	48	EG2=2.5e-05	E=8.8e-18	mismap=3.69e-11
198	chrM	11907	41	+	16571	63/1	1	41	-	50	41	EG2=0.0097	E=4.8e-15	mismap=4.95e-08
200	chrM	11765	50	+	16571	63/2	0	50	+	50	50	EG2=0.0061	E=3e-15	mismap=3.13e-08
204	chrM	13054	42	+	16571	64/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=4.69e-09
198	chrM	12879	44	+	16571	64/2	3	44	+	50	44	EG2=0.0097	E=4.8e-15	mismap=1.85e-08
167	chrM	1767	46	+	16571	65/2	2	46	-	50	46	EG2=12	E=7.8e-12	mismap=3.79e-05
228	chrM	16146	46	+	16571	67/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=2.23e-11
258	chrM	4154	47	+	16571	73/1	2	47	-	50	47	EG2=1.1e-08	E=1.6e-21	mismap=4.81e-10
210	chrM	12100	47	+	16571	74/1	0	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.03e-10
153	chrM	11892	48	+	16571	74/2	2	48	+	50	48	EG2=2.8e+02	E=2.2e-10	mismap=0.000421
228	chrM	4854	46	+	16571	75/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=5.56e-11
199	chrM	4585	44	+	16571	75/2	3	44	+	50	44	EG2=0.0077	E=3.8e-15	mismap=4.26e-08
276	chrM	5469	50	+	16571	78/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=7.79e-12
193	chrM	1833	46	+	16571	82/2	4	46	+	50	46	EG2=0.03	E=1.6e-14	mismap=0.0014
228	chrM	10231	50	+	16571	84/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=4.63e-07
203	chrM	11684	50	+	16571	86/2	0	50	+	50	50	EG2=0.0031	E=1.4e-15	mismap=4.42e-09
204	chrM	8530	42	+	16571	87/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=3.54e-09
175	chrM	8329	49	+	16571	87/2	1	49	+	50	49	EG2=1.9	E=1.2e-12	mismap=2.71e-06
210	chrM	1543	47	+	16571	89/1	0	47	+	50	47	EG2=0.00062	E=2.7e-16	mismap=2.86e-05
180	chrM	4414	50	+	16571	91/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=1.1e-06
183	chrM	4245	46	+	16571	91/2	4	46	+	50	46	EG2=0.3	E=1.7e-13	mismap=5.55e-07
157	chrM	7286	50	+	16571	95/1	0	50	-	50	50	EG2=1.1e+02	E=8.4e-11	mismap=0.000375
251	chrM	7140	50	+	16571	95/2	0	50	+	50	50	EG2=5.3e-08	E=1e-20	mismap=1.71e-13
180	chrM	9555	50	+	16571	98/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.7e-07
154	chrM	9350	26	+	16571	98/2	0	26	+	50	26	EG2=2.3e+02	E=1.7e-10	mismap=0.000335
# distance sample size: 35
# distance quartiles: 222 248 277
# estimated mean distance: 248
# estimated standard deviation of distance: 40.7715
# fraglen=248 sdev=40.7715 disjoint=0.01 genome=16571
156	chrM	11	50	+	16571	4/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000222
180	chrM	9112	50	+	16571	6/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.89e-07
220	chrM	8903	47	+	16571	6/2	3	47	+	50	47	EG2=6.3e-05	E=2.3e-17	mismap=9.34e-11
166	chrM	12726	38	+	16571	7/2	0	38	+	50	38	EG2=15	E=9.9e-12	mismap=2.31e-05
180	chrM	14235	50	+	16571	9/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=9.02e-07
253	chrM	14024	50	+	16571	9/2	0	50	+	50	50	EG2=3.4e-08	E=6e-21	mismap=4.26e-14
252	chrM	14182	50	+	16571	13/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
202	chrM	13997	49	+	16571	13/2	1	49	+	50	49	EG2=0.0039	E=1.8e-15	mismap=5.85e-09
210	chrM	4992	39	+	16571	14/2	2	39	+	50	39	EG2=0.00062	E=2.7e-16	mismap=2.86e-05
194	chrM	2765	47	+	16571	15/2	0	47	+	50	47	EG2=0.024	E=1.3e-14	mismap=0.00111
228	chrM	349	50	+	16571	16/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=2.13e-11
204	chrM	5387	50	+	16571	18/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.24e-09
226	chrM	5164	50	+	16571	18/2	0	50	+	50	50	EG2=1.6e-05	E=5.4e-18	mismap=2.75e-11
174	chrM	8013	41	+	16571	19/1	9	41	-	50	41	EG2=2.3	E=1.5e-12	mismap=3.82e-06
154	chrM	7826	50	+	16571	19/2	0	50	+	50	50	EG2=2.3e+02	E=1.7e-10	mismap=0.000372
210	chrM	4644	47	+	16571	20/1	3	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.12e-10
207	chrM	4452	50	+	16571	20/2	0	50	+	50	50	EG2=0.0012	E=5.5e-16	mismap=1.81e-09
156	chrM	607	46	+	16571	24/2	0	46	+	50	46	EG2=1.4e+02	E=1.1e-10	mismap=0.00026
156	chrM	9781	50	+	16571	25/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000269
199	chrM	9554	49	+	16571	25/2	1	49	+	50	49	EG2=0.0077	E=3.8e-15	mismap=1.42e-08
240	chrM	3458	48	+	16571	26/2	2	48	+	50	48	EG2=6.6e-07	E=1.7e-19	mismap=2.97e-08
156	chrM	16102	50	+	16571	31/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000595
206	chrM	15963	50	+	16571	31/2	0	50	+	50	50	EG2=0.0016	E=7e-16	mismap=6.34e-09
300	chrM	459	50	+	16571	34/2	0	50	+	50	50	EG2=7.3e-13	E=9.6e-27	mismap=0
204	chrM	15490	50	+	16571	35/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.31e-09
150	chrM	15318	49	+	16571	35/2	1	49	+	50	49	EG2=5.6e+02	E=4.4e-10	mismap=0.00101
186	chrM	9210	43	+	16571	36/1	4	43	-	50	43	EG2=0.15	E=8.5e-14	mismap=2.62e-07
176	chrM	9030	50	+	16571	36/2	0	50	+	50	50	EG2=1.5	E=9.2e-13	mismap=2.59e-06
276	chrM	10929	50	+	16571	39/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=0
205	chrM	10700	50	+	16571	39/2	0	50	+	50	50	EG2=0.002	E=8.9e-16	mismap=3.73e-09
168	chrM	15071	48	+	16571	43/1	0	48	-	50	48	EG2=9.2	E=6.2e-12	mismap=0.000191
281	chrM	14777	47	+	16571	43/2	0	47	+	50	47	EG2=5.6e-11	E=2.9e-24	mismap=0
264	chrM	3543	48	+	16571	44/1	0	48	-	50	48	EG2=2.7e-09	E=3.3e-22	mismap=1.42e-14
154	chrM	3288	27	+	16571	44/2	0	27	+	50	27	EG2=2.3e+02	E=1.7e-10	mismap=0.00082
204	chrM	358	50	+	16571	47/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.61e-09
215	chrM	130	50	+	16571	47/2	0	50	+	50	50	EG2=0.0002	E=7.9e-17	mismap=3.71e-10
222	chrM	419	41	+	16571	48/1	9	41	-	50	41	EG2=4e-05	E=1.4e-17	mismap=5.74e-11
156	chrM	207	50	+	16571	48/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.00021
213	chrM	13205	46	+	16571	55/1	4	46	-	50	46	EG2=0.00031	E=1.3e-16	mismap=5.32e-10
156	chrM	13027	50	+	16571	55/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000248
252	chrM	10892	50	+	16571	59/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
157	chrM	10691	34	+	16571	59/2	16	34	+	50	34	EG2=1.1e+02	E=8.4e-11	mismap=0.000167
168	chrM	6611	40	+	16571	60/1	0	40	-	50	40	EG2=9.2	E=6.2e-12	mismap=1.37e-05
224	chrM	6412	48	+	16571	60/2	0	48	+	50	48	EG2=2.5e-05	E=8.8e-18	mismap=3.69e-11
198	chrM	11907	41	+	16571	63/1	1	41	-	50	41	EG2=0.0097	E=4.8e-15	mismap=4.95e-08
200	chrM	11765	50	+	16571	63/2	0	50	+	50	50	EG2=0.0061	E=3e-15	mismap=3.13e-08
204	chrM	13054	42	+	16571	64/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=4.69e-09
198	chrM	12879	44	+	16571	64/2	3	44	+	50	44	EG2=0.0097	E=4.8e-15	mismap=1.85e-08
167	chrM	1767	46	+	16571	65/2	2	46	-	50	46	EG2=12	E=7.8e-12	mismap=3.79e-05
228	chrM	16146	46	+	16571	67/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=2.23e-11
258	chrM	4154	47	+	16571	73/1	2	47	-	50	47	EG2=1.1e-08	E=1.6e-21	mismap=4.81e-10
210	chrM	12100	47	+	16571	74/1	0	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.03e-10
153	chrM	11892	48	+	16571	74/2	2	48	+	50	48	EG2=2.8e+02	E=2.2e-10	mismap=0.000421
228	chrM	4854	46	+	16571	75/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=5.56e-11
199	chrM	4585	44	+	16571	75/2	3	44	+	50	44	EG2=0.0077	E=3.8e-15	mismap=4.26e-08
276	chrM	5469	50	+	16571	78/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=7.79e-12
193	chrM	1833	46	+	16571	82/2	4	46	+	50	46	EG2=0.03	E=1.6e-14	mismap=0.0014
228	chrM	10231	50	+	16571	84/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=4.63e-07
203	chrM	11684	50	+	16571	86/2	0	50	+	50	50	EG2=0.0031	E=1.4e-15	mismap=4.42e-09
204	chrM	8530	42	+	16571	87/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=3.54e-09
175	chrM	8329	49	+	16571	87/2	1	49	+	50	49	EG2=1.9	E=1.2e-12	mismap=2.71e-06
210	chrM	1543	47	+	16571	89/1	0	47	+	50	47	EG2=0.00062	E=2.7e-16	mismap=2.86e-05
180	chrM	4414	50	+	16571	91/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=1.1e-06
183	chrM	4245	46	+	16571	91/2	4	46	+	50	46	EG2=0.3	E=1.7e-13	mismap=5.55e-07
157	chrM	7286	50	+	16571	95/1	0	50	-	50	50	EG2=1.1e+02	E=8.4e-11	mismap=0.000375
251	chrM	7140	50	+	16571	95/2	0	50	+	50	50	EG2=5.3e-08	E=1e-20	mismap=1.71e-13
180	chrM	9555	50	+	16571	98/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.7e-07
154	chrM	9350	26	+	16571	98/2	0	26	+	50	26	EG2=2.3e+02	E=1.7e-10	mismap=0.000335
# distance sample size: 35
# distance quartiles: 222 248 277
# estimated mean distance: 248
# estimated standard deviation of distance: 40.7715
# fraglen=248 sdev=40.7715 disjoint=0.01 genome=16571
a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.000222
s chrM   11 50 + 16571 TATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTC
s 4/2   0 50 +    50 TATCACCCTATTAACCACTCACGAAAACTCTCCATACATTTAATATTTTC
q 4/2                ?BCC@?CCBCC@CBCACCCCCCCC@ABCCAABBCCBC@A=B@CBCACCC@

a score=180 EG2=0.59 E=3.6e-13 mismap=8.89e-07
s chrM   9112 50 + 16571 CTAATTCTACTGACTATCCTAGAAATCGCTGTCGCCTTAATCCAAGCCTA
s 6/1     0 50 -    50 CTAATTCTACTAACTATCCTAAAAATCGCTATCACCTTAATCCAAACCTA
q 6/1                  FFFDFECECEBCDCA:FAEFDCBD:EEBEAAECCEGGGGCFCDCFGFGGG

a score=220 EG2=6.3e-05 E=2.3e-17 mismap=9.34e-11
s chrM   8903 47 + 16571 CCCACTTCTTACCACAAGGCACACCTACACCCCTTATCCCCATACTA
s 6/2     3 47 +    50 CCCACTTCTTACCACAAAACACACCTACACCCCTTATCCCCCTACTA
q 6/2                  ?C<@BA;5AB:AB@B?B<AB>A<AB:@AB9;=@AA43A>@;,5=@@A

a score=166 EG2=15 E=9.9e-12 mismap=2.31e-05
s chrM   12726 38 + 16571 CTTAGTTACCGCTAACAACCTATTCCAACTGTTCATCG
s 7/2      0 38 +    50 CTTAATTACCACTAACAACCTATTCCAACTATTCATCG
q 7/2                   BCC@>CCBCBBCBCBB;B?6,:>@?BB<?:+4>89ABB

a score=180 EG2=0.59 E=3.6e-13 mismap=9.02e-07
s chrM   14235 50 + 16571 CATACAAAGCCCCCGCACCAATAGGATCCTCCCGAATCAACCCTGACCCC
s 9/1      0 50 -    50 CATACAAAACCCCCACACCAATAAAATCCTCCCGAATCAACCCTAACCCC
q 9/1                   FFFFFFFFFFFFFFFDFFFFEEEEEDEFFFFBFFFEEECDFFFFFFFBEF

a score=253 EG2=3.4e-08 E=6e-21 mismap=4.26e-14
s chrM   14024 50 + 16571 CTAAAACAATTTCACAGCACCAAATCTCCACCTCCATCATCACCTCAACC
s 9/2      0 50 +    50 CTAAAACAATTTCACAACACCAAAGCTCCACCTCCATCATCACCTCAACC
q 9/2                   B=9=B+BBBB>BBB=:BB@;<@@B:<A4@?A:?=BABAAAA>>@<A>@AA

a score=252 EG2=4.2e-08 E=7.7e-21 mismap=5.68e-14
s chrM   14182 50 + 16571 TACACCAACAAACAATGTTCAACCAGTAACCACTACTAATCAACGCCCAT
s 13/1     0 50 -    50 TACACCAACAAACAATATTCAACCAATAACCACTACTAATCAACGCCCAT
q 13/1                  D=D?BCDBBDCDCAEEE?DEFFFFFCCCC?CCCC6EEBEDDDDFFFDDFD

a score=202 EG2=0.0039 E=1.8e-15 mismap=5.85e-09
s chrM   13997 49 + 16571 ACCTAACCTGACTAGAAAAGCTATTACCTAAAACAATTTCACAGCACCA
s 13/2     1 49 +    50 ACCTAACCTAACTAAAAAAACTATTACCTAAAACAATTTCACAACACCA
q 13/2                  @BBAB?ACBA@@BAA>ABAB:>AA@9?@AA9=@BB;@>@A8><2?99A@

a score=210 EG2=0.00062 E=2.7e-16 mismap=2.86e-05
s chrM   4992 39 + 16571 CTACGCAAAATCTTAGCATACTCCTCAATTACCCACATA
s 14/2    2 39 +    50 CTACGCAAAATCTTAACATACTCCTCAATTACCCACATA
q 14/2                 CACCC<BBA5=CCCCCB><ACBC@C+BCBBBCC<@B9BB

a score=194 EG2=0.024 E=1.3e-14 mismap=0.00111
s chrM   2765 47 + 16571 AACAAACCCACAGGTCCTAAACTACCAAACCTGCATTAAAAATTTCG
s 15/2    0 47 +    50 AACAGACCCACAAATCCTAAACTACCAAACCTACATTAAAAATTTCG
q 15/2                 2@1>(?@)@C>(>CCC>)@1:15=>9@=ACCC@>A?*=B?8,>4@BC

a score=228 EG2=1e-05 E=3.3e-18 mismap=2.13e-11
s chrM   349 50 + 16571 CCAAACCCCAAAAACAAAGAACCCTAACACCAGCCTAACCAGATTTCAAA
s 16/1   0 50 -    50 CCAaaccccaaaaacaaaaaacccTAACACCAACCTAACCAAATTTCAAA
q 16/1                GGGGFGGGGGEBGGGFGGGGCFFDFGGFGGEGGFEGGGGEGGFGGGDFFB

a score=204 EG2=0.0025 E=1.1e-15 mismap=4.24e-09
s chrM   5387 50 + 16571 CATATCTAACAACGTAAAAATAAAATGACAGTTTGAACATACAAAACCCA
s 18/1    0 50 -    50 CATATCCAACAACGTAAAAATAAAATAacaatTtaaACATACAAAACCCA
q 18/1                 DDCDDFDEEEGEFFFGGFGGGGGGGFFFEEEEGGFGGGBGDGGGGGGGGG

a score=226 EG2=1.6e-05 E=5.4e-18 mismap=2.75e-11
s chrM   5164 50 + 16571 GCACCTGAAACAAGCTAACATGACTAACACCCTTAATTCCATCCACCCTC
s 18/2    0 50 +    50 GCACCTAAAACAAACTAACATAACTAACACCCTTAATTCCATCCACCCTC
q 18/2                 (>CCCCCCBCC=ABCBB=ACCBC3*@B@C=9CBC,CA=-C@CB=CA>9?/

a score=174 EG2=2.3 E=1.5e-12 mismap=3.82e-06
s chrM   8013 41 + 16571 TACTCCCGATTGAAGCCCCCATTCGTATAATAATTACATCA
s 19/1    9 41 -    50 TACTCCCAATTAAAACCCCCATTCGTATAATAATTACATCA
q 19/1                 FGGFGFGFGGGFCGFGFGGGGDGGGGFGGGGDGGGAGFGFF

a score=154 EG2=2.3e+02 E=1.7e-10 mismap=0.000372
s chrM   7826 50 + 16571 CTACGCATCCTTTACATAACAGACGAGGTCAACGATCCCTCCCTTACCAT
s 19/2    0 50 +    50 CTACACATACTTTACATAACAAACGAAATCAAGAATCCCTCCCTTACCAT
q 19/2                 =@9@B>C<';>)*8?'>67@1;(?<;A<@:@@(9A@@)(38>9+@<8?6@

a score=210 EG2=0.00062 E=2.7e-16 mismap=9.12e-10
s chrM   4644 47 + 16571 TATTTCCTCACGCAAGCAACCGCATCCATAATCCTTCTAATAGCTAT
s 20/1    3 47 -    50 TATTTCCTCACGCAAACAACCACATCCATAATCCTTCTAATAACTAT
q 20/1                 FFF?FFF?EEFFEFFFFEEDEEEEDEEEEEDEEEEE5FFFFFGDGGD

a score=207 EG2=0.0012 E=5.5e-16 mismap=1.81e-09
s chrM   4452 50 + 16571 TATACCCTTCCCGTACTAATTAATCCCCTGGCCCAACCCGTCATCTACTC
s 20/2    0 50 +    50 TATACCCTTCCCGTACTAATTAGTCCCCTAACCCAACCCACCATCTACTC
q 20/2                 )@,27@CBC:1?>C?-2AA?59'<)?B7A46)@>:*?8,5)<24?)5>?;

a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.00026
s chrM   607 46 + 16571 ACACTGAAAATGTTTAGACGGGCTCACATCACCCCATAAACAAATA
s 24/2   0 46 +    50 ACACTAAAAATATTTAAACGAACTCACATCACCCCATAAACAAATA
q 24/2                BACBCCCACCBCCCCCCCCCACBCCABABACBCCCCCCBBCCCBCB

a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.000269
s chrM   9781 50 + 16571 TCTACGGCTCAACATTTTTTGTAGCCACAGGCTTCCACGGACTTCACGTC
s 25/1    0 50 -    50 TCTACGACTCAACATTTTTTATAACCACAAACTTCCACGAACTTCACGTC
q 25/1                 @5CCCBEBEEEAEEBEBEE?GGGGGGGGGGFGGGGGGGGGDGFGGA=FEF

a score=199 EG2=0.0077 E=3.8e-15 mismap=1.42e-08
s chrM   9554 49 + 16571 GCCCCCAACAGGCATCACCCCGCTAAATCCCCTAGAAGTCCCACTCCTA
s 25/2    1 49 +    50 GCCCCCAACAGACATCACCCCACTAAATCCCCTACAAATCCCACTACTA
q 25/2                 B=BBCCC@25)=CB9B7;BC7>?A4A;5@BA+;6)()1BCA748@)=1?

a score=240 EG2=6.6e-07 E=1.7e-19 mismap=2.97e-08
s chrM   3458 48 + 16571 ACGCCATAAAACTCTTCACCAAAGAGCCCCTAAAACCCGCCACATCTA
s 26/2    2 48 +    50 ACTCCATAAAACTCTTCACCAaaaaacccctaaaaccCGCCACATCTA
q 26/2                 1@'(=CBBABA,AC664*>0?B@.9,BA<:=?>-B6.)*=>0?59;'3

a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.000595
s chrM   16102 50 + 16571 TACTGCCAGCCACCATGAATATTGTACGGTACCATAAATACTTGACCACC
s 31/1     0 50 -    50 TACTACCAACCACCATAAATATTATACGATAccataaataCTTAACCACC
q 31/1                  GGGGGGGGGGGGAGGGGGGEFFFFFEFFFFEEEEEEEEEDEEEEEECEEE

a score=206 EG2=0.0016 E=7e-16 mismap=6.34e-09
s chrM   15963 50 + 16571 AAAAGTCTTTAACTCCACCATTAGCACCCAAAGCTAAGATTCTAATTTAA
s 31/2     0 50 +    50 AAAAATCTTTAACTCCACCATTAACACCCAAAACTAAAATTCTAATTTAA
q 31/2                  BBBB@,AB>B=BB:@>AA8?B<>7@A@?A@A<@@BAA>?.9,9A?A<@?A

a score=300 EG2=7.3e-13 E=9.6e-27 mismap=0
s chrM   459 50 + 16571 CCTCCCACTCCCATACTACTAATCTCATCAATACAACCCCCGCCCATCCT
s 34/2   0 50 +    50 CCTCCCACTCCCATACTACTAATCTCATCAATACAACCCCCGCCCATCCT
q 34/2                BCABBBBCBBBCAACB;BCAB@BBBBBCABBA?AB>B@5@B@BA5BCB7?

a score=204 EG2=0.0025 E=1.1e-15 mismap=4.31e-09
s chrM   15490 50 + 16571 CCTCCTAGGCGACCCAGACAATTATACCCTAGCCAACCCCTTAAACACCC
s 35/1     0 50 -    50 CCTCCTAAACGACCCAAACAATTATACCCTAACCAACCCCTTAAACACCC
q 35/1                  AC?C<CEBEE@@?CC:A>C>EEEEEE?=DECC=CADECEEEEEEEDEEEE

a score=150 EG2=5.6e+02 E=4.4e-10 mismap=0.00101
s chrM   15318 49 + 16571 CCCTAGCAGCACTCCACCTCCTATTCTTGCACGAAACGGGATCAAACAA
s 35/2     1 49 +    50 CCCTAACAACACTCCACCTCCTATTCTTACACAAAACGAAATCAAACAA
q 35/2                  BC?CAAB@B=CCCCBBB@B@>@@BBACB?@B@CB@CBACC>A4BBACAA

a score=186 EG2=0.15 E=8.5e-14 mismap=2.62e-07
s chrM   9210 43 + 16571 ACCCACCAATCACATGCCTATCATATAGTAAAACCCAGCCCAT
s 36/1    4 43 -    50 ACCCACCAATCACATACCTATCAtataATAAAACCCAACCCAT
q 36/1                 C@DEEFBD7EEE?FFFEEE:EFFFEE??D?B9EB?D5CEEE?C

a score=176 EG2=1.5 E=9.2e-13 mismap=2.59e-06
s chrM   9030 50 + 16571 CCTACTCATGCACCTAATTGGAAGCGCCACCCTAGCAATATCAACCATTA
s 36/2    0 50 +    50 CCTACTCATACACCTAATTAAAAACGCCACCCTAACAATATCAACCATTA
q 36/2                 BB:@@BCB4BBCB?*@A>BB>99BCBBCC>53B<<(>)(<9(87:3B??C

a score=276 EG2=1.8e-10 E=1.2e-23 mismap=0
s chrM   10929 50 + 16571 CCTCCGACCCCCTAACAACCCCCCTCCTAATACTAACTACCTGACTCCTA
s 39/1     0 50 -    50 CCTCCGACCCCCTAACAACCCCCCTCCTAATACTAACTACCTAACTCCTA
q 39/1                  ?CBDD3;@4@B>>,>C7@>>.?=BB>7?4@C8>C3EFFFFFFFFFFDFFF

a score=205 EG2=0.002 E=8.9e-16 mismap=3.73e-09
s chrM   10700 50 + 16571 ACTAGTCTCAATCTCCAACACATATGGCCTAGACTACGTACATAACCTAA
s 39/2     0 50 +    50 ACTAATCTCAATCTCCAACACATATGACCTAAACTACATACATAACCTAA
q 39/2                  B?=BBBBB?B@BABBB=@BABBAAA@BBAA>;0===8>=;7.5:<;@@@A

a score=168 EG2=9.2 E=6.2e-12 mismap=0.000191
s chrM   15071 48 + 16571 TACTCAGAAACCTGAAACATCGGCATTATCCTCCTGCTTGCAACTATA
s 43/1     0 48 -    50 TACTCAAAAACCTAAAACATCGACATTATCCTCCTACTTACAACTATA
q 43/1                  FBGGGFGGGGGEEGGGFGEEGGGEGBDGEFGGFFGGGFGFGGGGGGGG

a score=281 EG2=5.6e-11 E=2.9e-24 mismap=0
s chrM   14777 47 + 16571 ATAAAATTAATTAACCACTCATTCATCGACCTCCCCACCCCATCCAA
s 43/2     0 47 +    50 ATAAAATTAATTAACCACTCATTCATCGACCTCCCCACCCCATCCAA
q 43/2                  BC@BCCCC=B=BBB@B@A@B?B@B9AAB>>BBBA):B<>??05:?7@

a score=264 EG2=2.7e-09 E=3.3e-22 mismap=1.42e-14
s chrM   3543 48 + 16571 CACCATCGCTCTTCTACTATGAACCCCCCTCCCCATACCCAACCCCCT
s 44/1    0 48 -    50 CACCATCGCTCTTCTACTATAAACCCCCCTCCCCATACCCAACCCCCT
q 44/1                 FFFFFFDFFDFFFDFDDDD:FFFFEEEBEEEFFFFGGGEGGGFGGGGG

a score=154 EG2=2.3e+02 E=1.7e-10 mismap=0.00082
s chrM   3288 27 + 16571 AATTCCTCTTCTTAACAACATACCCAT
s 44/2    0 27 +    50 AATTCCTCTTCTTAACAACATACCCAT
q 44/2                 >*>C+3(>:'<CC?B@-7'':0%'@-B

a score=204 EG2=0.0025 E=1.1e-15 mismap=4.61e-09
s chrM   358 50 + 16571 AAAAACAAAGAACCCTAACACCAGCCTAACCAGATTTCAAATTTTATCTT
s 47/1   0 50 -    50 AAAAACAAAAAAGCCTAACACCAACCTAACCAAATTTCAAATTTTATCTT
q 47/1                EE?EEBEEEEBCEEDBEE1EDDEEEEEB6A?=@CC;AAAADDCD@CB;;@

a score=215 EG2=0.0002 E=7.9e-17 mismap=3.71e-10
s chrM   130 50 + 16571 TCTTTGATTCCTGCCTCATTCTATTATTTATCGCACCTACGTTCAATATT
s 47/2   0 50 +    50 TCTTTAATTCCTACCTCATTCTATTATTTATCGCACCAACATTCAATATT
q 47/2                AC>=ABBAA@CB6BCCB9?BAB=CCB@<B;C<CA@C@,9CCBAB>?=BBA

a score=222 EG2=4e-05 E=1.4e-17 mismap=5.74e-11
s chrM   419 41 + 16571 CACTTTTAACAGTCACCCCCCAACTAACACATTATTTTCCC
s 48/1   9 41 -    50 CACTTTTAACAATCACCCCCCAACTAACACATTATTTTCCC
q 48/1                GGGGGFGGGGGGGFDGFFGGEGEGGGEEFFFDGDFGGGGEG

a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.00021
s chrM   207 50 + 16571 TTAATtaattaaTGCTTGTAGGACATAATAATAACAATTGAATGTCTGCA
s 48/2   0 50 +    50 TTAATTAATTAATACTTATAAAACATAAtaataaCAATTAAATATCTGCA
q 48/2                BCBBB@BBBB@BBB?BB>B?BBBAB@>A?B?<AAAB>B9BAB<@A@>ABA

a score=213 EG2=0.00031 E=1.3e-16 mismap=5.32e-10
s chrM   13205 46 + 16571 TCTGCGCCCTTACACAAAATGACATCAAAAAAATCGTAGCCTTCTC
s 55/1     4 46 -    50 TCTACGCCCTTACACAAAATAACATCAAAAAAATCATAGCCTTCTC
q 55/1                  0*3(7%:??47@@@B?BCC9EECCB6;EAEEGFGFGGEFFDGGFDG

a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.000248
s chrM   13027 50 + 16571 CCCTGACTCCCCTCAGCCATAGAAGGCCCCACCCCAGTCTCAGCCCTACT
s 55/2     0 50 +    50 CCCTAACTCCCCTCAACCATAAAAAACCCCACCCCAATCTCAGCCCTACT
q 55/2                  BBCBBBBC@C>B@B?@BBBBABAAB?A@B@BB>BBB?AACA@@BABBBBA

a score=252 EG2=4.2e-08 E=7.7e-21 mismap=5.68e-14
s chrM   10892 50 + 16571 ATCAACAACAACCTATTTAGCTGTTCCCCAACCTTTTCCTCCGACCCCCT
s 59/1     0 50 -    50 ATCAACAACAACCTATTTAACTATTCCCCAACCTTTTCCTCCGACCCCCT
q 59/1                  GEGFGDFFFFEBDBEAGGGFFFEFFGFGGGGGGGGDGGGGEGGGFGGGGG

a score=157 EG2=1.1e+02 E=8.4e-11 mismap=0.000167
s chrM   10691 34 + 16571 CCTAGCCCTACTAGTCTCAATCTCCAACACATAT
s 59/2    16 34 +    50 CCTAACCCTACTAATCTCAATCTCCAACACATAT
q 59/2                  ?A0*:BABB?B?=6><?<BA@5A99>024)?B>A

a score=168 EG2=9.2 E=6.2e-12 mismap=1.37e-05
s chrM   6611 40 + 16571 ATTTTTCGGTCACCCTGAAGTTTATATTCTTATCCTACCA
s 60/1    0 40 -    50 ATTTTTCGATCACCCTAAAATTTATATTCTTATCCTACCA
q 60/1                 FFFFDFFGFFGGGFGGGDGGGFGGGEEEEEEGGGGGGEGG

a score=224 EG2=2.5e-05 E=8.8e-18 mismap=3.69e-11
s chrM   6412 48 + 16571 ATATAAAACCCCCTGCCATAACCCAATACCAAACGCCCCTCTTCGTCT
s 60/2    0 48 +    50 ATATAAAACCCCCTACCGTAACCCAATACCAAACACCCCTCTTCGTCT
q 60/2                 B>B>=)=(@3:87/;7;(;:B8((:+@<)@7<B@:27(<ABAABAAB@

a score=198 EG2=0.0097 E=4.8e-15 mismap=4.95e-08
s chrM   11907 41 + 16571 TAACCACGTTCTCCTGATCAAATATCACTCTCCTACTTACA
s 63/1     1 41 -    50 TAACCACATTCTCCTAATCAAATATCACTCTCCTACTTACA
q 63/1                  C>?;@CAC:>@@>6?DDCDBD@>@CC@CAB1=;5DD:DDDD

a score=200 EG2=0.0061 E=3e-15 mismap=3.13e-08
s chrM   11765 50 + 16571 CGCACTCACAGTCGCATCATAATCCTCTCTCAAGGACTTCAAACTCTACT
s 63/2     0 50 +    50 CGCTCTCCCAATCGCATCATAATCCTCTCTCAAAAAGTTCAAACACTACT
q 63/2                  3)>(<@))@A*<*>(0(94<:0)<+BA:>:7>B8'?%66-?>B5).A6@0

a score=204 EG2=0.0025 E=1.1e-15 mismap=4.69e-09
s chrM   13054 42 + 16571 CCCACCCCAGTCTCAGCCCTACTCCACTCAAGCACTATAGTT
s 64/1     0 42 -    50 CCCACCCCAATCTCAGCCCTACTCCACTCAAACACTATAGTT
q 64/1                  GCDDGDF?GDDGDGG@FDGGAFFDGECCE@GDGDFEGDGBFG

a score=198 EG2=0.0097 E=4.8e-15 mismap=1.85e-08
s chrM   12879 44 + 16571 TTTCATCCTCGCCTTAGCATGATTTATCCTACACTCCAACTCAT
s 64/2     3 44 +    50 TTTCATCCTCACCTTAACATAATTTATCCTACACTCCAACTCAT
q 64/2                  CAC@CCACCCC=CCCC0CCCCCCACCCBCCCCCCC@CCBC>C@C

a score=167 EG2=12 E=7.8e-12 mismap=3.79e-05
s chrM   1767 46 + 16571 TGGCGCAATAGATATAGTACCGCAAGGGAAAGATGAAAAATTATAA
s 65/2    2 46 -    50 TGGCGTAACAGATATAGTATTGTAAGGGAAAGATGAAAAATTATAA
q 65/2                 ??@B@<2;+A:B=>BB/CB?BBBB?@A??>BBB<60:>;>B;B;AA

a score=228 EG2=1e-05 E=3.3e-18 mismap=2.23e-11
s chrM   16146 46 + 16571 ACCACCTGTAGTACATAAAAACCCAACCCACATCAaaccccccccc
s 67/1     4 46 -    50 ACCACCTATAATACATAAAAACCCAACCCACATCAAaccccccccC
q 67/1                  GGGGGFGGGEGGGGGBGGGFGGGGFGGGGGGGGGFGGGGGGFFFFE

a score=258 EG2=1.1e-08 E=1.6e-21 mismap=4.81e-10
s chrM   4154 47 + 16571 ACCAACTCATACACCTCCTATGAAAAAACTTCCTACCACTCACCCTA
s 73/1    2 47 -    50 ACCAACTCATACACCTCCTATAAAAAAACTTCCTACCACTCACCCTA
q 73/1                 BCDBBBBDBBDBCBCCCCDDBDBBCCCC0CCCAFFDGBDDDCDDDD?

a score=210 EG2=0.00062 E=2.7e-16 mismap=9.03e-10
s chrM   12100 47 + 16571 ATCCCTCAACCCCGACATCATTACCGGGTTTTCCTCTTGTAAATATA
s 74/1     0 47 -    50 ATCCCTCAACCCCGACATCATTACCGAATTTTCCTCTTATAAATATA
q 74/1                  FGGGGGGGFGGGGEGFGGGGGGF?GGEGFEGGGGGGGGGGGGGGGGE

a score=153 EG2=2.8e+02 E=2.2e-10 mismap=0.000421
s chrM   11892 48 + 16571 AACTCTCTGTGCTAGTAACCACGTTCTCCTGATCAAATATCACTCTCC
s 74/2     2 48 +    50 AACTCTCTATACTAATAACCACATTCTCCTAATCAAATATCACTGTCC
q 74/2                  C>*3B6)<CBA;=@3AC:B3BBB;A>)@A?CCBBCB:3>)*@B=+??C

a score=228 EG2=1e-05 E=3.3e-18 mismap=5.56e-11
s chrM   4854 46 + 16571 CTTCTTCTCACATGACAAAAACTAGCCCCCATCTCAATCATATACC
s 75/1    4 46 -    50 CTTCTTCTCACATAACAAAAACTAACCCCCATCTCAATCATATACC
q 75/1                 GEGGGGCCEBEEDEEAFGBFEFGGDGFGGGEBFGGFGGDDFFFEFF

a score=199 EG2=0.0077 E=3.8e-15 mismap=4.26e-08
s chrM   4585 44 + 16571 CTTTTATTCCAGTTCTAACCAAAAAAATAAACCCTCGTTCCACA
s 75/2    3 44 +    50 CTTTTATTCCAATTCTAACCAAAAAAATAAAACCTCATTCCACA
q 75/2                 ?CA==A=BCBCABAAA9BBBB;@B?B@B@7@.CA9AACBB==.B

a score=276 EG2=1.8e-10 E=1.2e-23 mismap=7.79e-12
s chrM   5469 50 + 16571 ACGCTACTCCTACCTATCTCCCCTTTTATACTAATAATCTTATAGAAATT
s 78/1    0 50 -    50 ACGCTACTCCTACCTATCTCCCCTTTTATACTAATAATCTTATAAAAATT
q 78/1                 EEEEEE?BAEEBEEEBE:BEEEDDEEE:DBEEEDEEEBEEBEEEEEBDEE

a score=193 EG2=0.03 E=1.6e-14 mismap=0.0014
s chrM   1833 46 + 16571 ACTAACCCCTATACCTTCTGCATAATGAATTAACTAGAAATAACTT
s 82/2    4 46 +    50 ACTAACCCCTATACCTTCTACAAAATAAATTAACTAAAAATAACTT
q 82/2                 B5BCBCBC?CBCBBAA<BAB@4'AC@ABCCBCCBACBA@B@BCB7<

a score=228 EG2=1e-05 E=3.3e-18 mismap=4.63e-07
s chrM   10231 50 + 16571 TAGCTATTACCTTCTTATTATTTGATCTAGAAATTGCCCTCCTTTTACCC
s 84/1     0 50 -    50 TAGCTATTACCTTCTTATTATTTAATCTAAAAATTACCCTCCTTTTACCC
q 84/1                  EEEEEEEEEEAEEEEEEEEEDBDB@EEEEEEEEEEEEEE=EEEBEEAE:D

a score=203 EG2=0.0031 E=1.4e-15 mismap=4.42e-09
s chrM   11684 50 + 16571 TTCACCGGCGCAGTCATTCTCATAATCGCCCACGGACTCACATCCTCATT
s 86/2     0 50 +    50 TTCACCGACACAATCATTCTCATAATCGCCCACGAACTCACATCCTCATT
q 86/2                  ,A;B<1@?BC>48C@B?C9BB.<?BC=:CB=BBC=BAB*):C>=='<CBB

a score=204 EG2=0.0025 E=1.1e-15 mismap=3.54e-09
s chrM   8530 42 + 16571 AACGAAAATCTGTTCGCTTCATTCATTGCCCCCACAATCCTA
s 87/1    0 42 -    50 AACGAAAATCTATTCGCTTCATTCATTACCCCCACAATCCTA
q 87/1                 EGDFEGFGGGGGDGGGFGGGFGGCGGGEGGEFBGGGGFGGGG

a score=175 EG2=1.9 E=1.2e-12 mismap=2.71e-06
s chrM   8329 49 + 16571 TTAAAGATTAAGAGAACCAACACCTCTTTACAGTGAAATGCCCCAACTA
s 87/2    1 49 +    50 TTAAAAATTAAAAGAACCAACACCTCTTTACAATAAAATACCCCAACTA
q 87/2                 ACC@=BBBC<C=C@BBA@C=AB@BB=?@BB?BBBBABACAA4A?BBB?A

a score=210 EG2=0.00062 E=2.7e-16 mismap=2.86e-05
s chrM   1543 47 + 16571 TTATATAGAGGAGACAAGTCGTAACATGGTAAGTGTACTGGAaagtg
s 89/1    0 47 +    50 TTATATAGAGGAGATAAGTCGTAATATGGTAAGTGTATTGGAAAGTG
q 89/1                 EEEBEEEEEAFFFFFEGGGGGGGAEGGADGFGFGAGGGGGGGGEGGA

a score=180 EG2=0.59 E=3.6e-13 mismap=1.1e-06
s chrM   4414 50 + 16571 TAAATAAGCTATCGGGCCCATACCCCGAAAATGTTGGTTATACCCTTCCC
s 91/1    0 50 -    50 TAAATAAACTATCGAACCCATACCCCGAAAATATTGATTATACCCTTCCC
q 91/1                 BEEEEEEEBEEBECACGGFGFAAFFGDDGEGGDGGEGEGGFGGGFGGGGG

a score=183 EG2=0.3 E=1.7e-13 mismap=5.55e-07
s chrM   4245 46 + 16571 CATTCCCCCTCAAACCTAAGAAATATGTCTGATAAAAGAGTTACTT
s 91/2    4 46 +    50 CATTCCCCCTCAAACCTAAGAAATATATCTAATAAAAAAATTACTT
q 91/2                 @6@>7@4:BABBB==B@BBB@(@59A@BAA@B9>B??32B<A<B=;

a score=157 EG2=1.1e+02 E=8.4e-11 mismap=0.000375
s chrM   7286 50 + 16571 TCTAACAGCAGTAATattaataattttcatgatTTGAGAAGCCTTCGCTT
s 95/1    0 50 -    50 TCTAACAACAAtaatattaataattttcataatttaAAAAACCTTCGCTT
q 95/1                 A06+0@@;>@C:@@@;FFFFGFGGFFFFFFFEFECGGGGGFGFGGGEGGD

a score=251 EG2=5.3e-08 E=1e-20 mismap=1.71e-13
s chrM   7140 50 + 16571 CATTTCACTATCATATTCATCGGCGTAAATCTAACTTTCTTCCCACAACA
s 95/2    0 50 +    50 CATTTCACTATCATATTCATCGACATAAATCTAACTTTCTTCCCACAACA
q 95/2                 B@BABA>BB=.=?*@@?2<A?6BA=*?@>>8:?0;=6>@A8?76@7<?@?

a score=180 EG2=0.59 E=3.6e-13 mismap=8.7e-07
s chrM   9555 50 + 16571 CCCCCAACAGGCATCACCCCGCTAAATCCCCTAGAAGTCCCACTCCTAAA
s 98/1    0 50 -    50 CCCCCAACAAACATCACCCCACTAAATCCCCTAAAAATCCCACTCCTAAA
q 98/1                 EDEEBEBCEBFDFEFDGGGGEGGGGFGEGGGBGGGGGGFGGGGGG?EEEE

a score=154 EG2=2.3e+02 E=1.7e-10 mismap=0.000335
s chrM   9350 26 + 16571 AACCAACACACTAACCATATACCAAT
s 98/2    0 26 +    50 AACCAACACACTAACCATATACCAAT
q 98/2                 9BBBA6=?8*<(;BAB@7B+AA<<A<

# distance sample size: 35
# distance quartiles: 222 248 277
# estimated mean distance: 248
# estimated standard deviation of distance: 40.7715
# fraglen=248 sdev=40.7715 disjoint=0.01 genome=16571
156	chrM	11	50	+	16571	4/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000222
180	chrM	9112	50	+	16571	6/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.89e-07
220	chrM	8903	47	+	16571	6/2	3	47	+	50	47	EG2=6.3e-05	E=2.3e-17	mismap=9.34e-11
166	chrM	12726	38	+	16571	7/2	0	38	+	50	38	EG2=15	E=9.9e-12	mismap=2.31e-05
180	chrM	14235	50	+	16571	9/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=9.02e-07
253	chrM	14024	50	+	16571	9/2	0	50	+	50	50	EG2=3.4e-08	E=6e-21	mismap=4.26e-14
252	chrM	14182	50	+	16571	13/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
202	chrM	13997	49	+	16571	13/2	1	49	+	50	49	EG2=0.0039	E=1.8e-15	mismap=5.85e-09
210	chrM	4992	39	+	16571	14/2	2	39	+	50	39	EG2=0.00062	E=2.7e-16	mismap=2.86e-05
194	chrM	2765	47	+	16571	15/2	0	47	+	50	47	EG2=0.024	E=1.3e-14	mismap=0.00111
228	chrM	349	50	+	16571	16/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=2.13e-11
204	chrM	5387	50	+	16571	18/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.24e-09
226	chrM	5164	50	+	16571	18/2	0	50	+	50	50	EG2=1.6e-05	E=5.4e-18	mismap=2.75e-11
174	chrM	8013	41	+	16571	19/1	9	41	-	50	41	EG2=2.3	E=1.5e-12	mismap=3.82e-06
154	chrM	7826	50	+	16571	19/2	0	50	+	50	50	EG2=2.3e+02	E=1.7e-10	mismap=0.000372
210	chrM	4644	47	+	16571	20/1	3	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.12e-10
207	chrM	4452	50	+	16571	20/2	0	50	+	50	50	EG2=0.0012	E=5.5e-16	mismap=1.81e-09
156	chrM	607	46	+	16571	24/2	0	46	+	50	46	EG2=1.4e+02	E=1.1e-10	mismap=0.00026
156	chrM	9781	50	+	16571	25/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000269
199	chrM	9554	49	+	16571	25/2	1	49	+	50	49	EG2=0.0077	E=3.8e-15	mismap=1.42e-08
240	chrM	3458	48	+	16571	26/2	2	48	+	50	48	EG2=6.6e-07	E=1.7e-19	mismap=2.97e-08
156	chrM	16102	50	+	16571	31/1	0	50	-	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000595
206	chrM	15963	50	+	16571	31/2	0	50	+	50	50	EG2=0.0016	E=7e-16	mismap=6.34e-09
300	chrM	459	50	+	16571	34/2	0	50	+	50	50	EG2=7.3e-13	E=9.6e-27	mismap=0
204	chrM	15490	50	+	16571	35/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.31e-09
150	chrM	15318	49	+	16571	35/2	1	49	+	50	49	EG2=5.6e+02	E=4.4e-10	mismap=0.00101
186	chrM	9210	43	+	16571	36/1	4	43	-	50	43	EG2=0.15	E=8.5e-14	mismap=2.62e-07
176	chrM	9030	50	+	16571	36/2	0	50	+	50	50	EG2=1.5	E=9.2e-13	mismap=2.59e-06
276	chrM	10929	50	+	16571	39/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=0
205	chrM	10700	50	+	16571	39/2	0	50	+	50	50	EG2=0.002	E=8.9e-16	mismap=3.73e-09
168	chrM	15071	48	+	16571	43/1	0	48	-	50	48	EG2=9.2	E=6.2e-12	mismap=0.000191
281	chrM	14777	47	+	16571	43/2	0	47	+	50	47	EG2=5.6e-11	E=2.9e-24	mismap=0
264	chrM	3543	48	+	16571	44/1	0	48	-	50	48	EG2=2.7e-09	E=3.3e-22	mismap=1.42e-14
154	chrM	3288	27	+	16571	44/2	0	27	+	50	27	EG2=2.3e+02	E=1.7e-10	mismap=0.00082
204	chrM	358	50	+	16571	47/1	0	50	-	50	50	EG2=0.0025	E=1.1e-15	mismap=4.61e-09
215	chrM	130	50	+	16571	47/2	0	50	+	50	50	EG2=0.0002	E=7.9e-17	mismap=3.71e-10
222	chrM	419	41	+	16571	48/1	9	41	-	50	41	EG2=4e-05	E=1.4e-17	mismap=5.74e-11
156	chrM	207	50	+	16571	48/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.00021
213	chrM	13205	46	+	16571	55/1	4	46	-	50	46	EG2=0.00031	E=1.3e-16	mismap=5.32e-10
156	chrM	13027	50	+	16571	55/2	0	50	+	50	50	EG2=1.4e+02	E=1.1e-10	mismap=0.000248
252	chrM	10892	50	+	16571	59/1	0	50	-	50	50	EG2=4.2e-08	E=7.7e-21	mismap=5.68e-14
157	chrM	10691	34	+	16571	59/2	16	34	+	50	34	EG2=1.1e+02	E=8.4e-11	mismap=0.000167
168	chrM	6611	40	+	16571	60/1	0	40	-	50	40	EG2=9.2	E=6.2e-12	mismap=1.37e-05
224	chrM	6412	48	+	16571	60/2	0	48	+	50	48	EG2=2.5e-05	E=8.8e-18	mismap=3.69e-11
198	chrM	11907	41	+	16571	63/1	1	41	-	50	41	EG2=0.0097	E=4.8e-15	mismap=4.95e-08
200	chrM	11765	50	+	16571	63/2	0	50	+	50	50	EG2=0.0061	E=3e-15	mismap=3.13e-08
204	chrM	13054	42	+	16571	64/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=4.69e-09
198	chrM	12879	44	+	16571	64/2	3	44	+	50	44	EG2=0.0097	E=4.8e-15	mismap=1.85e-08
167	chrM	1767	46	+	16571	65/2	2	46	-	50	46	EG2=12	E=7.8e-12	mismap=3.79e-05
228	chrM	16146	46	+	16571	67/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=2.23e-11
258	chrM	4154	47	+	16571	73/1	2	47	-	50	47	EG2=1.1e-08	E=1.6e-21	mismap=4.81e-10
210	chrM	12100	47	+	16571	74/1	0	47	-	50	47	EG2=0.00062	E=2.7e-16	mismap=9.03e-10
153	chrM	11892	48	+	16571	74/2	2	48	+	50	48	EG2=2.8e+02	E=2.2e-10	mismap=0.000421
228	chrM	4854	46	+	16571	75/1	4	46	-	50	46	EG2=1e-05	E=3.3e-18	mismap=5.56e-11
199	chrM	4585	44	+	16571	75/2	3	44	+	50	44	EG2=0.0077	E=3.8e-15	mismap=4.26e-08
276	chrM	5469	50	+	16571	78/1	0	50	-	50	50	EG2=1.8e-10	E=1.2e-23	mismap=7.79e-12
193	chrM	1833	46	+	16571	82/2	4	46	+	50	46	EG2=0.03	E=1.6e-14	mismap=0.0014
228	chrM	10231	50	+	16571	84/1	0	50	-	50	50	EG2=1e-05	E=3.3e-18	mismap=4.63e-07
203	chrM	11684	50	+	16571	86/2	0	50	+	50	50	EG2=0.0031	E=1.4e-15	mismap=4.42e-09
204	chrM	8530	42	+	16571	87/1	0	42	-	50	42	EG2=0.0025	E=1.1e-15	mismap=3.54e-09
175	chrM	8329	49	+	16571	87/2	1	49	+	50	49	EG2=1.9	E=1.2e-12	mismap=2.71e-06
210	chrM	1543	47	+	16571	89/1	0	47	+	50	47	EG2=0.00062	E=2.7e-16	mismap=2.86e-05
180	chrM	4414	50	+	16571	91/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=1.1e-06
183	chrM	4245	46	+	16571	91/2	4	46	+	50	46	EG2=0.3	E=1.7e-13	mismap=5.55e-07
157	chrM	7286	50	+	16571	95/1	0	50	-	50	50	EG2=1.1e+02	E=8.4e-11	mismap=0.000375
251	chrM	7140	50	+	16571	95/2	0	50	+	50	50	EG2=5.3e-08	E=1e-20	mismap=1.71e-13
180	chrM	9555	50	+	16571	98/1	0	50	-	50	50	EG2=0.59	E=3.6e-13	mismap=8.7e-07
154	chrM	9350	26	+	16571	98/2	0	26	+	50	26	EG2=2.3e+02	E=1.7e-10	mismap=0.000335