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	25.1%	783.5 / 59.9	362	intergenic (‑80/+76)	malZ/proY	alpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand): ref base A (195/74); new base C (0/90); total (195/167)
Fisher's exact test for biased strand distribution p-value = 2.22e-39
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 383 total aligned reads displayed.

GTAATGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGC  >  DH1/3448370‑3448670

  aaTGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggagg                                                                                                                                                      >  8:1346477/1‑151 (MQ=255)

  aaTGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggag                                                                                                                                                       <  4:840345/150‑1 (MQ=255)

     gaggagTTGATCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagg                                                                                                                                                   <  2:266392/151‑1 (MQ=255)

      aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagg                                                                                                                                                   >  2:1646796/1‑150 (MQ=255)

       gCAGTTGATCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGAc                                                                                                                                                 <  1:1068208/151‑1 (MQ=255)

       gCAGTTGATCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGAc                                                                                                                                                 >  1:1013544/1‑151 (MQ=255)

             gATCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGCTCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgggggg                                                                                                                                           >  8:129186/1‑150 (MQ=255)

             gATCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggggg                                                                                                                                            <  2:1279212/150‑1 (MQ=255)

             gATCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGc                                                                                                                                           >  7:828579/1‑151 (MQ=255)

              aTCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCg                                                                                                                                          >  3:960709/1‑151 (MQ=255)

              aTCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCg                                                                                                                                          >  8:1679687/1‑151 (MQ=255)

               tCTTTGCTTTGTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCg                                                                                                                                          <  8:1091131/150‑1 (MQ=255)

                    gctttgTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAAc                                                                                                                                    >  8:1014354/1‑151 (MQ=255)

                     ctttgTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATaagt                                                                                                                                   >  2:1142732/1‑149 (MQ=255)

                       ttgTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACttt                                                                                                                                 >  4:1204881/1‑151 (MQ=255)

                        tgTTTAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTc                                                                                                                                >  5:1436681/1‑151 (MQ=255)

                          tttAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                              >  3:703384/1‑151 (MQ=255)

                            tAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGa                                                                                                                             >  5:595785/1‑150 (MQ=255)

                            tAACAAATGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGa                                                                                                                             >  4:420162/1‑150 (MQ=255)

                             acatatgggggg‑CACCGGCAGGTGACATTCATTTAACATCATAATTCCCCTGTGATTAAAATGTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATGTGAATTATATTTGGGGCCGGCCGCCGGGGGCGGATAACTTTCGGAAt                                                                                                                           <  6:740325/144‑1 (MQ=255)

                              aCAACTGCGTG‑CACCGCCAGGAGCCATTAAATTAACATCATAATTCCCCTTTCATTACACTCTCATCACCTTGCCATCGCCTAATGCAGGATCACTCATCATCTGATGGATTCTTGCGGGAGGCGGACGGGGGCGGATAACTTTCGGAATa                                                                                                                          <  8:1208132/151‑1 (MQ=255)

                                ataTGGCGG‑CCCCGGGAGGTGCCATGCATTTAAAATCATAATTCCCCTGTGATTAAAATCTCTTCACCTTGCCATAGGCTAAAGGACGATCACTCCTCATCTGAAGGTTTTTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGc                                                                                                                        <  5:240402/149‑1 (MQ=255)

                                 aaTGGGGCCCACCGGCCGGTGCCATTAAATAAACAACATAATTCCCCTTTGATGAAAATCTGACCACCTTCCCAAAGCCTAATGCAGCATCCCTCATCATCTGATGTATATTTGGGGGAGGAGGCCGGGGCCGGATAACTTTCGGAATAGc                                                                                                                        <  4:1452780/151‑1 (MQ=255)

                                  aTGGGGG‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCg                                                                                                                      >  8:197648/1‑151 (MQ=255)

                                      ggg‑CCCCGGCACGTGCCATGCAATTAACATCATAATTCCCCTCTCAACAAAAGCGCACCACCTTCCCCTCGCCTAATCCAGGATCACTCCTCATCTGAAGGATCTTTGGGGGCGGCCGCCGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                  <  8:38877/151‑1 (MQ=255)

                                      ggg‑CACCGGCTGGTGCCATGCATTAAACATCTAATTTACCCTCTGATTAAACTCTGATCACCTGCCCATAGACTATTGCAGGATCCCTCATCATCTGACGGATTTTTGGGGGAGGCCGCCGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                  <  1:55636/151‑1 (MQ=255)

                                      ggg‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                  >  1:41271/1‑151 (MQ=255)

                                        g‑CACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                  >  7:1549095/1‑149 (MQ=255)

                                          cACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAgg                                                                                                               >  3:142755/1‑151 (MQ=255)

                                           cccggccggGGCCATGCATTTAACATCATAATTCCCCTGTGATTAAAATCTGCACCCCTTGCCATAGCCTAATGCCGCATCACTCATCATCTGAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGa                                                                                                              <  6:1418817/150‑1 (MQ=255)

                                              ggCTGGTGCCATTCATTAAACATCATAATTCCCCTCTCATGAAACTCTCATCCCCTTCCCCTAGCCTACTGCACCTTCACTCATCATCTTTTGGATTGGGGGGGGCGGCCGCCGCCCGCGGATAACTTTCGGAATAGCCGGATCAGGAttt                                                                                                           <  2:1395394/151‑1 (MQ=255)

                                              ggCAGGTGCCATTCATTTAAAATCAAAACTCCTCTCTTAATAAAATCTTATCACCTTGCCATCCCCAAAATCCCCATCACTCATCTTCTTATTGATATTTGGGGGCCGACGTCGCCGGCGGATAACTTTCGGAATAGCCGGATCAGGAttt                                                                                                           <  6:560004/151‑1 (MQ=255)

                                              ggCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAttt                                                                                                           >  8:1575116/1‑151 (MQ=255)

                                               gCAGGTGCCATCAATATACCACCATAATTCCCCTGAGATCAAAATCTGAGCACCTTGCCATAGCCTACTGCAGGATCACTCATCATCTGACGGATAGTGGGGGGCGGCGGCCGCCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtc                                                                                                          <  5:1620852/151‑1 (MQ=255)

                                                 aGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                        >  5:869378/1‑151 (MQ=255)

                                                 aGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGAGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                        >  5:751122/1‑151 (MQ=255)

                                                 aGGAGCCATGCATTTAACATCATAATTCCCCTCTCATCAAAATCTGACCACCATGCCATACCCTAATGCAGCATCATTAATCATCTGAACCATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                        <  1:762282/151‑1 (MQ=255)

                                                    tGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTcctcct                                                                                                     >  2:1520933/1‑151 (MQ=255)

                                                     gCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTg                                                                                                    >  8:757423/1‑151 (MQ=255)

                                                       cccgCATTTAACAACCAACATCCCCTGTGCTGAAAACCCGCACCCCTTTCCCTACCCTAAACCACGATCCCTACACATCTTAACGCTTTTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGaa                                                                                                  <  4:304453/148‑1 (MQ=255)

                                                          gCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCg                                                                                               >  8:456816/1‑151 (MQ=255)

                                                          ccATTTAACCTCATACTTCCCCTCTCATCAATATCGCATCCCCTAGACATAGCCTAATCCAGCACCCCTCATCATCTGACGGTTATTGGGGCGCGGCGGCCGCGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCg                                                                                               <  6:1065624/150‑1 (MQ=255)

                                                            aTTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGc                                                                                             >  8:70032/1‑151 (MQ=255)

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

                                                            aTTTAACATCATAATTCCCCTGTGATGAAAATCAGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATACCCGGATCAGGATTTCCTCCTGAACCGGc                                                                                             >  3:53685/1‑151 (MQ=255)

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

                                                             tttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                            >  5:972837/1‑151 (MQ=255)

                                                               tAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGCAGGAGGACGGGTGCGGATAACTTTCGGAATAGCCGGATAAGGATTTCCTCCTGAACCGGAtt                                                                                           >  3:434052/1‑150 (MQ=255)

                                                               aaacaacATAATCCCCCTGTCATCAAAATCCGACCACCTTCCCATAGCCTAATCCACGATCACTCATCATCGGACGGATATTTGGGGGCGGCCGCCGGGGGCGGATAACTTTCGGAATAGCCAGATCAGGATTTCCTCCTGAACCGGCtt                                                                                           <  2:932512/149‑1 (MQ=255)

                                                                 acacCATAAGCCCCCTGTCATCAAAAGCCGCCCCCCTTCCCATCGCCTCATCCACGATCACTCATCATCTGACCGATATTTGGGGCCGGCCGCCGGCCGCGGATAACTTTCGGAATAGCCGGAGCAGGATTTCCTCCTGAACCGGCTTACa                                                                                        <  6:1667078/151‑1 (MQ=255)

                                                                 aCATCATAATTCCCCTGTCAAGAACACCGGAGCCCCTTCCCATAGCCAACACCCGGACCCCTAAACACCTGAAGGATATTGGGGGGACGAGCCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCTTCCTGAACCGGCTTACa                                                                                        <  1:506226/151‑1 (MQ=255)

                                                                   atcatAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                         >  1:1539989/1‑148 (MQ=255)

                                                                          ttCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                 >  3:329253/1‑149 (MQ=255)

                                                                           cccccTGTGATCAAAATCCGCTCACCTTCCCATACCCTAATCCCGGATCACTCATCATCTGCAGGATATCTCGGCGCGGCGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAgg                                                                              <  4:553812/150‑1 (MQ=255)

                                                                             cccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAt                                                                            >  8:868610/1‑151 (MQ=255)

                                                                              ccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTCCAGAACGCAAGGAtt                                                                           >  7:225297/1‑151 (MQ=255)

                                                                                tgtgATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTAc                                                                         >  3:815430/1‑151 (MQ=255)

                                                                                 gtcatgaAAATCTGATCACCTTCCCATAGCCTAATGCAGGAGCACTCAGCATCTGAAGGATATTTGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACt                                                                        <  3:1167268/148‑1 (MQ=255)

                                                                                   gatcaAAATCTGAGCCCCTTGCCATACCCTAACGCACGAGCACTCATCATCTGCAGGATCTGTGGGGGAGGAGGCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTgg                                                                      <  7:36562/151‑1 (MQ=255)

                                                                                     agaAACTCTGAACACCAGCACATAGCATACTCCAGCAGCACTCATCATCCGACCGAGTGTGGGGGGCCGAGGCCGGCCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtt                                                                    <  1:16185/150‑1 (MQ=255)

                                                                                      caaaaTCTCCACACCATCCCATCGCCAAACGCCCGAGCACACATCATCTGAAGGATAGGTGGGCGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGttt                                                                   <  7:712048/150‑1 (MQ=255)

                                                                                       aaaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                  >  5:1048262/1‑151 (MQ=255)

                                                                                       aaaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTCCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                  >  7:766701/1‑151 (MQ=255)

                                                                                         aatcgcaTCACCTAGCCATAGCCTACCCCCGGATCACTCATCTTCTGACGGATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTca                                                                <  4:1285559/145‑1 (MQ=255)

                                                                                         aaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttta                                                                 >  6:998227/1‑149 (MQ=255)

                                                                                             tGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCTGAATAGCCGCATCAGGATTTCCTCCTGAACAGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                             >  4:1169947/1‑150 (MQ=255)

                                                                                               ccccccTTGCCATACCCAAAAGCAGCCTCCCTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagccagc                                                          <  1:1219210/147‑1 (MQ=255)

                                                                                                ccACCTTGCCATAGCCTAAGCCAGGAGCACTCATCATCTGAACGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                         <  2:506731/150‑1 (MQ=255)

                                                                                                 cACCTTGCCATAGCCTAAGGCAGGACCACCCATCATCTGACGGATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTg                                                        <  2:1234118/151‑1 (MQ=255)

                                                                                                  cccTTGCCATAGCCAAAGGCCGGAGCACTCATCATCTGTAGGATTTTTGGGGGCGGCGGCGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGa                                                       <  3:1501993/150‑1 (MQ=255)

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

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

                                                                                                      ttccATATCCTAATTCAGTCTCACTTATCATCTTAAGTTTTTTTTGGGGAGGTGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACg                                                     <  4:703736/147‑1 (MQ=255)

                                                                                                      tgccaaccccTAACGCAGGATCACTCATCATCTGAAGGATATTGGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACg                                                     <  4:1127000/141‑1 (MQ=255)

                                                                                                       gCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTAACAGAACGCAAGCATTACTGGTTTACAGCCAGCTGACGATc                                                  >  7:551099/1‑151 (MQ=255)

                                                                                                        ccATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCg                                                 >  7:488882/1‑151 (MQ=255)

                                                                                                         cATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGt                                                >  8:487552/1‑151 (MQ=255)

                                                                                                         cATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAGTACTGGTTTTCAGCCAGCTGACGATCGt                                                >  7:1207865/1‑151 (MQ=255)

                                                                                                             cccTAATGCAGGAGCACTCAGGCTCGGCCGCATAGGGGGGCGAGGCGCCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                            <  1:19136/150‑1 (MQ=255)

                                                                                                                   ggcaccagcACTCCTCATCTGAAGGATATGGGGGGGAGGCGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGttt                                      <  7:1418073/143‑1 (MQ=255)

                                                                                                                    ccagtcccACTCCTCATCTGAAGGATAGTTGGGGGAGCCGGCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTa                                     <  3:750267/144‑1 (MQ=255)

                                                                                                                     cAGGATCACTCATCATCTGAATGATATTCGGGGGCCCAGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTaa                                    <  1:659972/151‑1 (MQ=255)

                                                                                                                      cgcaccACTCCGCAGCTGAAGGCGATGGGGGGGAGGCGCACGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTaaa                                   <  3:1080901/146‑1 (MQ=255)

                                                                                                                       cgcgcccTCATCCTCTGAAGGATTTTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAAc                                  <  2:721090/145‑1 (MQ=255)

                                                                                                                        gATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCCGCCAGCTGACGATCGTGGCGGCGGTTAAAc                                  >  2:566290/1‑150 (MQ=255)

                                                                                                                            ccTCATCATGGGAAGGATATTGGGGGGCCGCGCCCGCGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                             <  1:1112971/150‑1 (MQ=255)

                                                                                                                            aCTCATCATCTTCATGATTTTTGGGGGCGGTGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                             <  6:1196568/151‑1 (MQ=255)

                                                                                                                               cctcagcTGAAGGATAGTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAg                            <  1:128531/147‑1 (MQ=255)

                                                                                                                                atcatcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGcc                          >  6:15403/1‑150 (MQ=255)

                                                                                                                                 acatcTGATGGTTTTTTGGGGGCGGAGGCCGCGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCaa                        <  4:1035779/150‑1 (MQ=255)

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

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

                                                                                                                                    acTGAAGGATATTTGGGGGAGGAGCACGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCa                     <  5:406936/150‑1 (MQ=255)

                                                                                                                                       gAAGGTTATTGGCGGGCGGCGGCCGGGGGCGGATAACTTTCGGACTAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCa                     <  5:1350878/148‑1 (MQ=255)

                                                                                                                                       gAAGGATCTTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACa                  <  6:935826/151‑1 (MQ=255)

                                                                                                                                       gAAGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAAc                   <  3:732176/150‑1 (MQ=255)

                                                                                                                                            gtttgggggggAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcaca             <  7:661961/145‑1 (MQ=255)

                                                                                                                                            atatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAACCATCCAGCCAATCAACAGcac              >  7:1266093/1‑150 (MQ=255)

                                                                                                                                                ttGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaaca          >  7:1132448/1‑150 (MQ=255)

                                                                                                                                                    gggAGCAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCa     <  6:370885/151‑1 (MQ=255)

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