BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	4,844,974	T→A	13.5%	intergenic (+50/+69)	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,974	0	T	A	13.5%	241.1 / 6.7	150	intergenic (+50/+69)	apt/recR	Adenine phosphoribosyltransferase/Recombination protein RecR
Reads supporting (aligned to +/- strand): ref base T (24/84); new base A (13/7); total (59/91)
Fisher's exact test for biased strand distribution p-value = 2.79e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

CCGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGC  >  ACB143/4844841‑4845109

 cGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACg                                                                                                                                              <  2:121890/136‑1 (MQ=255)

  ggAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑CCGGCTACg                                                                                                                                              >  2:216047/1‑135 (MQ=255)

   gAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGaa                                                                                                                                            <  7:374152/136‑1 (MQ=255)

     aCTGGAGGGGTGGGGGCGGTTGCACGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAg                                                                                                                                          <  4:329373/136‑1 (MQ=38)

      cTGGACGGGTCGCGGCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑c                                                                                                                                        <  1:52590/136‑1 (MQ=255)

         gACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑Cc‑gc                                                                                                                                    <  8:367774/136‑1 (MQ=255)

           cGGGTCGCGTGGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑Gt                                                                                                                                 <  6:337382/136‑1 (MQ=255)

           cGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑Gt                                                                                                                                 <  6:139926/136‑1 (MQ=255)

                cgcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑tt                                                                                                                           <  2:258034/136‑1 (MQ=255)

                 gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAAAAC‑GCGCCTACGCAAG‑CC‑GC‑GGATG‑ttt                                                                                                                          >  2:61322/1‑136 (MQ=255)

                    tCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑GCGCCAACGCAAG‑CC‑GC‑GGATG‑tttttt                                                                                                                       >  4:304887/1‑134 (MQ=39)

                     cgttggCAAGCGGCCGGGTTGCCTACCTTCTGCCTGACCGAGTTATCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTa                                                                                                                      <  8:303748/131‑1 (MQ=255)

                     cGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTa                                                                                                                      <  5:173232/136‑1 (MQ=255)

                     cGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTa                                                                                                                      <  7:72010/136‑1 (MQ=255)

                     cGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTa                                                                                                                      <  8:188273/136‑1 (MQ=255)

                      ggTTCCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTa                                                                                                                      <  7:93118/135‑1 (MQ=255)

                       gTTGCAAGCGGCCGGGGGGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAgg                                                                                                                    <  6:169223/136‑1 (MQ=255)

                           cAAGCGGCCGGGGTGCCTACCTTATGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑CCAACTACGAAAT‑AC‑CC‑GTATG‑TTTTTTTGGcg‑cg                                                                                                               >  1:299484/1‑136 (MQ=255)

                           aaaGCGGACGGGGTGCCTACATTCTGCCTTACCGAGTTTTCCCTCAGCGCGTATTGATTGGGCATTTGGGGCTGA‑ATTTCAGCACAAT‑‑CG‑CCGGTAAAACC‑CCGCCGAAGAAAA‑CC‑GC‑GGATT‑TTTTTTATgcg‑cc                                                                                                               >  2:263309/2‑135 (MQ=255)

                                ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATTAGCG‑GGGGTAAACCC‑‑CTCGAACAAAAG‑CC‑GC‑GAGTT‑TTTTTTTGGCg‑gggtgt                                                                                                           >  4:15870/1‑133 (MQ=255)

                                ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑GCAACGAACGAAT‑CC‑GC‑GGATG‑GTTTTTAGGCG‑CCGGgtt                                                                                                          >  1:378834/1‑136 (MQ=255)

                                ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑CCGCCGACCCAAG‑CC‑GC‑GGATG‑TTTTTTTGGGG‑CGGGgtt                                                                                                          >  4:141549/1‑136 (MQ=255)

                                ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GGATG‑TTTTGTAGGCG‑CGGGttt                                                                                                          >  5:195903/1‑136 (MQ=255)

                                    ggggTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑GCGCCGACGGAAG‑CC‑GC‑GGTTG‑TGTTGTAGGCG‑CGGGGTTA‑ccc                                                                                                     >  5:266157/1‑136 (MQ=38)

                                                     ccTGCCCGAGTTTTCCCGCACCTAGTATTGAGTGGGCTGTTGGGGCTGC‑TTTGCCGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACCAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAg                                                                                   <  6:229374/136‑1 (MQ=37)

                                                     ccTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑CCACCAGCGTAAG‑CC‑GC‑GGTTG‑TTTTTTATGCC‑CGGGGTTT‑ACCGCC‑CATGGGGCCCCAAg                                                                                   >  8:172932/1‑136 (MQ=255)

                                                      ctcaCCGAGTTTCCCCTCAGCGAGTATTGCGTGGGCTTTTGGGGCTGC‑TTTCCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTCCGAACG‑CC‑GC‑GTATG‑TTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGc                                                                                  <  4:84109/133‑1 (MQ=38)

                                                      cTGACCGAGTTTTCCCTCCCCTATTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACCAACG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGc                                                                                  <  1:44964/136‑1 (MQ=38)

                                                       tGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑CCTCCTACACAAA‑CC‑GC‑CGATG‑TTTTTTATGCG‑CGGGTTTT‑ACCCCC‑AAAGGgccgcccaccg                                                                                 >  6:96373/1‑131 (MQ=25)

                                                       tGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑AG‑CGGGTAAACAC‑GCGACAACCCAAC‑CC‑CC‑GGATT‑TTTTGTAGGCG‑GGGGTTTA‑ACCCCC‑AATGGGCCGCCAAACg                                                                                 >  5:68536/1‑136 (MQ=255)

                                                        gACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑GCTGCTACGAAAG‑CC‑GC‑GGATG‑TTTTGTTTGCG‑CGGGGTTT‑CCCGCG‑CAAGGGCCGCCaacagg                                                                                >  2:109702/1‑132 (MQ=37)

                                                        gACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAAACC‑CCGCCAACCAAAT‑CC‑GC‑CGATG‑TTTTTTAAGCG‑CGGGTTTA‑ACCGCC‑AAATGGCCGCCAACCgg                                                                                >  3:44475/1‑136 (MQ=255)

                                                          ccGAGTTTTCCCTCTCCGAGTATTGAGTGGGCTTGTGGGGCGGC‑TTGGCAGCCAATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAC‑CC‑GC‑GTATG‑TTTTGTAGGCC‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGcc                                                                              <  2:279331/136‑1 (MQ=37)

                                                           cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑CCGCCAACGAAAG‑CC‑GC‑GGATG‑TTTTGTAGGCG‑CGGGGTTA‑CCCGCG‑AATGGGCCGCCCAGCGGccc                                                                             >  7:95833/1‑136 (MQ=38)

                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑CCCCCTAACAAAG‑CC‑CC‑GTATT‑TTTTTTTGGCG‑CGGGTTTA‑CCCCCG‑AATGGGCCGCCAAACCGcccc                                                                            >  6:267263/1‑136 (MQ=37)

                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑CCTCCGAACCAAG‑CC‑CC‑GGTTGTTTTTTTTGCCC‑CGGGTTTA‑CCC‑CG‑AAAGGGGCCCCCACCGGcccc                                                                            >  3:165137/1‑136 (MQ=255)

                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑CCGCCGACGAAAG‑CC‑GC‑GGTTG‑TTTTTTTTGCG‑CGGGGTTT‑ACCGCGCATTGGGCCCCCCCCCGGccc                                                                             >  3:176967/1‑136 (MQ=255)

                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGcccc                                                                            >  2:309734/1‑136 (MQ=255)

                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCACATT‑‑CG‑CGGGTTAACCC‑GCTCCGACGAAAG‑CC‑GC‑GGATG‑TTTTGTAGGCG‑CGGGTTTT‑ACCGCG‑CATGGGGCGCCCAGCCGcccc                                                                            >  5:91497/1‑136 (MQ=255)

                                                                  tCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑CCTCCAACCAAAG‑CC‑CC‑GGATT‑TTTTTTACGCC‑CCGGTTTT‑ACCCCC‑AATGCGGCGCCCAGCCGCCCCAAAGa                                                                       >  2:220827/1‑135 (MQ=19)

                                                                     cTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTTGGCG‑CGGGTTTA‑CCCGCG‑AAAGGGCCGCCCAGCGGCCCCAATGAAt                                                                     >  5:303689/1‑134 (MQ=255)

                                                                          cGATTATTGCGGGGGCTTTTGGGGGTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAACCCC‑GAGCCTCCGAACG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑g                                                             <  8:253223/136‑1 (MQ=37)

                                                                          cGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCCCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑g                                                             <  7:299185/136‑1 (MQ=255)

                                                                           gAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑GCGGCTAAGAAAG‑CC‑GC‑GTATT‑TTTTTTATGCC‑CGGGGTTT‑CCCGCG‑GATTGGCCCCCAAGCGGCCCCCATGAATCAGAGC‑g                                                             >  8:80066/1‑135 (MQ=255)

                                                                           gAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GGATT‑TTTTGTTGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATAAGAGC‑g                                                             >  5:183453/1‑135 (MQ=255)

                                                                           gAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCGACTACGAAAG‑CC‑GC‑GGATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑g                                                             >  3:308942/1‑135 (MQ=255)

                                                                                  gtgggggCTTTTGGGGCTGC‑TTTCCTGCCAATT‑‑CG‑CGGTTAAACCC‑GCTCCTACTAAAG‑CC‑GG‑TTTTG‑TTTTTTAGGCG‑GGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATggg                                                      <  4:50663/131‑1 (MQ=255)

                                                                                    gggggCTTTTGCGCCTGC‑ATTGAAGACACTT‑‑CG‑CTGCTAACCCC‑CCTCCTACCAAAC‑CC‑GC‑TTAGG‑TTTTGTAGCGG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                   <  6:104991/134‑1 (MQ=19)

                                                                                    gTGGGCTTTTGGGGCGGC‑TGTGCCGCCCATT‑‑CG‑CGGGTAACCCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                   <  7:167221/136‑1 (MQ=38)

                                                                                     tGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CGCCGGGTAAACCC‑‑CGCCATCGCAAGCCC‑GG‑GTTTG‑TTTT‑TTGGCGCCGGGTTTG‑CCCCCC‑AAGGGGCC‑CCCACCGGGCCCAATTAATACGCGC‑GCAATTGGCCgg                                                  >  8:156446/1‑136 (MQ=255)

                                                                                     tGGGCTTTTGGGGCTGC‑TTTGCAGCCCATG‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                   <  3:50639/135‑1 (MQ=255)

                                                                                     gtgggttttGGGGGGGCATTTGCCCCCAATC‑‑CG‑C‑GGTAAACCC‑CCCCCAACCCAAC‑CCGGC‑GTTTT‑TTTTGAGGGCG‑GGGTTTAA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                   <  5:32174/131‑1 (MQ=255)

                                                                                     gggggTTTGGGGGCTGG‑TTTCCCGCCCTTG‑‑CG‑CGGGTAACCCC‑GCTCCTACCAACG‑CC‑GC‑TTATT‑TTTTGTAGGGG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCg                                                   <  7:233511/134‑1 (MQ=25)

                                                                                          ttttGGGGCTGC‑TTTCCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACCAAAC‑CC‑GC‑GTTTG‑TTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGccc                                             <  6:174153/136‑1 (MQ=37)

                                                                                           tttGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑GCGGCTTCGAAAG‑CC‑GC‑GGTTG‑TTTTGTAGGCG‑CGGGGTTT‑CCCGCG‑ACTGGGGCGGCCAGCGGGCCCCAAGGTTAAGTGC‑GGAATGGGGCGGCGcc                                              >  7:261184/1‑134 (MQ=16)

                                                                                           tttGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑GCGCCGTCGAAAC‑CC‑CC‑GGTTT‑TTTTGTTGGCG‑CGGGGTTT‑CCCGCG‑CATGGGGCGGCAAGCGGGCCCAATGGATAAGCGC‑GGAATTGGCCGGGGcccc                                            >  6:14700/1‑135 (MQ=255)

                                                                                           tttGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTGCTACGAAAG‑CC‑GC‑GGATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCt                                            >  3:293872/1‑136 (MQ=255)

                                                                                            ttGGGGCTGC‑TTGGCCGCCCATT‑‑CG‑CGGTGAACCCC‑GGGCCTACGAAAG‑CG‑CG‑GTATG‑TTTTGTAGGCG‑CGGTTTTT‑CCCGGG‑AATGGGCCGGCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCt                                            <  8:192966/135‑1 (MQ=255)

                                                                                            ttGGGGCGGC‑TTGGCAGCCCATT‑‑CG‑CGGGAAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCt                                            <  2:2707/135‑1 (MQ=255)

                                                                                             tGGGGCGGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GAGCCAACGAACG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTg                                           <  7:372988/135‑1 (MQ=38)

                                                                                                ggCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTTAAACC‑CCGCCGACGGAAG‑GC‑GCGGTTTG‑TTTTTTTGGGC‑GGGGTTTA‑CCGCGG‑ATTGGGCCCCCACCCGG‑CCCACAGAATGAGCGCGGAAATGGGCC‑GCGccccggaa                                        >  4:270630/1‑130 (MQ=255)

                                                                                                ggCTGC‑GTTCCCCCCCTTT‑‑CG‑GGGGAAACCCC‑GCTCCTACCAAGC‑CC‑GC‑GTATT‑TTTTGTCGCCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCaaa                                       <  2:158307/136‑1 (MQ=25)

                                                                                                 ggggc‑tTTGCAGCCCATT‑‑CG‑CGGGAAAACCC‑GCTCATACGAAAG‑CC‑GC‑GTTTG‑TTTTGAAGGCG‑CGGTTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  4:353100/133‑1 (MQ=255)

                                                                                                 gctgc‑tTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  1:309734/136‑1 (MQ=255)

                                                                                                 gctgc‑tTTCCAGCCCTTT‑‑CG‑CGGGTAACCCC‑GCTCCTCCGAAAG‑CC‑GC‑TTTTG‑TTTTGTCGGGG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  8:36434/136‑1 (MQ=38)

                                                                                                 gcggc‑tTTGCAGCCCATT‑‑CG‑CGGCTAACCCC‑GATCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  1:101953/133‑1 (MQ=255)

                                                                                                 gcggc‑tTTGCAGCCAATG‑‑GG‑CGGGAAAACCG‑GCTCCTACGAACG‑CG‑GC‑GTATG‑TTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  5:229131/133‑1 (MQ=38)

                                                                                                  gcggatttGCAGCCCATT‑‑CG‑CGGTAAACCCC‑GCTCCTCCGAACG‑CC‑GC‑GTATG‑TTTTGTCGGCG‑GGGTTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAg                                      <  7:154047/131‑1 (MQ=255)

                                                                                                     c‑tTTGCAGCCCATT‑‑CG‑CGGGTTAACCC‑GCGGCGACGGAAG‑CC‑GG‑GGTTG‑TTTTGTTGGCG‑CGGGGTTT‑CCCGCG‑AATGGGGCGGCCCGCGGGCCCCATGGATGAGGGG‑GGAATGGGGCCGCGCCCGGGAAAGGGAt                                  >  8:263892/1‑136 (MQ=255)

                                                                                                       tttGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCGGCTACGAAAG‑GC‑GC‑GGTTG‑GTTTGTTGGCG‑CGGGGTTT‑CCCGCG‑AATGGGGCGCCCAGCGGGCCCCATGGATCAGAGC‑GCAATGGGCCGGCGCCCGGGCAAGGGATg                                 >  6:374516/1‑136 (MQ=255)

                                                                                                       tttGCAGCCCATT‑‑CG‑CGGGAAACCCC‑GCTCCAACGAACG‑CC‑GC‑GTATG‑TTTTGTCGGCG‑CGGGTTTC‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  8:121688/136‑1 (MQ=255)

                                                                                                       tttGCAGCCCATT‑‑CG‑CGGGAAAACCC‑GCTCCTACGAACG‑CC‑GC‑GTAGG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  7:314865/136‑1 (MQ=255)

                                                                                                       tttCCAGCCCATT‑‑GG‑CGGGTAAACCC‑GCTCCTACTAACG‑CC‑GC‑GTATG‑TTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  5:386537/136‑1 (MQ=39)

                                                                                                       ttggCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCGCCTACTAAAG‑CC‑GC‑GTATG‑TTTTTTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCAGCAAAGGCATg                                 <  4:332060/133‑1 (MQ=38)

                                                                                                       gtcgctgCCCTTT‑‑GG‑CGGGTACACCC‑GCTCCAACGAAAG‑CC‑GC‑GTATT‑TTTTGTCGGCG‑GGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  4:210564/130‑1 (MQ=255)

                                                                                                       gcggCACCCCATG‑‑CG‑CGGTACAACCC‑GCGCATCCCAATG‑CC‑CC‑GTATG‑TTGTGACGGCG‑CGGTTATA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  8:333960/133‑1 (MQ=255)

                                                                                                        ttggcccccATT‑‑GG‑GGGGTAAACCCGGCTCCACCGAACG‑CC‑GC‑TTTTG‑TTTTGTAGGGG‑GGGTTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  3:197687/130‑1 (MQ=255)

                                                                                                        ttGCAGCCCATT‑‑CG‑CGGGGTAAACC‑GCGGCGTCGGAAG‑GC‑GG‑GGTTG‑TTTTGTTGGCG‑GGGGGTTT‑CCCGCG‑GATGGGGCGCCCAGCGGGCCCCATGGATGAGAGC‑GCAATGGGGCGGCGCCCGGGAAAgggaggg                                >  6:193721/1‑131 (MQ=255)

                                                                                                        atcccgCCCTTT‑‑GG‑CGGGTAACCCG‑GCCCATCCTAACG‑CC‑GC‑GTTTG‑TTTGGTGGGCG‑‑GGGTTTACCCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  8:12438/130‑1 (MQ=255)

                                                                                                         tGCAGCCCATT‑‑CG‑CGGGTAACCCC‑GCTCCTACTAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  8:232609/134‑1 (MQ=255)

                                                                                                         tGCAGCCCATT‑‑CG‑CGGGTAACCCC‑GCGCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  1:136977/134‑1 (MQ=255)

                                                                                                         tGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGAAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  7:18316/134‑1 (MQ=255)

                                                                                                         tGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCGCCTACGAAAG‑CC‑GC‑GTTTG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  2:8185/134‑1 (MQ=255)

                                                                                                         tGCAGCCAATG‑‑CG‑GGGGAACACCC‑GCTCCTCCTAAAG‑CC‑GG‑TTTTG‑TTTTTTAGCCG‑CGGGTTAA‑CCCGCT‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  4:343551/134‑1 (MQ=255)

                                                                                                         ggCAGCCCATT‑‑CG‑CGGGAACACCC‑GCGCCTACGAAAC‑CC‑GC‑GTATG‑TTTTGTCGGCG‑GGGTTTAA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATg                                 <  1:384515/133‑1 (MQ=255)

                                                                                                          gCAGCCCATT‑‑CG‑CGGGAAAACCC‑GCTCCTACGAACG‑CC‑GC‑GTATG‑TTTTGTCGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATggg                               <  1:49092/135‑1 (MQ=255)

                                                                                                            aGCCCATT‑‑CG‑CGGGTAAACCG‑GCTCCTACGAACG‑CC‑GC‑GTATG‑TTTTGTCGGCG‑CGGGTTAA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCa                            <  2:313617/136‑1 (MQ=255)

                                                                                                            aGCCCATG‑‑GG‑GGGGTAACCGC‑GCCCCTACTAAAG‑CC‑GC‑GTTTG‑TTTTGTACGCG‑CGGTTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATTGGCCGGCGCCCTGCAAAGGCATGGGCCa                            <  2:115279/136‑1 (MQ=255)

                                                                                                                      tcccgGGTACACCG‑GCTCCACCGACCG‑CC‑GG‑GTTTG‑TGTGGTCGGCG‑CGGGTATC‑CCGGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTa                       <  4:102381/131‑1 (MQ=255)

                                                                                                                      cg‑cgGGTAACCCC‑CCCCCTACTAAAG‑CC‑GC‑GTATG‑TTTGGTCGGCG‑CGGGTTAA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccg                    <  3:229514/136‑1 (MQ=255)

                                                                                                                          gggTAAACCC‑GCTCCTACTAACG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccgccg                 <  7:267054/136‑1 (MQ=255)

                                                                                                                           gttAAACCC‑GCTCCTACGAACG‑CC‑GC‑GTATG‑TTTTGAAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGt                <  1:220502/134‑1 (MQ=255)

                                                                                                                           ggTAACCCC‑GCTCCTAGGAATG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑GGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGt                <  5:365665/136‑1 (MQ=255)

                                                                                                                           ggTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGt                <  2:250873/136‑1 (MQ=255)

                                                                                                                             tAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcc              <  7:220406/136‑1 (MQ=255)

                                                                                                                             tAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcc              <  6:357226/136‑1 (MQ=255)

                                                                                                                              aaaCCC‑GCGCCTGCGAAAG‑GC‑GC‑GGATG‑TTTTGTAGGCG‑CGGGGTTT‑GCCGCG‑CATGGGGCGCCAAGCGGCCCCAATGAATCAGAGG‑GAAATGGGCCGGCGCCCTGCAAAGGCAGGGGCCCGGGTGCCGCCCTccc             >  2:68144/1‑135 (MQ=255)

                                                                                                                                          tACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATTAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGc  >  2:65320/1‑136 (MQ=255)

                                                                                                                                          tACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGc  >  1:246400/1‑136 (MQ=255)

                                                                                                                                          tACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CCGCTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCCATGAATAAGCGC‑GCAAAGGGCCGGCGCCCGGCAAAGGCATGGGCCCGGGTACCCACGTCCACCCTCAACAGc  >  4:112901/1‑136 (MQ=255)

CCGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGC‑TTTGCAGCCCATT‑‑CG‑CGGGTAAACCC‑GCTCCTACGAAAG‑CC‑GC‑GTATG‑TTTTGTAGGCG‑CGGGTTTA‑CCCGCG‑AATGGGCCGCCAAGCGGCCCCAATGAATCAGAGC‑GAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGC  >  ACB143/4844841‑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: `‑`