BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	3,448,526	0	A	C	25.0%	768.3 / 85.2	390	intergenic (‑86/+70)	malZ/proY	alpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand): ref base A (215/76); new base C (1/96); total (217/173)
Fisher's exact test for biased strand distribution p-value = 9.27e-42
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 416 total aligned reads displayed.

AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAACC  >  DH1/3448376‑3448676

 gCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGAc                                                                                                                                                       >  5:1220426/1‑151 (MQ=255)

   aGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgg                                                                                                                                                     <  3:655140/151‑1 (MQ=255)

   aGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgg                                                                                                                                                     <  6:435325/151‑1 (MQ=255)

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

        aTCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCg                                                                                                                                                >  7:921561/1‑151 (MQ=255)

             tgcttttTTTAACAACTGGGGCCCCCGGCCGGGGCCAGGAATTTAACATCATAATTCCCCTGTGATGAAAATCTTCTCCCCTTGCCATACCCTAATGCAGGCTCACTCATCCTCTGCCCGATATTTGGGGGAGGAGGCCGGGGGCGGATa                                                                                                                                            <  2:136242/150‑1 (MQ=255)

              gctttgTTTAACACAGGGGGGCCCCGGCAGGTGCCATTCATTTAACATAATAATTCCCCTGTTATTAAAATTTGATCACCTTGCAATAGCCTAATGCAGGATCACTCATCATCTGAAGGTTTTTTGGGGGAGGAGGCCGGGGGCGGATAAc                                                                                                                                          <  4:994542/151‑1 (MQ=255)

                tttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACtt                                                                                                                                        >  8:566732/1‑151 (MQ=255)

                tttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACtt                                                                                                                                        >  4:219956/1‑151 (MQ=255)

                  tgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTc                                                                                                                                      >  7:1441284/1‑151 (MQ=255)

                   gTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCg                                                                                                                                     >  8:695406/1‑151 (MQ=255)

                    ataaccAAATGGGGGCACCTGCAGCTTCAATTAATTTAACATCATAATTCCCCTGATATCAAAATCTCATCACCATGCCATCCCCAAATCCAGCATCACTCATCATCTCAAGCATATTTGCGGGAGGAGGCCGGGGGCGGATAACTTTCg                                                                                                                                     <  8:1343466/145‑1 (MQ=255)

                        aCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                                 >  1:945850/1‑150 (MQ=255)

                          aaaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGc                                                                                                                              >  5:595742/1‑151 (MQ=255)

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

                                  gCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                        >  3:463734/1‑149 (MQ=255)

                                       ggCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAttt                                                                                                                 >  7:1236649/1‑151 (MQ=255)

                                         cAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcc                                                                                                               >  4:501776/1‑151 (MQ=255)

                                         cAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcc                                                                                                               >  3:114143/1‑151 (MQ=255)

                                            gTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcctcc                                                                                                            >  7:640916/1‑151 (MQ=255)

                                             tGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTcctcct                                                                                                           >  6:624874/1‑151 (MQ=255)

                                             ggCCATGCAATTAACATCATAATTCCCCTGTGATCAAAATCTGACCCCCTTGCCATAGCCTAATCCAGGATCACTCATCATCTGAAGGTTTTTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTcctcct                                                                                                           <  3:1751288/150‑1 (MQ=255)

                                               ccATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTg                                                                                                          >  8:1916154/1‑150 (MQ=255)

                                                cATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGa                                                                                                         >  4:851113/1‑150 (MQ=255)

                                                cATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGa                                                                                                         >  8:1241961/1‑150 (MQ=255)

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

                                                 aTGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATAACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTAGGAATAGCCGGATCAGGATTTCCTCCTGAAc                                                                                                       >  7:561222/1‑151 (MQ=255)

                                                 aTGAATTTAACATCATAATTCCCCTTTTATGAAAATCTTATCACCTTGCAATAGCCTAATGCAGGATCACTCATCATCTGTTGTATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAc                                                                                                       <  8:1835134/151‑1 (MQ=255)

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

                                                      tttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                                  >  6:1437453/1‑151 (MQ=255)

                                                       taaaCATCAAAATTCCCCTGTTATTAAAATCTGATCACCTTGCCATAGCCTAATCCAGGATCACTCATCATCTTAAGGATTTTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                                  <  2:1566894/148‑1 (MQ=255)

                                                        tAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATACCTTTCGGAATAGCCGGATCAGGATTTCCGCCTGAACCGGcgt                                                                                                 >  4:1497710/1‑148 (MQ=255)

                                                         aaCATCATAATTCCCCTGTTATGAAAATCTGATCACCTTTCCATAGCCTAATGCAGGATCACTCATCATCTGATTGATTTTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                               <  2:754849/151‑1 (MQ=255)

                                                         aaCATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                               >  4:1283922/1‑151 (MQ=255)

                                                          aaataatAATTCCCCTGTTATGAAATTCTGATCACCTTGCCATAGCCTAATGCAGTATCACTAATCATCTTAAGGATTTTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACa                                                                                              <  2:608110/146‑1 (MQ=255)

                                                             acatAATTCCCGTGTGAATAAAATCTCCGCACCTTCCCCCACCCAAATGCCCGATCCCTCATCAACTGAAGGATATTTGGGGGCGCCCGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGaa                                                                                           <  3:1069002/150‑1 (MQ=255)

                                                             aaaaaaTACCCCTGTGAAGAAAATCTCATCACCTTGCCATACCCTAATGCAGTATCACTCATCTTCTGATGGATATTGGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGaa                                                                                           <  6:1321877/147‑1 (MQ=255)

                                                              catcaatCCCCTGAGAAGAAAATCTCACCACCTTCCCATACCCAAATCCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAAc                                                                                          <  5:311188/145‑1 (MQ=255)

                                                              catAATGCGCCTGTTATGAAAAGCTGAGAAACTTGCCATCGCCTCCTGCCCCAGCACTCAACAGCTGAACGCTATTTGGCGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAAc                                                                                          <  1:1066164/151‑1 (MQ=255)

                                                                aaaTTCCCCTGGCATAAAACGCTGAACCCCTTCCCATACCCAACTCCACCCTCACACATCATCTCCAGGATCTTTGGGGGCGCACGCCGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGc                                                                                        <  5:453680/150‑1 (MQ=255)

                                                                  aTTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                       >  4:446651/1‑150 (MQ=255)

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

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

                                                                          gttatgaAAATCTGATCCCCTTTCCATAGCCTAATGCAGTATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGTCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTa                                                                                <  2:510973/146‑1 (MQ=255)

                                                                          gttatgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGAGCCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACt                                                                              <  1:1439935/148‑1 (MQ=255)

                                                                          gggaggaAAATCTGCGCACCGTGCCATAGCCTAATGCAGCATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACt                                                                              <  1:129313/146‑1 (MQ=255)

                                                                           tgatgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCAACATCTGATGGATATTTGCGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTACTGAACCGGCTTACAGAACGCAAGGATTACTg                                                                             >  5:1816855/1‑151 (MQ=255)

                                                                             atgaAAATCTGATCACCTTGCCATAGCCTACTGCAGGATCACTCATCATCTGAAGTATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTgg                                                                            <  8:1553562/150‑1 (MQ=255)

                                                                             aacaaaaTCTGCTCACCTGGCCATCCCCTAAGGCAGGAGCACTCATCATCTGAAGCATATGTGGGGGAGCCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGt                                                                           <  7:1770823/148‑1 (MQ=255)

                                                                               gaAAACCTGATCCCCTACCCATACCCTACCGCAGGACCACTCATCATCTGAAGGATATTGGGGGGCGGAGGCCGGGCCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGttt                                                                         <  4:1129454/151‑1 (MQ=255)

                                                                               caaaaTCTTATCACCTTGCCATAGCCTAATGCAGTATCAGTCATCATCTGAATGATATTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGttt                                                                         <  2:948092/150‑1 (MQ=255)

                                                                                 aaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGCTTACTGGTTTTc                                                                       >  7:462367/1‑151 (MQ=255)

                                                                                 aaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                       >  7:978502/1‑151 (MQ=255)

                                                                                   aTCAGATCACCTTGCCATCGCCTAATCCAGGATCACTCATCATCTGCAGGATCGTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcag                                                                     <  7:1372703/151‑1 (MQ=255)

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

                                                                                    gcTGCTCCCCTTGCCATAGCCGAATGGAGGATCACCCATCATCTGAAGGATTGCTGGGCGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagc                                                                    <  7:1296727/150‑1 (MQ=255)

                                                                                     cTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                                   >  5:1115895/1‑151 (MQ=255)

                                                                                      tGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCGGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                                   >  5:817131/1‑150 (MQ=255)

                                                                                      tGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                                   >  5:526741/1‑150 (MQ=255)

                                                                                       gATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                                   >  7:1596369/1‑149 (MQ=255)

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

                                                                                          cACCTTTCCATAGCCTACTGCCGGCTCACTCATCATCTGAAGCATTTTTGGGGGAGGTGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAGTACTGGTTTTCAGCCAGCTg                                                              <  2:1848611/151‑1 (MQ=255)

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

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

                                                                                              ttGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACg                                                           >  7:1915677/1‑150 (MQ=255)

                                                                                                  catccccTAATGCAGGATCACTCATCATCTGAACGATATTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGt                                                      <  3:881057/146‑1 (MQ=255)

                                                                                                  caaacccTACTGCCGGATCACTCCTCACCTGAAGGATATTGGGGGCCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGt                                                      <  6:1566619/146‑1 (MQ=255)

                                                                                                   aTAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                                     >  2:296770/1‑151 (MQ=255)

                                                                                                     cgCCTAATCCTCGAGCACTCTTCATCTGAAGGTTATTTGGGGGAGCCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggc                                                   <  4:378718/150‑1 (MQ=255)

                                                                                                     aGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggc                                                   >  4:596718/1‑151 (MQ=255)

                                                                                                      ggcTAATCCAGGATCACTGATCAGCTTATGGTTTTTTGGGGGCGGAGGAGGGGGGCGGATAAGTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                                  <  6:192349/149‑1 (MQ=255)

                                                                                                      gcgaactGCAGGATCACGCCTCCTCTGAAGGATTGTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                                  <  7:683015/145‑1 (MQ=255)

                                                                                                      gCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                                  >  8:1707752/1‑151 (MQ=255)

                                                                                                      ccctactgccccaccACTCAACATCTGAAGGATATGGGGGCGACCCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                                  <  1:373610/137‑1 (MQ=255)

                                                                                                       ccTAATGCCGGATCACTCATCAGCTGACGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcgg                                                 <  8:1439891/151‑1 (MQ=255)

                                                                                                         tAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGgcggcg                                               >  6:1604015/1‑151 (MQ=255)

                                                                                                         tAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGgcggcg                                               >  2:228256/1‑151 (MQ=255)

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

                                                                                                             gccggATCCCTCATCCTCTGACGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTa                                           <  6:815344/148‑1 (MQ=255)

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

                                                                                                                    cACTCCTCATCTGAAGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                                    <  4:1591427/151‑1 (MQ=255)

                                                                                                                    cACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                                    >  1:546347/1‑151 (MQ=255)

                                                                                                                    cACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATc                                     <  1:1564346/150‑1 (MQ=255)

                                                                                                                     ccTCAGCCTCTGAAGGCTTGTGGGGGGCGGAGCCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                                   <  5:1534987/150‑1 (MQ=255)

                                                                                                                       gcatcatCTGAAGGATATGGGGGGGCGGCCGCCGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGc                                 <  3:1237909/150‑1 (MQ=255)

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

                                                                                                                            atcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATc                            >  1:193864/1‑151 (MQ=255)

                                                                                                                               tGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAAc                         >  1:1870755/1‑151 (MQ=255)

                                                                                                                                 ccgtatatTTGGGGGAGGCTGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAg                       <  6:1666352/147‑1 (MQ=255)

                                                                                                                                 aaGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAAc                         >  3:1534343/1‑149 (MQ=255)

                                                                                                                                    gATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcac                    >  8:1143520/1‑151 (MQ=255)

                                                                                                                                    gATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcac                    >  5:1300002/1‑151 (MQ=255)

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

                                                                                                                                        tttGGGGGAGGTGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaac                 <  8:1508253/150‑1 (MQ=255)

                                                                                                                                            ggggAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATcc            <  3:1097061/151‑1 (MQ=255)

                                                                                                                                              ggaggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGAGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAc          >  1:1776060/1‑151 (MQ=255)

                                                                                                                                                   cggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGaa     <  6:1166584/150‑1 (MQ=255)

                                                                                                                                                    ggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGaa     >  4:758685/1‑150 (MQ=255)

                                                                                                                                                      aCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAAcc  >  1:970722/1‑151 (MQ=255)

                                                                                                                                                      aCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAAcc  >  6:998077/1‑151 (MQ=255)

AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAACC  >  DH1/3448376‑3448676

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