[med-svn] [Git][med-team/last-align][upstream] New upstream version 1651
Charles Plessy (@plessy)
gitlab at salsa.debian.org
Fri Feb 6 23:08:18 GMT 2026
Charles Plessy pushed to branch upstream at Debian Med / last-align
Commits:
f7f81f83 by Charles Plessy at 2026-02-07T07:50:44+09:00
New upstream version 1651
- - - - -
12 changed files:
- bin/last-dotplot
- doc/last-cookbook.rst
- doc/last-evalues.rst
- doc/last-train.rst
- doc/lastal.rst
- src/AlignmentWrite.cc
- src/LastEvaluer.hh
- src/LastalArguments.cc
- src/lastal.cc
- src/makefile
- test/last-pair-test.out
- test/last-test.out
Changes:
=====================================
bin/last-dotplot
=====================================
@@ -977,10 +977,10 @@ def lastDotplot(opts, args):
logging.info("cutting...")
coverDictB1 = dict(mergedRangesPerSeq(coverDictB1))
coverDictB2 = dict(mergedRangesPerSeq(coverDictB2))
- cutRangesB1 = trimmed(seqRangesB1, coverDictB1, minAlignedBases,
- maxGapB1, 0, 0)
- cutRangesB2 = trimmed(seqRangesB2, coverDictB2, minAlignedBases,
- maxGapB2, 0, 0)
+ cutRangesB1 = list(trimmed(seqRangesB1, coverDictB1, minAlignedBases,
+ maxGapB1, 0, 0))
+ cutRangesB2 = list(trimmed(seqRangesB2, coverDictB2, minAlignedBases,
+ maxGapB2, 0, 0))
logging.info("sorting...")
sortOut = allSortedRanges(opts, alignments, alignmentsB,
@@ -1096,6 +1096,12 @@ def lastDotplot(opts, args):
box = lMargin, i - opts.border_pixels, rMarginBeg, i
im.paste(opts.border_color, box)
+ if cutRangesB1 and cutRangesB2: # no 2ndary-versus-2ndary alignments
+ x = rangePixBegs1[len(cutRanges1)] - opts.border_pixels
+ y = rangePixBegs2[len(cutRanges2)] - opts.border_pixels
+ box = x, y, rMarginBeg, bMarginBeg
+ im.paste(opts.margin_color, box)
+
base = args[0]
if base.endswith((".maf", ".psl", ".tab", ".tsv", ".txt")):
base = base[:-4]
=====================================
doc/last-cookbook.rst
=====================================
@@ -32,12 +32,12 @@ line-wrapping. For example, you can use::
Each alignment looks like this (MAF_ format)::
- a score=39 EG2=0.092 E=4.9e-11
+ a score=39 EG=0.092 E=4.9e-11
s humanMito 12264 71 + 16571 CAACTTTTAAAGGATAACAGCTA-TCCATTGGTCTTAGGCCccaa...
s fuguMito 11840 73 + 16447 CAGCTTTTGAAGGATAATAGCTAATCCGTTGGTCTTAGGAACCAA...
-Score, EG2, and E are explained at last-evalues_. Lines starting with
+Score, EG, and E are explained at last-evalues_. Lines starting with
"s" contain: the sequence name, the start coordinate of the alignment,
the number of bases spanned by the alignment, the strand, the sequence
length, and the aligned bases.
@@ -332,6 +332,16 @@ For whole genome dotplots, high sensitivity is not necessary.
``RY128`` is fine for e.g. human versus crocodile, though not for
human versus fish.
+Aligning closely-related non-ape genomes
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+You can use the human-chimp recipe above, except unfortunately the
+best ``-U`` value is unclear. If the genomes have similar huge tandem
+repeats, ``-U`` must be high enough to mask them. For some deer
+genomes, ``-U800`` was necessary. Tentative suggestion: try
+``-U1000`` if the genomes are > ~90% identical, and longer than than a
+few hundred million base-pairs.
+
Moar faster
-----------
@@ -355,7 +365,7 @@ alignment column with lastal_ option ``-j4``::
The output looks like this::
- a score=17 EG2=9.3e+09 E=5e-06
+ a score=17 EG=3.1e+05 E=5e-06
s seqX 0 57 + 57 TGAAGTTAAAGGTATATGAATTCCAATTCTTAACCCCCCTATTAAACGAATATCTTG
s seqY 0 51 + 51 TGAAGTTAGAGGTAT--GGTTTTGAGTAGT----CCTCCTATTTTTCGAATATCTTG
p %*.14442011.(%##"%$$$$###""!!!""""&'(*,340.,,.~~~~~~~~~~~
=====================================
doc/last-evalues.rst
=====================================
@@ -4,7 +4,7 @@ LAST E-values
It's useful to know whether a similarity is likely to occur by chance
between random sequences (shuffled letters). Let's see an example::
- a score=37 EG2=72 E=0.00001
+ a score=37 EG=12 E=0.00001
s chr3 9 23 + 939557 TCTGTGAGTTGAAGTTTTCGCCCTAG
s seqA 2 21 + 2500 TCTGGGAGTTGAAGGTT--GCCCCAG
@@ -13,27 +13,30 @@ related rather than random. (It doesn't reflect percent identity: a
long, weak similarity can have the same score as a short, strong
similarity.)
-The "E-value" is the expected number of similarities with score ≥ 37
-between random sequences. This depends on how long the random
-sequences are.
+``EG`` is "expected hits per giga". It's the expected number of
+similarities with score ≥ 37 per billion letters of query sequence.
+(Query sequences are those given to ``lastal``, reference/database
+sequences to ``lastdb``.) In other words, ``EG`` is the expected
+number of hits between:
-``EG2=`` is the E-value for two random sequences of length 1 billion
-each. (Expected alignments per square giga.)
+* a random sequence of length 1 billion, and
+* random sequences with the same lengths as all the reference sequences.
-The meaning of ``E=`` depends on whether you use lastal_ option ``-H``.
-``E=`` is the E-value for random sequences with the same lengths as
+The meaning of ``E`` depends on whether you use lastal_ option ``-H``.
+``E`` is the expected number of hits between random sequences with the
+same lengths as
* all the reference sequences, and all the query sequences (with ``-H``),
* all the reference sequences, and that one query sequence (without ``-H``).
-Experience suggests that ``EG2=`` isn't useful, ``E=`` with ``-H`` is
-useful, and ``E=`` without ``-H`` is misleading.
+If there is just one query, these are the same. If there's more than
+one query, ``E`` without ``-H`` is confusing, and best ignored!
Default threshold
-----------------
By default, LAST reports similarities that are expected by chance at
-most once per million query letters (for a given database).
+most once per million query letters (EG ≤ 1000).
Setting a threshold
-------------------
@@ -47,7 +50,7 @@ Alternatively, you can use ``-H``::
lastal -H100 humdb fuguMito.fa > myalns.maf
-This will get similarities with ``E=`` ≤ 100.
+This will get similarities with ``E`` ≤ 100.
Advanced issues
---------------
@@ -64,20 +67,14 @@ If your sequences have very different letter frequencies (e.g. very
AT-rich DNA), the E-values will be misleading. The best solution is
to use a suitable score matrix, such as AT77 or ATMAP.
-Calculation of EG2
-~~~~~~~~~~~~~~~~~~
-
-EG2 is calculated from the score like this::
+Calculation of EG and E
+~~~~~~~~~~~~~~~~~~~~~~~
- EG2 = K * exp(-lambda * score) * (1 billion) * (1 billion)
+``EG`` is calculated from the score like this::
-The parameters lambda and K are printed in the header of lastal's
-output.
+ EG = K * exp(-lambda * score) * strands * (1 billion) * length(n, score) * N / n
-Calculation of E
-~~~~~~~~~~~~~~~~
-
-``E=`` is calculated like this. With ``-H``::
+``E`` is calculated like this. With ``-H``::
E = K * exp(-lambda * score) * strands * area(m, n, score) * M/m * N/n
@@ -85,7 +82,8 @@ Without ``-H``::
E = K * exp(-lambda * score) * strands * area(q, n, score) * N/n
-* ``strands`` is either 1 or 2 (if both DNA strands are searched).
+* ``lambda`` and ``K`` are printed in the header of lastal's output.
+* ``strands`` is 2 if two DNA strands are searched, else it is 1.
* ``m`` is the length of the longest query sequence.
* ``n`` is the length of the longest sequence in the database.
* ``M`` is the total length of all query sequences.
@@ -93,15 +91,7 @@ Without ``-H``::
* ``q`` is the length of that one query sequence.
* ``area(x, y, score)`` is `slightly less than
<https://doi.org/10.1186/1756-0500-5-286>`_ ``x * y``.
-
-Effect of option -D
-~~~~~~~~~~~~~~~~~~~
-
-Option -D sets the minimum alignment score, by this formula::
-
- K * exp(-lambda * score) * strands * D * length(n, score) * N / n ≤ 1,
-
-where ``length(n, score)`` is slightly less than ``n``.
+* ``length(n, score)`` is slightly less than ``n``.
Bit score
~~~~~~~~~
@@ -109,7 +99,6 @@ Bit score
Some people like to use "bit scores"::
Bit score = (lambda * score - ln[K]) / ln[2]
- = log2( 1e18 / EG2 )
Limitations
~~~~~~~~~~~
=====================================
doc/last-train.rst
=====================================
@@ -128,8 +128,7 @@ Alignment options
-D LENGTH Query letters per random alignment. (See `here
<doc/last-evalues.rst>`_.)
--E EG2 Maximum expected alignments per square giga. (See `here
- <doc/last-evalues.rst>`_.)
+-E EG2 Maximum expected alignments per square giga.
-s NUMBER Which query strand to use: 0=reverse, 1=forward, 2=both.
If specified, this parameter is written in last-train's
output, so it will override lastal's default.
=====================================
doc/lastal.rst
=====================================
@@ -71,7 +71,7 @@ Cosmetic options
**MAF** format looks like this::
- a score=15 EG2=4.7e+04 E=2.6e-05
+ a score=15 EG=88 E=2.6e-05
s chr3 9 23 + 939557 TTTGGGAGTTGAAGTTTTCGCCC
s seqA 2 21 + 25 TTTGGGAGTTGAAGGTT--GCCC
@@ -84,7 +84,7 @@ Cosmetic options
The same alignment in **TAB** format looks like this::
- 15 chr3 9 23 + 939557 seqA 2 21 + 25 17,2:0,4 EG2=4.7e+04 E=2.6e-05
+ 15 chr3 9 23 + 939557 seqA 2 21 + 25 17,2:0,4 EG=88 E=2.6e-05
The ``17,2:0,4`` shows the sizes and offsets of gapless blocks in
the alignment. In this case, we have a block of size 17, then an
@@ -105,8 +105,10 @@ Cosmetic options
**BlastTab+** format is the same as BlastTab, with 3 extra
columns at the end: length of query sequence, length of
- reference sequence, and (raw) score. More columns might be
- added in future.
+ reference sequence, and (raw) score. For DNA-versus-protein
+ searches, there is a 4th extra column showing which frame in the
+ DNA the alignment starts in. More columns might be added in
+ future.
For backwards compatibility, a NAME of 0 means TAB and 1 means
MAF.
@@ -135,8 +137,13 @@ E-value options
* All the reference sequences, and all the query sequences (with ``-H``)
* All the reference sequences, and that one query sequence (without ``-H``)
+ If you use both ``-D`` and ``-H``, then alignments will be
+ reported according to ``-D``. The meaning of ``E=`` will be
+ according to ``-H``.
+
-E EXPECT
- Maximum EG2_ (expected alignments per square giga). Experience
+ Maximum expected alignments per square giga (between 2 random
+ single-strand sequences of length 1 billion each). Experience
suggests that ``-D`` and ``-H`` are more convenient than ``-E``.
Score options
@@ -687,7 +694,6 @@ computer for each volume.
.. _lastdb: doc/lastdb.rst
.. _last-train: doc/last-train.rst
.. _last-split: doc/last-split.rst
-.. _EG2:
.. _E-values: doc/last-evalues.rst
.. _tantan: https://gitlab.com/mcfrith/tantan
.. _gentle masking: https://doi.org/10.1371/journal.pone.0028819
=====================================
src/AlignmentWrite.cc
=====================================
@@ -199,9 +199,9 @@ static char* writeTags( const LastEvaluer& evaluer, double queryLength,
double epa = evaluer.evaluePerArea( score );
double area = evaluer.area( score, queryLength );
*out++ = separator;
- static const char aString[] = "EG2=";
+ static const char aString[] = "EG=";
out = copyString(out, aString, sizeof aString);
- out += std::sprintf(out, "%.2g", 1e18 * epa);
+ out += std::sprintf(out, "%.2g", evaluer.areaPerGiga(score) * epa);
*out++ = separator;
static const char eString[] = "E=";
out = copyString(out, eString, sizeof eString);
@@ -545,6 +545,7 @@ AlignmentText Alignment::writeBlastTab(const MultiSequence& seq1,
if (isExtraColumns) {
sc.set(scoreFormat(extras.fullScore >= 0), score);
s += s1.size() + s2.size() + sc.size() + 3;
+ if (translationType) s += 2 + (strand2 == '-');
}
char *text = new char[s + 1];
@@ -553,7 +554,15 @@ AlignmentText Alignment::writeBlastTab(const MultiSequence& seq1,
w << n2 << t << n1 << t << mp << t << as << t << mm << t << go << t
<< b2 << t << e2 << t << b1 << t << e1;
if (evaluer.isGood()) w << t << ev << t << bs;
- if (isExtraColumns) w << t << s2 << t << s1 << t << sc;
+ if (isExtraColumns) {
+ w << t << s2 << t << s1 << t << sc;
+ if (translationType) {
+ char frame = '1' + (alnBeg2 - seqStart2) % 3;
+ w << t;
+ if (strand2 == '-') w << '-';
+ w << frame;
+ }
+ }
w << '\n' << '\0';
return AlignmentText(seqNum2, alnBeg2, alnEnd2, size2, strand2, score,
=====================================
src/LastEvaluer.hh
=====================================
@@ -106,6 +106,11 @@ public:
return areaMultiplier * evaluer.area(score, q, databaseMaxSeqLen);
}
+ // Don't call this in the "bad" state or before calling setSearchSpace:
+ double areaPerGiga(double score) const {
+ return databaseLenMultiplier * evaluer.area(score, 1e9, databaseMaxSeqLen);
+ }
+
// Don't call this in the "bad" state:
double bitScore(double score) const { return evaluer.bitScore(score); }
=====================================
src/LastalArguments.cc
=====================================
@@ -698,6 +698,8 @@ To proceed without E-values, set a score threshold with option -e.");
ERR("can't do split DNA-protein alignment");
if (gapPairCost > 0)
ERR("can't do split alignment with option -c");
+ if (scoreType == 1) // the split score threshold would be too high
+ ERR("can't do split alignment with option -J 1");
if (verbosity > 1) splitOpts.verbose = true;
splitOpts.setUnspecifiedValues(minScoreGapped, temperature);
}
=====================================
src/lastal.cc
=====================================
@@ -1362,7 +1362,8 @@ static void scanOneVolume(unsigned volume, unsigned numOfThreadsLeft) {
static void readAllSequenceLengths(char **fileNames,
countT &totLen, countT &maxLen) {
for (char **i = fileNames; *i; ++i) {
- if (isSingleDash(*i)) err("need real files to read query sequences twice");
+ if (isSingleDash(*i))
+ err("option -H requires real files, to read the query sequences twice");
mcf::izstream f;
openOrThrow(f, *i);
readSequenceLengths(f, totLen, maxLen);
@@ -1618,6 +1619,8 @@ void writeHeader(countT numOfRefSeqs, countT refLetters, std::ostream &out) {
std::cout << "#\n";
}
}
+
+ out.flush();
}
void lastal(int argc, char **argv) {
=====================================
src/makefile
=====================================
@@ -95,7 +95,7 @@ ScoreMatrixData.hh: ../data/*.mat
../build/mat-inc.sh ../data/*.mat > $@
VERSION1 = git describe --dirty
-VERSION2 = echo ' (HEAD -> main, tag: 1645) ' | sed -e 's/.*tag: *//' -e 's/[,) ].*//'
+VERSION2 = echo ' (HEAD -> main, tag: 1651) ' | sed -e 's/.*tag: *//' -e 's/[,) ].*//'
VERSION = \"`test -e ../.git && $(VERSION1) || $(VERSION2)`\"
=====================================
test/last-pair-test.out
=====================================
@@ -38,342 +38,342 @@ Options:
# 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
+a score=156 EG=0.0047 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 at CBCACCCCCCCC@ABCCAABBCCBC at A=B at CBCACCC@
-a score=180 EG2=0.59 E=3.6e-13 mismap=8.89e-07
+a score=180 EG=2e-05 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
+a score=220 EG=2.1e-09 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 at B?B<AB>A<AB:@AB9;=@AA43A>@;,5=@@A
-a score=166 EG2=15 E=9.9e-12 mismap=2.31e-05
+a score=166 EG=0.00048 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
+a score=180 EG=2e-05 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
+a score=253 EG=1.1e-12 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
+a score=252 EG=1.4e-12 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
+a score=202 EG=1.3e-07 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 at 9?@AA9=@BB;@>@A8><2?99A@
-a score=210 EG2=0.00062 E=2.7e-16 mismap=2.86e-05
+a score=210 EG=2.1e-08 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 at C+BCBBBCC<@B9BB
-a score=194 EG2=0.024 E=1.3e-14 mismap=0.00111
+a score=194 EG=8e-07 E=1.3e-14 mismap=0.00111
s chrM 2765 47 + 16571 AACAAACCCACAGGTCCTAAACTACCAAACCTGCATTAAAAATTTCG
s 15/2 0 47 + 50 AACAGACCCACAAATCCTAAACTACCAAACCTACATTAAAAATTTCG
q 15/2 2 at 1>(?@)@C>(>CCC>)@1:15=>9@=ACCC@>A?*=B?8,>4 at BC
-a score=228 EG2=1e-05 E=3.3e-18 mismap=2.13e-11
+a score=228 EG=3.4e-10 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
+a score=204 EG=8.1e-08 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
+a score=226 EG=5.3e-10 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 at C=9CBC,CA=-C at CB=CA>9?/
-a score=174 EG2=2.3 E=1.5e-12 mismap=3.82e-06
+a score=174 EG=7.7e-05 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
+a score=154 EG=0.0075 E=1.7e-10 mismap=0.000372
s chrM 7826 50 + 16571 CTACGCATCCTTTACATAACAGACGAGGTCAACGATCCCTCCCTTACCAT
s 19/2 0 50 + 50 CTACACATACTTTACATAACAAACGAAATCAAGAATCCCTCCCTTACCAT
q 19/2 =@9 at B>C<';>)*8?'>67 at 1;(?<;A<@:@@(9A@@)(38>9+@<8?6@
-a score=210 EG2=0.00062 E=2.7e-16 mismap=9.12e-10
+a score=210 EG=2.1e-08 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
+a score=207 EG=4.1e-08 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 at CBC:1?>C?-2AA?59'<)?B7A46)@>:*?8,5)<24?)5>?;
-a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.00026
+a score=156 EG=0.0047 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
+a score=156 EG=0.0047 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
+a score=199 EG=2.5e-07 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 at 25)=CB9B7;BC7>?A4A;5 at BA+;6)()1BCA748@)=1?
-a score=240 EG2=6.6e-07 E=1.7e-19 mismap=2.97e-08
+a score=240 EG=2.2e-11 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 at .9,BA<:=?>-B6.)*=>0?59;'3
-a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.000595
+a score=156 EG=0.0047 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
+a score=206 EG=5.1e-08 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 at A@?A at A<@@BAA>?.9,9A?A<@?A
-a score=300 EG2=7.3e-13 E=9.6e-27 mismap=0
+a score=300 EG=2.4e-17 E=9.6e-27 mismap=0
s chrM 459 50 + 16571 CCTCCCACTCCCATACTACTAATCTCATCAATACAACCCCCGCCCATCCT
s 34/2 0 50 + 50 CCTCCCACTCCCATACTACTAATCTCATCAATACAACCCCCGCCCATCCT
q 34/2 BCABBBBCBBBCAACB;BCAB at BBBBBCABBA?AB>B at 5@B at BA5BCB7?
-a score=204 EG2=0.0025 E=1.1e-15 mismap=4.31e-09
+a score=204 EG=8.1e-08 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
+a score=150 EG=0.019 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 at B=CCCCBBB at B@>@@BBACB?@B at CB@CBACC>A4BBACAA
-a score=186 EG2=0.15 E=8.5e-14 mismap=2.62e-07
+a score=186 EG=5e-06 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 at DEEFBD7EEE?FFFEEE:EFFFEE??D?B9EB?D5CEEE?C
-a score=176 EG2=1.5 E=9.2e-13 mismap=2.59e-06
+a score=176 EG=4.9e-05 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
+a score=276 EG=5.8e-15 E=1.2e-23 mismap=0
s chrM 10929 50 + 16571 CCTCCGACCCCCTAACAACCCCCCTCCTAATACTAACTACCTGACTCCTA
s 39/1 0 50 - 50 CCTCCGACCCCCTAACAACCCCCCTCCTAATACTAACTACCTAACTCCTA
q 39/1 ?CBDD3;@4 at B>>,>C7@>>.?=BB>7?4 at C8>C3EFFFFFFFFFFDFFF
-a score=205 EG2=0.002 E=8.9e-16 mismap=3.73e-09
+a score=205 EG=6.5e-08 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 at BABBB=@BABBAAA at BBAA>;0===8>=;7.5:<;@@@A
-a score=168 EG2=9.2 E=6.2e-12 mismap=0.000191
+a score=168 EG=0.0003 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
+a score=281 EG=1.9e-15 E=2.9e-24 mismap=0
s chrM 14777 47 + 16571 ATAAAATTAATTAACCACTCATTCATCGACCTCCCCACCCCATCCAA
s 43/2 0 47 + 50 ATAAAATTAATTAACCACTCATTCATCGACCTCCCCACCCCATCCAA
q 43/2 BC at BCCCC=B=BBB at B@A at B?B at B9AAB>>BBBA):B<>??05:?7@
-a score=264 EG2=2.7e-09 E=3.3e-22 mismap=1.42e-14
+a score=264 EG=9e-14 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
+a score=154 EG=0.0075 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 at -7'':0%'@-B
-a score=204 EG2=0.0025 E=1.1e-15 mismap=4.61e-09
+a score=204 EG=8.1e-08 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 at CB;;@
-a score=215 EG2=0.0002 E=7.9e-17 mismap=3.71e-10
+a score=215 EG=6.6e-09 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 at CB6BCCB9?BAB=CCB@<B;C<CA at C@,9CCBAB>?=BBA
-a score=222 EG2=4e-05 E=1.4e-17 mismap=5.74e-11
+a score=222 EG=1.3e-09 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
+a score=156 EG=0.0047 E=1.1e-10 mismap=0.00021
s chrM 207 50 + 16571 TTAATtaattaaTGCTTGTAGGACATAATAATAACAATTGAATGTCTGCA
s 48/2 0 50 + 50 TTAATTAATTAATACTTATAAAACATAAtaataaCAATTAAATATCTGCA
q 48/2 BCBBB at 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
+a score=213 EG=1e-08 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
+a score=156 EG=0.0047 E=1.1e-10 mismap=0.000248
s chrM 13027 50 + 16571 CCCTGACTCCCCTCAGCCATAGAAGGCCCCACCCCAGTCTCAGCCCTACT
s 55/2 0 50 + 50 CCCTAACTCCCCTCAACCATAAAAAACCCCACCCCAATCTCAGCCCTACT
q 55/2 BBCBBBBC at C>B at B?@BBBBABAAB?A at B@BB>BBB?AACA@@BABBBBA
-a score=252 EG2=4.2e-08 E=7.7e-21 mismap=5.68e-14
+a score=252 EG=1.4e-12 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
+a score=157 EG=0.0038 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 at 5A99>024)?B>A
-a score=168 EG2=9.2 E=6.2e-12 mismap=1.37e-05
+a score=168 EG=0.0003 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
+a score=224 EG=8.4e-10 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
+a score=198 EG=3.2e-07 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 at CAB1=;5DD:DDDD
-a score=200 EG2=0.0061 E=3e-15 mismap=3.13e-08
+a score=200 EG=2e-07 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 at 0
-a score=204 EG2=0.0025 E=1.1e-15 mismap=4.69e-09
+a score=204 EG=8.1e-08 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 at FDGGAFFDGECCE@GDGDFEGDGBFG
-a score=198 EG2=0.0097 E=4.8e-15 mismap=1.85e-08
+a score=198 EG=3.2e-07 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 at CCACCCC=CCCC0CCCCCCACCCBCCCCCCC at CCBC>C at C
-a score=167 EG2=12 E=7.8e-12 mismap=3.79e-05
+a score=167 EG=0.00038 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
+a score=228 EG=3.4e-10 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
+a score=258 EG=3.5e-13 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
+a score=210 EG=2.1e-08 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
+a score=153 EG=0.0094 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
+a score=228 EG=3.4e-10 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
+a score=199 EG=2.5e-07 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 at B@7 at .CA9AACBB==.B
-a score=276 EG2=1.8e-10 E=1.2e-23 mismap=7.79e-12
+a score=276 EG=5.8e-15 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
+a score=193 EG=1e-06 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 at 4'AC at ABCCBCCBACBA@B at BCB7<
-a score=228 EG2=1e-05 E=3.3e-18 mismap=4.63e-07
+a score=228 EG=3.4e-10 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 at EEEEEEEEEEEEEE=EEEBEEAE:D
-a score=203 EG2=0.0031 E=1.4e-15 mismap=4.42e-09
+a score=203 EG=1e-07 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 at B?C9BB.<?BC=:CB=BBC=BAB*):C>=='<CBB
-a score=204 EG2=0.0025 E=1.1e-15 mismap=3.54e-09
+a score=204 EG=8.1e-08 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
+a score=175 EG=6.1e-05 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 at BBA@C=AB at BB=?@BB?BBBBABACAA4A?BBB?A
-a score=210 EG2=0.00062 E=2.7e-16 mismap=2.86e-05
+a score=210 EG=2.1e-08 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
+a score=180 EG=2e-05 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
+a score=183 EG=9.9e-06 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 at 4:BABBB==B at BBB@(@59A at BAA@B9>B??32B<A<B=;
-a score=157 EG2=1.1e+02 E=8.4e-11 mismap=0.000375
+a score=157 EG=0.0038 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
+a score=251 EG=1.8e-12 E=1e-20 mismap=1.71e-13
s chrM 7140 50 + 16571 CATTTCACTATCATATTCATCGGCGTAAATCTAACTTTCTTCCCACAACA
s 95/2 0 50 + 50 CATTTCACTATCATATTCATCGACATAAATCTAACTTTCTTCCCACAACA
q 95/2 B at BABA>BB=.=?*@@?2<A?6BA=*?@>>8:?0;=6>@A8?76 at 7<?@?
-a score=180 EG2=0.59 E=3.6e-13 mismap=8.7e-07
+a score=180 EG=2e-05 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
+a score=154 EG=0.0075 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 at 7B+AA<<A<
@@ -383,771 +383,771 @@ q 98/2 9BBBA6=?8*<(;BAB at 7B+AA<<A<
# 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
+156 chrM 11 50 + 16571 4/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000222
+180 chrM 9112 50 + 16571 6/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.89e-07
+220 chrM 8903 47 + 16571 6/2 3 47 + 50 47 EG=2.1e-09 E=2.3e-17 mismap=9.34e-11
+166 chrM 12726 38 + 16571 7/2 0 38 + 50 38 EG=0.00048 E=9.9e-12 mismap=2.31e-05
+180 chrM 14235 50 + 16571 9/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=9.02e-07
+253 chrM 14024 50 + 16571 9/2 0 50 + 50 50 EG=1.1e-12 E=6e-21 mismap=4.26e-14
+252 chrM 14182 50 + 16571 13/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+202 chrM 13997 49 + 16571 13/2 1 49 + 50 49 EG=1.3e-07 E=1.8e-15 mismap=5.85e-09
+210 chrM 4992 39 + 16571 14/2 2 39 + 50 39 EG=2.1e-08 E=2.7e-16 mismap=2.86e-05
+228 chrM 349 50 + 16571 16/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=2.13e-11
+204 chrM 5387 50 + 16571 18/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.24e-09
+226 chrM 5164 50 + 16571 18/2 0 50 + 50 50 EG=5.3e-10 E=5.4e-18 mismap=2.75e-11
+174 chrM 8013 41 + 16571 19/1 9 41 - 50 41 EG=7.7e-05 E=1.5e-12 mismap=3.82e-06
+154 chrM 7826 50 + 16571 19/2 0 50 + 50 50 EG=0.0075 E=1.7e-10 mismap=0.000372
+210 chrM 4644 47 + 16571 20/1 3 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.12e-10
+207 chrM 4452 50 + 16571 20/2 0 50 + 50 50 EG=4.1e-08 E=5.5e-16 mismap=1.81e-09
+156 chrM 607 46 + 16571 24/2 0 46 + 50 46 EG=0.0047 E=1.1e-10 mismap=0.00026
+156 chrM 9781 50 + 16571 25/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000269
+199 chrM 9554 49 + 16571 25/2 1 49 + 50 49 EG=2.5e-07 E=3.8e-15 mismap=1.42e-08
+240 chrM 3458 48 + 16571 26/2 2 48 + 50 48 EG=2.2e-11 E=1.7e-19 mismap=2.97e-08
+156 chrM 16102 50 + 16571 31/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000595
+206 chrM 15963 50 + 16571 31/2 0 50 + 50 50 EG=5.1e-08 E=7e-16 mismap=6.34e-09
+300 chrM 459 50 + 16571 34/2 0 50 + 50 50 EG=2.4e-17 E=9.6e-27 mismap=0
+204 chrM 15490 50 + 16571 35/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.31e-09
+186 chrM 9210 43 + 16571 36/1 4 43 - 50 43 EG=5e-06 E=8.5e-14 mismap=2.62e-07
+176 chrM 9030 50 + 16571 36/2 0 50 + 50 50 EG=4.9e-05 E=9.2e-13 mismap=2.59e-06
+276 chrM 10929 50 + 16571 39/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=0
+205 chrM 10700 50 + 16571 39/2 0 50 + 50 50 EG=6.5e-08 E=8.9e-16 mismap=3.73e-09
+168 chrM 15071 48 + 16571 43/1 0 48 - 50 48 EG=0.0003 E=6.2e-12 mismap=0.000191
+281 chrM 14777 47 + 16571 43/2 0 47 + 50 47 EG=1.9e-15 E=2.9e-24 mismap=0
+264 chrM 3543 48 + 16571 44/1 0 48 - 50 48 EG=9e-14 E=3.3e-22 mismap=1.42e-14
+154 chrM 3288 27 + 16571 44/2 0 27 + 50 27 EG=0.0075 E=1.7e-10 mismap=0.00082
+204 chrM 358 50 + 16571 47/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.61e-09
+215 chrM 130 50 + 16571 47/2 0 50 + 50 50 EG=6.6e-09 E=7.9e-17 mismap=3.71e-10
+222 chrM 419 41 + 16571 48/1 9 41 - 50 41 EG=1.3e-09 E=1.4e-17 mismap=5.74e-11
+156 chrM 207 50 + 16571 48/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.00021
+213 chrM 13205 46 + 16571 55/1 4 46 - 50 46 EG=1e-08 E=1.3e-16 mismap=5.32e-10
+156 chrM 13027 50 + 16571 55/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000248
+252 chrM 10892 50 + 16571 59/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+157 chrM 10691 34 + 16571 59/2 16 34 + 50 34 EG=0.0038 E=8.4e-11 mismap=0.000167
+168 chrM 6611 40 + 16571 60/1 0 40 - 50 40 EG=0.0003 E=6.2e-12 mismap=1.37e-05
+224 chrM 6412 48 + 16571 60/2 0 48 + 50 48 EG=8.4e-10 E=8.8e-18 mismap=3.69e-11
+198 chrM 11907 41 + 16571 63/1 1 41 - 50 41 EG=3.2e-07 E=4.8e-15 mismap=4.95e-08
+200 chrM 11765 50 + 16571 63/2 0 50 + 50 50 EG=2e-07 E=3e-15 mismap=3.13e-08
+204 chrM 13054 42 + 16571 64/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=4.69e-09
+198 chrM 12879 44 + 16571 64/2 3 44 + 50 44 EG=3.2e-07 E=4.8e-15 mismap=1.85e-08
+167 chrM 1767 46 + 16571 65/2 2 46 - 50 46 EG=0.00038 E=7.8e-12 mismap=3.79e-05
+228 chrM 16146 46 + 16571 67/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=2.23e-11
+258 chrM 4154 47 + 16571 73/1 2 47 - 50 47 EG=3.5e-13 E=1.6e-21 mismap=4.81e-10
+210 chrM 12100 47 + 16571 74/1 0 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.03e-10
+153 chrM 11892 48 + 16571 74/2 2 48 + 50 48 EG=0.0094 E=2.2e-10 mismap=0.000421
+228 chrM 4854 46 + 16571 75/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=5.56e-11
+199 chrM 4585 44 + 16571 75/2 3 44 + 50 44 EG=2.5e-07 E=3.8e-15 mismap=4.26e-08
+276 chrM 5469 50 + 16571 78/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=7.79e-12
+228 chrM 10231 50 + 16571 84/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=4.63e-07
+203 chrM 11684 50 + 16571 86/2 0 50 + 50 50 EG=1e-07 E=1.4e-15 mismap=4.42e-09
+204 chrM 8530 42 + 16571 87/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=3.54e-09
+175 chrM 8329 49 + 16571 87/2 1 49 + 50 49 EG=6.1e-05 E=1.2e-12 mismap=2.71e-06
+210 chrM 1543 47 + 16571 89/1 0 47 + 50 47 EG=2.1e-08 E=2.7e-16 mismap=2.86e-05
+180 chrM 4414 50 + 16571 91/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=1.1e-06
+183 chrM 4245 46 + 16571 91/2 4 46 + 50 46 EG=9.9e-06 E=1.7e-13 mismap=5.55e-07
+157 chrM 7286 50 + 16571 95/1 0 50 - 50 50 EG=0.0038 E=8.4e-11 mismap=0.000375
+251 chrM 7140 50 + 16571 95/2 0 50 + 50 50 EG=1.8e-12 E=1e-20 mismap=1.71e-13
+180 chrM 9555 50 + 16571 98/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.7e-07
+154 chrM 9350 26 + 16571 98/2 0 26 + 50 26 EG=0.0075 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
+156 chrM 11 50 + 16571 4/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000213
+180 chrM 9112 50 + 16571 6/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=9.4e-07
+220 chrM 8903 47 + 16571 6/2 3 47 + 50 47 EG=2.1e-09 E=2.3e-17 mismap=9.88e-11
+166 chrM 12726 38 + 16571 7/2 0 38 + 50 38 EG=0.00048 E=9.9e-12 mismap=2.49e-05
+180 chrM 14235 50 + 16571 9/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=9.6e-07
+253 chrM 14024 50 + 16571 9/2 0 50 + 50 50 EG=1.1e-12 E=6e-21 mismap=5.68e-14
+252 chrM 14182 50 + 16571 13/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+202 chrM 13997 49 + 16571 13/2 1 49 + 50 49 EG=1.3e-07 E=1.8e-15 mismap=5.63e-09
+210 chrM 4992 39 + 16571 14/2 2 39 + 50 39 EG=2.1e-08 E=2.7e-16 mismap=1.44e-05
+194 chrM 2765 47 + 16571 15/2 0 47 + 50 47 EG=8e-07 E=1.3e-14 mismap=0.00056
+228 chrM 349 50 + 16571 16/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=1.97e-11
+204 chrM 5387 50 + 16571 18/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.66e-09
+226 chrM 5164 50 + 16571 18/2 0 50 + 50 50 EG=5.3e-10 E=5.4e-18 mismap=3.02e-11
+174 chrM 8013 41 + 16571 19/1 9 41 - 50 41 EG=7.7e-05 E=1.5e-12 mismap=3.6e-06
+154 chrM 7826 50 + 16571 19/2 0 50 + 50 50 EG=0.0075 E=1.7e-10 mismap=0.000351
+210 chrM 4644 47 + 16571 20/1 3 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=8.91e-10
+207 chrM 4452 50 + 16571 20/2 0 50 + 50 50 EG=4.1e-08 E=5.5e-16 mismap=1.77e-09
+156 chrM 607 46 + 16571 24/2 0 46 + 50 46 EG=0.0047 E=1.1e-10 mismap=0.000241
+156 chrM 9781 50 + 16571 25/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000297
+199 chrM 9554 49 + 16571 25/2 1 49 + 50 49 EG=2.5e-07 E=3.8e-15 mismap=1.57e-08
+240 chrM 3458 48 + 16571 26/2 2 48 + 50 48 EG=2.2e-11 E=1.7e-19 mismap=1.49e-08
+156 chrM 16102 50 + 16571 31/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000635
+206 chrM 15963 50 + 16571 31/2 0 50 + 50 50 EG=5.1e-08 E=7e-16 mismap=6.77e-09
+300 chrM 459 50 + 16571 34/2 0 50 + 50 50 EG=2.4e-17 E=9.6e-27 mismap=0
+204 chrM 15490 50 + 16571 35/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.01e-09
+150 chrM 15318 49 + 16571 35/2 1 49 + 50 49 EG=0.019 E=4.4e-10 mismap=0.00094
+186 chrM 9210 43 + 16571 36/1 4 43 - 50 43 EG=5e-06 E=8.5e-14 mismap=2.44e-07
+176 chrM 9030 50 + 16571 36/2 0 50 + 50 50 EG=4.9e-05 E=9.2e-13 mismap=2.41e-06
+276 chrM 10929 50 + 16571 39/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=0
+205 chrM 10700 50 + 16571 39/2 0 50 + 50 50 EG=6.5e-08 E=8.9e-16 mismap=4.13e-09
+168 chrM 15071 48 + 16571 43/1 0 48 - 50 48 EG=0.0003 E=6.2e-12 mismap=0.000127
+281 chrM 14777 47 + 16571 43/2 0 47 + 50 47 EG=1.9e-15 E=2.9e-24 mismap=0
+264 chrM 3543 48 + 16571 44/1 0 48 - 50 48 EG=9e-14 E=3.3e-22 mismap=1.42e-14
+154 chrM 3288 27 + 16571 44/2 0 27 + 50 27 EG=0.0075 E=1.7e-10 mismap=0.000861
+204 chrM 358 50 + 16571 47/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=5.09e-09
+215 chrM 130 50 + 16571 47/2 0 50 + 50 50 EG=6.6e-09 E=7.9e-17 mismap=4.1e-10
+222 chrM 419 41 + 16571 48/1 9 41 - 50 41 EG=1.3e-09 E=1.4e-17 mismap=5.94e-11
+156 chrM 207 50 + 16571 48/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000218
+213 chrM 13205 46 + 16571 55/1 4 46 - 50 46 EG=1e-08 E=1.3e-16 mismap=4.97e-10
+156 chrM 13027 50 + 16571 55/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000232
+252 chrM 10892 50 + 16571 59/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+157 chrM 10691 34 + 16571 59/2 16 34 + 50 34 EG=0.0038 E=8.4e-11 mismap=0.000171
+168 chrM 6611 40 + 16571 60/1 0 40 - 50 40 EG=0.0003 E=6.2e-12 mismap=1.34e-05
+224 chrM 6412 48 + 16571 60/2 0 48 + 50 48 EG=8.4e-10 E=8.8e-18 mismap=3.61e-11
+198 chrM 11907 41 + 16571 63/1 1 41 - 50 41 EG=3.2e-07 E=4.8e-15 mismap=5.78e-08
+200 chrM 11765 50 + 16571 63/2 0 50 + 50 50 EG=2e-07 E=3e-15 mismap=3.65e-08
+204 chrM 13054 42 + 16571 64/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=4.34e-09
+198 chrM 12879 44 + 16571 64/2 3 44 + 50 44 EG=3.2e-07 E=4.8e-15 mismap=1.72e-08
+167 chrM 1767 46 + 16571 65/2 2 46 - 50 46 EG=0.00038 E=7.8e-12 mismap=4.04e-05
+228 chrM 16146 46 + 16571 67/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=2.46e-11
+258 chrM 4154 47 + 16571 73/1 2 47 - 50 47 EG=3.5e-13 E=1.6e-21 mismap=2.42e-10
+210 chrM 12100 47 + 16571 74/1 0 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.41e-10
+153 chrM 11892 48 + 16571 74/2 2 48 + 50 48 EG=0.0094 E=2.2e-10 mismap=0.000439
+228 chrM 4854 46 + 16571 75/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=5.37e-11
+199 chrM 4585 44 + 16571 75/2 3 44 + 50 44 EG=2.5e-07 E=3.8e-15 mismap=4.11e-08
+276 chrM 5469 50 + 16571 78/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=3.92e-12
+193 chrM 1833 46 + 16571 82/2 4 46 + 50 46 EG=1e-06 E=1.6e-14 mismap=0.000704
+228 chrM 10231 50 + 16571 84/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=2.33e-07
+203 chrM 11684 50 + 16571 86/2 0 50 + 50 50 EG=1e-07 E=1.4e-15 mismap=4.48e-09
+204 chrM 8530 42 + 16571 87/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=3.52e-09
+175 chrM 8329 49 + 16571 87/2 1 49 + 50 49 EG=6.1e-05 E=1.2e-12 mismap=2.7e-06
+210 chrM 1543 47 + 16571 89/1 0 47 + 50 47 EG=2.1e-08 E=2.7e-16 mismap=1.44e-05
+180 chrM 4414 50 + 16571 91/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=1.02e-06
+183 chrM 4245 46 + 16571 91/2 4 46 + 50 46 EG=9.9e-06 E=1.7e-13 mismap=5.14e-07
+157 chrM 7286 50 + 16571 95/1 0 50 - 50 50 EG=0.0038 E=8.4e-11 mismap=0.000372
+251 chrM 7140 50 + 16571 95/2 0 50 + 50 50 EG=1.8e-12 E=1e-20 mismap=1.71e-13
+180 chrM 9555 50 + 16571 98/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=9.06e-07
+154 chrM 9350 26 + 16571 98/2 0 26 + 50 26 EG=0.0075 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
+156 chrM 11 50 + 16571 4/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000222
+180 chrM 9112 50 + 16571 6/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.89e-07
+220 chrM 8903 47 + 16571 6/2 3 47 + 50 47 EG=2.1e-09 E=2.3e-17 mismap=9.34e-11
+166 chrM 12726 38 + 16571 7/2 0 38 + 50 38 EG=0.00048 E=9.9e-12 mismap=2.31e-05
+180 chrM 14235 50 + 16571 9/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=9.02e-07
+253 chrM 14024 50 + 16571 9/2 0 50 + 50 50 EG=1.1e-12 E=6e-21 mismap=4.26e-14
+252 chrM 14182 50 + 16571 13/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+202 chrM 13997 49 + 16571 13/2 1 49 + 50 49 EG=1.3e-07 E=1.8e-15 mismap=5.85e-09
+210 chrM 4992 39 + 16571 14/2 2 39 + 50 39 EG=2.1e-08 E=2.7e-16 mismap=2.86e-05
+194 chrM 2765 47 + 16571 15/2 0 47 + 50 47 EG=8e-07 E=1.3e-14 mismap=0.00111
+228 chrM 349 50 + 16571 16/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=2.13e-11
+204 chrM 5387 50 + 16571 18/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.24e-09
+226 chrM 5164 50 + 16571 18/2 0 50 + 50 50 EG=5.3e-10 E=5.4e-18 mismap=2.75e-11
+174 chrM 8013 41 + 16571 19/1 9 41 - 50 41 EG=7.7e-05 E=1.5e-12 mismap=3.82e-06
+154 chrM 7826 50 + 16571 19/2 0 50 + 50 50 EG=0.0075 E=1.7e-10 mismap=0.000372
+210 chrM 4644 47 + 16571 20/1 3 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.12e-10
+207 chrM 4452 50 + 16571 20/2 0 50 + 50 50 EG=4.1e-08 E=5.5e-16 mismap=1.81e-09
+156 chrM 607 46 + 16571 24/2 0 46 + 50 46 EG=0.0047 E=1.1e-10 mismap=0.00026
+156 chrM 9781 50 + 16571 25/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000269
+199 chrM 9554 49 + 16571 25/2 1 49 + 50 49 EG=2.5e-07 E=3.8e-15 mismap=1.42e-08
+240 chrM 3458 48 + 16571 26/2 2 48 + 50 48 EG=2.2e-11 E=1.7e-19 mismap=2.97e-08
+156 chrM 16102 50 + 16571 31/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000595
+206 chrM 15963 50 + 16571 31/2 0 50 + 50 50 EG=5.1e-08 E=7e-16 mismap=6.34e-09
+300 chrM 459 50 + 16571 34/2 0 50 + 50 50 EG=2.4e-17 E=9.6e-27 mismap=0
+204 chrM 15490 50 + 16571 35/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.31e-09
+150 chrM 15318 49 + 16571 35/2 1 49 + 50 49 EG=0.019 E=4.4e-10 mismap=0.00101
+186 chrM 9210 43 + 16571 36/1 4 43 - 50 43 EG=5e-06 E=8.5e-14 mismap=2.62e-07
+176 chrM 9030 50 + 16571 36/2 0 50 + 50 50 EG=4.9e-05 E=9.2e-13 mismap=2.59e-06
+276 chrM 10929 50 + 16571 39/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=0
+205 chrM 10700 50 + 16571 39/2 0 50 + 50 50 EG=6.5e-08 E=8.9e-16 mismap=3.73e-09
+168 chrM 15071 48 + 16571 43/1 0 48 - 50 48 EG=0.0003 E=6.2e-12 mismap=0.000191
+281 chrM 14777 47 + 16571 43/2 0 47 + 50 47 EG=1.9e-15 E=2.9e-24 mismap=0
+264 chrM 3543 48 + 16571 44/1 0 48 - 50 48 EG=9e-14 E=3.3e-22 mismap=1.42e-14
+154 chrM 3288 27 + 16571 44/2 0 27 + 50 27 EG=0.0075 E=1.7e-10 mismap=0.00082
+204 chrM 358 50 + 16571 47/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.61e-09
+215 chrM 130 50 + 16571 47/2 0 50 + 50 50 EG=6.6e-09 E=7.9e-17 mismap=3.71e-10
+222 chrM 419 41 + 16571 48/1 9 41 - 50 41 EG=1.3e-09 E=1.4e-17 mismap=5.74e-11
+156 chrM 207 50 + 16571 48/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.00021
+213 chrM 13205 46 + 16571 55/1 4 46 - 50 46 EG=1e-08 E=1.3e-16 mismap=5.32e-10
+156 chrM 13027 50 + 16571 55/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000248
+252 chrM 10892 50 + 16571 59/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+157 chrM 10691 34 + 16571 59/2 16 34 + 50 34 EG=0.0038 E=8.4e-11 mismap=0.000167
+168 chrM 6611 40 + 16571 60/1 0 40 - 50 40 EG=0.0003 E=6.2e-12 mismap=1.37e-05
+224 chrM 6412 48 + 16571 60/2 0 48 + 50 48 EG=8.4e-10 E=8.8e-18 mismap=3.69e-11
+198 chrM 11907 41 + 16571 63/1 1 41 - 50 41 EG=3.2e-07 E=4.8e-15 mismap=4.95e-08
+200 chrM 11765 50 + 16571 63/2 0 50 + 50 50 EG=2e-07 E=3e-15 mismap=3.13e-08
+204 chrM 13054 42 + 16571 64/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=4.69e-09
+198 chrM 12879 44 + 16571 64/2 3 44 + 50 44 EG=3.2e-07 E=4.8e-15 mismap=1.85e-08
+167 chrM 1767 46 + 16571 65/2 2 46 - 50 46 EG=0.00038 E=7.8e-12 mismap=3.79e-05
+228 chrM 16146 46 + 16571 67/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=2.23e-11
+258 chrM 4154 47 + 16571 73/1 2 47 - 50 47 EG=3.5e-13 E=1.6e-21 mismap=4.81e-10
+210 chrM 12100 47 + 16571 74/1 0 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.03e-10
+153 chrM 11892 48 + 16571 74/2 2 48 + 50 48 EG=0.0094 E=2.2e-10 mismap=0.000421
+228 chrM 4854 46 + 16571 75/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=5.56e-11
+199 chrM 4585 44 + 16571 75/2 3 44 + 50 44 EG=2.5e-07 E=3.8e-15 mismap=4.26e-08
+276 chrM 5469 50 + 16571 78/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=7.79e-12
+193 chrM 1833 46 + 16571 82/2 4 46 + 50 46 EG=1e-06 E=1.6e-14 mismap=0.0014
+228 chrM 10231 50 + 16571 84/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=4.63e-07
+203 chrM 11684 50 + 16571 86/2 0 50 + 50 50 EG=1e-07 E=1.4e-15 mismap=4.42e-09
+204 chrM 8530 42 + 16571 87/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=3.54e-09
+175 chrM 8329 49 + 16571 87/2 1 49 + 50 49 EG=6.1e-05 E=1.2e-12 mismap=2.71e-06
+210 chrM 1543 47 + 16571 89/1 0 47 + 50 47 EG=2.1e-08 E=2.7e-16 mismap=2.86e-05
+180 chrM 4414 50 + 16571 91/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=1.1e-06
+183 chrM 4245 46 + 16571 91/2 4 46 + 50 46 EG=9.9e-06 E=1.7e-13 mismap=5.55e-07
+157 chrM 7286 50 + 16571 95/1 0 50 - 50 50 EG=0.0038 E=8.4e-11 mismap=0.000375
+251 chrM 7140 50 + 16571 95/2 0 50 + 50 50 EG=1.8e-12 E=1e-20 mismap=1.71e-13
+180 chrM 9555 50 + 16571 98/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.7e-07
+154 chrM 9350 26 + 16571 98/2 0 26 + 50 26 EG=0.0075 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
+156 chrM 11 50 + 16571 4/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000222
+180 chrM 9112 50 + 16571 6/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.89e-07
+220 chrM 8903 47 + 16571 6/2 3 47 + 50 47 EG=2.1e-09 E=2.3e-17 mismap=9.34e-11
+166 chrM 12726 38 + 16571 7/2 0 38 + 50 38 EG=0.00048 E=9.9e-12 mismap=2.31e-05
+180 chrM 14235 50 + 16571 9/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=9.02e-07
+253 chrM 14024 50 + 16571 9/2 0 50 + 50 50 EG=1.1e-12 E=6e-21 mismap=4.26e-14
+252 chrM 14182 50 + 16571 13/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+202 chrM 13997 49 + 16571 13/2 1 49 + 50 49 EG=1.3e-07 E=1.8e-15 mismap=5.85e-09
+210 chrM 4992 39 + 16571 14/2 2 39 + 50 39 EG=2.1e-08 E=2.7e-16 mismap=0.000288
+228 chrM 349 50 + 16571 16/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=2.13e-11
+204 chrM 5387 50 + 16571 18/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.24e-09
+226 chrM 5164 50 + 16571 18/2 0 50 + 50 50 EG=5.3e-10 E=5.4e-18 mismap=2.75e-11
+174 chrM 8013 41 + 16571 19/1 9 41 - 50 41 EG=7.7e-05 E=1.5e-12 mismap=3.82e-06
+154 chrM 7826 50 + 16571 19/2 0 50 + 50 50 EG=0.0075 E=1.7e-10 mismap=0.000372
+210 chrM 4644 47 + 16571 20/1 3 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.12e-10
+207 chrM 4452 50 + 16571 20/2 0 50 + 50 50 EG=4.1e-08 E=5.5e-16 mismap=1.81e-09
+156 chrM 607 46 + 16571 24/2 0 46 + 50 46 EG=0.0047 E=1.1e-10 mismap=0.00026
+156 chrM 9781 50 + 16571 25/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000269
+199 chrM 9554 49 + 16571 25/2 1 49 + 50 49 EG=2.5e-07 E=3.8e-15 mismap=1.42e-08
+240 chrM 3458 48 + 16571 26/2 2 48 + 50 48 EG=2.2e-11 E=1.7e-19 mismap=2.99e-07
+156 chrM 16102 50 + 16571 31/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000595
+206 chrM 15963 50 + 16571 31/2 0 50 + 50 50 EG=5.1e-08 E=7e-16 mismap=6.34e-09
+300 chrM 459 50 + 16571 34/2 0 50 + 50 50 EG=2.4e-17 E=9.6e-27 mismap=0
+204 chrM 15490 50 + 16571 35/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.31e-09
+150 chrM 15318 49 + 16571 35/2 1 49 + 50 49 EG=0.019 E=4.4e-10 mismap=0.00101
+186 chrM 9210 43 + 16571 36/1 4 43 - 50 43 EG=5e-06 E=8.5e-14 mismap=2.62e-07
+176 chrM 9030 50 + 16571 36/2 0 50 + 50 50 EG=4.9e-05 E=9.2e-13 mismap=2.59e-06
+276 chrM 10929 50 + 16571 39/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=0
+205 chrM 10700 50 + 16571 39/2 0 50 + 50 50 EG=6.5e-08 E=8.9e-16 mismap=3.73e-09
+168 chrM 15071 48 + 16571 43/1 0 48 - 50 48 EG=0.0003 E=6.2e-12 mismap=0.000191
+281 chrM 14777 47 + 16571 43/2 0 47 + 50 47 EG=1.9e-15 E=2.9e-24 mismap=0
+264 chrM 3543 48 + 16571 44/1 0 48 - 50 48 EG=9e-14 E=3.3e-22 mismap=1.42e-14
+154 chrM 3288 27 + 16571 44/2 0 27 + 50 27 EG=0.0075 E=1.7e-10 mismap=0.00082
+204 chrM 358 50 + 16571 47/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.61e-09
+215 chrM 130 50 + 16571 47/2 0 50 + 50 50 EG=6.6e-09 E=7.9e-17 mismap=3.71e-10
+222 chrM 419 41 + 16571 48/1 9 41 - 50 41 EG=1.3e-09 E=1.4e-17 mismap=5.74e-11
+156 chrM 207 50 + 16571 48/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.00021
+213 chrM 13205 46 + 16571 55/1 4 46 - 50 46 EG=1e-08 E=1.3e-16 mismap=5.32e-10
+156 chrM 13027 50 + 16571 55/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000248
+252 chrM 10892 50 + 16571 59/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+157 chrM 10691 34 + 16571 59/2 16 34 + 50 34 EG=0.0038 E=8.4e-11 mismap=0.000167
+168 chrM 6611 40 + 16571 60/1 0 40 - 50 40 EG=0.0003 E=6.2e-12 mismap=1.37e-05
+224 chrM 6412 48 + 16571 60/2 0 48 + 50 48 EG=8.4e-10 E=8.8e-18 mismap=3.69e-11
+198 chrM 11907 41 + 16571 63/1 1 41 - 50 41 EG=3.2e-07 E=4.8e-15 mismap=4.95e-08
+200 chrM 11765 50 + 16571 63/2 0 50 + 50 50 EG=2e-07 E=3e-15 mismap=3.13e-08
+204 chrM 13054 42 + 16571 64/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=4.69e-09
+198 chrM 12879 44 + 16571 64/2 3 44 + 50 44 EG=3.2e-07 E=4.8e-15 mismap=1.86e-08
+167 chrM 1767 46 + 16571 65/2 2 46 - 50 46 EG=0.00038 E=7.8e-12 mismap=3.8e-05
+228 chrM 16146 46 + 16571 67/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=2.23e-11
+258 chrM 4154 47 + 16571 73/1 2 47 - 50 47 EG=3.5e-13 E=1.6e-21 mismap=4.85e-09
+210 chrM 12100 47 + 16571 74/1 0 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.03e-10
+153 chrM 11892 48 + 16571 74/2 2 48 + 50 48 EG=0.0094 E=2.2e-10 mismap=0.000421
+228 chrM 4854 46 + 16571 75/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=5.56e-11
+199 chrM 4585 44 + 16571 75/2 3 44 + 50 44 EG=2.5e-07 E=3.8e-15 mismap=4.26e-08
+276 chrM 5469 50 + 16571 78/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=7.86e-11
+228 chrM 10231 50 + 16571 84/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=4.67e-06
+203 chrM 11684 50 + 16571 86/2 0 50 + 50 50 EG=1e-07 E=1.4e-15 mismap=4.42e-09
+204 chrM 8530 42 + 16571 87/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=3.54e-09
+175 chrM 8329 49 + 16571 87/2 1 49 + 50 49 EG=6.1e-05 E=1.2e-12 mismap=2.71e-06
+210 chrM 1543 47 + 16571 89/1 0 47 + 50 47 EG=2.1e-08 E=2.7e-16 mismap=0.000288
+180 chrM 4414 50 + 16571 91/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=1.1e-06
+183 chrM 4245 46 + 16571 91/2 4 46 + 50 46 EG=9.9e-06 E=1.7e-13 mismap=5.55e-07
+157 chrM 7286 50 + 16571 95/1 0 50 - 50 50 EG=0.0038 E=8.4e-11 mismap=0.000375
+251 chrM 7140 50 + 16571 95/2 0 50 + 50 50 EG=1.8e-12 E=1e-20 mismap=1.71e-13
+180 chrM 9555 50 + 16571 98/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.7e-07
+154 chrM 9350 26 + 16571 98/2 0 26 + 50 26 EG=0.0075 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
+156 chrM 11 50 + 16571 4/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000222
+180 chrM 9112 50 + 16571 6/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.89e-07
+220 chrM 8903 47 + 16571 6/2 3 47 + 50 47 EG=2.1e-09 E=2.3e-17 mismap=9.34e-11
+166 chrM 12726 38 + 16571 7/2 0 38 + 50 38 EG=0.00048 E=9.9e-12 mismap=2.31e-05
+180 chrM 14235 50 + 16571 9/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=9.02e-07
+253 chrM 14024 50 + 16571 9/2 0 50 + 50 50 EG=1.1e-12 E=6e-21 mismap=4.26e-14
+252 chrM 14182 50 + 16571 13/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+202 chrM 13997 49 + 16571 13/2 1 49 + 50 49 EG=1.3e-07 E=1.8e-15 mismap=5.85e-09
+210 chrM 4992 39 + 16571 14/2 2 39 + 50 39 EG=2.1e-08 E=2.7e-16 mismap=2.86e-05
+194 chrM 2765 47 + 16571 15/2 0 47 + 50 47 EG=8e-07 E=1.3e-14 mismap=0.00111
+228 chrM 349 50 + 16571 16/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=2.13e-11
+204 chrM 5387 50 + 16571 18/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.24e-09
+226 chrM 5164 50 + 16571 18/2 0 50 + 50 50 EG=5.3e-10 E=5.4e-18 mismap=2.75e-11
+174 chrM 8013 41 + 16571 19/1 9 41 - 50 41 EG=7.7e-05 E=1.5e-12 mismap=3.82e-06
+154 chrM 7826 50 + 16571 19/2 0 50 + 50 50 EG=0.0075 E=1.7e-10 mismap=0.000372
+210 chrM 4644 47 + 16571 20/1 3 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.12e-10
+207 chrM 4452 50 + 16571 20/2 0 50 + 50 50 EG=4.1e-08 E=5.5e-16 mismap=1.81e-09
+156 chrM 607 46 + 16571 24/2 0 46 + 50 46 EG=0.0047 E=1.1e-10 mismap=0.00026
+156 chrM 9781 50 + 16571 25/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000269
+199 chrM 9554 49 + 16571 25/2 1 49 + 50 49 EG=2.5e-07 E=3.8e-15 mismap=1.42e-08
+240 chrM 3458 48 + 16571 26/2 2 48 + 50 48 EG=2.2e-11 E=1.7e-19 mismap=2.97e-08
+156 chrM 16102 50 + 16571 31/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000595
+206 chrM 15963 50 + 16571 31/2 0 50 + 50 50 EG=5.1e-08 E=7e-16 mismap=6.34e-09
+300 chrM 459 50 + 16571 34/2 0 50 + 50 50 EG=2.4e-17 E=9.6e-27 mismap=0
+204 chrM 15490 50 + 16571 35/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.31e-09
+150 chrM 15318 49 + 16571 35/2 1 49 + 50 49 EG=0.019 E=4.4e-10 mismap=0.00101
+186 chrM 9210 43 + 16571 36/1 4 43 - 50 43 EG=5e-06 E=8.5e-14 mismap=2.62e-07
+176 chrM 9030 50 + 16571 36/2 0 50 + 50 50 EG=4.9e-05 E=9.2e-13 mismap=2.59e-06
+276 chrM 10929 50 + 16571 39/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=0
+205 chrM 10700 50 + 16571 39/2 0 50 + 50 50 EG=6.5e-08 E=8.9e-16 mismap=3.73e-09
+168 chrM 15071 48 + 16571 43/1 0 48 - 50 48 EG=0.0003 E=6.2e-12 mismap=0.000191
+281 chrM 14777 47 + 16571 43/2 0 47 + 50 47 EG=1.9e-15 E=2.9e-24 mismap=0
+264 chrM 3543 48 + 16571 44/1 0 48 - 50 48 EG=9e-14 E=3.3e-22 mismap=1.42e-14
+154 chrM 3288 27 + 16571 44/2 0 27 + 50 27 EG=0.0075 E=1.7e-10 mismap=0.00082
+204 chrM 358 50 + 16571 47/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.61e-09
+215 chrM 130 50 + 16571 47/2 0 50 + 50 50 EG=6.6e-09 E=7.9e-17 mismap=3.71e-10
+222 chrM 419 41 + 16571 48/1 9 41 - 50 41 EG=1.3e-09 E=1.4e-17 mismap=5.74e-11
+156 chrM 207 50 + 16571 48/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.00021
+213 chrM 13205 46 + 16571 55/1 4 46 - 50 46 EG=1e-08 E=1.3e-16 mismap=5.32e-10
+156 chrM 13027 50 + 16571 55/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000248
+252 chrM 10892 50 + 16571 59/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+157 chrM 10691 34 + 16571 59/2 16 34 + 50 34 EG=0.0038 E=8.4e-11 mismap=0.000167
+168 chrM 6611 40 + 16571 60/1 0 40 - 50 40 EG=0.0003 E=6.2e-12 mismap=1.37e-05
+224 chrM 6412 48 + 16571 60/2 0 48 + 50 48 EG=8.4e-10 E=8.8e-18 mismap=3.69e-11
+198 chrM 11907 41 + 16571 63/1 1 41 - 50 41 EG=3.2e-07 E=4.8e-15 mismap=4.95e-08
+200 chrM 11765 50 + 16571 63/2 0 50 + 50 50 EG=2e-07 E=3e-15 mismap=3.13e-08
+204 chrM 13054 42 + 16571 64/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=4.69e-09
+198 chrM 12879 44 + 16571 64/2 3 44 + 50 44 EG=3.2e-07 E=4.8e-15 mismap=1.85e-08
+167 chrM 1767 46 + 16571 65/2 2 46 - 50 46 EG=0.00038 E=7.8e-12 mismap=3.79e-05
+228 chrM 16146 46 + 16571 67/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=2.23e-11
+258 chrM 4154 47 + 16571 73/1 2 47 - 50 47 EG=3.5e-13 E=1.6e-21 mismap=4.81e-10
+210 chrM 12100 47 + 16571 74/1 0 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.03e-10
+153 chrM 11892 48 + 16571 74/2 2 48 + 50 48 EG=0.0094 E=2.2e-10 mismap=0.000421
+228 chrM 4854 46 + 16571 75/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=5.56e-11
+199 chrM 4585 44 + 16571 75/2 3 44 + 50 44 EG=2.5e-07 E=3.8e-15 mismap=4.26e-08
+276 chrM 5469 50 + 16571 78/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=7.79e-12
+193 chrM 1833 46 + 16571 82/2 4 46 + 50 46 EG=1e-06 E=1.6e-14 mismap=0.0014
+228 chrM 10231 50 + 16571 84/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=4.63e-07
+203 chrM 11684 50 + 16571 86/2 0 50 + 50 50 EG=1e-07 E=1.4e-15 mismap=4.42e-09
+204 chrM 8530 42 + 16571 87/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=3.54e-09
+175 chrM 8329 49 + 16571 87/2 1 49 + 50 49 EG=6.1e-05 E=1.2e-12 mismap=2.71e-06
+210 chrM 1543 47 + 16571 89/1 0 47 + 50 47 EG=2.1e-08 E=2.7e-16 mismap=2.86e-05
+180 chrM 4414 50 + 16571 91/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=1.1e-06
+183 chrM 4245 46 + 16571 91/2 4 46 + 50 46 EG=9.9e-06 E=1.7e-13 mismap=5.55e-07
+157 chrM 7286 50 + 16571 95/1 0 50 - 50 50 EG=0.0038 E=8.4e-11 mismap=0.000375
+251 chrM 7140 50 + 16571 95/2 0 50 + 50 50 EG=1.8e-12 E=1e-20 mismap=1.71e-13
+180 chrM 9555 50 + 16571 98/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.7e-07
+154 chrM 9350 26 + 16571 98/2 0 26 + 50 26 EG=0.0075 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
+156 chrM 11 50 + 16571 4/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000222
+180 chrM 9112 50 + 16571 6/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.89e-07
+220 chrM 8903 47 + 16571 6/2 3 47 + 50 47 EG=2.1e-09 E=2.3e-17 mismap=9.34e-11
+166 chrM 12726 38 + 16571 7/2 0 38 + 50 38 EG=0.00048 E=9.9e-12 mismap=2.31e-05
+180 chrM 14235 50 + 16571 9/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=9.02e-07
+253 chrM 14024 50 + 16571 9/2 0 50 + 50 50 EG=1.1e-12 E=6e-21 mismap=4.26e-14
+252 chrM 14182 50 + 16571 13/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+202 chrM 13997 49 + 16571 13/2 1 49 + 50 49 EG=1.3e-07 E=1.8e-15 mismap=5.85e-09
+210 chrM 4992 39 + 16571 14/2 2 39 + 50 39 EG=2.1e-08 E=2.7e-16 mismap=2.86e-05
+194 chrM 2765 47 + 16571 15/2 0 47 + 50 47 EG=8e-07 E=1.3e-14 mismap=0.00111
+228 chrM 349 50 + 16571 16/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=2.13e-11
+204 chrM 5387 50 + 16571 18/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.24e-09
+226 chrM 5164 50 + 16571 18/2 0 50 + 50 50 EG=5.3e-10 E=5.4e-18 mismap=2.75e-11
+174 chrM 8013 41 + 16571 19/1 9 41 - 50 41 EG=7.7e-05 E=1.5e-12 mismap=3.82e-06
+154 chrM 7826 50 + 16571 19/2 0 50 + 50 50 EG=0.0075 E=1.7e-10 mismap=0.000372
+210 chrM 4644 47 + 16571 20/1 3 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.12e-10
+207 chrM 4452 50 + 16571 20/2 0 50 + 50 50 EG=4.1e-08 E=5.5e-16 mismap=1.81e-09
+156 chrM 607 46 + 16571 24/2 0 46 + 50 46 EG=0.0047 E=1.1e-10 mismap=0.00026
+156 chrM 9781 50 + 16571 25/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000269
+199 chrM 9554 49 + 16571 25/2 1 49 + 50 49 EG=2.5e-07 E=3.8e-15 mismap=1.42e-08
+240 chrM 3458 48 + 16571 26/2 2 48 + 50 48 EG=2.2e-11 E=1.7e-19 mismap=2.97e-08
+156 chrM 16102 50 + 16571 31/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000595
+206 chrM 15963 50 + 16571 31/2 0 50 + 50 50 EG=5.1e-08 E=7e-16 mismap=6.34e-09
+300 chrM 459 50 + 16571 34/2 0 50 + 50 50 EG=2.4e-17 E=9.6e-27 mismap=0
+204 chrM 15490 50 + 16571 35/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.31e-09
+150 chrM 15318 49 + 16571 35/2 1 49 + 50 49 EG=0.019 E=4.4e-10 mismap=0.00101
+186 chrM 9210 43 + 16571 36/1 4 43 - 50 43 EG=5e-06 E=8.5e-14 mismap=2.62e-07
+176 chrM 9030 50 + 16571 36/2 0 50 + 50 50 EG=4.9e-05 E=9.2e-13 mismap=2.59e-06
+276 chrM 10929 50 + 16571 39/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=0
+205 chrM 10700 50 + 16571 39/2 0 50 + 50 50 EG=6.5e-08 E=8.9e-16 mismap=3.73e-09
+168 chrM 15071 48 + 16571 43/1 0 48 - 50 48 EG=0.0003 E=6.2e-12 mismap=0.000191
+281 chrM 14777 47 + 16571 43/2 0 47 + 50 47 EG=1.9e-15 E=2.9e-24 mismap=0
+264 chrM 3543 48 + 16571 44/1 0 48 - 50 48 EG=9e-14 E=3.3e-22 mismap=1.42e-14
+154 chrM 3288 27 + 16571 44/2 0 27 + 50 27 EG=0.0075 E=1.7e-10 mismap=0.00082
+204 chrM 358 50 + 16571 47/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.61e-09
+215 chrM 130 50 + 16571 47/2 0 50 + 50 50 EG=6.6e-09 E=7.9e-17 mismap=3.71e-10
+222 chrM 419 41 + 16571 48/1 9 41 - 50 41 EG=1.3e-09 E=1.4e-17 mismap=5.74e-11
+156 chrM 207 50 + 16571 48/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.00021
+213 chrM 13205 46 + 16571 55/1 4 46 - 50 46 EG=1e-08 E=1.3e-16 mismap=5.32e-10
+156 chrM 13027 50 + 16571 55/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000248
+252 chrM 10892 50 + 16571 59/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+157 chrM 10691 34 + 16571 59/2 16 34 + 50 34 EG=0.0038 E=8.4e-11 mismap=0.000167
+168 chrM 6611 40 + 16571 60/1 0 40 - 50 40 EG=0.0003 E=6.2e-12 mismap=1.37e-05
+224 chrM 6412 48 + 16571 60/2 0 48 + 50 48 EG=8.4e-10 E=8.8e-18 mismap=3.69e-11
+198 chrM 11907 41 + 16571 63/1 1 41 - 50 41 EG=3.2e-07 E=4.8e-15 mismap=4.95e-08
+200 chrM 11765 50 + 16571 63/2 0 50 + 50 50 EG=2e-07 E=3e-15 mismap=3.13e-08
+204 chrM 13054 42 + 16571 64/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=4.69e-09
+198 chrM 12879 44 + 16571 64/2 3 44 + 50 44 EG=3.2e-07 E=4.8e-15 mismap=1.85e-08
+167 chrM 1767 46 + 16571 65/2 2 46 - 50 46 EG=0.00038 E=7.8e-12 mismap=3.79e-05
+228 chrM 16146 46 + 16571 67/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=2.23e-11
+258 chrM 4154 47 + 16571 73/1 2 47 - 50 47 EG=3.5e-13 E=1.6e-21 mismap=4.81e-10
+210 chrM 12100 47 + 16571 74/1 0 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.03e-10
+153 chrM 11892 48 + 16571 74/2 2 48 + 50 48 EG=0.0094 E=2.2e-10 mismap=0.000421
+228 chrM 4854 46 + 16571 75/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=5.56e-11
+199 chrM 4585 44 + 16571 75/2 3 44 + 50 44 EG=2.5e-07 E=3.8e-15 mismap=4.26e-08
+276 chrM 5469 50 + 16571 78/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=7.79e-12
+193 chrM 1833 46 + 16571 82/2 4 46 + 50 46 EG=1e-06 E=1.6e-14 mismap=0.0014
+228 chrM 10231 50 + 16571 84/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=4.63e-07
+203 chrM 11684 50 + 16571 86/2 0 50 + 50 50 EG=1e-07 E=1.4e-15 mismap=4.42e-09
+204 chrM 8530 42 + 16571 87/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=3.54e-09
+175 chrM 8329 49 + 16571 87/2 1 49 + 50 49 EG=6.1e-05 E=1.2e-12 mismap=2.71e-06
+210 chrM 1543 47 + 16571 89/1 0 47 + 50 47 EG=2.1e-08 E=2.7e-16 mismap=2.86e-05
+180 chrM 4414 50 + 16571 91/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=1.1e-06
+183 chrM 4245 46 + 16571 91/2 4 46 + 50 46 EG=9.9e-06 E=1.7e-13 mismap=5.55e-07
+157 chrM 7286 50 + 16571 95/1 0 50 - 50 50 EG=0.0038 E=8.4e-11 mismap=0.000375
+251 chrM 7140 50 + 16571 95/2 0 50 + 50 50 EG=1.8e-12 E=1e-20 mismap=1.71e-13
+180 chrM 9555 50 + 16571 98/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.7e-07
+154 chrM 9350 26 + 16571 98/2 0 26 + 50 26 EG=0.0075 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
+a score=156 EG=0.0047 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 at CBCACCCCCCCC@ABCCAABBCCBC at A=B at CBCACCC@
-a score=180 EG2=0.59 E=3.6e-13 mismap=8.89e-07
+a score=180 EG=2e-05 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
+a score=220 EG=2.1e-09 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 at B?B<AB>A<AB:@AB9;=@AA43A>@;,5=@@A
-a score=166 EG2=15 E=9.9e-12 mismap=2.31e-05
+a score=166 EG=0.00048 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
+a score=180 EG=2e-05 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
+a score=253 EG=1.1e-12 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
+a score=252 EG=1.4e-12 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
+a score=202 EG=1.3e-07 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 at 9?@AA9=@BB;@>@A8><2?99A@
-a score=210 EG2=0.00062 E=2.7e-16 mismap=2.86e-05
+a score=210 EG=2.1e-08 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 at C+BCBBBCC<@B9BB
-a score=194 EG2=0.024 E=1.3e-14 mismap=0.00111
+a score=194 EG=8e-07 E=1.3e-14 mismap=0.00111
s chrM 2765 47 + 16571 AACAAACCCACAGGTCCTAAACTACCAAACCTGCATTAAAAATTTCG
s 15/2 0 47 + 50 AACAGACCCACAAATCCTAAACTACCAAACCTACATTAAAAATTTCG
q 15/2 2 at 1>(?@)@C>(>CCC>)@1:15=>9@=ACCC@>A?*=B?8,>4 at BC
-a score=228 EG2=1e-05 E=3.3e-18 mismap=2.13e-11
+a score=228 EG=3.4e-10 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
+a score=204 EG=8.1e-08 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
+a score=226 EG=5.3e-10 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 at C=9CBC,CA=-C at CB=CA>9?/
-a score=174 EG2=2.3 E=1.5e-12 mismap=3.82e-06
+a score=174 EG=7.7e-05 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
+a score=154 EG=0.0075 E=1.7e-10 mismap=0.000372
s chrM 7826 50 + 16571 CTACGCATCCTTTACATAACAGACGAGGTCAACGATCCCTCCCTTACCAT
s 19/2 0 50 + 50 CTACACATACTTTACATAACAAACGAAATCAAGAATCCCTCCCTTACCAT
q 19/2 =@9 at B>C<';>)*8?'>67 at 1;(?<;A<@:@@(9A@@)(38>9+@<8?6@
-a score=210 EG2=0.00062 E=2.7e-16 mismap=9.12e-10
+a score=210 EG=2.1e-08 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
+a score=207 EG=4.1e-08 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 at CBC:1?>C?-2AA?59'<)?B7A46)@>:*?8,5)<24?)5>?;
-a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.00026
+a score=156 EG=0.0047 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
+a score=156 EG=0.0047 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
+a score=199 EG=2.5e-07 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 at 25)=CB9B7;BC7>?A4A;5 at BA+;6)()1BCA748@)=1?
-a score=240 EG2=6.6e-07 E=1.7e-19 mismap=2.97e-08
+a score=240 EG=2.2e-11 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 at .9,BA<:=?>-B6.)*=>0?59;'3
-a score=156 EG2=1.4e+02 E=1.1e-10 mismap=0.000595
+a score=156 EG=0.0047 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
+a score=206 EG=5.1e-08 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 at A@?A at A<@@BAA>?.9,9A?A<@?A
-a score=300 EG2=7.3e-13 E=9.6e-27 mismap=0
+a score=300 EG=2.4e-17 E=9.6e-27 mismap=0
s chrM 459 50 + 16571 CCTCCCACTCCCATACTACTAATCTCATCAATACAACCCCCGCCCATCCT
s 34/2 0 50 + 50 CCTCCCACTCCCATACTACTAATCTCATCAATACAACCCCCGCCCATCCT
q 34/2 BCABBBBCBBBCAACB;BCAB at BBBBBCABBA?AB>B at 5@B at BA5BCB7?
-a score=204 EG2=0.0025 E=1.1e-15 mismap=4.31e-09
+a score=204 EG=8.1e-08 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
+a score=150 EG=0.019 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 at B=CCCCBBB at B@>@@BBACB?@B at CB@CBACC>A4BBACAA
-a score=186 EG2=0.15 E=8.5e-14 mismap=2.62e-07
+a score=186 EG=5e-06 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 at DEEFBD7EEE?FFFEEE:EFFFEE??D?B9EB?D5CEEE?C
-a score=176 EG2=1.5 E=9.2e-13 mismap=2.59e-06
+a score=176 EG=4.9e-05 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
+a score=276 EG=5.8e-15 E=1.2e-23 mismap=0
s chrM 10929 50 + 16571 CCTCCGACCCCCTAACAACCCCCCTCCTAATACTAACTACCTGACTCCTA
s 39/1 0 50 - 50 CCTCCGACCCCCTAACAACCCCCCTCCTAATACTAACTACCTAACTCCTA
q 39/1 ?CBDD3;@4 at B>>,>C7@>>.?=BB>7?4 at C8>C3EFFFFFFFFFFDFFF
-a score=205 EG2=0.002 E=8.9e-16 mismap=3.73e-09
+a score=205 EG=6.5e-08 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 at BABBB=@BABBAAA at BBAA>;0===8>=;7.5:<;@@@A
-a score=168 EG2=9.2 E=6.2e-12 mismap=0.000191
+a score=168 EG=0.0003 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
+a score=281 EG=1.9e-15 E=2.9e-24 mismap=0
s chrM 14777 47 + 16571 ATAAAATTAATTAACCACTCATTCATCGACCTCCCCACCCCATCCAA
s 43/2 0 47 + 50 ATAAAATTAATTAACCACTCATTCATCGACCTCCCCACCCCATCCAA
q 43/2 BC at BCCCC=B=BBB at B@A at B?B at B9AAB>>BBBA):B<>??05:?7@
-a score=264 EG2=2.7e-09 E=3.3e-22 mismap=1.42e-14
+a score=264 EG=9e-14 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
+a score=154 EG=0.0075 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 at -7'':0%'@-B
-a score=204 EG2=0.0025 E=1.1e-15 mismap=4.61e-09
+a score=204 EG=8.1e-08 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 at CB;;@
-a score=215 EG2=0.0002 E=7.9e-17 mismap=3.71e-10
+a score=215 EG=6.6e-09 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 at CB6BCCB9?BAB=CCB@<B;C<CA at C@,9CCBAB>?=BBA
-a score=222 EG2=4e-05 E=1.4e-17 mismap=5.74e-11
+a score=222 EG=1.3e-09 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
+a score=156 EG=0.0047 E=1.1e-10 mismap=0.00021
s chrM 207 50 + 16571 TTAATtaattaaTGCTTGTAGGACATAATAATAACAATTGAATGTCTGCA
s 48/2 0 50 + 50 TTAATTAATTAATACTTATAAAACATAAtaataaCAATTAAATATCTGCA
q 48/2 BCBBB at 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
+a score=213 EG=1e-08 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
+a score=156 EG=0.0047 E=1.1e-10 mismap=0.000248
s chrM 13027 50 + 16571 CCCTGACTCCCCTCAGCCATAGAAGGCCCCACCCCAGTCTCAGCCCTACT
s 55/2 0 50 + 50 CCCTAACTCCCCTCAACCATAAAAAACCCCACCCCAATCTCAGCCCTACT
q 55/2 BBCBBBBC at C>B at B?@BBBBABAAB?A at B@BB>BBB?AACA@@BABBBBA
-a score=252 EG2=4.2e-08 E=7.7e-21 mismap=5.68e-14
+a score=252 EG=1.4e-12 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
+a score=157 EG=0.0038 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 at 5A99>024)?B>A
-a score=168 EG2=9.2 E=6.2e-12 mismap=1.37e-05
+a score=168 EG=0.0003 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
+a score=224 EG=8.4e-10 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
+a score=198 EG=3.2e-07 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 at CAB1=;5DD:DDDD
-a score=200 EG2=0.0061 E=3e-15 mismap=3.13e-08
+a score=200 EG=2e-07 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 at 0
-a score=204 EG2=0.0025 E=1.1e-15 mismap=4.69e-09
+a score=204 EG=8.1e-08 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 at FDGGAFFDGECCE@GDGDFEGDGBFG
-a score=198 EG2=0.0097 E=4.8e-15 mismap=1.85e-08
+a score=198 EG=3.2e-07 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 at CCACCCC=CCCC0CCCCCCACCCBCCCCCCC at CCBC>C at C
-a score=167 EG2=12 E=7.8e-12 mismap=3.79e-05
+a score=167 EG=0.00038 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
+a score=228 EG=3.4e-10 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
+a score=258 EG=3.5e-13 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
+a score=210 EG=2.1e-08 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
+a score=153 EG=0.0094 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
+a score=228 EG=3.4e-10 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
+a score=199 EG=2.5e-07 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 at B@7 at .CA9AACBB==.B
-a score=276 EG2=1.8e-10 E=1.2e-23 mismap=7.79e-12
+a score=276 EG=5.8e-15 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
+a score=193 EG=1e-06 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 at 4'AC at ABCCBCCBACBA@B at BCB7<
-a score=228 EG2=1e-05 E=3.3e-18 mismap=4.63e-07
+a score=228 EG=3.4e-10 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 at EEEEEEEEEEEEEE=EEEBEEAE:D
-a score=203 EG2=0.0031 E=1.4e-15 mismap=4.42e-09
+a score=203 EG=1e-07 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 at B?C9BB.<?BC=:CB=BBC=BAB*):C>=='<CBB
-a score=204 EG2=0.0025 E=1.1e-15 mismap=3.54e-09
+a score=204 EG=8.1e-08 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
+a score=175 EG=6.1e-05 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 at BBA@C=AB at BB=?@BB?BBBBABACAA4A?BBB?A
-a score=210 EG2=0.00062 E=2.7e-16 mismap=2.86e-05
+a score=210 EG=2.1e-08 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
+a score=180 EG=2e-05 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
+a score=183 EG=9.9e-06 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 at 4:BABBB==B at BBB@(@59A at BAA@B9>B??32B<A<B=;
-a score=157 EG2=1.1e+02 E=8.4e-11 mismap=0.000375
+a score=157 EG=0.0038 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
+a score=251 EG=1.8e-12 E=1e-20 mismap=1.71e-13
s chrM 7140 50 + 16571 CATTTCACTATCATATTCATCGGCGTAAATCTAACTTTCTTCCCACAACA
s 95/2 0 50 + 50 CATTTCACTATCATATTCATCGACATAAATCTAACTTTCTTCCCACAACA
q 95/2 B at BABA>BB=.=?*@@?2<A?6BA=*?@>>8:?0;=6>@A8?76 at 7<?@?
-a score=180 EG2=0.59 E=3.6e-13 mismap=8.7e-07
+a score=180 EG=2e-05 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
+a score=154 EG=0.0075 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 at 7B+AA<<A<
@@ -1157,71 +1157,71 @@ q 98/2 9BBBA6=?8*<(;BAB at 7B+AA<<A<
# 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
+156 chrM 11 50 + 16571 4/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000222
+180 chrM 9112 50 + 16571 6/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.89e-07
+220 chrM 8903 47 + 16571 6/2 3 47 + 50 47 EG=2.1e-09 E=2.3e-17 mismap=9.34e-11
+166 chrM 12726 38 + 16571 7/2 0 38 + 50 38 EG=0.00048 E=9.9e-12 mismap=2.31e-05
+180 chrM 14235 50 + 16571 9/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=9.02e-07
+253 chrM 14024 50 + 16571 9/2 0 50 + 50 50 EG=1.1e-12 E=6e-21 mismap=4.26e-14
+252 chrM 14182 50 + 16571 13/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+202 chrM 13997 49 + 16571 13/2 1 49 + 50 49 EG=1.3e-07 E=1.8e-15 mismap=5.85e-09
+210 chrM 4992 39 + 16571 14/2 2 39 + 50 39 EG=2.1e-08 E=2.7e-16 mismap=2.86e-05
+194 chrM 2765 47 + 16571 15/2 0 47 + 50 47 EG=8e-07 E=1.3e-14 mismap=0.00111
+228 chrM 349 50 + 16571 16/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=2.13e-11
+204 chrM 5387 50 + 16571 18/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.24e-09
+226 chrM 5164 50 + 16571 18/2 0 50 + 50 50 EG=5.3e-10 E=5.4e-18 mismap=2.75e-11
+174 chrM 8013 41 + 16571 19/1 9 41 - 50 41 EG=7.7e-05 E=1.5e-12 mismap=3.82e-06
+154 chrM 7826 50 + 16571 19/2 0 50 + 50 50 EG=0.0075 E=1.7e-10 mismap=0.000372
+210 chrM 4644 47 + 16571 20/1 3 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.12e-10
+207 chrM 4452 50 + 16571 20/2 0 50 + 50 50 EG=4.1e-08 E=5.5e-16 mismap=1.81e-09
+156 chrM 607 46 + 16571 24/2 0 46 + 50 46 EG=0.0047 E=1.1e-10 mismap=0.00026
+156 chrM 9781 50 + 16571 25/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000269
+199 chrM 9554 49 + 16571 25/2 1 49 + 50 49 EG=2.5e-07 E=3.8e-15 mismap=1.42e-08
+240 chrM 3458 48 + 16571 26/2 2 48 + 50 48 EG=2.2e-11 E=1.7e-19 mismap=2.97e-08
+156 chrM 16102 50 + 16571 31/1 0 50 - 50 50 EG=0.0047 E=1.1e-10 mismap=0.000595
+206 chrM 15963 50 + 16571 31/2 0 50 + 50 50 EG=5.1e-08 E=7e-16 mismap=6.34e-09
+300 chrM 459 50 + 16571 34/2 0 50 + 50 50 EG=2.4e-17 E=9.6e-27 mismap=0
+204 chrM 15490 50 + 16571 35/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.31e-09
+150 chrM 15318 49 + 16571 35/2 1 49 + 50 49 EG=0.019 E=4.4e-10 mismap=0.00101
+186 chrM 9210 43 + 16571 36/1 4 43 - 50 43 EG=5e-06 E=8.5e-14 mismap=2.62e-07
+176 chrM 9030 50 + 16571 36/2 0 50 + 50 50 EG=4.9e-05 E=9.2e-13 mismap=2.59e-06
+276 chrM 10929 50 + 16571 39/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=0
+205 chrM 10700 50 + 16571 39/2 0 50 + 50 50 EG=6.5e-08 E=8.9e-16 mismap=3.73e-09
+168 chrM 15071 48 + 16571 43/1 0 48 - 50 48 EG=0.0003 E=6.2e-12 mismap=0.000191
+281 chrM 14777 47 + 16571 43/2 0 47 + 50 47 EG=1.9e-15 E=2.9e-24 mismap=0
+264 chrM 3543 48 + 16571 44/1 0 48 - 50 48 EG=9e-14 E=3.3e-22 mismap=1.42e-14
+154 chrM 3288 27 + 16571 44/2 0 27 + 50 27 EG=0.0075 E=1.7e-10 mismap=0.00082
+204 chrM 358 50 + 16571 47/1 0 50 - 50 50 EG=8.1e-08 E=1.1e-15 mismap=4.61e-09
+215 chrM 130 50 + 16571 47/2 0 50 + 50 50 EG=6.6e-09 E=7.9e-17 mismap=3.71e-10
+222 chrM 419 41 + 16571 48/1 9 41 - 50 41 EG=1.3e-09 E=1.4e-17 mismap=5.74e-11
+156 chrM 207 50 + 16571 48/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.00021
+213 chrM 13205 46 + 16571 55/1 4 46 - 50 46 EG=1e-08 E=1.3e-16 mismap=5.32e-10
+156 chrM 13027 50 + 16571 55/2 0 50 + 50 50 EG=0.0047 E=1.1e-10 mismap=0.000248
+252 chrM 10892 50 + 16571 59/1 0 50 - 50 50 EG=1.4e-12 E=7.7e-21 mismap=5.68e-14
+157 chrM 10691 34 + 16571 59/2 16 34 + 50 34 EG=0.0038 E=8.4e-11 mismap=0.000167
+168 chrM 6611 40 + 16571 60/1 0 40 - 50 40 EG=0.0003 E=6.2e-12 mismap=1.37e-05
+224 chrM 6412 48 + 16571 60/2 0 48 + 50 48 EG=8.4e-10 E=8.8e-18 mismap=3.69e-11
+198 chrM 11907 41 + 16571 63/1 1 41 - 50 41 EG=3.2e-07 E=4.8e-15 mismap=4.95e-08
+200 chrM 11765 50 + 16571 63/2 0 50 + 50 50 EG=2e-07 E=3e-15 mismap=3.13e-08
+204 chrM 13054 42 + 16571 64/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=4.69e-09
+198 chrM 12879 44 + 16571 64/2 3 44 + 50 44 EG=3.2e-07 E=4.8e-15 mismap=1.85e-08
+167 chrM 1767 46 + 16571 65/2 2 46 - 50 46 EG=0.00038 E=7.8e-12 mismap=3.79e-05
+228 chrM 16146 46 + 16571 67/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=2.23e-11
+258 chrM 4154 47 + 16571 73/1 2 47 - 50 47 EG=3.5e-13 E=1.6e-21 mismap=4.81e-10
+210 chrM 12100 47 + 16571 74/1 0 47 - 50 47 EG=2.1e-08 E=2.7e-16 mismap=9.03e-10
+153 chrM 11892 48 + 16571 74/2 2 48 + 50 48 EG=0.0094 E=2.2e-10 mismap=0.000421
+228 chrM 4854 46 + 16571 75/1 4 46 - 50 46 EG=3.4e-10 E=3.3e-18 mismap=5.56e-11
+199 chrM 4585 44 + 16571 75/2 3 44 + 50 44 EG=2.5e-07 E=3.8e-15 mismap=4.26e-08
+276 chrM 5469 50 + 16571 78/1 0 50 - 50 50 EG=5.8e-15 E=1.2e-23 mismap=7.79e-12
+193 chrM 1833 46 + 16571 82/2 4 46 + 50 46 EG=1e-06 E=1.6e-14 mismap=0.0014
+228 chrM 10231 50 + 16571 84/1 0 50 - 50 50 EG=3.4e-10 E=3.3e-18 mismap=4.63e-07
+203 chrM 11684 50 + 16571 86/2 0 50 + 50 50 EG=1e-07 E=1.4e-15 mismap=4.42e-09
+204 chrM 8530 42 + 16571 87/1 0 42 - 50 42 EG=8.1e-08 E=1.1e-15 mismap=3.54e-09
+175 chrM 8329 49 + 16571 87/2 1 49 + 50 49 EG=6.1e-05 E=1.2e-12 mismap=2.71e-06
+210 chrM 1543 47 + 16571 89/1 0 47 + 50 47 EG=2.1e-08 E=2.7e-16 mismap=2.86e-05
+180 chrM 4414 50 + 16571 91/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=1.1e-06
+183 chrM 4245 46 + 16571 91/2 4 46 + 50 46 EG=9.9e-06 E=1.7e-13 mismap=5.55e-07
+157 chrM 7286 50 + 16571 95/1 0 50 - 50 50 EG=0.0038 E=8.4e-11 mismap=0.000375
+251 chrM 7140 50 + 16571 95/2 0 50 + 50 50 EG=1.8e-12 E=1e-20 mismap=1.71e-13
+180 chrM 9555 50 + 16571 98/1 0 50 - 50 50 EG=2e-05 E=3.6e-13 mismap=8.7e-07
+154 chrM 9350 26 + 16571 98/2 0 26 + 50 26 EG=0.0075 E=1.7e-10 mismap=0.000335
=====================================
test/last-test.out
=====================================
The diff for this file was not included because it is too large.
View it on GitLab: https://salsa.debian.org/med-team/last-align/-/commit/f7f81f83f2d317b0af14c0c070fd4ef5b357775d
--
View it on GitLab: https://salsa.debian.org/med-team/last-align/-/commit/f7f81f83f2d317b0af14c0c070fd4ef5b357775d
You're receiving this email because of your account on salsa.debian.org.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://alioth-lists.debian.net/pipermail/debian-med-commit/attachments/20260206/c5da700a/attachment-0001.htm>
More information about the debian-med-commit
mailing list