BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	4,844,977	G→C	17.1%	intergenic (+53/+66)	apt → / ← recR	Adenine phosphoribosyltransferase/Recombination protein RecR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	4,844,977	0	G	C	17.1%	262.4 / 16.4	151	intergenic (+53/+66)	apt/recR	Adenine phosphoribosyltransferase/Recombination protein RecR
Reads supporting (aligned to +/- strand): ref base G (45/67); new base C (13/12); total (61/90)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

CGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGC  >  ACB143/4844842‑4845109

cGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACg                                                                                                                                            <  2:121890/136‑1 (MQ=255)

    aCTGGAGGGGTGGGGGCGGTTGCACGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAg                                                                                                                                        <  4:329373/136‑1 (MQ=38)

     cTGGACGGGTCGCGGCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑c                                                                                                                                      <  1:52590/136‑1 (MQ=255)

        gACGGGTCGCGTCGGTTGGAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑Ccgc                                                                                                                                   <  2:357320/136‑1 (MQ=255)

        gACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑Ccgc                                                                                                                                   <  8:367774/136‑1 (MQ=255)

          cGGGTCGCGTGGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑Gt                                                                                                                                <  6:337382/136‑1 (MQ=255)

          cGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑Gt                                                                                                                                <  6:139926/136‑1 (MQ=255)

            ggTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAAACCGCGCCGACGGAAG‑CCGC‑GTAt                                                                                                                              >  7:253223/1‑136 (MQ=255)

               cgcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGtt                                                                                                                           <  2:258034/136‑1 (MQ=255)

                gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTTGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGGTTAAACCCCTACATAGAAAAG‑CCGG‑GTGTGtt                                                                                                                           >  3:214100/1‑136 (MQ=255)

                gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAAACCGCTCCAACGCAAG‑CCGC‑GGTTGttt                                                                                                                          >  2:81295/1‑136 (MQ=39)

                gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAAAACGCGCCTACGCAAG‑CCGC‑GGATGttt                                                                                                                          >  2:61322/1‑136 (MQ=255)

                gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCAGCGAAGAAAG‑CCGC‑CTATGttt                                                                                                                          >  6:156764/1‑136 (MQ=39)

                   tCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAACCCGCGCCAACGCAAG‑CCGC‑GGATGtttttt                                                                                                                       >  4:304887/1‑134 (MQ=39)

                    cgttggCAAGCGGCCGGGTTGCCTACCTTCTGCCTGACCGAGTTATCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTa                                                                                                                      <  8:303748/131‑1 (MQ=255)

                    cGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTa                                                                                                                      <  8:188273/136‑1 (MQ=255)

                    cGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTa                                                                                                                      <  7:72010/136‑1 (MQ=255)

                    cGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTa                                                                                                                      <  5:173232/136‑1 (MQ=255)

                     ggTTCCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTa                                                                                                                      <  7:93118/135‑1 (MQ=255)

                      gTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAgg                                                                                                                    <  4:138655/136‑1 (MQ=255)

                          cAAGCGGCCGGGGTGCCTACCTTATGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCCCAACTACGAAAT‑ACCC‑GTATGTTTTTTTGGcg‑cg                                                                                                               >  1:299484/1‑136 (MQ=255)

                               ggCCGGGGTGCCTACCTTCTGCCTGCCCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAACCCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGttt                                                                                                          <  3:54352/136‑1 (MQ=255)

                               ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCTGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGGTTAAACCCCCGCACAGGAAAC‑CCCC‑GTTTTTTTTTTAGGGG‑GGGGtt                                                                                                           >  6:135253/1‑136 (MQ=255)

                               ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATTAGCG‑GGGG‑TAAACCC‑CTCGAACAAAAG‑CCGC‑GAGTTTTTTTTTGGCg‑gggtgt                                                                                                           >  4:15870/1‑133 (MQ=255)

                               ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAAACCGCAACGAACGAAT‑CCGC‑GGATGGTTTTTAGGCG‑CCGGgtt                                                                                                          >  1:378834/1‑136 (MQ=255)

                                                  tGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAAACCCCACCAACCAAAA‑CCGCGGTTTGTTTTTTGGGGC‑GGGGTTTA‑CCCCGC‑AATGGgccccc                                                                                      >  7:168780/1‑136 (MQ=255)

                                                    ccTGCCCGAGTTTTCCCGCACCTAGTATTGAGTGGGCTGTTGGGGC‑TGCTTTGCCGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACCAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAg                                                                                   <  6:229374/136‑1 (MQ=37)

                                                     cTGACCGAGTTTTCCCTCCCCTATTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACCAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGc                                                                                  <  1:44964/136‑1 (MQ=38)

                                                      tGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑AG‑CGGG‑TAAACACGCGACAACCCAAC‑CCCC‑GGATTTTTTGTAGGCG‑GGGGTTTA‑ACCCCC‑AATGGGCCGCCAAACg                                                                                 >  5:68536/1‑136 (MQ=255)

                                                          cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCCCGCCAACGAAAG‑CCGC‑GGATGTTTTGTAGGCG‑CGGGGTTA‑CCCGCG‑AATGGGCCGCCCAGCGGccc                                                                             >  7:95833/1‑136 (MQ=38)

                                                           gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAACCCGCTCCGACGCAAC‑CCGC‑GGATGGTTTTTAGGCG‑CGGGGTTT‑CCCGCG‑AATGGGGCCCCCAGCGGcccc                                                                            >  5:392701/1‑136 (MQ=37)

                                                           gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAACCCGCGCCAAACCAAG‑CCCC‑CGATGTTTTTTTTGCG‑CGGGTTTT‑CCCCCC‑AAATGGCCCCCCACCGGgccc                                                                            >  5:302619/1‑136 (MQ=16)

                                                           gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAACCCCCTCCGACGAAAG‑CCGC‑GTTTGTTTTTTAAGCG‑CGGGGTTA‑ACCGCG‑AAAGGGCCGCCCAGCGGccc                                                                             >  1:48927/1‑135 (MQ=37)

                                                           gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAACCCCCGCCAACGAAAA‑CCGC‑GGATTTTTTTTATGCG‑CGGGGTTT‑ACCGCG‑AATTGGCCGCCCAGCGGgcc                                                                             >  5:363973/1‑135 (MQ=255)

                                                           gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAAACCCCGCCGACGAAAG‑CCGC‑GGTTGTTTTTTTTGCG‑CGGGGTTT‑ACCGCGCATTGGGCCCCCCCCCGGccc                                                                             >  3:176967/1‑136 (MQ=255)

                                                           gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGcccc                                                                            >  2:309734/1‑136 (MQ=255)

                                                           gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAAACCGCTCCAACCAAAC‑CCGC‑GGATGTTTTTTAGGCG‑CGGGTTTA‑CCCCCG‑AATGGGCCGCCCACCCGcccc                                                                            >  8:188238/1‑136 (MQ=37)

                                                           gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCACATT‑‑CG‑CGGG‑TTAACCCGCTCCGACGAAAG‑CCGC‑GGATGTTTTGTAGGCG‑CGGGTTTT‑ACCGCG‑CATGGGGCGCCCAGCCGcccc                                                                            >  5:91497/1‑136 (MQ=255)

                                                                ttCCCTCAGCGAGTATGGAGTGGGCTTTTGGGGC‑GGCTTTGCAGCCCATT‑‑CG‑CGGG‑AAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTGGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATg                                                                        <  5:103297/135‑1 (MQ=38)

                                                                 tCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCCCTCCAACCAAAG‑CCCC‑GGATTTTTTTTACGCC‑CCGGTTTT‑ACCCCC‑AATGCGGCGCCCAGCCGCCCCAAAGa                                                                       >  2:220827/1‑135 (MQ=19)

                                                                    cTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTTGGCG‑CGGGTTTA‑CCCGCG‑AAAGGGCCGCCCAGCGGCCCCAATGAAt                                                                     >  5:303689/1‑134 (MQ=255)

                                                                         gtaGTTTTGAGTGGGCTTTGGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAACCCCGCTCCTACGAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑g                                                             <  6:264861/134‑1 (MQ=38)

                                                                         cGATTATTGCGGGGGCTTTTGGGGG‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAACCCCGAGCCTCCGAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑g                                                             <  8:253223/136‑1 (MQ=37)

                                                                         cGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCCCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑g                                                             <  7:299185/136‑1 (MQ=255)

                                                                          gAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAAACCGCGGCTAAGAAAG‑CCGC‑GTATTTTTTTTATGCC‑CGGGGTTT‑CCCGCG‑GATTGGCCCCCAAGCGGCCCCCATGAATCAGAGC‑g                                                             >  8:80066/1‑135 (MQ=255)

                                                                          gAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAAACCCCTGCGTCGCAAG‑CCGC‑GGTTGTTTTGTAGGCG‑CGGGGTTT‑CCCGCG‑AATGGGCCCCCCACCGGCCCCAATGAATAAGcgc‑g                                                             >  7:35531/1‑135 (MQ=255)

                                                                          gAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GGATTTTTTGTTGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATAAGAGC‑g                                                             >  5:183453/1‑135 (MQ=255)

                                                                          gAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTAACTAAG‑GCGC‑GTATGTTTTGTAGGAG‑CCGGTTTA‑ACCGCC‑AATTGGGCGCCAAACCGCCCCAAAGAATCAGAGC‑g                                                             >  7:378865/1‑135 (MQ=25)

                                                                          gAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCGACTACGAAAG‑CCGC‑GGATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑g                                                             >  3:308942/1‑135 (MQ=255)

                                                                              aTTGAGTGCGGTTTTGGGGC‑GTCTTCCCA‑CCAATTG‑CC‑CGGT‑AAAACCCGCTCAAACTAAAG‑CGGC‑TTTTGTTTTGTCGCCG‑GGGTTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATg                                                        <  4:395567/136‑1 (MQ=255)

                                                                                 gtgggggCTTTTGGGGC‑TGCTTTCCTGCCAATT‑‑CG‑CGGT‑TAAACCCGCTCCTACTAAAG‑CCGG‑TTTTGTTTTTTAGGCG‑GGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATggg                                                      <  4:50663/131‑1 (MQ=255)

                                                                                    tGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CGCCGGG‑TAAACCC‑CGCCATCGCAAGCCCGG‑GTTTGTTTT‑TTGGCGCCGGGTTTG‑CCCCCC‑AAGGGGCC‑CCCACCGGGCCCAATTAATACGCGC‑GCAATTGGCCgg                                                  >  8:156446/1‑136 (MQ=255)

                                                                                    tGGGCTTTTGGGGC‑TGCTTTGCAGCCCATG‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                   <  3:50639/135‑1 (MQ=255)

                                                                                    ggggCTTTGGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                   <  2:263323/134‑1 (MQ=255)

                                                                                     tgggttttGGGGCGTTCTTTGCTGCCCTTC‑‑CG‑CGGG‑TAACCCCGCTCTTCCGAACG‑CCGC‑TTTTGTTTTGTAGGGG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCgg                                                  <  8:227611/132‑1 (MQ=25)

                                                                                      ggCTTTTGGGGC‑TGCTTTCCAGCCAATT‑‑CG‑CGGG‑AAACCCCGCTCCTACCAACC‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGgc                                                 <  4:340065/135‑1 (MQ=38)

                                                                                         ttttGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCGCCTACGAAAG‑CCGC‑GGATGTTTTTTAGGCG‑CGGGGTTA‑CCCGCG‑AATGGGCCGCCCAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGccc                                             >  7:105491/1‑136 (MQ=38)

                                                                                         ttttGGGGC‑TGCTTTCCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACCAAAC‑CCGC‑GTTTGTTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGccc                                             <  6:174153/136‑1 (MQ=37)

                                                                                          tttGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTTACCCCCGGCGTCGAGAG‑CCGG‑GGTTGTTTTGTTGGGG‑CGGGGTTT‑GCCGCG‑AATTGGCCCCCCACCGGGCCCCAAGAAAGAGAGC‑CCAAAGGGGCCGCGcccg                                            >  6:325399/1‑135 (MQ=255)

                                                                                          tttGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAACCCGCGGCTTCGAAAG‑CCGC‑GGTTGTTTTGTAGGCG‑CGGGGTTT‑CCCGCG‑ACTGGGGCGGCCAGCGGGCCCCAAGGTTAAGTGC‑GGAATGGGGCGGCGcc                                              >  7:261184/1‑134 (MQ=16)

                                                                                           ttGGGGC‑GGCTTGGCAGCCCATT‑‑CG‑CGGG‑AAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCt                                            <  2:2707/135‑1 (MQ=255)

                                                                                            tGGGGC‑GGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGAGCCAACGAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTg                                           <  7:372988/135‑1 (MQ=38)

                                                                                               ggC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAACCCGCTGCTACGAAAG‑CCGC‑GGTTGTTTTGTAGGCG‑CGGGGTTA‑CCCGCG‑CATGGGCCGCCCAGCGGCCCCAATGGATCAGAGC‑GATATGGGCCGGCGCCCGGCa                                         >  3:304530/1‑134 (MQ=37)

                                                                                               ggC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TTAAACCCCGCCGACGGAAG‑GCGCGGTTTGTTTTTTTGGGC‑GGGGTTTA‑CCGCGG‑ATTGGGCCCCCACCCGG‑CCCACAGAATGAGCGCGGAAATGGGCC‑GCGccccggaa                                        >  4:270630/1‑130 (MQ=255)

                                                                                               ggC‑TGCGTTCCCCCCCTTT‑‑CG‑GGGG‑AAACCCCGCTCCTACCAAGC‑CCGC‑GTATTTTTTGTCGCCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCaaa                                       <  2:158307/136‑1 (MQ=25)

                                                                                                gc‑tgctTTGCAGCCCTTT‑‑CG‑CGGG‑TAACCCCGCTCCTACGAAAC‑CCGC‑GTATGTTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  6:370130/136‑1 (MQ=255)

                                                                                                gc‑tgctTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  1:309734/136‑1 (MQ=255)

                                                                                                gc‑tgctTTGCAGCCCATT‑‑CG‑CGGG‑AAAACCCGCGCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGTTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  7:279688/136‑1 (MQ=255)

                                                                                                gc‑tgctTTCCAGCCCTTT‑‑CG‑CGGG‑TAACCCCGCTCCTCCGAAAG‑CCGC‑TTTTGTTTTGTCGGGG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  8:36434/136‑1 (MQ=38)

                                                                                                 gcggctTTGCAGCCCATT‑‑CG‑CGGC‑TAACCCCGATCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  1:101953/133‑1 (MQ=255)

                                                                                                 gcggctTTGCAGCCAATG‑‑GG‑CGGG‑AAAACCGGCTCCTACGAACG‑CGGC‑GTATGTTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  5:229131/133‑1 (MQ=38)

                                                                                                 gcggatttGCAGCCCATT‑‑CG‑CGGT‑AAACCCCGCTCCTCCGAACG‑CCGC‑GTATGTTTTGTCGGCG‑GGGTTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  7:154047/131‑1 (MQ=255)

                                                                                                 c‑tgctTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACTAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAgg                                     <  1:55620/136‑1 (MQ=255)

                                                                                                     ctTTGCAGCCCATT‑‑CG‑CGGG‑TTAACCCGCGGCGACGGAAG‑CCGG‑GGTTGTTTTGTTGGCG‑CGGGGTTT‑CCCGCG‑AATGGGGCGGCCCGCGGGCCCCATGGATGAGGGG‑GGAATGGGGCCGCGCCCGGGAAAGGGAt                                  >  8:263892/1‑136 (MQ=255)

                                                                                                      tttGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCGGCTACGAAAG‑GCGC‑GGTTGGTTTGTTGGCG‑CGGGGTTT‑CCCGCG‑AATGGGGCGCCCAGCGGGCCCCATGGATCAGAGC‑GCAATGGGCCGGCGCCCGGGCAAGGGATg                                 >  6:374516/1‑136 (MQ=255)

                                                                                                      tttGCAGCCCATT‑‑CG‑CGGG‑AAACCCCGCTCCAACGAACG‑CCGC‑GTATGTTTTGTCGGCG‑CGGGTTTC‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  8:121688/136‑1 (MQ=255)

                                                                                                      gtcgctgCCCTTT‑‑GG‑CGGG‑TACACCCGCTCCAACGAAAG‑CCGC‑GTATTTTTTGTCGGCG‑GGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  4:210564/130‑1 (MQ=255)

                                                                                                      gcggCACCCCATG‑‑CG‑CGGT‑ACAACCCGCGCATCCCAATG‑CCCC‑GTATGTTGTGACGGCG‑CGGTTATA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  8:333960/133‑1 (MQ=255)

                                                                                                       atcccgCCCTTT‑‑GG‑CGGG‑TAACCCGGCCCATCCTAACG‑CCGC‑GTTTGTTTGGTGGGCG‑‑GGGTTTACCCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  8:12438/130‑1 (MQ=255)

                                                                                                        tGCAGCCCATT‑‑CG‑CGGG‑TAACCCCGCTCCTACTAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  8:232609/134‑1 (MQ=255)

                                                                                                        tGCAGCCCATT‑‑CG‑CGGG‑TAACCCCGCGCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  1:136977/134‑1 (MQ=255)

                                                                                                        tGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGAAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  7:18316/134‑1 (MQ=255)

                                                                                                        tGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCGCCTACGAAAG‑CCGC‑GTTTGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  2:8185/134‑1 (MQ=255)

                                                                                                        tGCAGCCCATT‑‑CG‑CGGG‑AAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATTGGCCGGCGCCCTGCAAAGGCATg                                 <  2:256620/134‑1 (MQ=255)

                                                                                                        tGCAGCCAATG‑‑CG‑GGGG‑AACACCCGCTCCTCCTAAAG‑CCGG‑TTTTGTTTTTTAGCCG‑CGGGTTAA‑CCCGCT‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  4:343551/134‑1 (MQ=255)

                                                                                                        tGCACCCCATT‑‑CG‑CGGG‑CAAACCCGCTCCTACCAAAG‑CCGC‑TTATGTTTTGACGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  5:206888/134‑1 (MQ=39)

                                                                                                        ggCAGCCCATT‑‑CG‑CGGG‑AACACCCGCGCCTACGAAAC‑CCGC‑GTATGTTTTGTCGGCG‑GGGTTTAA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  1:384515/133‑1 (MQ=255)

                                                                                                           aGCCCATG‑‑GG‑GGGG‑TAACCGCGCCCCTACTAAAG‑CCGC‑GTTTGTTTTGTACGCG‑CGGTTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATTGGCCGGCGCCCTGCAAAGGCATGGGCCa                            <  2:115279/136‑1 (MQ=255)

                                                                                                                  t‑‑CG‑GGGC‑TAAACCCGCTCCTACGAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAc                      <  6:178916/135‑1 (MQ=255)

                                                                                                                     cg‑cgGG‑TAACCCCCCCCCTACTAAAG‑CCGC‑GTATGTTTGGTCGGCG‑CGGGTTAA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccg                    <  3:229514/136‑1 (MQ=255)

                                                                                                                      g‑cgGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccgc                   <  2:287538/136‑1 (MQ=255)

                                                                                                                          gg‑TAACCCCGCTCCTAGGAATG‑CCGC‑GTATGTTTTGTAGGCG‑GGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGt                <  5:365665/136‑1 (MQ=255)

                                                                                                                           gttAAACCCGCTCCTACGAACG‑CCGC‑GTATGTTTTGAAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGt                <  1:220502/134‑1 (MQ=255)

                                                                                                                             tAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcc              <  7:220406/136‑1 (MQ=255)

                                                                                                                                    gcgcgtACGAACG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAtct       <  2:23337/131‑1 (MQ=255)

                                                                                                                                         tACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATTAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGc  >  2:65320/1‑136 (MQ=255)

                                                                                                                                         tACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGc  >  1:246400/1‑136 (MQ=255)

                                                                                                                                         tACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CCGCTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCCATGAATAAGCGC‑GCAAAGGGCCGGCGCCCGGCAAAGGCATGGGCCCGGGTACCCACGTCCACCCTCAACAGc  >  4:112901/1‑136 (MQ=255)

                                                                                                                                         tACGAAAG‑CCGC‑GTATGTTTTGTAGGAG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCCCCATCTCCAGc  >  3:170879/1‑136 (MQ=255)

                                                                                                                                           cGAACG‑CCGC‑GTAGTTTTTGTAGGCG‑GGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGc  <  5:95840/134‑1 (MQ=255)

CGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGC‑TGCTTTGCAGCCCATT‑‑CG‑CGGG‑TAAACCCGCTCCTACGAAAG‑CCGC‑GTATGTTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGC  >  ACB143/4844842‑4845109

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`