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	30.2%	588.0 / 100.3	362	intergenic (‑86/+70)	malZ/proY	alpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand): ref base A (192/57); new base C (2/106); total (195/167)
Fisher's exact test for biased strand distribution p-value = 1.40e-45
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 380 total aligned reads displayed.

AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAAC  >  DH1/3448376‑3448675

aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                       <  8:481429/151‑1 (MQ=255)

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

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

  cAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACg                                                                                                                                                     <  6:1196918/151‑1 (MQ=255)

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

   aGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgg                                                                                                                                                    >  3:1454716/1‑151 (MQ=255)

   aGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgg                                                                                                                                                    >  3:1643644/1‑151 (MQ=255)

    gTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggg                                                                                                                                                   <  2:1484126/151‑1 (MQ=255)

     ttGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggg                                                                                                                                                   <  2:498266/150‑1 (MQ=255)

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

              gctttgTTTAACAACTGGGGGCCCCGGCCGGGGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCCCCTTGCCATAGCATAATGCCCGCTCACTCATCATCTGAATGTTCTTTGGGGGCGGCGGCCGGGGGCGGATaa                                                                                                                                          <  6:946547/150‑1 (MQ=255)

              gctttgTTCCACCAATGGGGGCACCGGCAGGTGCCATGCATTAAACAACATAATTCCCCTGTCAAGACAATCTGATCACCTTGCCTGAGCCTACTGCCCGCTCACCCTTCCTCTGACGGATATGGGGGGGCGGCGGCCGGGGCCGGATaa                                                                                                                                          <  6:240444/150‑1 (MQ=255)

                   gTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCg                                                                                                                                    >  7:244999/1‑151 (MQ=255)

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

                     ttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                                   >  3:1056763/1‑150 (MQ=255)

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

                         aagctGGGGGCACGGGCAGGTGACATTAATTTAACATAATAATTCCCCTGTTATTAAAATCTGAGCACCTTGCCATAGCCTAATGCAGGATCACTCATCTTCTGAATGATCTTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAg                                                                                                                              <  1:41581/147‑1 (MQ=255)

                           acggggggCACCTGCAGGTGCAATGAATTTAACATCATAATTCCCCTGTGATGAAAATCTTCTCACCTTGCCATAGCCTAATGCAGGATCACTCTTCATCTGACGTATATTTGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGcc                                                                                                                            <  7:16360/148‑1 (MQ=255)

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

                             tGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCgg                                                                                                                          >  2:171509/1‑151 (MQ=255)

                              gggtgctccgtatgttGCCATGCATTTAAAATAATATTTGCGTTCTTTTGAAAATCTGATCCCATTTCCATACCCTAATGCAGGTTCACTAATTTTTTTATTGATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGa                                                                                                                         <  2:1529317/136‑1 (MQ=255)

                               ggggCACCGGCCGGTGCCATGCATATAACATAATAATTCCACTGTGATTAAAATCTGACCACCTTCCCATAGCCTAATGCAGTATCACTCAACATCTGAAGGATTTTGGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGAt                                                                                                                        <  5:1383912/151‑1 (MQ=255)

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

                                ggcgaCCGGCAGGTGCCATGAATTTAACATCATAATTCCCCTGTTATGAAAATCTGATCACCTTGCCATACCCTAATGCAGGATCCCTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                       <  6:1121554/147‑1 (MQ=255)

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

                                    aCCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAggt                                                                                                                   >  8:1299477/1‑150 (MQ=255)

                                    aCCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAgg                                                                                                                    >  4:113772/1‑150 (MQ=255)

                                       ggCCGGTGCCATGAATTTAACATCATAATTCCCCTGTGATCAAAATCTGATCACCTGCCCCTAGCCTAATCCCCGATCACTCATCATCTGAAGTATATTGGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAttt                                                                                                                <  1:375459/151‑1 (MQ=255)

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

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

                                              gcaattcATTTAACATAATAATTCCCCTGTGATCAAAATCTCACCCCCTTCCCATACCCTAATGCAGCATCACTCATCATCTGAAGGTTTTTTGGGGCCGGCGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTg                                                                                                         <  7:1078283/145‑1 (MQ=255)

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

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

                                                   tcaattAACATCATAATTCCCATGTTATGAAAATCTAATCCCCTTTACATATCCTAATGCAGGAGCACTACTCATCTGAAGGTTATTTGGTGGCGTAGGAGGGGGGCTTATAACTTTCGGAATAGCCGGATCATTATTTCCTCCTTAAcc                                                                                                     <  6:706164/146‑1 (MQ=255)

                                                   gCATTTAACATCAAAATGCCCCTGTGATGAAACTCTGATCACCTTGCCATAGCCTACTGCAGGATCACTCATCATCTGAAGGATATTGGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCg                                                                                                    <  7:784143/151‑1 (MQ=255)

                                                   ccATTTAACACCATACTTCCCCCCTGAACAAAATCACATCCCCTGGCCATCCCCTAATGCAGGCCCACTCAACCTCTCAAGGCTCTGCGGGGCCGGCGGCCGCCCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCg                                                                                                    <  4:1354519/150‑1 (MQ=255)

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

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

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

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

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

                                                      tttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGACCCGGCt                                                                                                 >  7:970605/1‑151 (MQ=255)

                                                      tttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                                 >  8:1098902/1‑151 (MQ=255)

                                                       ttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                                >  5:1629274/1‑151 (MQ=255)

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

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

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

                                                        aaaCATCATAATCCCCCTGTGATCAAAATCTGAGCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTa                                                                                               <  2:1240816/150‑1 (MQ=255)

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

                                                                  aTTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGACCCGGCTTACAGAACGCa                                                                                      >  8:1062244/1‑150 (MQ=255)

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

                                                                     cgccTGTGCTGAAAATCTGCTCACCTTGCCCTCGCCTCCTGCCGGATCACTCCTCATATGCAGTATATGTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAgg                                                                                   <  2:1320315/148‑1 (MQ=255)

                                                                     ccccTGTGATCAAAATCTGATCACCTTCCCATCGCCCAATGCAGGATCACTCATCCTCTGCAGCATCTGTGGCGCCCGCCCCCCCGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGa                                                                                  <  8:1222765/151‑1 (MQ=255)

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

                                                                         tgtgATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTAc                                                                              >  5:462950/1‑151 (MQ=255)

                                                                            gatgaAAATCTCATCACCTTGCCATAGCATAATCCAGCATCACTCCTAATCTTAATGATATTTGGGGGTGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGTATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTgg                                                                           <  5:295207/151‑1 (MQ=255)

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

                                                                                aaaaTCTGATCCCCTTGCCATAGCATAATGCACGATCACTCATCATCTGAAGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                       <  7:1226339/151‑1 (MQ=255)

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

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

                                                                                      ggctCACCTTGCCATCGCCTACTGCAGGAGCACTCATCATCTGAAGGATATTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcca                                                                 <  1:509140/148‑1 (MQ=255)

                                                                                       gATCACCTTGCCATAGCCTAATGCACTAGCACTCCTAATAACAATTATAGTTGGGGGAGCCTGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcag                                                                    <  6:82371/147‑1 (MQ=255)

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

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

                                                                                                gCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAGAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                       >  6:1614928/1‑151 (MQ=255)

                                                                                                gCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                       >  7:1184692/1‑151 (MQ=255)

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

                                                                                                 ccATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGCACGCAAGGATTACTGGTTTTCCGCCAGCTGACGATCg                                                      >  1:629638/1‑151 (MQ=255)

                                                                                                  cATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGt                                                     >  4:124772/1‑151 (MQ=255)

                                                                                                       ccTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcgg                                                >  2:1096211/1‑151 (MQ=255)

                                                                                                            tGCAGGATCTCTGTTCTTCTGTAGGATATTTTGGGTAGGCGGACTGGCGCGTATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGtt                                            <  4:1056921/150‑1 (MQ=255)

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

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

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

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

                                                                                                                     aCTCATCATCTGACGGATATGGGGGGGCGGCGGACGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                                  <  7:732981/151‑1 (MQ=255)

                                                                                                                     aCTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGGTTTCCGCCAGCTGACGATCGTGGCGGCGTTTAAACATccc                                  >  4:1509786/1‑150 (MQ=255)

                                                                                                                       tcatcaTCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGc                                >  7:1396828/1‑151 (MQ=255)

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

                                                                                                                          tcatcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCa                              >  6:1431382/1‑150 (MQ=255)

                                                                                                                          gcatcTGAAGGATATTGGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCaa                             <  7:494429/150‑1 (MQ=255)

                                                                                                                          gcagctcaaGGATATGTGCGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCaa                             <  7:429823/144‑1 (MQ=255)

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

                                                                                                                            atcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATACCCGGATCTGCATTGCGTCGTGAACCGGCTTACAGAACCCAAGGATAACTGGTTTACAGCCAGCTGACGATCGTGGCGGGGTTTAAAAATCCAGCGaagc                           >  4:565780/1‑149 (MQ=255)

                                                                                                                               tGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCAAATaaa                         >  1:851714/1‑150 (MQ=255)

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

                                                                                                                                 aaGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAg                      <  6:1469532/151‑1 (MQ=255)

                                                                                                                                    gATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcac                   <  5:1554649/151‑1 (MQ=255)

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

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

                                                                                                                                           gggggAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAAt             <  3:1109158/150‑1 (MQ=255)

                                                                                                                                             gggAGGTGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCa          <  4:1455429/151‑1 (MQ=255)

                                                                                                                                              ggaggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAc         >  3:820750/1‑151 (MQ=255)

                                                                                                                                                cggaggcCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcgc       <  3:1254125/150‑1 (MQ=255)

                                                                                                                                                aggaggaCTGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcgc       <  5:619515/151‑1 (MQ=255)

                                                                                                                                                  gcggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGa     <  8:1531081/149‑1 (MQ=255)

                                                                                                                                                  gaggcCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGa     <  6:40294/151‑1 (MQ=255)

                                                                                                                                                  gaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcgcg      >  5:919012/1‑150 (MQ=255)

                                                                                                                                                    ggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGaaa   <  8:473942/151‑1 (MQ=255)

                                                                                                                                                     gaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTTGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAAc  <  4:921238/151‑1 (MQ=255)

AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAAC  >  DH1/3448376‑3448675

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