BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	3,448,520	0	A	C	21.7%	745.5 / 48.5	318	intergenic (‑80/+76)	malZ/proY	alpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand): ref base A (183/62); new base C (0/68); total (183/135)
Fisher's exact test for biased strand distribution p-value = 1.08e-32
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 335 total aligned reads displayed.

GTAATGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGC  >  DH1/3448370‑3448670

   aTGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggagga                                                                                                                                                     <  1:284535/151‑1 (MQ=255)

   aTGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggagg                                                                                                                                                      >  3:280582/1‑150 (MQ=255)

      aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                  <  3:1054640/151‑1 (MQ=255)

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

       gCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGAc                                                                                                                                                 <  2:510149/151‑1 (MQ=255)

       gCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGAc                                                                                                                                                 <  6:352614/151‑1 (MQ=255)

       gCACTTGATCGTTGCGTTGTTTAACAAATTGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                  <  6:1177464/150‑1 (MQ=255)

          gTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggg                                                                                                                                              <  5:111457/151‑1 (MQ=255)

          gTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggg                                                                                                                                              <  3:809085/151‑1 (MQ=255)

           ttGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgggg                                                                                                                                             >  3:171001/1‑151 (MQ=255)

              aTCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATAATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGAAGGATCACTCATCATCTGAAGTATATTTGGGGGAGGAGGACGGGGGc                                                                                                                                           <  7:865079/150‑1 (MQ=255)

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

                     ctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACt                                                                                                                                   >  4:811368/1‑151 (MQ=255)

                     ctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACt                                                                                                                                   >  4:59525/1‑151 (MQ=255)

                        tgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTc                                                                                                                                >  3:865485/1‑151 (MQ=255)

                        tgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCAATTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTc                                                                                                                                >  6:823747/1‑151 (MQ=255)

                          tttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                              >  7:1274283/1‑151 (MQ=255)

                          tttAAAAAATGGGGGCACCGGCAGGTTCCATTCATTTAACATAATAATTCCCTTGTTATTAAAATCTTATCACGTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAATTATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCg                                                                                                                               <  5:985627/150‑1 (MQ=255)

                           ttAAAAAATGGGGGCACCGTCAGGTGCCATTAATTAAACATCATATTTCCCCTGTTATCAAAATCTGATCAGCTTGCCATAGCCTATTGCATGATCACTGATCATCTGAAGGATATTGGGGGGCGGAGGCCGGGGGCGGATAACTTTCgg                                                                                                                              <  7:480379/150‑1 (MQ=255)

                            tAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGa                                                                                                                             >  2:263373/1‑150 (MQ=255)

                             aaCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                           >  6:242134/1‑151 (MQ=255)

                              aCAAAGGGGGCCCCCGGCAGTTGCCATGAATTTAACATCATAATTCCCCTGTGATCAAAATCTGATCACCTTTCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGTTATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATa                                                                                                                          <  8:389879/151‑1 (MQ=255)

                                aaaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGc                                                                                                                        >  8:759674/1‑151 (MQ=255)

                                 aaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGcc                                                                                                                       >  8:143883/1‑151 (MQ=255)

                                   tGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCgg                                                                                                                     >  8:590946/1‑151 (MQ=255)

                                    gggggCCCGGGCAGGGGCCATGCATTTAACATCATAATTCCCGTGTGATGAAATTCTGATCACCTTGCCATAGCCTATTGCAGGATCACTCATCATGTGAAGTATTTTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCgg                                                                                                                     <  1:1133334/150‑1 (MQ=255)

                                    gggcccccctgCAGGTGCCATGAATTTAACATCATAATTCCCCTGTTATTAAAATCTCATCACCTTCCCATACCCTAATGCCGGACCACTCATCATCTGATGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGa                                                                                                                    <  4:1143349/141‑1 (MQ=255)

                                      gggCACCGGCAGGTGCCATGCATTAAACATCATAAGTCCCCTGTCATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATCTTGGGGGGCCGCGCCCGCCGCCGGATAACTTTCGGAATAGCCGGATc                                                                                                                  <  7:989180/151‑1 (MQ=255)

                                       ggCACCGGCAGGTGCCATTCATTTAAAATCATAATTCCCCTGTTATTAAAATCTGATCACCTTGCCATAGCCTAATCCAGGACCACTCATCCTCTTAAGTATATTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCa                                                                                                                 <  7:323510/151‑1 (MQ=255)

                                           ccGGCAGGTGCCATCCATTAAACATCATAAATCCCCTTTGATTAAATTCTGATCCCCTTGCCATAGCCTAATGCAGTATCACTCATCATCTGAAGTATATTTGGGGGCGCCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAt                                                                                                             <  2:824779/151‑1 (MQ=255)

                                            cctcagcttccattaaTTAAACATCATAATTCCCCTCTCACCAATATCTCATCACCTTGCCATACCCTAATCCACCATCACTCATCATCTGAAGCATATTTGGGGGCGCCGGCCGGCCGCGGATAACTTTCGGAATAGCCGGATCAGGAtt                                                                                                            <  6:203490/136‑1 (MQ=255)

                                            cGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTANTGCANGANCACTCATCANCNNAAGGATANTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAtt                                                                                                            >  4:702966/1‑151 (MQ=255)

                                                agggtccATGCATTTAACATCATAATTCCCCTTTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGCAGCACTCATCATCTGATGTATATTGGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                        <  6:550768/146‑1 (MQ=255)

                                                    gCCATCCATTAAACATCATAATGCCCCTGTCATGAAAATCTGATCCCCTTCCCATAGCCTCATGCAGGATCACTCATCATCTGCCGGCTCGTGGCGCGCCCCGCCCGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTg                                                                                                    <  7:1003862/151‑1 (MQ=255)

                                                    cccATGCATTTAACATCATAATTCCCCTGTTACTAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGTTTTTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTcctcct                                                                                                     <  2:900908/149‑1 (MQ=255)

                                                     ccATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTg                                                                                                    >  7:791813/1‑150 (MQ=255)

                                                      cATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGa                                                                                                   >  5:675486/1‑150 (MQ=255)

                                                       aTGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAc                                                                                                 >  1:802352/1‑151 (MQ=255)

                                                       aTGCATTTAACATCATAATGCCCCTGTTATCAAAATCTGATCACCTTGCCATACCCTAATGCCGGATCACTCATCATCTGAATGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAc                                                                                                 <  4:300620/151‑1 (MQ=255)

                                                        tGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAGTTCCTCCTGAAcc                                                                                                >  2:300096/1‑151 (MQ=255)

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

                                                              tAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                           >  5:578256/1‑150 (MQ=255)

                                                              tAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                           >  6:778219/1‑150 (MQ=255)

                                                              tAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                           >  2:1062868/1‑150 (MQ=255)

                                                                   tcatAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGa                                                                                      >  7:71372/1‑150 (MQ=255)

                                                                    catACTTCCCCTGTGATGAAAATCTGACCCCCTTGCCATAGCCTACTGGAGTACCACTCATAATTTGAAGTATATTTGGGGGTGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAAc                                                                                    <  7:204248/151‑1 (MQ=255)

                                                                       aaTTCCCCTGTGATGAAACTCTGATCACCTTCCCATAGCCTACTGCAGCAGCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                 <  2:1282203/151‑1 (MQ=255)

                                                                       aaTACCCCTGTCATAAAAAGATCAACAACAAGCCATAGCCAAATCCAGCATCACTCATCATCTGAAGTATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                 <  6:700458/151‑1 (MQ=255)

                                                                        aTTACCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAACGGCTTACAGAACGCa                                                                                 >  8:1261380/1‑150 (MQ=255)

                                                                         ttCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGTGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTACCTCCTGAACCGGCTTACATAACGCa                                                                                 >  1:518328/1‑149 (MQ=255)

                                                                         ttCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                 >  5:304937/1‑149 (MQ=255)

                                                                         gcccccTGTGATCAAAAGCTGACCACCTTCCCATACCCTAATGCAGGAGCACTCATCATCTGAAGGATATTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAg                                                                               <  8:686557/149‑1 (MQ=255)

                                                                             ccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAtt                                                                           >  3:714674/1‑151 (MQ=255)

                                                                                 tgatgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTg                                                                       >  5:406946/1‑151 (MQ=255)

                                                                                  gatgaAAAGCTGATCACCTTCCCATCGCCAAATGCAGCATCCCTCATCATCTGACGGATATTTGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTg                                                                       <  3:598895/150‑1 (MQ=255)

                                                                                  catgaAAATCTGAGCCCCTTCCAATAGCCTAATGCAGCATCACTCATCATCTGAAGGATATTTGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTgg                                                                      <  4:402485/150‑1 (MQ=255)

                                                                                      aaaaTCTGATCACCTGGCCATAGCCTAATGCAGGATCACTCATCATCTGACGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                  <  7:1084316/151‑1 (MQ=255)

                                                                                       aaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                 >  4:689741/1‑151 (MQ=255)

                                                                                          tCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagc                                                              >  1:366627/1‑151 (MQ=255)

                                                                                             gaacACCTTGCCATACCCCACTCCAGGATCACTCATCATCTGAAGGATTTTGGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcca                                                            <  2:532436/147‑1 (MQ=255)

                                                                                             gATCACCTTGCCATACCCTAATGCAGGATCACTCATCATCTCAAGTATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagccag                                                           <  8:1108629/151‑1 (MQ=255)

                                                                                             gATCACCTTCCCATAGCCTAATGCAGGATCACTCATCATCTGATGGATTTTGGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                             <  4:1248001/149‑1 (MQ=255)

                                                                                               tCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                         >  1:1285506/1‑151 (MQ=255)

                                                                                                  ccTTGCCATAGCCTAATGCAGCACCACTCATCATCTGCAGGATATTGGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGAc                                                      <  3:994284/151‑1 (MQ=255)

                                                                                                  ccTTCCCATACCCTAATGCAGCCGCACTCATCCTCTGACCGATCTTTGGGGGAGGCCGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGAc                                                      <  4:171251/151‑1 (MQ=255)

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

                                                                                                     tGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGAt                                                   >  4:910867/1‑151 (MQ=255)

                                                                                                      gCCATAGCCTAATGCAGCATCACTCATCATCTGAAGGATATTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                  <  3:36156/151‑1 (MQ=255)

                                                                                                       ccATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCATCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGATGACGATCg                                                 >  5:245553/1‑151 (MQ=255)

                                                                                                         aTAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCAGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                               >  6:153505/1‑151 (MQ=255)

                                                                                                         aTAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCACGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                               >  5:1182991/1‑151 (MQ=255)

                                                                                                            gCCTAAGGCAGGATCACTCATCCTCAGAAGGATTTTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                            <  3:562168/151‑1 (MQ=255)

                                                                                                            cccTAATGCCGGAGCCCTCATCATCTGAAGGATATTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                            <  3:1201602/150‑1 (MQ=255)

                                                                                                                actgcagccgccctcctcAGCGGAAGGAGATTGGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcggc                                          <  1:1209541/137‑1 (MQ=255)

                                                                                                                aaTGCACCATCACCCATCTTCTGTAGGATCTTGGGGGGAGGCGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcggc                                          <  6:1119559/149‑1 (MQ=255)

                                                                                                                aaTCCACGAACACTCATCATCTCAAGGATATTTGGGGGAGCAGGACGGCCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGt                                        <  7:1375479/151‑1 (MQ=255)

                                                                                                                 aTGCAGGATCTCTCATCCTGGGTAGCATAGTGGGGGGCGGCGCCCGCGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGtt                                       <  5:1324509/151‑1 (MQ=255)

                                                                                                                   ccAGGACCACGCAACATCTGAAGGAGATGGGGGCGCCGCGCCCCGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTa                                     <  3:417575/150‑1 (MQ=255)

                                                                                                                    caccaccACTCCTGATCTGCCGGATCGGTGGGGGCGGCGCCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTaa                                    <  8:1159097/145‑1 (MQ=255)

                                                                                                                      ggATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAAc                                  >  1:1352915/1‑151 (MQ=255)

                                                                                                                       gATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAAc                                  >  5:1132232/1‑150 (MQ=255)

                                                                                                                           ccgccgcatcTGCAGGATATTGGGGGGCGGAGGCCGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                             <  5:548717/145‑1 (MQ=255)

                                                                                                                            cTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                              >  2:684170/1‑149 (MQ=255)

                                                                                                                             acatcatCTGACGGCTCTTTGGCGGACGCCGACGGCCCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGc                           <  2:7397/150‑1 (MQ=255)

                                                                                                                              catcttCTGAAGGATATTTGGGGGAGGAGGAGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGc                           <  4:982442/150‑1 (MQ=255)

                                                                                                                                 catcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGAGGATCGTGGCGGCGTTTATGCATCCAGCGAAt                       >  8:297569/1‑151 (MQ=255)

                                                                                                                                 catcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAAt                       >  5:908446/1‑151 (MQ=255)

                                                                                                                                   gcTGAAGGATAGTGGGGGGCGCCGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCa                     <  3:912346/150‑1 (MQ=255)

                                                                                                                                    cTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCa                     >  1:1089285/1‑150 (MQ=255)

                                                                                                                                      gAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAAc                   >  3:1138844/1‑150 (MQ=255)

                                                                                                                                        aGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGc                >  5:1368477/1‑151 (MQ=255)

                                                                                                                                        aGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGc                >  5:81359/1‑151 (MQ=255)

                                                                                                                                         ggATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGc                >  3:650020/1‑150 (MQ=255)

                                                                                                                                            tatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGccca             >  3:112479/1‑150 (MQ=255)

                                                                                                                                             atTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaac           >  1:1059331/1‑151 (MQ=255)

                                                                                                                                             atTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaac           >  5:193912/1‑151 (MQ=255)

                                                                                                                                              tttGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGCCGGCGTTTAAAGATCCAGCCAATCAACAGcacaccc          >  6:435912/1‑148 (MQ=255)

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

                                                                                                                                                  ggggCGGAGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATcc      <  7:1112175/151‑1 (MQ=255)

                                                                                                                                                    ggaggcggcCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAc    <  2:1295641/151‑1 (MQ=255)

GTAATGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGC  >  DH1/3448370‑3448670

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