BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	613,597	0	G	C	31.0%	240.4 / 50.9	209	intergenic (+458/+565)	garK/yhaC	glycerate 2‑kinase/pentapeptide repeats‑related protein
Reads supporting (aligned to +/- strand): ref base G (60/60); new base C (2/52); total (69/140)
Fisher's exact test for biased strand distribution p-value = 1.63e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

CCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTG  >  DH1/613449‑613725

    aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑Gat                                                                                                                            <  8:302131/151‑1 (MQ=255)

      tGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGATA‑Gaacg                                                                                                                          >  7:1277464/1‑148 (MQ=255)

       gatgcggggTATTCTGCTTGAGCGAGTTATCGATTGCTGGACTAGATTATTGACTGTACACTAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTTCTTTTGGAGGGGCAGAAA‑Gat                                                                                                                            <  2:1325382/142‑1 (MQ=14)

        aaCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGTCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑Gatgaat                                                                                                                        <  4:500622/151‑1 (MQ=255)

        aaCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑Gatgaat                                                                                                                        <  8:682649/151‑1 (MQ=255)

        aaCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTAGATCGGAAGATA‑Gaagtgt                                                                                                                        >  8:1024594/1‑146 (MQ=25)

         aCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAAACGGCTTATCTGTCAGGTTCACGTGATTTACGTAAAAACCCGCTTCGGCGGGTTGTTGCTTTTGGAGGGGCAGAAC‑Gatgcatg                                                                                                                       >  5:845669/1‑151 (MQ=25)

          ccgggttattgtttttGAGGCAGTTATCGATTGATTGCCTAGATCACTGTCTATACACCACAGAACCCGGCTTATCGGTCAGTTTCACCTGCTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑Gatgaatga                                                                                                                      <  7:223548/137‑1 (MQ=14)

          cccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑Gatgaatga                                                                                                                      <  7:919994/151‑1 (MQ=255)

           cgggtttttgggctTGAGGCGGGGATCGATTGCTGGCCTAGCCAAGAGACTCTCCCCCACAGAACCCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAA‑GATGAATGAc                                                                                                                     <  1:1330607/140‑1 (MQ=14)

           ccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGAc                                                                                                                     <  2:1019769/151‑1 (MQ=255)

             gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTg                                                                                                                   >  1:459826/1‑151 (MQ=255)

             gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTg                                                                                                                   <  7:424530/151‑1 (MQ=255)

             gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAACCCCAATATCAGTCAGATAAACCAGAATTACGTAAAAACCATCTTCTGAGAATTTTTGCTTTTGTATGGGCAGAAA‑GATTAATAACt                                                                                                                    >  8:374604/1‑150 (MQ=255)

              ggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGAATTACGTAAAAACCACCTTCAGAGGAATATTGATTTTTTTTGTGGAGAAA‑GATGAATTACTGt                                                                                                                  >  6:1070417/1‑151 (MQ=25)

                 aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTccc                                                                                                               >  7:835350/1‑150 (MQ=255)

                 aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACACAAACCGGCTAATACGTCAGTTTCACCTGATTTACGTAAAAACCACAGTCGGCTGGATTTAGCTTTTGGAGGGGAACAAA‑GATGAATTACTGTcc                                                                                                                >  5:184279/1‑150 (MQ=25)

                  ggCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAA‑GATGAATGAcagaacgc                                                                                                              >  4:1171737/1‑144 (MQ=255)

                   gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCacg                                                                                                             <  7:1364507/151‑1 (MQ=255)

                   gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCacg                                                                                                             <  4:16550/151‑1 (MQ=255)

                   gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACACAACCCGCATTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGGTTTTCCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCccg                                                                                                             >  1:1403264/1‑151 (MQ=21)

                      gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCacgac                                                                                                           <  6:1295420/150‑1 (MQ=255)

                      gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCacgac                                                                                                           <  2:213868/150‑1 (MQ=255)

                      gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACACGACATATCTGTAAGTATAACCTGATTTAAGTAAAAAAAAGCTTCGTCGGCTTATTAATTTTTAAATAGCAGAAA‑GATGAATGACTGTCCacgagc                                                                                                          >  6:1091930/1‑149 (MQ=255)

                        ttGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAAACAGCTTATGCGCGCTGTTGTTTTGGGAGGTGAAGAGA‑GAGGAATGACTATCCacgcggct                                                                                                        >  5:1247578/1‑146 (MQ=21)

                         tGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACAAGATACGAAGAGATTATGGTATTGAATGGGGAGAAA‑GATTAATGAATATCAACGAGGca                                                                                                        >  6:322208/1‑149 (MQ=19)

                          gAGCCATTTAGCTATTGCTGGCCTAGATGAATGACTGTCCCCGACATAACCCGGCTTATCGGTCAGTTTCACCTGATTTAGTTAAAAACCCGCTTCGGCGGTTTTTTGCTTTTGGAGGGGCATAAA‑TATTAATGACTGTCCACGACGCtat                                                                                                      <  8:712747/151‑1 (MQ=18)

                              cAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCCGCTTATAGGTCAGTTTCACATGATATACGTAAAAACCAGCTTACGCGGGTTTTTCCGTTTGGTGGGGCAGTAA‑GATGAATCACTGTCCACGACCCTATAccc                                                                                                  >  7:774787/1‑151 (MQ=21)

                               aGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCa                                                                                                 <  8:465290/151‑1 (MQ=255)

                               aGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGAATAAAGGTCAGTTTCACCTGATTTACGTAAAAACCCCCTTCGGCGGGGTTTTGCTGTTGTAGGGGGAGAAA‑GATGAATGACTGTCCACGCCGCGATccccc                                                                                                 >  5:42557/1‑150 (MQ=21)

                                gTGAGCGATGGCTGGCCTAGATGAATGACTGGCCACGACAGAACCCGGCTTATCGGGCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCaa                                                                                                <  5:688113/151‑1 (MQ=255)

                                 tGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGATTATAGCTCAGTTTCGCATACTTTACCTAAAAACACACTTCCTAGCCTTTTTTCTTTAGTATGGGGAGAAA‑GATGAATGACTGTCGACGCCGCTATACCCaa                                                                                                >  5:806988/1‑150 (MQ=255)

                                  gAGCTATTGCTGGCCTAGATTAATGACTGCCCACGACACAACCCGGCTTATCGGTCAGTTTCACCTGATTTACTTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGCATGTGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCaaa                                                                                               <  6:245824/150‑1 (MQ=18)

                                  gAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACACAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACAGAGGTGAACCGGTATGTTCTTTTCTCGCGTCAGAAA‑GATGAATGACTGTTGACGGCTCTTTACCCaaa                                                                                               >  3:321178/1‑150 (MQ=16)

                                   aGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCTCCTTATTTACGTAAAAACACGCTTCGGCGGGTTTTTGCTGTTGTTGCGGCAGAAA‑GATGAATTACTGTCCCCGACGCTATAcccgca                                                                                               >  1:521000/1‑146 (MQ=21)

                                    gCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGGGTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACCCTATACCCaaa                                                                                               >  1:1145740/1‑148 (MQ=255)

                                      gATTGCTGGCCTAGATGAATCACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAaaagaaa                                                                                          <  4:650740/151‑1 (MQ=255)

                                       aTTGCTGGCCTAGATGAATGACTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGGGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAaaagaaag                                                                                         <  5:119156/151‑1 (MQ=255)

                                       aTTGCTGGCCTAGATGAATGACTGCCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTTATTTACTTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGGCCACGACGCTATACCCAaaagaaag                                                                                         <  7:592682/151‑1 (MQ=255)

                                        ttGGTGGCGTAGATGAATGAGTGTCCACGACAGAACCCGGCTTATGGGTCAGTTTCACCTGATTTACGTAAAACCCGTCTTGGGCGGTTTTTTGCTTTTGGAGGGGCAGAAA‑TATGATTGACTGTCCACGACGCTATACCCAaaagaaag                                                                                         <  3:1369666/150‑1 (MQ=17)

                                          ggtGGCCTAGATTATTCCCTGCCCACGACATAACCCGGGTTCGGGGGCATTTTCAGGTTTTTTCCGTAAAACCGCGTTTCGGGGGTTTTTTGCTTTTGGCGGGGCACAAA‑TATGAATCCATGCCCACGACGCTATACCCAAAAGAAAGCg                                                                                       <  3:1448409/148‑1 (MQ=255)

                                            tGGCCTAGATGAATGACTGTCCACGACACAAAACCCAAAATAAATCAGATACAAATCAATTACTTAAAAAAAAGCTTCTTAAGGTTTTTGATTATGTAGGTGCAGAAA‑GAATAATTAATATAAACGACACTATACCCAAAAGAAAGCggat                                                                                    >  6:828091/1‑149 (MQ=255)

                                             ggcgttgATGAATGACTGTCCACGACAGAACCCGGCTTATGGGTGTGTTTAACCTGATTTACGTAAAAACCGGGTTCGGCGGGTTTTTGCTTTTGGAGGGGAAGAAA‑GATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCtt                                                                                   <  1:326030/145‑1 (MQ=17)

                                             ggCCTATATTAATCACTCCCCACGACAGAACCCGGCTTAGCGTTCAGTTTCACCTTATTTACTTAAAAACCCTCTTCGGCGGTTTTTTGATTTTTTATGTTCTTAAA‑TATCAATTACTGTCCACGACGCTATACCCAAAAGAAAGCGGCt                                                                                    <  4:1211963/150‑1 (MQ=12)

                                             ggCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCGACGACGCTCTACCCAAAAGAAAGCGGCtt                                                                                   >  5:569674/1‑151 (MQ=255)

                                             ggCCTAGATGAATGACTCTCCACCACAGAACCCGGCTTATCGGTCATTTTCACCTGATTTACTTAAAAACCATTTTCGGCGGTTTTTTGCTTTTGTATGTTCTGAAA‑TATGAATTACTGTCCACGACGCTATACCCAAAAGAAAGCGGCt                                                                                    <  2:92958/150‑1 (MQ=17)

                                               cgtAGATGATTGAGTGGCCACGACAGACCCGGGCTTCTCTGTCAGTTTCCCCTGATTTACTTAAAAACCCGCTTCGGGGGGTTTTTGCTTTTGGAGGGGCACAAA‑GATTAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTAt                                                                                 <  3:747973/149‑1 (MQ=14)

                                                cTAGATGAATGACTGTCCACGACAGAAAACGCCTAATAGGTCAGTTTCACCTGATTTAAGTAAAAAACCGCTTAGGAGGGTTTTTACTTTTGGAGGGGCAGAAA‑GATGAATTAATCTCCACGACACTATACCCAaaagaaaatcgcataac                                                                                >  8:641450/1‑141 (MQ=14)

                                                  atattaatgaGTGTCCACGACAGAACCCGGCTTATCGGTCATTTTCACGTGATTTACTTAAAAACCCGGTTGGGGGTTTTTTTGATTTTGGAGGGGCACAAA‑TATTAATGACTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                                <  7:1053733/144‑1 (MQ=17)

                                                  agtttattgaATGTCCACCACAGAACACGGATTATGGGTCTTTTTCACCTGATTTACGTAAAAACGCGCTTCGGCGGGTTTTTGCTTTTGTTGGGGCCTAAC‑GATGATTTACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                                <  3:130710/142‑1 (MQ=14)

                                                  aGATTAATGACTTTCCACGAAAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACTTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCACAAA‑GATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCgg                                                                              <  8:226186/151‑1 (MQ=14)

                                                    attaatgaCTGCCCACTACAGAACCAGGCTTATCGGTCAGTTTCACGTGATTTACGTAAAAACCGGCTGCGGCGGGTTTTTGTTTTTGTAGGGGCACAAA‑TATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTc                                                                            <  6:819417/148‑1 (MQ=14)

                                                      tattgaCTTCCCACTACACAACCCGGTTTTTGGGTATTTTTCAGTTGATTTACTTAATAACCGTCTTCTGCGGTTTTTTGTTTTTGTTGGGGCACAAA‑TTTTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                           <  8:551861/147‑1 (MQ=2)

                                                      gaatgaCTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                           <  5:198017/150‑1 (MQ=255)

                                                      gaatgaCTGCCCACCACAGAACCCGGCTTATCGGTCAGTTGCACCTGATTTCCGTAAAAACCCGCTTCGGGGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                           <  3:1402694/150‑1 (MQ=17)

                                                       aatgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAAGCAACTTCGGACGTGTTTTCATTATTGAGGGGCATCAA‑GATGGTTTACTGTCAACGACTCTATAAAAAAAAGTAAGCGGGTTATCGGTCAgg                                                                         >  7:192172/1‑150 (MQ=12)

                                                          gaCTGTCCAGTACAGAACCCGGGTTATGGGGAAGTTTCACCTGATTTACTTAAAACCGGGCTGCTGCGGGTTTTTGGTTTTGTAGGGGCAGAAA‑GATCAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTc                                                                      <  1:515006/151‑1 (MQ=14)

                                                           aCTGGCCCCGACAGAACCCGGCTTATCGGGAAGTTTCACGTGATTTACGTAACAACGCGCTTCGGCGGTTTTTTGGTTTTGGTGGGGCACAAC‑GATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCa                                                                     <  3:1183918/151‑1 (MQ=14)

                                                             tGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTcc                                                                                                                >  7:519941/1‑106 (MQ=255)

                                                                    accAGAACCCGGCTTATGGGACATTTTCACTTCATTTCCGTAAAAACGCTCTTCGGCGGGTTTTTGTTTTTGTAGGGGCACAAT‑TTTGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTa                                                            <  2:1343841/149‑1 (MQ=14)

                                                                      cacaaCCCGCCTTATCGGTCTCTTTCACCTGATTTACTTAAAACCCGGCTTCGGCGGGTTTTTGCTTTTGGCGGGGCAGAAA‑TTTGAAGGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACg                                                          <  4:1354434/148‑1 (MQ=14)

                                                                      cAGAACCCGGCTTATCGGTCAGTTTAACCTGAATAACGTAAAAACCACCTTAAGCGAAATATTGATAATATAAGTGAAGAAA‑TATGAATTACACTAAACAAAAAAAAAAAAAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTgg                                                                >  2:1088179/1‑145 (MQ=255)

                                                                       acaccccGGCTTATCGGGCAGTTTCACCTGATTTCCGTAACACCCCGCGGCGGCGGGTTTTTGGTTTTGCTGGGGCACAAC‑GATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGt                                                         <  7:763252/147‑1 (MQ=11)

                                                                        caaCCCGGCTTATCGGGCAGTTGCACCTGATTTACTTAAAAACCCGCTGCGGCGGGTTTTTGCTTTTGCAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTa                                                        <  1:877957/150‑1 (MQ=255)

                                                                             cGGCTTATCGGTCATTATCACATGATTTACGTAAAAACAAACTTAGGATAGTTTTTGCTTTTGGAGGGTCAGAAA‑GATAAATGAATGTAAAAGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                           >  8:1162554/1‑127 (MQ=2)

                                                                                  tATCGGTCAGTTACACCTGATTTACGTAAAAACACGCTTCGGAGTGTTTTTGCTTTTGGAGGGGCATAAA‑GATGAATGACTGTCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCAGCTGGAGAACGTaatctcgggat                                              >  8:513567/1‑142 (MQ=11)

                                                                                  caTCGGCCACTTGCACCTGATTTCCGTAACCACCGGCTTCGGCGGGTTTTTGGTTTTGGTGGGGCCTAAA‑TATGATTGCCTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCg                                                <  2:333471/148‑1 (MQ=11)

                                                                                       gTCAGTTTCAGCTGATTTCCTTAAAAACGCGCTGCGGCGGGTTTTTGCTTTTGGCGGGGCACAAA‑CTTGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcggc                                         <  3:267479/151‑1 (MQ=11)

                                                                                       gTCAGTTTCACCTGATTTACGTAAAAACCCGATTCGGCGGGTTTATGATTTTTGAGTGCCATAAA‑TAAAAATAACTGAACACCAACCTATACCCAAAAAAAATATACTTATCGGTCACTATAACATGTTGTACTTaaaaaaaaaaaaaggc                                         >  4:577451/1‑140 (MQ=255)

                                                                                        tCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGATA‑GAGGAATGACTGTCCCCGACGCGATAGCCAAAAGAAAGCGGCTGATCGGTCAGTTTAACCTGCTGTACGTAATAACCCGCTTcggcg                                        >  3:25451/1‑151 (MQ=11)

                                                                                        tCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTACACGACCCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTGTGACCTGGTGTACGTAATACCCCggtaggggg                                        >  5:1326989/1‑143 (MQ=11)

                                                                                         aattttCACGTCATTTCCTTAAAAACCGTCTGCTGCGGGTTTTTGTTTTTGTAGGTGCACAAC‑GATGAATGCCTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcggcgg                                       <  2:778248/148‑1 (MQ=11)

                                                                                          aGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTTCTTTTGGAGGGGCAGAAA‑AAAGAATGACAGTCCTCGAGGGTAAAACAGAAAAAAATAGTCTTATAGTTCAGTTTCGCCTTATTTAATTAAAAAAAATCTTCGGCggg                                      >  4:617267/1‑151 (MQ=2)

                                                                                             ttCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGCTTTACGTAATAACCCGTTGCGGCGGGttt                                   >  5:339487/1‑151 (MQ=255)

                                                                                              gcACCTCATTACCGAAAAACCCCGCTGCCGCGGGGATTGGCTTTTGGAGGGGCCGACC‑GAGGAATCACGCGCCACGACGCGATCCCCAACAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGtttt                                  <  6:529537/150‑1 (MQ=2)

                                                                                               cACCTGATTTACGTAAACACCGTCTTCGGCGGGTTTTTGGTTTTGGCGGGGCCCAAC‑CATGAATGACTCTCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGtttt                                  <  1:820773/150‑1 (MQ=11)

                                                                                               cACCTCATATACGAAAAACCCCGCTTCGGCGGGTTTTTGCTTTTGGCCGGGCACAAC‑CCTCAAGGACGGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttttt                                 <  3:1019980/151‑1 (MQ=11)

                                                                                                 ccTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGATA‑GCTGGAGGACAGGGCACGACGCGATTCACAAGAGAACGGGGCGTATGGGGGAGTGTCCGGTGGTTTACGTAAAACCCCGGGTGGGCGGGttttggg                               >  7:1148258/1‑148 (MQ=2)

                                                                                                  cTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCCGTCAGTTTCACCTGGTTTACGTAATAACCCGATGCGGCGGGTTtgtgat                              >  7:667934/1‑149 (MQ=11)

                                                                                                  cTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGAAGAAA‑GATGAATGACTGTCCACGACCATATAACCAAAAGAAAGAGGCTTATCGTTCAGTTTAACATGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCt                              >  4:436767/1‑151 (MQ=11)

                                                                                                   tGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGTTGTACGTAAAAAACCGGTCCGGCGGGTTTGTGAtt                             >  3:1354304/1‑151 (MQ=11)

                                                                                                   tGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGATTTACGTAAAAAACCGGTGCGGCGGGTTTGTGCtt                             >  1:1008280/1‑151 (MQ=11)

                                                                                                    gATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTTGATGGGCAGATA‑GATGAATGGCTGTCCAGGACGCTATAGCCAAAAGGAAGCGGGTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGGTTGGGCGGGGTTTTGGttt                            >  5:1153038/1‑151 (MQ=11)

                                                                                                     aTTTACGTAAAACCCCCGATCGGCGGGTTTTTGGTTTTGCCGGCGCACAAA‑CATGAATGACTGGCCACGACGCTATACCCAAAAGAAACCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                           <  3:529972/151‑1 (MQ=11)

                                                                                                         acttAAAAACCGTCTTCTGCGGGTTTTTGTTTTTGTTGGTGCATAAA‑TATGATTCCCGGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGa                        <  6:1305798/147‑1 (MQ=11)

                                                                                                         aCGAAAAAACCCGCTTCGCCGGGTTTTTGCTTTTGCAGGGGCACAAA‑CATGAATGACTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAg                       <  1:591850/151‑1 (MQ=11)

                                                                                                          cGTAAAACCGCTTTTGGGCGGGTTTTTGTTTTTGGTGGTGCATAAT‑TATGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAgg                      <  2:461406/151‑1 (MQ=11)

                                                                                                           gtaaacccggtctgcggcggGTTTTTGCTTTTGCCCGGGCCCAAC‑CTTGAAGGACTGGCCACGACGCTCTACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAggg                     <  5:1423567/137‑1 (MQ=11)

                                                                                                           gaaacccccccCTTCGCCGGGTTTTTGCTTTTGCCGGCGCCCAAA‑CATGAATGACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAggg                     <  5:750546/144‑1 (MQ=11)

                                                                                                            tAAAAACGCTCTGCTGCGGGTTTTTGTTTTTGTTGGTGCACAAT‑TATGATTGACTTTCGACGAGGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAg                       <  4:1240425/148‑1 (MQ=11)

                                                                                                                accgggCTTCTGCGGGTTTTTGTTTTTGTCGGTGGTTAAT‑TATGATTGACTTTCCTCGTCGCTATACCCTAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                <  4:803854/146‑1 (MQ=11)

                                                                                                                  ccccCTGCGGGGGGTTTTTGCTCTTGGCGCCCCCAAAC‑CATCAATGACTCGACAAGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa              <  7:1337873/151‑1 (MQ=11)

                                                                                                                  cccGCTGCGGCGGGTTTTGGCTTTTGCCGGCGCACAAG‑TATCATTCCCTCTCCACGACGCTATACCCTAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa              <  7:188357/151‑1 (MQ=11)

                                                                                                                      gtgaggcggGTTTTTGTTTTTGGTGGCGTATCAC‑TATGATGGACTCGACAAGACTCTATACCCTAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAATAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatg          <  1:981587/147‑1 (MQ=11)

                                                                                                                       ttCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGAAGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTTACCTGGTTTACGTAAAAAACCGGTTCGGCGGGTTTCTGCTTGTTGAgtgggctttagatct         >  1:280309/1‑140 (MQ=11)

                                                                                                                          ggcggGTTTTTGTTTTTGTCGGAG‑AGAAATTATGAATGACTTTCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatg      <  8:17640/151‑1 (MQ=11)

                                                                                                                          cgcggGTTTTTGTTTTTGCCGGTGCCCACT‑GTTGATTGACTTGCCACGACGCTCTACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatg      <  3:1228844/150‑1 (MQ=11)

                                                                                                                            cggGTTTTTGTTTTTGCTGGGGCCCTAT‑CATGAATGACTCTCCACTACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGAc    <  4:759454/151‑1 (MQ=11)

                                                                                                                             gggTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGATTTACGTAATAACCCGATGCGGCGGGTTTGAGCTTTTGGAGGGCCATaaataacaaactct   >  8:437173/1‑140 (MQ=11)

                                                                                                                              ggTTATTGTTTTTGTAGGGGCTTAAA‑TATGATTGACTTTCTAATACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg  <  8:1215362/151‑1 (MQ=11)

CCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAA‑GATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTG  >  DH1/613449‑613725

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: `‑`