BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,304,547	C→T	4.5%	intergenic (+156/+559)	eco → / ← mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,304,547	0	C	T	4.5%	284.1 / 2.9	105	intergenic (+156/+559)	eco/mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand): ref base C (58/41); new base T (3/2); total (61/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

Only 100 of 106 total aligned reads displayed.

TGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCC  >  NC_000913/2304312‑2304784

tgtATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTTCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCAc                                                                                                                                                                                                                                               >  2:1087829/1‑238 (MQ=255)

tgtATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCAc                                                                                                                                                                                                                                               >  2:1331046/1‑238 (MQ=255)

    taCGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccgc                                                                                                                                                                                                                                          >  2:1394599/1‑239 (MQ=255)

    taCGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCAcgccgc                                                                                                                                                                                                                                          >  2:1598489/1‑239 (MQ=255)

     aCGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCAcgccgc                                                                                                                                                                                                                                          >  1:276163/1‑238 (MQ=255)

       gCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAt                                                                                                                                                                                                                                        >  1:623189/1‑238 (MQ=255)

         cAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCg                                                                                                                                                                                                                                     >  2:286383/1‑239 (MQ=255)

          aGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCg                                                                                                                                                                                                                                     >  2:961712/1‑238 (MQ=255)

          aGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCACCCg                                                                                                                                                                                                                                     >  2:1067930/1‑238 (MQ=255)

            aCAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCcaca                                                                                                                                                                                                                                  >  2:607661/1‑239 (MQ=255)

            aCAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACa                                                                                                                                                                                                                                  >  2:1383611/1‑239 (MQ=255)

            aCAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAAGCGGTCACGCCGCGTCCGACATCCAACGCCCGAGCCGGTTGCCCGA‑TGCGACGCACGCGCGGCTTATCAGGCCCAGACCGCTGACTAGTGCTCCTCCCCGGAGGCCGGACAAGGCCGGCCCGCCGCACCCgc                                                                                                                                                                                                                                    >  1:49820/1‑236 (MQ=18)

                 gTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCCTCCGACATC‑Ca                                                                                                                                                                                                                             >  2:1348138/1‑238 (MQ=255)

                 gTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑Ca                                                                                                                                                                                                                             >  2:1173321/1‑238 (MQ=255)

                 gTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑Ca                                                                                                                                                                                                                             >  2:1100902/1‑238 (MQ=255)

                     aTGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAAccc                                                                                                                                                                                                                         >  2:622399/1‑236 (MQ=255)

                     aTGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGc                                                                                                                                                                                                                         >  2:1249383/1‑238 (MQ=255)

                     aTGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGc                                                                                                                                                                                                                         >  2:1521150/1‑238 (MQ=255)

                          aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAg                                                                                                                                                                                                                    >  2:1118949/1‑238 (MQ=255)

                          aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAg                                                                                                                                                                                                                    >  2:135615/1‑238 (MQ=255)

                          aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAg                                                                                                                                                                                                                    >  1:1575870/1‑238 (MQ=255)

                          aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAg                                                                                                                                                                                                                    >  2:866023/1‑238 (MQ=255)

                          aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAg                                                                                                                                                                                                                    >  2:360470/1‑238 (MQ=255)

                          aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGGCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAg                                                                                                                                                                                                                    >  2:636733/1‑238 (MQ=255)

                              aCCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCACACAAGGCGGGCACGCCGCATCCGACATC‑CAACGCCCGCCCCggg                                                                                                                                                                                                               >  2:1198492/1‑238 (MQ=255)

                              aCCGCGTCAGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACGGCCGTGAAGGGCGGCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCGGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCgg                                                                                                                                                                                                                <  2:132249/238‑1 (MQ=255)

                               ccGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGtt                                                                                                                                                                                                              >  1:67081/1‑239 (MQ=255)

                                    tCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTg                                                                                                                                                                                                         >  2:1013381/1‑239 (MQ=255)

                                       ggAAGGCGGAAGAGAAAATTGACAACGCGTTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATg                                                                                                                                                                                                      <  2:168197/239‑1 (MQ=255)

                                       ggAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGAt                                                                                                                                                                                                       >  2:439111/1‑238 (MQ=255)

                                       ggAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGAt                                                                                                                                                                                                       >  2:1457334/1‑238 (MQ=255)

                                        gAAGGAGGAAAAGAAAATGGAAAACGGGGTAGTGCGCCAAACGGCCGGGAAGTGCGGCCCCCCGTGGGTCAGACAAGGCGGCCACGCCGCATCCGACACCCAACGCCCGAGCCGGTTGCCTGA‑GGGGACGCTGGCGCGGCTTAACAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATg                                                                                                                                                                                                      <  1:439111/238‑1 (MQ=25)

                                             cGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCCGACAAGGCGGTCACGCGGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCCGGCCTACACCGCTGTGCAGTGCGCCACCCCGCCGGTCAGACAAGGCGGTCACGTCGCATCCGCCACC‑CAACGCCCGACCCGGTGGACTGATGCGACGc                                                                                                                                                                                                >  1:430435/1‑239 (MQ=17)

                                                 agagAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGGGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                                                             <  2:1575870/238‑1 (MQ=255)

                                                       aaTTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGCGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTCGGTCACACAAGGCGGTCACGCCGCATCCGCCATC‑CAACGCCCGCGCCGGTTGCCTGATGCGACGCtagacctcca                                                                                                                                                                                      >  1:1417635/1‑230 (MQ=17)

                                                          tGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                   <  2:788407/239‑1 (MQ=255)

                                                          ggACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                   <  2:812874/238‑1 (MQ=255)

                                                              tggcctcgaccgggCTAACCGGCCGCGAAGTGCTCCACCTCGTAGGACGGATGTGGCGGTTACGCCGCATCCGACATCTACCGCCCGAGCGGGTTGCGTGA‑TGCGACGCTGGCGCCTCTGATCAGGCCTACACCGCTGTGAAGTGCTCCCCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAg                                                                                                                                                                                <  1:286383/225‑1 (MQ=11)

                                                               aCGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAgg                                                                                                                                                                               >  2:1272706/1‑238 (MQ=255)

                                                                     tgactCGATAAACGGCAGTTACGTGCGCGCCCCAGTAGGTAAGAAAAGGCGGTCACGCCGCGTCCGACACCCGACGCCCGAGCCGGTGGCCTGA‑TGCGACGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGGGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                         <  1:1100902/234‑1 (MQ=14)

                                                                          cGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCATACAAGGCGGTCACGCCGCATCCGCCATC‑CAACGCCCGAGCCGGTTGCCGGATGCGACGCTAGCGCGACTTATCACGCCTACACCGCtg                                                                                                                                                                   >  2:1722813/1‑239 (MQ=255)

                                                                          cGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAcc                                                                                                                                                                       >  2:1163485/1‑235 (MQ=255)

                                                                          cGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGCTGCCTGATTCGACGCTGGCGCGTCTTATCCGGCCTACAcc                                                                                                                                                                       >  1:1465246/1‑235 (MQ=255)

                                                                              aaagcgccggtaaggccGGCACCCCGTAGGTCGGACAGGGCGGCCACCGCTCAGCCCAGAACCCACCCCCGGCGCGGTTGCGGGA‑TGCGACGCGGGCGCGTCGTAGCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGaa                                                                                                                                                               <  1:1163485/223‑1 (MQ=14)

                                                                              aaaCGGCGGTGAAGGGCGGCACCCCGTGGGCAAGACAGGGCGGGACAGCGGCATCGGCGTTCCAACGCCCGAGCGGGTGGCGTGA‑GGCGACGCTGGCGCGTCTGATCAGGGCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGaa                                                                                                                                                               <  1:866023/239‑1 (MQ=14)

                                                                                     ggTGAAGGGCGGCACCCCGCGGGTCAGACAAGGAGGCCACGCCGCGTCCGCGATCCAACGCCCGAGCCGGTGGCCGGA‑TGCGACGTTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTcc                                                                                                                                                        <  1:1322628/238‑1 (MQ=14)

                                                                                     cGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTc                                                                                                                                                         >  2:1477003/1‑238 (MQ=40)

                                                                                       tGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCa                                                                                                                                                       >  2:630072/1‑238 (MQ=40)

                                                                                          aGTGCGGCACCCCGTGGGGCAGGCAGGGCGGCCACGCCGCAGCCGACAGCCAACGCCCGAGCCGGTTGCCTGA‑TGGGACGCTGGCGCGTCTGATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcccc                                                                                                                                                   <  1:1485908/239‑1 (MQ=14)

                                                                                              cGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCTGAGCCGGTTGCCTGATGCGAACCTGGCGCGTCTTATCAGACCTACACCGCTGTGAAGTGCTCCACCCCGTAg                                                                                                                                               >  2:385171/1‑239 (MQ=255)

                                                                                                     ccGGAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAt                                                                                                                                        <  2:1465246/239‑1 (MQ=255)

                                                                                                      cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCGGACAAGGCGGTTACGCCGCATCCGACATC‑CAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAACCGCCGTGAAGTGCTCCACCCCGTAGGTCGGCTa                                                                                                                                       >  1:592545/1‑239 (MQ=255)

                                                                                                      cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                                                                                                                                       >  2:1479457/1‑239 (MQ=255)

                                                                                                       gTAGGTCGGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGa                                                                                                                                                                                                                                    <  2:869148/146‑1 (MQ=255)

                                                                                                       gTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                      <  2:623189/239‑1 (MQ=255)

                                                                                                             cgggaacgtgggcgaagcagcATCCGAGAACCAACGGTCGTGCC‑GATGGCTGAGTGCGACGCTGGCGCGTCTGATGTGGCCTAAAACGCGGTGAAGTGCGGCACCCGGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGt                                                                                                                                <  1:1578717/220‑1 (MQ=14)

                                                                                                              aGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATACGGCGGtt                                                                                                                               >  2:1081287/1‑239 (MQ=255)

                                                                                                              aGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGtt                                                                                                                               >  2:389767/1‑239 (MQ=255)

                                                                                                              aGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGAGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGAGCGGCACCCCGTAGGGCAGATAAGGCGGTCACCCCGCATCCGACATC‑CAACGCCCGCGGCGGTTGCCCAATGGTACCCTTGGCCGTCTTATCAGGCCTACTCCGCTGTGAACGGCTGCACCCCGTACGTCTGACAAGGCgtg                                                                                                                                >  1:1099212/1‑236 (MQ=11)

                                                                                                                aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGATAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCCCTGTGGAGAGCCCCACCCCGAAGGGCGGATAAGGCGGTTAc                                                                                                                             >  1:1140715/1‑239 (MQ=11)

                                                                                                                 cAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGc                                                                                                                                                                                                                         >  1:869148/1‑146 (MQ=255)

                                                                                                                 cAACGCGGTCACGCGTCATCCAACAGCCATCGCCCGCGCCGATCGCCGGA‑TGGGCCGCTGGGGCGTCTTACGCGGCCTACACCGCTGTGAGGGGCGGCTCCCCGTCGGTCAGACAAGGCGGTCCCGCCGCATCCGTCATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCCCCCCGTAGGTCGGATAAGGCGGTTAcg                                                                                                                            <  1:242772/239‑1 (MQ=14)

                                                                                                                  aaGGCGGTCACGCCGCATCCGACAGCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgc                                                                                                                           <  2:276163/239‑1 (MQ=255)

                                                                                                                      cgcacgcgccgtATGCGACACCCGAGGCCGGAGCCGTTGGCTGGA‑TGCGACGCTGACACGGCTTATGTGGCCTACACCGCGGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATACGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg                                                                                                                         <  1:60340/231‑1 (MQ=14)

                                                                                                                            cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACGCCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGAc                                                                                                                 <  2:1417635/239‑1 (MQ=255)

                                                                                                                               cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCa                                                                                                                                                                                                                                                <  2:557095/110‑1 (MQ=255)

                                                                                                                               cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCa                                                                                                                                                                                                                                                >  1:557095/1‑110 (MQ=255)

                                                                                                                                             aaggCCCGGGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCg                                                                                               <  2:757344/236‑1 (MQ=255)

                                                                                                                                                 cccGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                   <  1:1422265/152‑1 (MQ=255)

                                                                                                                                                     aGCCGGTTGGCTAA‑GGCGACGCGGGCGCGTGTTATAGGGCGTACACCGCTGTGAAGGCCGGCACCCCGTAGGTCGGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCt                                                                                        <  1:795433/238‑1 (MQ=11)

                                                                                                                                                         ggtgggctggatgcaaaccTGGCGCGCCTGATCAGGCCTACCCCGCGGTGAATTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACGTC‑CACCGCCCGAGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACA‑TCCAACGCCCGAGCCGGTTGCCTGATg                                                                                    <  1:1272706/221‑1 (MQ=11)

                                                                                                                                                           ttGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                         >  2:1422265/1‑152 (MQ=255)

                                                                                                                                                             gCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACg                                                                               <  2:430435/239‑1 (MQ=40)

                                                                                                                                                                   aggggACGCGGGCGAGTCGGAGCAGCCCTACCCCGCTGTGAAGGGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTgg                                                                           <  1:688457/234‑1 (MQ=14)

                                                                                                                                                                            tGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                 <  2:67081/239‑1 (MQ=40)

                                                                                                                                                                            aggcagctcgtcgcagtccgaacccGCGGTTAAGTGCGCCACCCCGTCGGCCATTCAAGGCGGTCACGCTGCATCCGACATC‑CAACGCCCGAGCCGTTTGCCTGATGCGACGCTGGCGCGTCTTATCTGGCCTACACCGCTGTGAAGTGCTCCCCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                 <  1:372774/216‑1 (MQ=11)

                                                                                                                                                                             ggCGCGTCTTATCATGCCTACACCGCTGTGAGGTGCGGCACCCCGTAGTACGGACAAGGCGGCCACGCCGCATCCGCCATC‑GAACGCCCGAGCCGGTTGCCTGATGCGACGTTGGCGCGTTTTATCAGGCCTACACGGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                 <  1:487737/238‑1 (MQ=11)

                                                                                                                                                                             ggCGCGCCGTATGAGGCCTACGCCGCGGGGAAGTGCGCCCCCGCGTGGGTCAGACAGGGCGGGCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                 <  1:135615/238‑1 (MQ=11)

                                                                                                                                                                                   tCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCTGCGCAGTGCTCCACCCCGTAGGTCGGATCAGGCGGTTACGCCCCCTCCTACATTCTACCACCC‑CGCCCGTTTCCTGGCGCGGCGCCGTCGCGCCTTCTCAGccc                                                           >  1:1101384/1‑238 (MQ=14)

                                                                                                                                                                                      tATCAGGGCTACACCGCGGGTAAGTGCGGCTCCCCGTAGGGCGGCGAAGGCGGCCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                       <  1:1457334/239‑1 (MQ=11)

                                                                                                                                                                                      gaTCAGGCCGACACCGCTGTGAAGTGCGGCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACAGC‑CAACGCGCGTGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                       <  1:1383611/238‑1 (MQ=255)

                                                                                                                                                                                        tCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                      >  2:158292/1‑238 (MQ=40)

                                                                                                                                                                                                 cacCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCCGATGCGACGCCGGCGCGGCTTATCCGGCCTACACCGCCGTGGAg                                            >  2:1719996/1‑239 (MQ=255)

                                                                                                                                                                                                     ggcgtgAACGGCGGCACCCGGCAGGGTGGAAGAGGCGGGCACGGCGCATCCGCCAGC‑GACCGCCCGAGCCGGTTGCCGGATGCGACGCTGGCGCGGCTTATGAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCt                                        <  1:563882/236‑1 (MQ=11)

                                                                                                                                                                                                      cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTTACGCCGCATCCGACATC‑TAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcccc                                                                                                                                                   >  1:1330936/1‑132 (MQ=12)

                                                                                                                                                                                                      cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTTACGCCGCATCCGACATC‑TAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcccc                                                                                                                                                   <  2:1330936/132‑1 (MQ=12)

                                                                                                                                                                                                      cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                                                                                                                                              <  2:620715/137‑1 (MQ=37)

                                                                                                                                                                                                      cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACC‑TCTAATGCCCAACCCGGTTTCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTggc                                        >  2:1059535/1‑236 (MQ=255)

                                                                                                                                                                                                      cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCCCATCCGCCA‑TCTAACGCCCAAGCCGGGTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCt                                        >  2:36107/1‑238 (MQ=255)

                                                                                                                                                                                                            aaggcagtcccccGTAGAGCAGACAGGGCGGTCACCCCGCATCCGA‑ATCACACCGCCCGGGCGGGTTGCCGGATGGGCCGCTGGCGCGTCTTATCCGGCCTACACCGCTGTGAAGTGCTGCGCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcccc                                 <  1:302073/230‑1 (MQ=11)

                                                                                                                                                                                                                cGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTGCTCCACCCCGTAg                             >  2:1035817/1‑239 (MQ=38)

                                                                                                                                                                                                                 gccACCCCGGAGGACAGACATGGCGGCCAGGCCGCATCCGACACC‑CAACGCCCGAGCCGGTTGCCGGATGCGCCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                            <  1:1545495/237‑1 (MQ=11)

                                                                                                                                                                                                                        cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATCAGGCGGATACGCCCCATCCCCCC‑TCTTCCCCCCCCGGCGGGTGCCTGATGCGCCACCCACTCGTCTTCTCCGGCCTCCACCCCCGTGCAGCGCTCCCCCCCGTAgctcgcct                      >  1:1410192/1‑231 (MQ=11)

                                                                                                                                                                                                                        cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACC‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATa                     >  2:828792/1‑239 (MQ=37)

                                                                                                                                                                                                                        cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCCACC‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCCGGCCTCCACCCCTGTGAAGGTCTCCACCCCGTAGGTTGGAtt                     >  1:1607557/1‑238 (MQ=255)

                                                                                                                                                                                                                        cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACCCCGCATCCGACC‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACACCGTAGGCTGGTTa                     >  2:792546/1‑239 (MQ=255)

                                                                                                                                                                                                                                aGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGAc                                                                                                                 <  2:84514/140‑1 (MQ=38)

                                                                                                                                                                                                                                aGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGAc                                                                                                                 >  1:84514/1‑140 (MQ=38)

                                                                                                                                                                                                                                       gCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcc                                                                                                                                                     >  1:652589/1‑97 (MQ=37)

                                                                                                                                                                                                                                            cACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCCTCCGCCA‑TCTAACCCCCAAGCCGGGTGCCTGGTGCCGCGCCCGCGCGGCTTATCAGGCCTACACCCCTGGGGAGTGCCCCACCCCGAGGGTTGGATAAGgggggtcaccccgatcc  >  1:1621834/1‑225 (MQ=11)

TGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGA‑TGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATC‑CAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACA‑TCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCC  >  NC_000913/2304312‑2304784

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`