BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	1,165,444	0	A	C	26.2%	388.7 / 34.6	169	T71P (ACG→CCG)	yfiF	methyltransferase
Reads supporting (aligned to +/- strand): ref base A (16/108); new base C (44/0); total (60/109)
Fisher's exact test for biased strand distribution p-value = 2.37e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

TCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTCC  >  DH1/1165300‑1165594

tCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑Gtt                                                                                                                                                       <  4:741093/150‑1 (MQ=255)

  tGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑Ggttgc                                                                                                                                                    >  1:410691/1‑149 (MQ=255)

   gATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑Gttt                                                                                                                                                      <  1:203714/148‑1 (MQ=255)

    aTAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCcg                                                                                                                                                  <  5:350047/151‑1 (MQ=255)

       aaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGGGCCGG‑GGTTGCcgccc                                                                                                                                               >  6:765735/1‑149 (MQ=255)

       aaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GGTTCCcgcgg                                                                                                                                               >  3:189449/1‑150 (MQ=255)

       aaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GGTTCCcgcg‑c                                                                                                                                              >  1:574053/1‑151 (MQ=255)

       aaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCcgcg‑c                                                                                                                                              <  8:648367/151‑1 (MQ=255)

         cGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCC‑GGTTGCCgcg‑cgc                                                                                                                                            >  3:724218/1‑151 (MQ=255)

            aCCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGGCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGGCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCgg                                                                                                                                         <  5:804311/151‑1 (MQ=255)

             cccACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAATCAGAGTCAGCCCCGTGACCGCAAGTCGGAAGATTCGCCGT‑GGGGGAAG‑GGTTGCCGCG‑CGCCggg                                                                                                                                        >  6:1065570/1‑151 (MQ=255)

              ccACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACTGTCGCCGTCGCCCCGCACGCGATGACAAACAGAGTCAGCCCCGTGACAGCAAGTGGGAAGATTAGCCGT‑GGGGCACA‑GTTTCCAgcg‑ggcgccgt                                                                                                                                       >  6:894536/1‑146 (MQ=255)

                 caACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGCCAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGGCCC‑GGTTCCCGCG‑CGGCggggg                                                                                                                                      >  8:899742/1‑147 (MQ=255)

                 caACCCGCGTACAGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCCGGACCGCAAGTGCGAAGCTTCGCCGA‑GGGCCCCG‑GGGTCCcgcg‑cgggcggggcc                                                                                                                                    >  6:811621/1‑142 (MQ=255)

                   aCCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGtgatga                                                                                                                                  <  3:84877/151‑1 (MQ=255)

                      cgcgTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GGTTGCCGCC‑CGCCGGGGGGTGAGAc                                                                                                                               >  1:371122/1‑151 (MQ=255)

                      cgcgTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCC‑GGTTGCCGCG‑CGCCggggggtgggcc                                                                                                                               >  5:184808/1‑142 (MQ=255)

                           aCCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGGCCCG‑GGTT‑CCGCG‑CGCCGCGGGGTGAGCCC‑CCgg                                                                                                                         >  7:19068/1‑151 (MQ=255)

                             cGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGGCCG‑GTTTGCCGCG‑CGCCGGGGGATGAGCCG‑CCggc                                                                                                                        >  1:145822/1‑149 (MQ=255)

                             cGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GGTTCCCGCG‑CGCCCGGGGCTGAGACC‑CCggc                                                                                                                        >  7:992768/1‑149 (MQ=255)

                             cGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGaa                                                                                                                       <  3:307861/151‑1 (MQ=255)

                              gggAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGaaa                                                                                                                      <  3:50154/151‑1 (MQ=255)

                                gAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAg                                                                                                                    <  7:407260/151‑1 (MQ=255)

                                gAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGGCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAg                                                                                                                    <  5:603152/151‑1 (MQ=255)

                                  aaGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTGCCCCG‑CGCCGGGTGCTGAGCCG‑CCGGAAAAGGc                                                                                                                  >  2:89416/1‑151 (MQ=255)

                                   aGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGcc                                                                                                                 <  3:1055015/151‑1 (MQ=255)

                                     ggggCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGGCCGCGTTTCCCCCC‑CCCCGGG‑GGTGAGCCG‑CCGCGAAAGGCCgc                                                                                                               >  5:292315/1‑150 (MQ=255)

                                      gggcgggcGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAt                                                                                                              <  4:287785/151‑1 (MQ=255)

                                         cgggagTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GGTTGCCGCG‑CGCCGGGGGGTGAGCCG‑CCGGGAAAgggccggccc                                                                                                           >  6:1021137/1‑143 (MQ=255)

                                          gggcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACg                                                                                                         <  5:711227/151‑1 (MQ=255)

                                          gggcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACg                                                                                                         <  8:734448/151‑1 (MQ=255)

                                              gTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCC‑GGTTACCGCG‑CGCCGGGGGGTGCGGCG‑CCGGCAAAGGccccgcacgggtg                                                                                                      >  1:85761/1‑140 (MQ=255)

                                              gTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTgg                                                                                                     <  5:981326/151‑1 (MQ=255)

                                                ccAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGGGGCCG‑GTTTACCGCG‑CGCCGGGGGGTGAGTCG‑CCGGCAAAGGCCGCG‑CACGGTgggc                                                                                                   >  5:1007453/1‑149 (MQ=255)

                                                ccAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GGTTCCCGCG‑CGCCGGGGGTTGGGCCG‑CCGGCAAAGGCCGCT‑CACCGTGGcc                                                                                                   >  3:578002/1‑150 (MQ=255)

                                                ccAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCa                                                                                                   <  6:118949/151‑1 (MQ=255)

                                                    gAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAg                                                                                               <  1:91533/151‑1 (MQ=255)

                                                       aaTCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAACGGCCGAT‑CACGGTGGCATCAGTgg                                                                                            <  5:46741/151‑1 (MQ=255)

                                                       aaTCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATAAGTgg                                                                                            <  5:1066105/151‑1 (MQ=255)

                                                          ctcGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GGTTGCCGCG‑CGCCGGGGGGTGGGCCG‑CCGGAAAAGGCCGCG‑CACCGTGGGCTCAGTggga                                                                                          >  5:443838/1‑148 (MQ=255)

                                                            cGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTa                                                                                          <  2:1053685/148‑1 (MQ=255)

                                                            cGTGCTGACGCTGGCCGGCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGGCGGGTGATGAGACG‑CCGGAAAAGTCCGAT‑CACGGTGGCATCAGTGGTa                                                                                          <  3:112887/148‑1 (MQ=255)

                                                             gTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCC‑GGTTGCCCCG‑CGCCCGGGGGTGTGACGCCCGGAAAAGG‑CGCT‑CACGGGGGGCGCAGTgggtaa                                                                                        >  1:399197/1‑145 (MQ=255)

                                                             gTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAg                                                                                      <  3:976714/151‑1 (MQ=255)

                                                                  gACGGTGGCCGTCGCCCCGCCGTCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGttt                                                                                   <  3:338450/149‑1 (MQ=255)

                                                                  gACGGTGGCCGTCGCCCCGCCCGCGATGACAAAGAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCTCG‑GTTTCCCGCG‑CGCCGGGAGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGtttt                                                                                  <  6:215520/150‑1 (MQ=255)

                                                                  gACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTa                                                                                 <  4:167079/151‑1 (MQ=255)

                                                                   aCGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTAt                                                                                <  8:80856/151‑1 (MQ=255)

                                                                    gggTGGCCGTCGCCGCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCTGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTAtt                                                                               <  7:679506/150‑1 (MQ=255)

                                                                    cGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTAtt                                                                               <  8:88530/151‑1 (MQ=255)

                                                                      gTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGGGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGa                                                                             <  4:587377/151‑1 (MQ=255)

                                                                       gggCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGAt                                                                            <  1:516044/150‑1 (MQ=255)

                                                                        ggCCGTCGCCCCGCCCGCGATGACAATCAGAGTCAGCGCGGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCGCGCG‑CGCCGGGTGATGAGACG‑CCGGAATAGGGCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGAt                                                                            <  2:69122/150‑1 (MQ=255)

                                                                        ggCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGATG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGAt                                                                            <  8:1049176/150‑1 (MQ=255)

                                                                            gTCGCCCCGCCAGCTATTACAAACAGAGTCAGCCACGTGACAGCAAGTGGGAAGATTCGCCGT‑GGCGGACG‑TTTTCCCGGG‑CGCCGTGTGATGAGACG‑CCGGAAAAGGCCGAT‑AACGGTTGCATCTGTGGTAAAAGTTTTATTGATCCGGa                                                                       <  3:1033775/151‑1 (MQ=255)

                                                                             gcgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGaa                                                                      <  6:224631/150‑1 (MQ=255)

                                                                              cgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGaa                                                                      <  1:491499/150‑1 (MQ=255)

                                                                              cgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAg                                                                     <  6:589020/151‑1 (MQ=255)

                                                                              cgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAg                                                                     <  3:174549/151‑1 (MQ=255)

                                                                              cgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAg                                                                     <  8:302670/151‑1 (MQ=255)

                                                                                ccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GGTTCCCGCG‑CGCCGGGGGATGAGCCG‑CCGGAAAACGGCGCG‑CCCGGTGGGCTCAGTGGGTAAAGTTTTTTTGATACGGCAgtg                                                                   >  5:1082290/1‑151 (MQ=255)

                                                                                ccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAgtg                                                                   <  6:955331/151‑1 (MQ=255)

                                                                                  ccgcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GGTTACCGCC‑CGCCCGGGGATGAGACG‑CCGGAAAAGGGCCAA‑CAACGTGGCATACGTGGTAAAAATTTTATTGGTCCGGCAGTGtt                                                                 >  6:270823/1‑151 (MQ=255)

                                                                                  ccgcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGtt                                                                 <  4:190554/151‑1 (MQ=255)

                                                                                    gcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGc                                                               <  7:859452/151‑1 (MQ=255)

                                                                                    gcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGc                                                               <  4:812729/151‑1 (MQ=255)

                                                                                    gcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGc                                                               <  6:548626/151‑1 (MQ=255)

                                                                                     cccgcGATGACAAACAGAGTCAGCCCCGTGAGCGGATGTGGGAAGATTCGCCGT‑GGCGCACG‑GGTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTTTTGATCCGGAAGTGTTGcg                                                              <  4:37109/151‑1 (MQ=255)

                                                                                         ggATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGcgtcgt                                                          <  8:935213/150‑1 (MQ=255)

                                                                                           aTGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCAGAT‑CACTGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCa                                                        <  7:262603/151‑1 (MQ=255)

                                                                                            tGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACC‑GTTTGCCCCG‑CGCCGGGGGATGAGGAG‑CCGGCAAAGGCCGCT‑GACGGGGGCCTGAGGGGGTAAAGTTTTTTTGATCCCGCAGTGTTgcggcggc                                                         >  3:444400/1‑147 (MQ=255)

                                                                                             gACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAg                                                       <  1:16655/150‑1 (MQ=255)

                                                                                                     gagTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑Ggaa                                             <  5:558715/151‑1 (MQ=255)

                                                                                                       gTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTGCCGCC‑CGCCCGGGGATGAGACC‑CCGGAAAAGGCCGCT‑CACCGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGGAGTGTTTCGTCGTCAGCGTGC‑ggg                                              >  4:553937/1‑147 (MQ=255)

                                                                                                          aGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAAc                                        <  2:339200/151‑1 (MQ=255)

                                                                                                             cccGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCACACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCcg                                     >  5:532722/1‑151 (MQ=255)

                                                                                                               cGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCcgcg                                   <  4:866532/151‑1 (MQ=255)

                                                                                                               cGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCcgcg                                   <  3:100639/151‑1 (MQ=255)

                                                                                                                gTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GGTTGCCCCG‑CGCCGGGGGGTGAGACC‑CCGGCAAAGGGCGCT‑CACCGTGGGCTCAGTGGTTAAAATTTTTTTTATACGGGAGTGTTGGGGGGGCAGCGGGC‑GGGAGACACCcgcgg                                  >  7:538594/1‑150 (MQ=255)

                                                                                                                   aCCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTa                               <  2:1040364/151‑1 (MQ=255)

                                                                                                                   aCCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTa                               <  2:1035448/151‑1 (MQ=255)

                                                                                                                   aCCGCAAGTGGGAAGATTCGCCGT‑GGCCCCCG‑GTTTCCCCCGCCGCCGGG‑GAGTAGGCG‑GCGGGAAAGGGCGATACAC‑GTGGCATGAGTGGGTAAAGGTTTTTTTATACCGACGGGTTGGGTCGTGCGCGTGCGGGACGAAAccccggctc                                >  3:972316/1‑145 (MQ=255)

                                                                                                                    ccGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTAc                              <  2:28570/151‑1 (MQ=255)

                                                                                                                       cAAGTGGGAAGATTCGCCGT‑GGCGCCCG‑GTTTCCCGCG‑CGCCGGGTGATGAGCCG‑CCGGGAAAGGCCGCT‑CACGGTGGCCTGAGTGGGAAAACTTTTATTGATCCGGAAGTGTTGGGTCGTCACCGTGC‑GGGAGAAACCCGCGGCTAAGGc                           >  3:607174/1‑151 (MQ=255)

                                                                                                                        aaGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCg                          <  4:187186/151‑1 (MQ=255)

                                                                                                                          gTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGaa                        <  1:724272/151‑1 (MQ=255)

                                                                                                                            gggAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGaaaa                      <  8:509102/151‑1 (MQ=255)

                                                                                                                            gggAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGaaa                       <  2:675559/150‑1 (MQ=255)

                                                                                                                               aaGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGGGATGAGACG‑CCGGAAAAGGCCGAT‑CAAGGGGGCATCAGTGGTATAACTTTTATTGATACGGAAGTGTTTGGTCGTGAGCGTGG‑GGGAGAATCCCGCGTCTACGGCGCGAATGc                   >  2:558836/1‑151 (MQ=255)

                                                                                                                                 gATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCAt                 <  7:59472/151‑1 (MQ=255)

                                                                                                                                 gATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCAt                 <  6:1033167/151‑1 (MQ=255)

                                                                                                                                  aTTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATg                <  5:270003/151‑1 (MQ=255)

                                                                                                                                     cGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCa             <  7:353040/151‑1 (MQ=255)

                                                                                                                                       ccGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAgg           >  4:529668/1‑151 (MQ=255)

                                                                                                                                         gT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCa         <  7:451499/151‑1 (MQ=255)

                                                                                                                                             gcgcACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGTCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACtt      <  7:553228/151‑1 (MQ=255)

                                                                                                                                                cACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTc   >  8:459018/1‑151 (MQ=255)

                                                                                                                                                cACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTc   >  2:505740/1‑151 (MQ=255)

                                                                                                                                                 aCG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTcc  <  1:118377/151‑1 (MQ=255)

                                                                                                                                                 aCG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTcc  >  5:367448/1‑151 (MQ=255)

TCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGT‑GGCGCACG‑GTTTCCCGCG‑CGCCGGGTGATGAGACG‑CCGGAAAAGGCCGAT‑CACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGC‑GGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTCC  >  DH1/1165300‑1165594

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