BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	613,601	0	G	T	24.1%	288.1 / 16.9	194	intergenic (+462/+561)	garK/yhaC	glycerate 2‑kinase/pentapeptide repeats‑related protein
Reads supporting (aligned to +/- strand): ref base G (51/78); new base T (7/34); total (61/133)
Fisher's exact test for biased strand distribution p-value = 8.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

GTACTGAACCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACC  >  DH1/613451‑613741

  aCTGAACCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAGCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGat                                                                                                                                            <  5:1531549/151‑1 (MQ=255)

   cTGAACCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACACGGATTATCGGTCAGTATCGCAAGATTTAAGTAAAAACCAG‑CTTAGGCGGGTTTTTCGTTTTGGAGGGGCAGAAAGatt                                                                                                                                           >  5:1198989/1‑150 (MQ=255)

     gAACCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaa                                                                                                                                         <  5:765038/151‑1 (MQ=255)

     gAACCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatg                                                                                                                                           <  3:1050381/149‑1 (MQ=255)

     gAACCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTCCTTTTGGAGGTGCAGAAAGatca                                                                                                                                          >  4:519546/1‑148 (MQ=255)

     gAACCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACACGCCTTATCGATCAGTTTCACCTGATTTACGTAAAAACCAG‑CTGCGGCGGGTTTTTGGTTTTGGAGGGGCAGAAAGatga                                                                                                                                          >  1:1492815/1‑150 (MQ=37)

     gAACCCGGG‑GAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatga                                                                                                                                          <  4:1369330/150‑1 (MQ=255)

      gaCCCGGGTTATGCTTCTTGAGCCAGTGAGCGATTGCTGGCCTAGATCACAGACTCTCCACCACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGatg                                                                                                                                           <  3:46153/148‑1 (MQ=18)

      aaCCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGGTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                                                                                                                                        <  7:696403/151‑1 (MQ=255)

      aaCCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                                                                                                                                        <  2:511011/151‑1 (MQ=255)

      aaCCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                                                                                                                                        <  1:868804/151‑1 (MQ=255)

      aaCCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTCCTTTTGGAGGGGCAGAAAGatgaa                                                                                                                                         >  4:254016/1‑150 (MQ=255)

       aCCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatg                                                                                                                                       <  3:1241658/151‑1 (MQ=255)

       aCCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCATCGGCATATAGATTTTGAAGGTGAAGAAATaagaat                                                                                                                                        >  4:105995/1‑150 (MQ=37)

         ccGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGAc                                                                                                                                     >  5:729642/1‑151 (MQ=255)

          cGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACt                                                                                                                                    <  1:1431531/151‑1 (MQ=255)

          cGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAAACCG‑CTTCGGCGGGTTTTTCCTTTTGTAGGGACAGAAGAATGAATGAc                                                                                                                                     >  6:319189/1‑150 (MQ=25)

           ggg‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                                                                                                                                   <  7:141275/151‑1 (MQ=255)

           ggg‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                                                                                                                                   <  5:1095799/151‑1 (MQ=255)

                 ggCTGCTTGAGCCAGTGAGCGATTGCGGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCac                                                                                                                              <  7:994390/151‑1 (MQ=255)

                  gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAGCCGGCTTATAGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTAGGCGGGTTTTTGCTTTTGGAGGGGCTGAAAGATGAATGAATGTCCacg                                                                                                                             >  5:1319323/1‑151 (MQ=255)

                  gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacg                                                                                                                             <  3:182743/151‑1 (MQ=255)

                   ctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACACAAAACCCATTATCAATCACTTTAACCAGAATTACATAAAAACACT‑CTTCCTCGAGTTTTTGATTTTGGAGGTGCAGAAAGATGAATAACtataaacg                                                                                                                             >  2:535243/1‑143 (MQ=255)

                     gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgacg                                                                                                                          <  2:795952/151‑1 (MQ=255)

                     gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                                           <  2:1020381/150‑1 (MQ=255)

                     gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacga                                                                                                                            <  7:1227848/149‑1 (MQ=255)

                       ttGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACACGGCTTATCCGTCAGTTTCAGCTGATTAACGTAAAAACCCC‑CTTCGGCCCCTTTTTGGTATTGCAGGTGCAGAAAGATGAATGACTGTCCACGACCCt                                                                                                                        >  1:1045228/1‑151 (MQ=21)

                        ggAGCCAGGGAGCCATGGCTGGCCTAGCTGAATGACGGTCCACGACAGAACCCGGCTTATCGGGCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGAGGAATGACTGTCCACGACGCta                                                                                                                       <  2:932509/150‑1 (MQ=18)

                            ccAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACTCTATAcc                                                                                                                   >  7:1501260/1‑151 (MQ=255)

                            ccAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAAACCGCTAATAAGACAGTTTAAAATGATTTAAGTAAAAAACAA‑CTTCTAAATGTTTTTGCTATTGGATTGACAGAAATATGAATGACTGTACACGACACTATAcc                                                                                                                   >  8:1323907/1‑151 (MQ=255)

                             cAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTTATGTAACTAAAAACCCG‑AATCGCCTGGTTTTTGCTATTTGAGTGGCTCAAAGATGACTGCCTGTCCTAGACACtaacac                                                                                                                   >  6:822029/1‑146 (MQ=19)

                             cAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATAGGTCAGTTTAACCTGATTTACGTAAAAACCAG‑CTTAGGCGGGTGTTTGGTTTTGGAGGGGCAGATAGATGAATGACTGTCCACGACGCTATAccc                                                                                                                  >  5:236888/1‑151 (MQ=21)

                                 gAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAAACCCCAAATCAATCAGATTCAAATGAATTACGTAAAAACAAG‑CTTCTGCCGCTTTTTGATTTTGTAGGTGCAGAAAGATGAATTAATATAAAAGACACTATACCCaaa                                                                                                               >  2:822779/1‑150 (MQ=12)

                                    cGATGGCTGGCCTAGATGAATGACTGGCCACGACAGAACCCGGCTTAGCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaagaa                                                                                                           <  4:197084/151‑1 (MQ=255)

                                      aTTGCTGGCGTAGATCAATGACTGTCCCCCACAGACCCCGGCTTAGCGGCCTGTTTCACCTGATTTACGTAAAACCCGG‑CTTCGGGGGGTTTTTGGTTTTGGAGGGGCAGAAAGATCAATGACTGGCCACGACGCTATACCCAaaagaaa                                                                                                          <  7:1303324/150‑1 (MQ=255)

                                      aTTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACGGGATTTACGTAAAAACGCG‑CTTCGGCGGGGTTTTGCTTTTGGAGTGTCAGAAAGATGAATTACTGTCCACGACCCTATACCCAaaagaa                                                                                                           >  3:585860/1‑149 (MQ=18)

                                       ttgggggCATAGATTAATGAGTGTCCACCACATAACCCGGCTTATGGGTCAGTTTCACATGATTTACGTAAAAACCCG‑CTTCGGGGGTTTTTTGCTTTTGGAGGGGAATAAATATGAATGACTGTCCACGACGCTATACCCAaaagaaag                                                                                                         <  7:1065076/145‑1 (MQ=17)

                                         gCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCg                                                                                                       <  5:849382/150‑1 (MQ=255)

                                          cTGGCCTAGATGAATGACTGTCCACGACAGAAAACGCATTATAGCTCCGTTTCACATGATTTAAGTAAAAAACAC‑CTTCGGCGGGTTTTTGGTTTTGGAGGGGCAGAAAGATGAATGAATGTCCACGACGCTATACCCCACAGAAAGCggg                                                                                                     >  7:344998/1‑150 (MQ=17)

                                             gCCTAGATGAATGACTGCCCACGACAGAACCCGGGTTAGCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGGTTTTGGAGGGGCAGAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTa                                                                                                  <  7:259232/151‑1 (MQ=255)

                                              ccTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGTTTTTTGCTTTTGTCGGGGCACAAATATCAATGACTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTAt                                                                                                 <  8:206981/151‑1 (MQ=18)

                                                 aGATGAATGACTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCATAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCgg                                                                                              <  6:647398/151‑1 (MQ=255)

                                                 aGATGAATGACTGACCACGACAGACCCCGGGTTATCGGTCAGTTTCACCTGATTTACGTAAAACCCCG‑CTTCGGCGGGTTTTTGCTTTTGGCGGGGCAGAAAGATGAATGACTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTATCg                                                                                               <  7:924350/150‑1 (MQ=17)

                                                  gATTAATGACTGACCACGACAGAACCCGCGTTATGGGTCTGTTTCACCTGATTTACTTAAAAACCCG‑CTGGGGGGGTTTTTTGGTTTTGGTGGGGCACAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGt                                                                                             <  5:500809/151‑1 (MQ=17)

                                                  gATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTAGGAAAACACCAC‑CTGAGACCGCGATTGTCTAATGGAGGGCCAGAATGCTGAATGCGTGTCGAAGACCCTATACCCAAAAGAAAGCggggggggggg                                                                                             >  7:125930/1‑142 (MQ=12)

                                                  gATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTAAGTAAAAACCCG‑CTGCGCAGGGGTTTTGAGTTTGGAGGGGCAGAAAGATGAATTACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGt                                                                                             >  7:93612/1‑151 (MQ=17)

                                                   tttaatgaCTGTCTACGACAGAACCCGGCTTATCGGTCAGTTTCACCTTATTTAGGTAAAAACCCG‑CTTCGGGGGTTTTTTGGTTTTGGAGGGGCAGAAATATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGt                                                                                             <  5:681335/147‑1 (MQ=17)

                                                   atcaatgaCTGTCCACGCCACAACCCGGCTTTTCGGTCAGTTTCACGTCATTTACGTAAAAACCCG‑CTTGGGCGGGTTTTTGCTTTTGGCGGGGCACAAACTTCAATCACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGt                                                                                             <  3:847212/147‑1 (MQ=14)

                                                    tgaatgaCTGTCCACGACAGAAACCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCT‑CTTCGTCGGGATTTTGCTTTTTGAGGGTCATAAACATGAATGAATAACAAAGAATATATACCAAAAAGAAAGCGGCTTATCGGTc                                                                                            >  4:1201592/1‑150 (MQ=14)

                                                    taactgaCTCTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGa                                                                                                                                             <  3:333176/97‑1 (MQ=255)

                                                       atgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACCCTATACCCAaaggaaaaccggagaggggtc                                                                                            >  4:1532065/1‑133 (MQ=17)

                                                         gaCTGTCCACGACAGAACCCGGCTTATGGGTAATTTTCACCTGATTTACTTAAAAACGCG‑CTTCGGCGGGTTTTTGTTTTTGTTGGGGCTTAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTc                                                                                      <  6:1341448/151‑1 (MQ=17)

                                                         gaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGATGAATGACTGTCCa                                                                                                                               >  4:1148355/1‑110 (MQ=255)

                                                         gaCTCTCCACTACAGAACCCGGCTTATCGGTAATTTTCACCTGATTTACTTAAAAACCCG‑CTTCGGCGGTTTTTTGTTTTTGGTGGGGCATAAATATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                                       <  6:1057421/150‑1 (MQ=14)

                                                         gaCTCTCCACCACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGATGAATGACTGTCCa                                                                                                                               <  3:1138905/110‑1 (MQ=255)

                                                          aCTGTACACAACAGAACCCTGCTTAGCGGTAATTTTCCCCTGATTTACTTAAAACCCGG‑CTGAGGCGGGTTTTTGCTTTTGAAGCCGCACAAAGATGAATGACTTTCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTc                                                                                      <  1:72638/150‑1 (MQ=11)

                                                           cTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTAACCTGATTTACGTAATAACCAC‑CGTAGACCGTATTTTACTTTAAGAGTGTCACAAAAATTAATGAAATACGACGACGCAATAACCAAAAAAAAGAGGCTTATCGGTCAgggtg                                                                                      >  2:1453167/1‑145 (MQ=255)

                                                           cTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAAACCT‑CTTAAGCGCGTATTTTATTTAAAAAGGGCAGAAATATGAATAACTGTCCAAGAACCTATAAACAAAAAAAAAAAGCTTATCGGTCAGTTTGAc                                                                                    >  2:387722/1‑151 (MQ=255)

                                                            tGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAATATGAAGGCCTGTGCACTCCGCTCTACACATACGAAAGCGGCTTATCGGTCAGgttgtc                                                                                    >  5:960228/1‑147 (MQ=14)

                                                              cccACGACACAACCCGGCTTCTCGGGCAGTTTCACCTGATTAAGGAAAAACCCCG‑CTGCGGGGGTTTTTTGGTTTTGGAGGGGCACAAAGATCAATCACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTg                                                                                 <  7:537048/150‑1 (MQ=14)

                                                               cctctactgAACCCTGATTATGGGTAATTTTCTCGTCATTTA‑GTTAAAACCGGTCTGCTGCGGGTTTTTGTTTTTGTTGGTGGTCAATGATGAATGACTTTCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTgg                                                                                <  2:286551/143‑1 (MQ=2)

                                                                caccacACAACCCGGCTTATGGTTAATTTTCACTTGATTTACTTAAAAACGGT‑CTGCTGCGGTTTTTTGTTTTTGTTGGTGTTTAATTATGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGt                                                                               <  2:1226206/151‑1 (MQ=9)

                                                                   caCAGAACCCGGCTTGTCGGGCAGTTTCACCTGATTTAGGTAAAACCCCT‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAACGTTGATTGCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTa                                                                            <  1:1129911/150‑1 (MQ=14)

                                                                   aaccgAACCCGGCTTATCGGACAGTTGCCCGTGATATCCTAAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGCAGGGGCACAAAGATCAATGACTCGAAACGACGCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTa                                                                            <  3:596681/147‑1 (MQ=14)

                                                                        aaCCCGGCTTCTGGGCCATTTGCGCGTGATTTACGTAACACCGCG‑CTGCTGCGGGTTTTTGCTTTTGTGGGGGCCCAACGATGAATGACTCGCCACGACGCTTTACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaa                                                                       <  4:1332634/151‑1 (MQ=9)

                                                                         aCCCGGCTTATGGGTCATTTTCACCTTATTTAGGTAAAAACGCC‑CTTCGGCGGGTTTTTGCTTTTGTAGGGGCCCAAACATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaa                                                                      <  1:1201165/151‑1 (MQ=14)

                                                                           ccGGCTTATGGCGCAGTTGCACGTCATTAAGGTAAAACCGCT‑CTGCGGCGGGTTTTTGCTTTTGTGGGCGCGCAACTTTGAATGACTTTCCACGACTCTATACCCAAAAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaaa                                                                    <  3:533672/151‑1 (MQ=2)

                                                                           ccGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACCCAAAAGAAAGc                                                                                                        >  3:886376/1‑115 (MQ=11)

                                                                            cGGCTTAGCGGGAATTTGCACGTCATTAACTAAACAACGGG‑CTTCTGCGGTTTTTTGGTTTTGGAGGGGCACAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAAc                                                                   <  1:267184/151‑1 (MQ=11)

                                                                              gCTTATGGGTATTTTTCAGGTTATTTACG‑AAAAACCGGCCTGCGGGGGTTTTTTGGTTTTGCTGCCGGATAAATTTGATTGAGTTGCCAAGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAAccc                                                                 <  7:460031/151‑1 (MQ=2)

                                                                                 tatggcgaactttCTCCTGATTTACTTAATACCGGT‑CTTCGGCGGGTTTTTGTTTTTGTTGGTGGCCTATCATGAATGACTTGCCAAGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGc                                                               <  2:1340191/140‑1 (MQ=2)

                                                                                   aggttaattttCACGTTATTTCGTTAAAAACCGT‑CTTCTGCGGGTTTTTGCTTTTGTAGGTGCTTAATTATGATTGACTTTCCACGACTCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTc                                                            <  4:727061/143‑1 (MQ=2)

                                                                                    ggggcattttCCCCTCATTAACGTAACACCCCG‑CTTCGGCGGGTTTTTGCTTTTGTCGGGGCAGAAAGATGAATGCCTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcg                                                           <  1:1287900/144‑1 (MQ=11)

                                                                                    ggggaattttCAGGTGATTTACTTAACACCGGT‑CTGCGGCGGGTTTTTTTTTTTGGGGGTTCTTCACTATGATGGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcg                                                           <  8:52709/144‑1 (MQ=2)

                                                                                    cGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑ATACGGCGGGTTTAAGATTTTTGAGTGCCATAAATAAAAAACACTTTCCACGACACAAAAAACAAAAGCAAACAACATATCGCAAAATATCACATGATATACGTaaaaaaaaaata                                                             >  8:1105325/1‑142 (MQ=255)

                                                                                        cAGTTTCACCTGATCTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGATGTACGTaaaaaacagatacga                                                          >  8:701471/1‑138 (MQ=11)

                                                                                        atgTTTCCCTTGATTT‑CCTAAATACCCGTCTTCTGCGGGTTTTTGTTTTTGTTGGTGCTTAATTTTGATTGACTTTACACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcggcg                                                        <  8:323271/148‑1 (MQ=2)

                                                                                           tttgcggtGATTTACTTAAAACGCGT‑CTTAGGGGGGTTTTTGTTTTTGTTGGCGCCGAATTTTGAATGACTTGCCACGACGCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGt                                                     <  6:21158/143‑1 (MQ=2)

                                                                                            ttCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGATTTACGTAAAAACCCGGTACGGCGGGttt                                                   >  7:647124/1‑151 (MQ=255)

                                                                                            gaccccTGA‑TTACGTTAAAACCCGTCTCGGGGGGGTTTTGGCTTTGGGCGGGGCGTAATAATGATGCACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGtt                                                    <  5:296887/146‑1 (MQ=2)

                                                                                             gcccttcaTTTACTTAAAACCCCT‑CTGCGGCGGGTTTTTGTTTTGGTCGGGGCTCAACTATCATTGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttt                                                   <  8:618390/143‑1 (MQ=11)

                                                                                              cACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACGTGGTGTACGTAAAAAACCGGGTCGGCGGGtttct                                                 >  3:1451097/1‑149 (MQ=11)

                                                                                              cACCTCATTTCCGTAAAACCCCG‑CTGCGGCGGGTTTTTGCTTTTGGCGGGGCACAAACATGAAGGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttttt                                                 <  1:1066453/151‑1 (MQ=11)

                                                                                                   gATTTCCGTAAAACCCCG‑CGTCGGCGGGTTTTTGGTTTTGCCGGCGCCCAAAGATGAAGGCCTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCttt                                            <  3:1468790/151‑1 (MQ=11)

                                                                                                   gATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGGGGCGTAGCGGTCAGTTTCACGTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCttt                                            >  7:600102/1‑151 (MQ=255)

                                                                                                         cGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAAACCGATGCGGCGGGTTTTAGATTTTTGAgt                                      >  6:443156/1‑150 (MQ=11)

                                                                                                             aaaaCCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCCCGACCACGTAACCGAACCACTGGGGCGTAGCGGGCAGTTTAAGCTGGTTTACTTAAAAAACCGCGTCCGCCGGTATTTCGGGTTGGAgggggc                                  >  7:1409028/1‑149 (MQ=2)

                                                                                                                aaccgccTTCTGCGGGTTTTTGCTTTTGTACGTGACGACATATCATTGACTTTCCCAGACGCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                                <  8:1042186/145‑1 (MQ=11)

                                                                                                                 cccggctgccgcggGTTTTTGATTTTGCAGGCGCCTAAAGATGAATGACTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                               <  1:169147/141‑1 (MQ=11)

                                                                                                                 cccgccTTCTGCGGGTTTTTGTTTTTGTTCGCGCTATAACTTGATTGACTTGCGTCTACGCTATCCCCTAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                               <  8:391223/146‑1 (MQ=2)

                                                                                                                  ccttctttcgcggGTTTTTGTTTTTGTTCGTGTATAATTTTGATTGACTTTCAAATAGGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa                              <  6:980946/142‑1 (MQ=11)

                                                                                                                   cG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAACCCGATACGGCGGGTTTATGATTTTTGAGGGCCatat                               >  8:560235/1‑145 (MQ=11)

                                                                                                                        tCGCCGGGTTTTTGCTTTTGGCGCGGCCCAACGATGAAGGACTGCCCACGACGCTATACCCAACAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatg                          <  4:405982/149‑1 (MQ=11)

                                                                                                                         ccgcggGTTTTTGGTTTTGTCGGGGCGAACAATTGATTGACTTGCCAAGACGCTATCCCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                       <  1:1392386/149‑1 (MQ=11)

                                                                                                                          ggcggGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTACAACACAATATACACATATAAAAAAAGCTTATCGCTCAGTTTCAACTGGTTTAAGTAAAAACAATCTTATAAATAATTTTGCTTTTGGAGGGGCAGAAAGatgaatg                      >  2:1496282/1‑151 (MQ=2)

                                                                                                                             gggTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTTACCTGATTTACGTAATAAACCGATCCGGCGGGTTTCAGATTTTTGAGTGCCAGatataaataacac                    >  2:206286/1‑138 (MQ=14)

                                                                                                                               gTTTTTGCTTTTGGAGGGGCAGAAAGATAAATCACAATCCACTACTCTATACAAAAAATAAAGCTCATTATCTATCACTTTCACCTGGTATAATAAAAAAAAAAAATCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAgatggatgggagg                 >  8:947033/1‑146 (MQ=2)

                                                                                                                                      gttttGGCGGGGCAGAAAGATGAATGACTCGCCAAGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGAAGAATGACTGTCCACGac          <  7:1198695/150‑1 (MQ=255)

                                                                                                                                         ttGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCt                                                                                                    >  3:780108/1‑58 (MQ=35)

                                                                                                                                         ttGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAAACCGATACGGCGGGTTTAAGCTTTTTGAGTGGCAGAAAGAAAAATCAATGTCAAAGacaat       >  8:1386660/1‑149 (MQ=11)

GTACTGAACCCGGG‑TAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACC  >  DH1/613451‑613741

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