BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	3,448,523	0	A	C	22.4%	955.7 / 71.2	430	intergenic (‑83/+73)	malZ/proY	alpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand): ref base A (246/83); new base C (0/95); total (246/184)
Fisher's exact test for biased strand distribution p-value = 4.20e-45
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 452 total aligned reads displayed.

ATGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGA  >  DH1/3448373‑3448672

aTGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACTGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggagga                                                                                                                                                       <  1:40009/151‑1 (MQ=255)

  gAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagg                                                                                                                                                     >  3:469962/1‑151 (MQ=255)

  gAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggag                                                                                                                                                      <  4:1196379/150‑1 (MQ=255)

  gAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggag                                                                                                                                                      <  8:291263/150‑1 (MQ=255)

   aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagg                                                                                                                                                     >  8:1708455/1‑150 (MQ=255)

   aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGAAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                    <  5:721393/151‑1 (MQ=255)

    gCAGTTGATCTTTGCTTTGTTTATCAAATGGGGGCACCGGCAGGTGGCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCTCTCATCATCTGAAGGATATTTGGGGGAGGAGGAc                                                                                                                                                   <  5:249320/151‑1 (MQ=255)

    gCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAACGATATTTGGGGGAGGAGGAc                                                                                                                                                   >  6:1146730/1‑151 (MQ=255)

     cAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACg                                                                                                                                                  <  8:1510441/151‑1 (MQ=255)

     cAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACg                                                                                                                                                  <  4:1896137/151‑1 (MQ=255)

     cAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACg                                                                                                                                                  >  1:668470/1‑151 (MQ=255)

      aGTTGATCTTTGCTTTGTTTAACAAATCGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGCATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgg                                                                                                                                                 >  4:950119/1‑151 (MQ=255)

        ttGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgggg                                                                                                                                               >  7:1475423/1‑151 (MQ=255)

         tGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggggg                                                                                                                                              >  6:1851553/1‑151 (MQ=255)

            tCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCgg                                                                                                                                           >  7:1842279/1‑151 (MQ=255)

                       tttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                                >  6:581002/1‑151 (MQ=255)

                       tttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                                >  3:74837/1‑151 (MQ=255)

                       ttaaacaacTGGCCGCACCGGCCGGTGCCATGCAATAAACATCATAATTCCCCTGTCATTAAAATCTGATCACCTTCCCATCGCCTAATGCAGGATCACTCATCATCTGAAGGTTATTGGGGGGCGGAGGCCGGGGCCGGATAACTTTCg                                                                                                                                 <  4:176052/147‑1 (MQ=255)

                        ttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                                >  5:448720/1‑150 (MQ=255)

                         tAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTCATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGa                                                                                                                               >  3:338985/1‑150 (MQ=255)

                          aaCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                             >  4:1642843/1‑151 (MQ=255)

                          aaCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                             >  8:1618606/1‑151 (MQ=255)

                            aaattGGGGGCACCGGCAGGTTCCATTAAATAAAAATAATAATTCCCCTGTTATCAAAATCTGATCACCTTCCCATAGCCTAATGCAGTATCACTCATCATCTGAAGTATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAg                                                                                                                           <  8:1187490/147‑1 (MQ=255)

                              aaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGcc                                                                                                                         >  5:544094/1‑151 (MQ=255)

                               aTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCg                                                                                                                        >  6:355197/1‑151 (MQ=255)

                                 gggggCACCGGCAGGTTCCATGAATTTAAAATCATAAATCCCCTGTTAATAAAATCTGATCACCTTCCCATAGCCTAATTCACGATAACTCATCATCTTAAGTATTTTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCgg                                                                                                                       <  7:525970/150‑1 (MQ=255)

                                 gggcgCCCCGGCCGGTGCCATGCATTAAACATCAAAATTCCCCTGTCATCAAAATCTGATCACCTTGCCATAGCCTAATCCAGGATCACTCATCATCTGAAGTATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGa                                                                                                                      <  8:346943/151‑1 (MQ=255)

                                        ccGCCAGGTGCCATTCATTTAAAAACATAATTCCCCTTTTATTAAAATCTGACCACCTTGCCATACCCTAATGCACCATCACTCATCATCTGAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAt                                                                                                               <  6:878200/151‑1 (MQ=255)

                                           gcatcttcaattcATTTAAAATAATAATTCCCCTTTAATGAAAATCTGATCACCTACCCATAGCCTAATCCAGGATCACTCATCATCTCAAGGATAGTGGGGGCCGCCGGCCGCCCGCGGATAACTTTCGGAATAGCCGGATCAGGATTtc                                                                                                            <  5:937584/139‑1 (MQ=255)

                                               gTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATACCTTTCGGAATAGCCGGATCAGGATTtcctcc                                                                                                        >  2:440417/1‑151 (MQ=255)

                                                     tGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAcc                                                                                                  >  7:795628/1‑151 (MQ=255)

                                                       cATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCgg                                                                                                >  2:1001735/1‑151 (MQ=255)

                                                        aTTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGc                                                                                               >  2:1697127/1‑151 (MQ=255)

                                                        aTTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGc                                                                                               >  7:1618848/1‑151 (MQ=255)

                                                          ttAACATCATAATTCCCCTGTCATCAAAATCTGATCACCTTGCCATAGCCTAACGCAGCATCACTCATCATCTGAAGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                              <  3:1783711/150‑1 (MQ=255)

                                                            aaCATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                           >  7:1962866/1‑151 (MQ=255)

                                                            aaCATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                           <  4:402644/151‑1 (MQ=255)

                                                             cctatctAATTCCCCTTTCATCAAAATCTGATCACCTTCCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGCCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                           <  6:40324/144‑1 (MQ=255)

                                                             aCATCATACATCCCCTGGGATGAAAATCTGACCACCTGGCCCAACCCAACTGCCGGCGCACCCACCACCTGCAGGATCTTGGGGGGAGCCGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACa                                                                                          <  3:1464042/151‑1 (MQ=255)

                                                             aCATCATAATTCCCCTGTGATGAAAATCTGATCACATTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCCGCGTAc                                                                                           >  8:1869907/1‑150 (MQ=255)

                                                              catcatAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                           >  1:2011176/1‑149 (MQ=255)

                                                              catcatAAATCCCCTGTCAGCAAAATCTGATCACTTTCCCATAGCCTAATGCAGGATCACTCATCATCTGACGTATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAg                                                                                         <  6:760134/151‑1 (MQ=255)

                                                                    aaTTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGc                                                                                    >  2:1560825/1‑150 (MQ=255)

                                                                     agcccccTGGCATGAAAAGCGGACCACCTCGCCATCGCCTCATCCCGCACCACTCCACCTCTGAAGCATATTGGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCaa                                                                                  <  3:1054539/148‑1 (MQ=255)

                                                                        ccccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAgg                                                                                >  5:856022/1‑150 (MQ=255)

                                                                        ccccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGCACCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAgg                                                                                <  7:684191/150‑1 (MQ=255)

                                                                         cccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACTCAAGGAt                                                                              >  7:1744661/1‑151 (MQ=255)

                                                                         cccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAt                                                                              >  4:1251938/1‑151 (MQ=255)

                                                                           cTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCGTCCTGAACCGGCTTACAGAACGCAAGGAtt                                                                             >  8:1466578/1‑150 (MQ=255)

                                                                           cTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACAGGCTTACAGAACGCAAGGAtt                                                                             >  4:2000332/1‑150 (MQ=255)

                                                                            tgtgATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTa                                                                            >  3:1585203/1‑150 (MQ=255)

                                                                              tgatgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGTATTTTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACt                                                                          <  2:668324/150‑1 (MQ=255)

                                                                                atcaaaaTCTGATCACCTTCCCATAGCCTACTCCAGGATCCCTCCTCATCTGCAGGATCTTTGGGGGCGCCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGt                                                                       <  4:612335/148‑1 (MQ=255)

                                                                                  gaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGttt                                                                     >  4:1387835/1‑151 (MQ=255)

                                                                                  aaaaaTCTGATCAGCTTGCCATAGCCTAATGCCGGATCACTCATCATCTGACGGTTATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGttt                                                                     <  5:1749769/150‑1 (MQ=255)

                                                                                   aaaaTCTTATCACCTTGCCATTCCATACTGCACTATCACTCATCATCTTTAGGATTTTTGGGGGCGGCGGAGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGttt                                                                     <  2:847578/150‑1 (MQ=255)

                                                                                          gATCACCTTGCCATAGCCTAATGCAGTATCACTCATCATGTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcag                                                                 <  2:76055/147‑1 (MQ=255)

                                                                                          gATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                               >  3:1224175/1‑149 (MQ=255)

                                                                                           aTCACCTTCCCATCCCCTACTGCTCGATCACTCATCATCTTAAGGTTTTTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagccag                                                             <  4:1033495/150‑1 (MQ=255)

                                                                                            tCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTAAGCCAGCt                                                           >  7:1217043/1‑151 (MQ=255)

                                                                                             cACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTg                                                          >  6:45173/1‑151 (MQ=255)

                                                                                              aCCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTg                                                          >  3:941907/1‑150 (MQ=255)

                                                                                                cTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACg                                                       >  2:1457434/1‑151 (MQ=255)

                                                                                                cTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACg                                                       >  1:804593/1‑151 (MQ=255)

                                                                                                cTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACg                                                       >  4:516551/1‑151 (MQ=255)

                                                                                                   cccatccccTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                    <  7:947336/144‑1 (MQ=255)

                                                                                                    ccATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCg                                                   >  4:1430600/1‑151 (MQ=255)

                                                                                                    ccATACCCTAAAGCAGGAGCACACCTCATCTGAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCg                                                   <  6:1501059/151‑1 (MQ=255)

                                                                                                      aTAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGACAGCTGACGATCGTg                                                 >  1:1597602/1‑151 (MQ=255)

                                                                                                       aacccTAATCCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGACGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTgg                                                <  3:1882858/148‑1 (MQ=255)

                                                                                                        aGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggc                                               >  8:824416/1‑151 (MQ=255)

                                                                                                         ccatAATGCAGGATAACTCATCATCTAAATGATTTTGGGGGGCGGAGGACGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                              <  5:1911767/148‑1 (MQ=255)

                                                                                                          ccaaaTCCAGGATCACTCATCCTCTGAAGGATTGTTGGGCGAGGCGCCCGGGGGCGGATAACTTTCGGAATAGCCGTATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcgg                                             <  5:1188011/148‑1 (MQ=255)

                                                                                                          ccTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcgg                                             >  3:1976965/1‑151 (MQ=255)

                                                                                                           cTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcggc                                            >  1:1315775/1‑151 (MQ=255)

                                                                                                            tAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGgcggcg                                           >  5:337442/1‑151 (MQ=255)

                                                                                                              aTGCAGGATCACTCATCATCTGAAGTATATTGGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGtt                                         <  2:2030951/151‑1 (MQ=255)

                                                                                                                gcacgagccctccacatcTGAAGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTa                                       <  5:1489098/137‑1 (MQ=255)

                                                                                                                 cAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTa                                       >  4:262851/1‑150 (MQ=255)

                                                                                                                 cAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCC‑‑ATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAAc                                    >  2:1900403/1‑151 (MQ=255)

                                                                                                                   ggATCACTCATCATCTGAAGGATATTTGGGCGAGGAGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAAc                                    <  4:74838/151‑1 (MQ=255)

                                                                                                                         cTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                                >  7:412466/1‑149 (MQ=255)

                                                                                                                          gcagcatcTGACGGATTTTTGGGGGCGGCGGCCGGCGGCGGAGAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGGTTAAACATCCAGc                             <  8:535166/147‑1 (MQ=255)

                                                                                                                           aatcatcTGAAGTATAGTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGc                             <  2:1262999/149‑1 (MQ=255)

                                                                                                                            atcatcTTATGGATCTTTTGGGGACGTGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTATAAATATCCAGCCa                           <  8:1060147/151‑1 (MQ=255)

                                                                                                                                tcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGACCCGGCTTACAGCACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAACCATCCAGCCAATc                        >  6:1626244/1‑150 (MQ=255)

                                                                                                                                   gAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCa                       <  1:956802/148‑1 (MQ=255)

                                                                                                                                    aaGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAg                   <  1:148371/151‑1 (MQ=255)

                                                                                                                                     aGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGc                  <  8:1377540/151‑1 (MQ=255)

                                                                                                                                     aGGATATTGGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGc                  <  3:562873/151‑1 (MQ=255)

                                                                                                                                       gaaatTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcac                <  1:103356/148‑1 (MQ=255)

                                                                                                                                        atatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcac                >  1:1749400/1‑150 (MQ=255)

                                                                                                                                        atatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTCGTATTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcac                >  8:1953846/1‑150 (MQ=255)

                                                                                                                                         tatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcaca               >  1:1228501/1‑150 (MQ=255)

                                                                                                                                          ttttGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaa              <  8:887257/149‑1 (MQ=255)

                                                                                                                                           tttGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaac             >  7:1323743/1‑150 (MQ=255)

                                                                                                                                            ttGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaaca            <  7:1692763/150‑1 (MQ=255)

                                                                                                                                             tGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCAcaacaa           <  2:974398/150‑1 (MQ=255)

                                                                                                                                               ggggAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATcc        >  2:333352/1‑151 (MQ=255)

                                                                                                                                                  gaggcggaCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcg     <  3:1257439/151‑1 (MQ=255)

ATGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGA  >  DH1/3448373‑3448672

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