BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	613,597	0	G	C	30.4%	168.5 / 43.2	167	intergenic (+458/+565)	garK/yhaC	glycerate 2‑kinase/pentapeptide repeats‑related protein
Reads supporting (aligned to +/- strand): ref base G (38/58); new base C (2/40); total (44/123)
Fisher's exact test for biased strand distribution p-value = 1.13e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

GCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGAC  >  DH1/613447‑613723

  ccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                              <  6:730416/151‑1 (MQ=255)

  ccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACACAAAACGCAATATAGCTCAGATTCACCTCATTTACGTAAAA‑ACACT‑CTTACGAGGGTTTTTTCTTTTGGAGGGCCAGa                                                                                                                               >  4:508914/1‑150 (MQ=255)

   cGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa                                                                                                                             <  6:611162/151‑1 (MQ=255)

     tACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTAAGTAAAA‑AAAAG‑CTTCGGCGGGTTATTGCTTTTGGATCGgaagaa                                                                                                                              >  8:1244217/1‑148 (MQ=37)

      aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGat                                                                                                                          <  2:148406/151‑1 (MQ=255)

      aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGa                                                                                                                           <  1:174979/150‑1 (MQ=255)

         gAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatga                                                                                                                        <  1:623695/150‑1 (MQ=255)

          aaCCCGTGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTCATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                                                                                                                      <  6:880997/151‑1 (MQ=255)

          aaCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGatgaat                                                                                                                      >  4:296234/1‑151 (MQ=255)

          aaCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCAAGACAGAACCCGGAATCTAGGTCACTATCAGCTGTTTTTCGTAAAA‑ACCCT‑CTTATGCGCCTTTTTACTATAGGAAGGACAGTAAGctgatt                                                                                                                      >  5:774144/1‑149 (MQ=255)

            cccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatga                                                                                                                    <  2:1156419/151‑1 (MQ=255)

              cGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaacctgattcact                                                                                                                  >  7:1292890/1‑140 (MQ=255)

               gggTAGGCTGCTTGAGGCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAGCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                                                                                                                 <  7:232691/151‑1 (MQ=255)

                  tAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGATTATCCGTCAGTTTCACCTGATTTACGTAAAA‑ACAAG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGACATAAAAatgacttagtgta                                                                                                               >  8:1218910/1‑141 (MQ=25)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACACGGATAATCCCTCAGATTCACATGAATTACGTAAAA‑ACACC‑ATTCGCAGCGTTTTTGCTTTTGTAGGGGCATAAAGATGAATGAATGTcc                                                                                                              >  3:838620/1‑150 (MQ=255)

                     gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacg                                                                                                           <  2:1272311/151‑1 (MQ=255)

                     gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGCCTTATCGGTCAGTTTCGCCTGATTTACCTAAAA‑CCCCG‑CTTCGGCTGGTGTTTGCTTTTGGAGGGGCAGAAAGCTGAATGCCTGGCCccg                                                                                                           >  1:510176/1‑151 (MQ=25)

                      ctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACACAAAAAGGAATAAAGCTCAGATTCAACTGAAATACGTAAAA‑ACAAT‑CTTAAGCTCGATTTTGATTAAGGAGGATAAGAAAGATGAATGACTGTCCacg                                                                                                           >  8:537533/1‑150 (MQ=255)

                        gctTGAGCGAGTGAGCTTTTGTTGGCCTAGTTGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACTTAAAA‑CCCCG‑CTTCGGCGGGTTTTTTCTTTTGGCGGGTCATAAAGATGAATGACTGTCCacgacg                                                                                                        <  5:932107/151‑1 (MQ=18)

                        gctTGAGCCAGTGTGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                         <  5:543130/150‑1 (MQ=255)

                        gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑TCCCA‑CGTCGCCGCGTATGAGCTCATGTAGGGGCTGAACGAGGAATGACTGTGCacgacg                                                                                                        >  7:682644/1‑151 (MQ=25)

                        gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgacg                                                                                                        <  3:712655/151‑1 (MQ=255)

                          ttGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCCGAAAATCCGTCAGTTTAACGTGATTTACGTAAAA‑ACACG‑ATTCGGCTGGAATTAGCTTTTGGAATGGCAGAAAGatgaatgaaaaaacaagaccct                                                                                                      >  1:783534/1‑137 (MQ=19)

                          ttGAGCAAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCTGTCAGTTTCAACTGAATTAACTAAAA‑ACACA‑CTTAGGCAGATATTTCCTGATGTAAGTAAAGAAAGATGAAGGAAAATCCacgaaaa                                                                                                       >  8:46704/1‑147 (MQ=255)

                           gccgCCAGGGTGCCATGGCTGGCCAAGCTCAATGACTGCCCACCCCAGACCCCGGCTTAGCGGGCAGTTTCACGTGATTTACGTAAAA‑ACCCG‑CTTCGGGGGGTTTTTGCTTTTGGAGGGGCATAAATATGAATTACTGTCCACGACGCt                                                                                                      <  8:1026823/147‑1 (MQ=14)

                            gAGCCAGTGTGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTTATTTACTTAAAA‑ACCCG‑CTTCGGGGGGTTTTTGTTTTTGTAGGGGCATAAATATGAATTACTGTCCACGACGCtat                                                                                                    <  2:959300/151‑1 (MQ=21)

                               ccAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAGAAGCAGTAAAGCTCAGAGTACACAGATAAACTTAAAA‑AACCG‑CATCGGCGGGTGGTTGCTTTTGGAGGGTAAGAAAGATGAATGCCTGTCCACGCCCCtatccg                                                                                                 >  1:1042829/1‑148 (MQ=255)

                                 ggTGAGCGATTGCTGGCCTAGATTAATGACTGGCCACCACAGAACCCGGCTTATCGGTCAGTTTCACCTTATTTACGTAAAA‑ACCCG‑CTGCGGGGGTTTTTTGCTTTTGGAGGGGCACAAAGATGAATGACTGTCCACGACGCTATACCCa                                                                                               <  1:617128/150‑1 (MQ=18)

                                 aGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGAGATCAAAGAAGAATCAGTGTCGGAAACACTATAccc                                                                                                >  2:597169/1‑150 (MQ=18)

                                 aGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACACAAAACCCAAAATAAATCACTTAAACCTGAATTAAGTAAAA‑AAACA‑CTTCGGCGAGTTTTTGCTATTGGAGGTGCAGAAAGATGAAAAAATGTCCACGACACTATAccc                                                                                                >  6:1159383/1‑150 (MQ=12)

                                    gtgCGATTGCTGGCGTATATTATTGACTGTCCCCGACAGAACCCGGGTTTGCGTTCCGTTTCACCTTATTTACTAAAAA‑ACCGG‑CTTCGGCGGTTTTTTGCTTTTGGAGGGGCACAAATATGAATGACTGTCCACGACGCTATACCCaaa                                                                                             <  5:1186265/148‑1 (MQ=17)

                                    gcgCGATTGCTGGCCTAGATGAATGACTGGCCACGACAGAACCCGGGTTATCGGGAAGTTTCACCTGATTTACTTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCACAAAGATCAATGACTGTCCACGACGCTATACCCaaaa                                                                                            <  5:761450/149‑1 (MQ=18)

                                    gAGCTATGGCGGGCCTAGTTTAATGACTGCCCACCACACAACCGGGCTTATCGGTAAGTTTCACCTGATTTAGTTAAAA‑ACCCG‑CTGCGGGGGGTTTTTGCTTTGGGCGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCaaa                                                                                             <  6:253813/150‑1 (MQ=17)

                                    gAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGAATCACCTCATTTACCTAAAA‑ACCAG‑CAGCGGCGGGTTTTTGCTGTTCGAGGGGCAGAAAGATGAATCACTGTCCCCGACCCTATCCCCcaa                                                                                             >  1:915799/1‑150 (MQ=21)

                                     cgCGATTGCGGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCC‑CTTCGGGGGGTTTTTGCTTTTGGCGGGGCACAAACATGAATGACTGCCCACGACGCTATACCCAaaag                                                                                           <  3:932621/150‑1 (MQ=18)

                                      gCGATTGCTGGCCTAGATGACTGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGatga                                                                                                                        <  7:860307/121‑1 (MQ=255)

                                        gATTGCTGGCCTAGATGAATGACTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaagaaa                                                                                        <  5:420759/151‑1 (MQ=255)

                                         atgggtGGCGTAGTTGAATGAGTGTCCACGACTGAACCGGGGTTTTCGGTCTGTGTCCCGTTATTTCCGTAAAA‑CCCGG‑CTGCGGGGGGTTTTTGGTTTTGGCGGGGCAGAAAGTTGATTGAGTGTCCACGACGCTATACCCAaaagaaa                                                                                        <  5:1032891/145‑1 (MQ=255)

                                         aTTGCTGGCCTAGATGAATGACTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaagaaag                                                                                       <  5:629999/151‑1 (MQ=255)

                                            tcTGGCCAACAGGAATGACTCGCCACCAAAGAACCCGCCTTATCGTTCATTTTCACCTGATTTACTTAAAA‑CCCCG‑CTTGGGCGGGTTTTTGCTTTTGGCGGGCCACAAACATCAATGACTGTCCACGACGCTATACCCAAAAGAAAGCg                                                                                     <  6:348333/149‑1 (MQ=255)

                                            gCTGGCCTATATTAATGACTGACCACGACATAACCGGGCTTTTGGTTCTTTTGAACTTTTTTTACGTAAAA‑ACCCG‑GTTCGGCGGGTTTTTTCTTTTGTCGGGGCAAAAATATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCg                                                                                     <  1:405218/150‑1 (MQ=255)

                                             cTGGCCTAGATGAATGACTGTCCACGACAGAAACCCCCTTATAGATCAGTTTCACCTCATTTACGTAAAA‑ACCCT‑CTTAGGATGGTTTTTACTATTGGAGGGGCAGAAAGATTAATAAATGTCCACGACCCTATACCCAAAAGAAATCgg                                                                                    >  8:175192/1‑150 (MQ=17)

                                               ggCCTATATGAATGACTGTCCACGACATAACCCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAA‑ACCCG‑GTTCGGCGGTTTTTTGCTTTTGGAGGGGCACAAAGATTAATCACTGTCCACGACGCTATACCCAAAAGAAAGCGGCtt                                                                                 <  3:910491/151‑1 (MQ=17)

                                               ggCCTAGTTGAATGCCTGCCCACGACAGAACCCGGCTTATCGGGCAGTTTCACCTTATTTACGTAAAA‑CCCCG‑CTTCGGGGGTTTTTTGGTTTTGGAGGGGCACAAAGATCAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCtt                                                                                 <  5:284073/151‑1 (MQ=17)

                                               ggCCTAGATCAATGCCTGTCCACGACAGAACCCGGCTTATCGCTCAGTTTCACCTGATTTAGTTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGTAGGGGCATAAAGATTAATTACTGTCCACGACGCTATACCCAAAAGAAAGCGGCtt                                                                                 <  2:728853/151‑1 (MQ=17)

                                                 ccTAGATGAATGACTGTCCACGACACAAAACCCAAAAACATACAAATACACAGGTTTTACGTAAAA‑AACAG‑CATCGAAGGGTATATCAGTTTGGAGGGGCAAAAAGATAAAAAAATGTCCACGACACTATACCCaaaaagatgcggagacg                                                                               >  4:765807/1‑138 (MQ=2)

                                                 ccTAGATGAAGGACTGCCCACCACAGAACCCGGCTTAGCGGTCAGTTTCCCCTGATTTACGTAAAA‑ACCCG‑CTTGGGCGGTTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTAt                                                                               <  3:1185612/151‑1 (MQ=17)

                                                  cTAGATGAATGACTGTCCACGACAGAAACCCGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACAAC‑ATTCGGCTGGTTTTTGCTTTTGGAGGTGCAGAAAGATGAATGACTGTCCACAACTCTCTACACAAAAGAAAACGGATTAt                                                                               >  8:868929/1‑150 (MQ=17)

                                                    aGATTAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑CCCCG‑CTTCGGCGGTTTTTTGCTTTTGTAGGGGCAGAAAGATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                              <  6:562926/149‑1 (MQ=255)

                                                    aGATGAATGAGTGGCCACGACAGAACCCGGCTTATCGTTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTGGGCGGGTTTTTGCTTTTGGTGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                              <  4:447958/149‑1 (MQ=255)

                                                    aGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑CCCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                              <  3:561038/149‑1 (MQ=255)

                                                    aGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCgg                                                                            <  7:814773/151‑1 (MQ=255)

                                                    aGATGAATGACTGCCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACTTAAAA‑ACCCG‑CTTCGGGGGGTTTTTGCTTTTGTAGGGGAAGAATGATGATTTACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCg                                                                             <  2:509465/150‑1 (MQ=17)

                                                       tgaatgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGATTAACCTGATTTACTTAAAA‑ACCCG‑CTTCGCAGGGTTTTTACTTTTGGAGGGGCAGAAAGATGAATGACTGTAAAAGACGCTATACCCAAAAGAAAGCGGCTTATCgg                                                                            >  4:662392/1‑148 (MQ=17)

                                                           ggccgctCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACTTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGTTCAAGGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                     <  1:1004162/145‑1 (MQ=14)

                                                            gtgaggcgccaaCAGAACCCGGATTATGGGCAAGTTTCCCGTGATTTAGTTAACA‑CCGGG‑CTGAGGCGGGTTTTTATTTTTGCCGGCGCCTAAACATGAAGGACTGGCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                     <  3:655963/139‑1 (MQ=2)

                                                            gaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑CCCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTc                                                                    <  7:431198/151‑1 (MQ=255)

                                                            gaCGGTCCACGACAGAACCCGGCTTATGGGTCATTTTCACCTGATTTACTTAAAA‑ACCGG‑CTGCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                     <  5:669825/150‑1 (MQ=14)

                                                             aCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCAGCtccgaggacccctcc                                                                    >  7:172183/1‑136 (MQ=17)

                                                             aCTGCCCCCGACAGAACCCGGGTTATCGGTCATTTTCACCTGATTTACGAAAAA‑ACCGG‑CTTCGGCGGTTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCa                                                                   <  3:599175/151‑1 (MQ=14)

                                                              gcatgctaaaCAGAACGCGGATTAT‑GGGCAATTTGGCCTTCTTTACGTAATA‑CCCGGTCTGATGCGGGTTTTTGTGTTTGTTGGTGTTTAATTATGATTGACTTTCCACGACACTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCAc                                                                  <  4:615484/142‑1 (MQ=2)

                                                               tGTCCACGAAAGAACCCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTTCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCa                                                                                                             <  7:778312/107‑1 (MQ=255)

                                                               gggccACGACAGAACCCGGATTATGGGTCATTTTCACCTGATTTAGTTAAAA‑ACCGG‑CTGCGGCGGGTTTTTGTTTTTGGAGGGGCAGAAAGATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCAcc                                                                 <  3:397070/148‑1 (MQ=14)

                                                                cgccaccaCACAACCCGGCTTATGGGGCAGTTGCACGTCATTTACTTAAAA‑CCGGG‑CTGCGGCGGGTTTTTGGTTTTGTTGGGGCGCAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCt                                                                <  3:1005263/149‑1 (MQ=14)

                                                                caccaccaCAGAACCCGGCTTATCGGTAAGTTACACCTGATTTACGTAAAA‑CCCGG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCCCAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCt                                                                <  7:1189660/149‑1 (MQ=14)

                                                                   caccacaccaCCCGGATTATGGCTCAGTTTCACGTCATTAACTTAAAA‑ACCCG‑CTGCGGCGGGTTTTTGCTTTTGTTGGGGCTCAAAGATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGt                                                             <  4:1073426/142‑1 (MQ=11)

                                                                       aCATAACCCGGCTTATCGCCCAGTTTCCCCTGATTAAGTTAAAA‑CCCGG‑CTTCGGCGGGTTTTTGCTTTTGTTGGCGCACAAAGATGATTGACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTAc                                                         <  8:381561/151‑1 (MQ=14)

                                                                         ccaaCCCGGCTGCTCGGCCAGTTGCACCTCATTCCCGTAACC‑CCCGC‑CTGCGGCGGGTTTTTGGGTTTGCCGGGGCCCAACCCTGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGt                                                       <  4:688931/149‑1 (MQ=9)

                                                                          caaCCCGGCTTATCGGTCAGTTGCACCTGATTTACGTAAAA‑CCCGG‑CTTCGGCGGGTTTTTGCTTTTGTAGGCGCACAATCATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTa                                                      <  5:118985/150‑1 (MQ=14)

                                                                             cccGGCGTATGGGGAAGTTGCACGTGATTTACTTAAAC‑ACCCG‑CTGCGGCGGGTTTTTGGTTTTGGCGGGGCCCAACTATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaa                                                   <  5:551785/151‑1 (MQ=11)

                                                                              ccgggttttggggccgTTTCCCCTGATTTACTTAACA‑CCCGC‑CTGCGGGGGGTTTTTGTTTTTGTTGGGGCTAAAATTTTAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaaa                                                  <  5:237834/136‑1 (MQ=2)

                                                                              ccGGCTTATCGGCAATTTTCCCCTGATTAACTAAAAC‑CCCCG‑CTTCGGCGGGTTTTTGGTTTTGCGGCGGCACAACTATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaaa                                                  <  1:210435/151‑1 (MQ=11)

                                                                              ccGGATTATGGGCAATTTTCCCCTCATTT‑CGTAAAACCCCCG‑CTGCTGCGTGTTTTTGTTTTTGCGGGGGCATAATGATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaaa                                                  <  1:848450/151‑1 (MQ=2)

                                                                              ccGGATCATGGGGCAGTTTCCCCTGATTAACGTAAAA‑CCCCG‑CTGCCGCGGGTTTTTGCTTTTGGCGGGGCAGAAACATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaaa                                                  <  1:1120785/151‑1 (MQ=11)

                                                                                  cTTATCGGTCAGTTTCACCTGAATAACGTAAAA‑ACCAC‑CTTGGGCGGGATTTTGATTATCGATGGGGAGATATATGAATGACTCTCAACGAGAAAAAAAAAAAAAGAAAGCGGCTTATCGGTCAGTGTCGCCTGGtgggcggaagagcggg                                              >  8:1258918/1‑136 (MQ=2)

                                                                                     atgggcaattttCACCTGATTTACTTAACA‑CCCGT‑CTGCGCCGGGTTTTTGGTTTTGGTGGCGCCTAACCATGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCt                                            <  6:787762/141‑1 (MQ=11)

                                                                                     aTCGGGCAGTTGCAGCTCATTTAGGTAAAC‑CCCGG‑CTTCGGCGGGTTTTTGGTTTTGGCGGCGCAGAACGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCtt                                           <  3:462632/151‑1 (MQ=11)

                                                                                       cgggaaGTTTCACCTGATTTACGTAAAA‑ACCGC‑CTTCGGCGGGTTTTTGCTTTTGTCGGGGCACAACCATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTc                                          <  3:505323/145‑1 (MQ=11)

                                                                                       cGGGCATTTGCACGTGTTTTCGGTAAAA‑ACCCG‑ATGCGGCGGGTTTTTGTTTTTGGTGGGGCAGAAAGTTGAATGCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGGCAGTTTCACCTGGTTTACGTAAAAACCCGCTTc                                          <  1:859528/150‑1 (MQ=11)

                                                                                        ggTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCAGGACGGTAGAGCCCAAGTCAGACGTCCTCTCCCTAATGCTCAGCTGGTTTACGTAAAAACCCGGttat                                         >  1:355565/1‑148 (MQ=2)

                                                                                        ggTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGATTTTTGAGTGGCATAAAGATCAATTACTGTACAAGAAGCTATACACAAAAGAAAGATGATTATAGGTCAGTTTTACCTGATTTAATTaaaaaaacgatccgg                                        >  6:964300/1‑141 (MQ=2)

                                                                                              tttCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGATATTTACTTTTTGATGTTGAGAAAGATGATTGCATTTAGACTACAATTAACAGGAAAGAATAGGGATTATATTTAAGTTTCAGGTGGTTTACGTATAAGCCGGCTTCGGGGGGtg                                  >  2:58982/1‑150 (MQ=2)

                                                                                               ttcccgtGATTTCCTTAATA‑ACCGC‑CTGCTGCGGGTTTTTGTTTTTGTCGGTGGCGACAGTTGATTGCCTTGCCACGACGCTATACACAAAAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGtt                                  <  6:739144/144‑1 (MQ=2)

                                                                                               ttCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGGTGTACGTAATAAACCGGTGCGGCGGGttt                                 >  5:451657/1‑151 (MQ=11)

                                                                                                gcACCTGATTAACGAAAAA‑ACCGC‑CTGCGGGGGGTTTTTGGTTTTGCCGCGGCACAACCATGAATGACTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGtttt                                <  5:1096523/150‑1 (MQ=11)

                                                                                                gcACCGGATTAACTTAAAC‑CCCGG‑CTGCGGCGGGTTTTTGGTTTTGTCGGCGCCCAAACATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttt                                 <  3:405060/149‑1 (MQ=11)

                                                                                                      caTTAACGTAAAA‑ACCCG‑CTGCGGCGGGTTTTTGGTTTTGCCGGCGCCCAACCATGAAGGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCt                            <  3:840058/148‑1 (MQ=11)

                                                                                                          accttAAAA‑CCCCG‑CTGCGGCGGGTTTTTGTTTTTGTTGCGGATGAATTTTGAATGACTTGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGa                      <  2:1009197/147‑1 (MQ=11)

                                                                                                          accgaaaaa‑cCCCG‑CTGCGGGGGGTTTTTGGTTTTGCCGGGGCACAACCATGACGGCCTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGa                      <  7:1132964/146‑1 (MQ=11)

                                                                                                           aCGTAAAA‑CCCCG‑CTGCGGGGGTTTTTTGCTTTTGCGGGCGCTAAACCTTGATGGCCTGCCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGa                      <  7:287930/150‑1 (MQ=11)

                                                                                                                  caatccttcttatgcggGTTTTTGCTTTTGCACGCGATGACTTATGACTGAGTGTCCTTCACGCTTTACACTAAAGAACGCGGCTTATCGTTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGc                 <  2:39300/137‑1 (MQ=2)

                                                                                                                    aaccgtcTTCTGCGGGTTTTTGTTTTTGTTCGTGCAATATATTTATTGACTTTACTCTACTCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa              <  4:1140469/145‑1 (MQ=2)

                                                                                                                      ccgcctgcgccggGTTTTTGTTTTTGCCCGGGCCAAACCATGAATGCCTCGCCACGACGCTATACCTAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaaa            <  7:990618/141‑1 (MQ=11)

                                                                                                                         tctgccgcggGTTTTTGTTTTTGCTGGTGATAAACTATGCTTGACTTGCCAAGACTCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAg           <  2:385468/143‑1 (MQ=11)

                                                                                                                           ttCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTTACCTGTTTTACGTAAAAACCCGGTACGGCGGGTTTGTGCTTGTGGAGTGGgagatatagct       >  3:489760/1‑144 (MQ=11)

                                                                                                                            gcgggggGTTTTTGCTTTTGGAGCGGCACAACTATGAATGACTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaa      <  7:338622/150‑1 (MQ=11)

                                                                                                                            gcggcggGTTTTTGTTTTTGTCGGCGTCCAATATTGAATGCCTATACAACACTCTATACCCTAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatga       <  7:1074543/149‑1 (MQ=11)

                                                                                                                             cggcggGTTTTTGGTTTTGTCGCTGCCGAAATATGAATGACTTTACACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat     <  4:967719/151‑1 (MQ=11)

                                                                                                                              ggcggGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCATGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatg    <  1:411970/151‑1 (MQ=25)

                                                                                                                                cggGTTTTTGGTTTTGGTGGGGGACACAAATGATTGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGAc  <  3:20856/151‑1 (MQ=11)

GCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAA‑ACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGAC  >  DH1/613447‑613723

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