BRESEQ :: Evidence

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

Only 100 of 205 total aligned reads displayed.

GCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACC  >  DH1/613447‑613741

gCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAg                                                                                                                                                  <  4:1107372/151‑1 (MQ=255)

   cGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                                                >  1:506077/1‑150 (MQ=255)

     tACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAACCCCCTTATCAATCAGATTCACCAGATATAAATAAAAACAAG‑ATTCGACAGGTTATAAAATTTGGAGGAGCAGaa                                                                                                                                                >  6:1348031/1‑148 (MQ=255)

      aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGat                                                                                                                                            <  5:1202674/151‑1 (MQ=255)

      aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGa                                                                                                                                             <  8:1357497/150‑1 (MQ=255)

      aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGa                                                                                                                                             <  8:1101010/150‑1 (MQ=255)

        tGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACATGATTTACGTAAAAACCAA‑CTTATTCGTGTTTTTTCTTTTTGAGGGGCAGAAAGatg                                                                                                                                           >  4:678245/1‑150 (MQ=37)

         gAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaa                                                                                                                                         <  7:862930/151‑1 (MQ=255)

            gcgggttattgtgctTGATCCAGAGATAGATTGATTGACTAGATCACTCCCTATCCACTACAGAACGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                                                                                                                                        <  2:505980/138‑1 (MQ=11)

            cccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGattcatcg                                                                                                                                      >  5:1107112/1‑145 (MQ=255)

             ccGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGAc                                                                                                                                     <  5:1241447/151‑1 (MQ=255)

              cGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACt                                                                                                                                    <  3:1314465/151‑1 (MQ=255)

              cGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTAACCTTATTTACTAAAAAAACCG‑CTTCGGCTGCTTTTTGCTATTGTATGTGAAGAAATATGAATGACt                                                                                                                                    >  4:20770/1‑151 (MQ=37)

               gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                                                                                                                                   <  7:92378/151‑1 (MQ=255)

               gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                                                                                                                                   <  2:266571/151‑1 (MQ=255)

               gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTACGGCGGGTTTTTGCTTTTGGAGGGGCAGTAAGatttgaggctt                                                                                                                                   >  1:1260197/1‑142 (MQ=255)

               gggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAAACCG‑CTTCGGCGGGTTTTTGCTTGTGGAGGGGCAGTAAGaggtgagccat                                                                                                                                   >  5:238267/1‑141 (MQ=25)

                ggtattgtgtttGAGGCAGTCACAGATTGCTTGACTTGAGAAGTGACTCTACACCAAAGAACCCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGATGAATGACTGt                                                                                                                                  <  7:1263195/141‑1 (MQ=2)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTcc                                                                                                                                <  6:1048456/150‑1 (MQ=255)

                   aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTcc                                                                                                                                <  5:924524/150‑1 (MQ=255)

                   aGGCGGCGTGAGCCAGTTACCGAATGCTGGCCTAGATGAATGCCTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCa                                                                                                                               <  1:938183/151‑1 (MQ=18)

                    ggCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGCCTGTCCac                                                                                                                              <  5:592574/151‑1 (MQ=255)

                    ggCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCac                                                                                                                              <  6:193425/151‑1 (MQ=255)

                     gctgctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacg                                                                                                                             <  7:458442/151‑1 (MQ=255)

                          ttGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaataaatctgaacaacga                                                                                                                         >  6:665621/1‑135 (MQ=255)

                          ttGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGatga                                                                                                                                          >  5:881465/1‑133 (MQ=255)

                           tgcgCCAGTGCGCGATGGCTGGCCTAGATCAATGACTGCCCACCACACAACCCGGCTTATCGGCCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCt                                                                                                                        <  4:187643/147‑1 (MQ=18)

                              gCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATAc                                                                                                                    <  5:203547/151‑1 (MQ=255)

                               ccAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTTAGTAAAAAACCG‑CTTCGGAGGGGTATAGTTTTTGGATTGGATGAAGTATGACTGAGTGTGCAGGACAATATAc                                                                                                                    >  3:116564/1‑150 (MQ=19)

                                    gAGCGATTGCTGCCCTAGATCAATGACTGACCACGACAGAACCCGGCTTATCGGGCAGTTTCCCCTTATTTACGTAAAAACCCG‑CTTCGGGGGTTTTTTGCTTTTGGAGGGGCAGAAAGATCACTGACTGTCCACGACGCTATACCCaaa                                                                                                               <  3:287796/150‑1 (MQ=18)

                                      gctaTTGCTGGCCTAGATTAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACGCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaag                                                                                                             <  6:438177/147‑1 (MQ=255)

                                       cGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGGTTGAGCTAAAAACACG‑AGTCGGCTCCAGTGTGCAGTTTGCGGGGCAGATGCATGAATGACTGTCCAAGAAACTATACCCAaaaga                                                                                                            >  7:1024116/1‑150 (MQ=12)

                                          ttGCTGGCCTAGATGAATGACTGTCCACGACACAAACCCGCTTATCACTCAGATTCACCTGATTAACGTAAAAACCCG‑CTACCGAGTGTTATTCATTTTAGAGGTCCAGAAAGATGAATTAATGTAAAAAAAAATATACCCAAAAGAAAGc                                                                                                        >  8:1237001/1‑151 (MQ=12)

                                           tGCTGGCCTAGTTCAATGACTCTCCACGACACAACCCGGCTTAGCGGTCAGTTTCACCTGATTTACGTAAAACCCGG‑CTGCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGc                                                                                                        <  2:698379/150‑1 (MQ=17)

                                            gggggCCTTGATTAATGACTGTCCACCACACAACCCGGCTTATGGGTCAGTTTCACCTTATTTACTTAAAACCCGG‑CTGGGGGTGTTTTTTGTTTTTGGTGGGGCAGAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCg                                                                                                       <  3:1359458/147‑1 (MQ=17)

                                            gCTGGCCTAGATTACTGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACTTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATCAATGACTGTCCACGACGCTATACCCAAAAGAAAGCg                                                                                                       <  2:1323366/150‑1 (MQ=255)

                                            gCTGGCCTACATCAATCACTGTCCACCACAGAACCCGGCTTATCGGTCAGTTTCACCTCATTTACGTAAAACCCCG‑CTGCGGCGGGTTTTTGCTTTTGGAGGGGCACAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCgg                                                                                                      <  4:426770/151‑1 (MQ=17)

                                             cTGGCCTAGATGAATGACTGTCCACGACAGAACACGGATTATCGGTCAGTTTAACCTGATTTAAGTAAAAACCCT‑CTTCGCAGGGGTTTTGATTTTGGATGGGCAGAAAGATGAAAGCCTGTCAAAGACACTATACCCAAAAGAACGCggg                                                                                                     >  5:821736/1‑150 (MQ=17)

                                              gggCCTAGATTATTGACTGTCCACGACATAACCGTGCTTTTCGGTCAGTTTCACGTTATTTACGTAAAACCCCG‑CTTGGGCGTTTTTTTGTTTTTGGAGGGGCATAAAGATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCt                                                                                                    <  8:611116/150‑1 (MQ=17)

                                               ggCATAGATGAAGGACTGTCCACGACACAACCCGCCTTATCGGTCAGTTTCACATGATTTACGTAAAAAGCCG‑CTTCGGGGGTTTTTTGATTTTTTAGGGGCAGAATTATTAATTACTGTCCACGACGCTATACCCAAAAGAAAGCGGc                                                                                                     <  4:80731/149‑1 (MQ=17)

                                                 ccTAGATGAATCACTGCCCCCCACAGAACCCGGCTTAGCGGTCAGTTTCACCTGATTTACGTAAAACCCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTAt                                                                                                 <  4:871155/151‑1 (MQ=17)

                                                  ggtctcccgacggtttCCACAACAGAACCGGGCTTATCGGTCAGTTTCACCTTTTTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACCCAAAAGAAAGCGGCTTAt                                                                                                 <  3:493603/135‑1 (MQ=1)

                                                    agattattgccggTCCACGCCAGAACCCGGGTTATCGGTCTGTTTCACCTGATTTACTTAAAAACCCG‑GTTCGGCGGGTTTTTGTTTTTGGTGGGGCAGAAAGATTAATTACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                                                <  4:1269265/137‑1 (MQ=14)

                                                    aGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAACCCCG‑CTTCGGCGGTTTTTTGCTTTTGGAGGGGCAGAAAGATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                                                <  1:439148/149‑1 (MQ=255)

                                                    aGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCAG‑CTGCGGACGGTTTTTGCGTTTGGAGGTGCAGAAAGATGAATGAATGTCAACGACACTATACCCAAAAGCCAGCGGCTTATggg                                                                                              >  3:775443/1‑151 (MQ=17)

                                                      cctcatcgcTGTCCCCTACAGAACGCGGATTATGGGTAAGTTGCCCGTGATTAAGTTAAACCCGGG‑CTGCGGCGGGTTTTTGGTTTTGTTGGGGCTCAACTATGAATGACTGCCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTc                                                                                            <  3:1291242/143‑1 (MQ=2)

                                                      atgaatgaCTGGCCACCAAAGAACCCGGATTATCGGTCATTTTCACCTTATTTACTTAAAAACCCG‑CTTCGGGGGGTTTTTGCTTTTGGCGGGGCATAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTc                                                                                            <  7:217022/151‑1 (MQ=14)

                                                       tgaatgcCTGGCCACGACAGAACCCGGCTTATCGGGCAGTTTCACCTGATTTACGTAAAACCCCG‑CTTCGGCGGGTTTTTGCTTTTGGTGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTc                                                                                            <  4:1334332/150‑1 (MQ=255)

                                                       tgaatgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCGCCTGATATACGTAAAAACCCG‑CTTCGGCCGGTTTTTGCTTTTGTGGGGGCAGAAGGATGAATGACTGTCCACGACCCTCTCCCCAAAAGAACGCGGCtgggcggcgc                                                                                           >  3:930266/1‑142 (MQ=17)

                                                        gaatgaCTGGCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAg                                                                                          <  1:794132/151‑1 (MQ=255)

                                                        gaatcaatttCCTCGACAGAACCCGGCTTTGCGGTCTTTTTTACGTTATTTACTTAAACACGCG‑CTGCGGCGGGTTTTTGCTTTTGGTGGGGGACAATTTTGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                                           <  2:122813/141‑1 (MQ=9)

                                                        cctctcctttcccccaCAGAACCCGGCTTATCGGTCAGTTTTTTTTTTTTTATTTAAAATCCCT‑CTTCGGTGTTTTTTTGCTTTTTTAGGGTTTTAATTATGAAGGATTTTCCACTACACTATACCCAAAAGAAAGCGGCTTATCGGTCAg                                                                                          <  2:1025720/136‑1 (MQ=2)

                                                        cccccccaaaccccaacacAACCCGCCTAATCGGCCAGTTTAACTTTTTTTACTTATACGCCGC‑CTTCGGCGTTTTTTTACTTTTTTGCGTCCATTTACATAAACGACGCTCCACGAATCTATACCCAAAAGAAAGCGGCTTATCGGTCAg                                                                                          <  8:27453/136‑1 (MQ=2)

                                                         aatgaCTGGCCACGACAGAACCCGGCTTTTGGGTCATTTTCACCTGATTTACGTAAAAACCGG‑CTGCGGCGGGTTTTTGCTTTTGTAGGGGCACAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAg                                                                                          <  4:637621/150‑1 (MQ=17)

                                                          cccccaaccccccaCAGAACCCGGCTTATCGGTCAGTTTCACCTTTTTTTTTAAATATCCGG‑CTTCGTCGGTTTTTTACTTTTGGTTTGCCACAATTTTCATTCACTTGGCTCTACAGTTTACACTAAAGAAAGCGGCTTATCGGTCAGtt                                                                                        <  6:601997/138‑1 (MQ=2)

                                                          atgaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCAAAATAAAGTAAAAACCCG‑AGTAGGAGGGATATTGATTATGGAGCGGAAGAAACATTAAAGACTGAACAAAAAGCTATACACAAAAGAAAGCGGCTGGTCGGTCAggg                                                                                        >  4:238794/1‑149 (MQ=255)

                                                            gaCTGGCCACCACAGAACCCGGCTTATCGGTCTGTTTCACTTGATTTACTTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGTGGGGCAGAAATATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                                       <  8:1299959/150‑1 (MQ=14)

                                                             acggcccaccaCAGACCCCGGCTTAGCGGGCAGTTTCACCTGATTTACTTAAAACCCCC‑CTGCGGCGGGTTTTTGCTTTTGGAGGGGCACAACGATCAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCa                                                                                     <  2:206904/146‑1 (MQ=14)

                                                             aCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CGTCGCCGGGTTTTTGCTTTTGGATGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAATAGAAAGCGGCTTATCGGTCAGtgagc                                                                                     >  1:1025588/1‑147 (MQ=255)

                                                               tGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTCCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACACAAAAGAAAGCGGCTTATAGATCAgaataacc                                                                                   >  4:493636/1‑144 (MQ=14)

                                                                 tacACGACAGAACCCGGGTATTCGGCACGTTACAGCTGTTTTACGTAAAACGCGT‑CGTCGGCGGGTTTTTTTTTTTGTTGGCGCTCAAACATGAATGACTGGCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTTACCt                                                                                  <  3:20797/148‑1 (MQ=2)

                                                                 taaaCGACATAACCCGGCTTATCTGTAAGTTTCAGTTTATTTAGTTAAAACCGGG‑CTTCGGCGGGTTTTTGTTTTTGTTGGTGCTTAATTATGATTGACTTGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCt                                                                                  <  8:334831/147‑1 (MQ=9)

                                                                  ccACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTc                                                                                                                                 >  7:1250244/1‑102 (MQ=255)

                                                                   cagtactcaaCCCGGCTTATCGGTCTTTTTCACCTGTTTTCGTTAACAAGCTT‑CTGCGGGGGGTTTTTGTTTTTGTTGGGGGTTTATTTTGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGt                                                                               <  6:54101/143‑1 (MQ=2)

                                                                      gACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCAG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATTACTGTAAAAGACACTATACCCAAAAGAAAGCGGCTTATCGGTCaagatcggaagaggtt                                                                            >  8:692309/1‑136 (MQ=11)

                                                                      cacatcaCCCGGATTTTCGGACAGTTTAACGTCTTTTACGTAAAACCGCT‑CTTGGGCGGTTTTTTGGTTTTGTAGGGGCTCAAAGATGATTGACTGTCCACGTCGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTa                                                                            <  5:162580/145‑1 (MQ=2)

                                                                      cacacAACCCGGCTTATCGGTCAGTTTCACCTGATTTACTTAAACACCCG‑CTGCGGCGGGTTTTTGTTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTa                                                                            <  3:76949/150‑1 (MQ=255)

                                                                         aGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTACTTTTGGAGGGGCAGa                                                                                                                                                 >  6:411990/1‑79 (MQ=255)

                                                                                    tATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCCGCGCGTGTTTCCTGGTGTAGCGGCAGAAAGATGAATCACTCTCCACGACCGTATAAAAAAAAGAAAGCGGCGTATCGGGCAGGTTGCCCTGGGGTACGTaaggtcgggcg                                                              >  1:580274/1‑142 (MQ=2)

                                                                                      acGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTAAGGCGGGTTTATGCAGCTGGAGCGGCAGAAACAAGAATGACTGTCCACGACCCTATACACAAAAGAAAGCCCATTATCGGTCAGTTTTACCTGTTGTACGAAATAACCCggtcc                                                            >  1:550543/2‑147 (MQ=255)

                                                                                       gggTCATTTTCACGTGATTTAGTTAAACACGCT‑CTTCTGCGGGTTTTTGTTTTTGTTGGTGCTTAAATATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcg                                                           <  4:567387/150‑1 (MQ=11)

                                                                                       cGGTAATTTTCCCCGCATTTAGGTAAAAACCCT‑CTTCGGCGGGTTTTTGTTTTTGGTGGGGCCCAACGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcg                                                           <  4:520071/151‑1 (MQ=11)

                                                                                          tCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAACGACGATGAGAGGTCAGTTCGACCTTCAATACGTAAAAACCCGCTTcgccg                                                        >  3:426934/1‑151 (MQ=11)

                                                                                           aatttgcACGTGATTTACTTAAAACCCGG‑CTTCGGCGGGTTTTTGCTTTTGGCGGCGCACAAAGATGAATGCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCGGTTTCACCTGGTTTACGTAAAAACCCGCTTcggcgg                                                       <  3:33427/145‑1 (MQ=11)

                                                                                            aGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTTACCTgatgt                                                                             >  1:860706/1‑124 (MQ=9)

                                                                                            aGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGGTGTACGTAATAACCCGCGTGgggggg                                                      >  3:1117532/1‑151 (MQ=11)

                                                                                            aGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAGAACCGCTGCgggggg                                                      >  7:584209/1‑151 (MQ=255)

                                                                                               ttcccgtGATTTCCGAAAACCCGCG‑CTGCGGCGGGTTTTTGCTTTTGGCGGGGCCCAACCATGAATGCCTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGtt                                                    <  1:192133/144‑1 (MQ=11)

                                                                                                tctcttGATTTACTTAAAACCGGT‑CTTCGGGGGGTTTTTGTTTTTGTTGGTGCTTAAATTTGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttt                                                   <  8:142166/145‑1 (MQ=11)

                                                                                                tCACCTGATTTACGTAAAAACCGG‑CGGCGGCGGGTTTTTGGTTTTGGTGGGGCCCAAAGATTATTGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttt                                                   <  5:536734/150‑1 (MQ=11)

                                                                                                 cACCTGATTAACGAAACAACCCC‑CTGCGCCGGGTTTTTGGTTTTGCCGGCGCCCAACCATCCATCCCTCCCCACCACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttttt                                                 <  3:961055/151‑1 (MQ=2)

                                                                                                 cACCTCATTTACGAAAAAACCCC‑CTTCGGCGGGTTTTTGCTTTTGGCGGCGCACAAACATCAAGCCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGttttt                                                 <  5:150640/151‑1 (MQ=11)

                                                                                                    cTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCgg                                                                                                      >  5:658312/1‑95 (MQ=35)

                                                                                                      gATTTTCTTAATAACCCG‑CTGCTGCGGTTTTTTGTTTTTGGCGGCGGCTTACGATGCAGCCCTCTACCCGACGCTATACCCAAAACAACCCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCt                                              <  5:1296411/149‑1 (MQ=2)

                                                                                                      gATTTACTTAAAACCCTG‑CTGCTGCGGGTTTTTGTTTTTGTTGGTGCTATAATATGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtt                                             <  8:480414/150‑1 (MQ=11)

                                                                                                      gATTAACGAAAACCCCCC‑CTGCGGCGGGTTTTTGCTTTTGGCGGCGCCCAACCATCAATGCCTGCCCCCGACGCTATACCCAAAAGAAACCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCttt                                            <  5:605812/151‑1 (MQ=11)

                                                                                                          tccGTAAACCCCCG‑CTGCGGCGGGTTTTTGGTTTTGGGCGGGCAGAACCATGAAGGCCTGGCCACGACGCTATACCCAAAAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTgg                                         <  7:351595/148‑1 (MQ=11)

                                                                                                          tACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAAACCGATAcggcgagatattaatatagt                                         >  2:986137/1‑135 (MQ=11)

                                                                                                           acttAAAACCGCG‑CTTCGGCGGGTTTTTGTTTTTGGAGGGGCAGAAATTTGATTGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAg                                       <  6:348378/148‑1 (MQ=11)

                                                                                                           aCGAAAAACCCCG‑CTGCGGCGGGTTTTTGCTTTTGCCGGGGCAGAAACATGAATGCCTGGCCACGACGCTATACCCAAAAGAACGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGa                                        <  5:681222/150‑1 (MQ=11)

                                                                                                            gtcaaaaaCCGG‑CTTCGGCGGGTTTTTGGTTTTGCCGGGGCCGAACCATGAAGGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAgg                                      <  3:584703/148‑1 (MQ=11)

                                                                                                             gTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATAGGTCTGTTTCGCCGGGGTATCGTATAAACGCGCGTCGGAGGGttttggatgttgtcggg                                     >  7:1101203/1‑138 (MQ=11)

                                                                                                              taaaaccggg‑CTGCGGCGGTTTTTGGGTTTTGTAGGCGCTGAAATTTTATTGCCTTTCCACGACGCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAgg                                      <  8:1060297/140‑1 (MQ=11)

                                                                                                               aaaaaCCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAACCCGATGCGGCGGGTTTGTGATTGTGGAgg                                      >  8:253767/1‑148 (MQ=11)

                                                                                                                acaaCCCGCCTGCGGCGGGTTTTGGGTTTTGTCGCGGCTTAATCATGATGGACTTGCCACGACTCTATACCCCAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGc                                   <  6:743448/149‑1 (MQ=11)

                                                                                                                     ccG‑CTGAGGCGTGTTTTGGTTTTTGCCGGCGCCCAAACTTGAATGACTGTCCACCACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAg                             <  5:1026943/151‑1 (MQ=11)

                                                                                                                           tCGGCGTGTTTTTGTTTTTGCCGGGGCCTAACCATCAATGACTAGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatga                         <  3:640681/150‑1 (MQ=11)

                                                                                                                           tCGGCGGGTTTTTGTTTTTGCGGGCGCCTAATCTTGAATGACTTGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatga                         <  8:745409/150‑1 (MQ=11)

                                                                                                                            cggcggGTTTTTGCTTTTGCACCGGCACTACGATGATTGCCTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaat                       <  3:529907/151‑1 (MQ=11)

                                                                                                                               gggttttttgtttttgGGGGGGTATAAGCTTGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGAc                    <  7:700088/147‑1 (MQ=11)

GCCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCG‑CTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATACC  >  DH1/613447‑613741

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