BRESEQ :: Evidence

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

Only 100 of 394 total aligned reads displayed.

ATGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGA  >  DH1/3448373‑3448672

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

  gAGCAGTTGATCTTTTCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCAGCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagg                                                                                                                                                     <  1:1160158/151‑1 (MQ=255)

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

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

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

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

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

       gTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggg                                                                                                                                                >  2:1138542/1‑151 (MQ=255)

          gATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGc                                                                                                                                             >  2:1674688/1‑151 (MQ=255)

               ttgctttgTTAAACAAAGGGGGGCCCCGGCAGGTGCCATGAATTAAACATAATAATTCCCCTGTTATTAAAATCTGATCACCTTCCCATAGCCTAATCCAGTATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGACGGGGGCGGAt                                                                                                                                         <  6:1419566/150‑1 (MQ=255)

                 gctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAAc                                                                                                                                      >  4:1580556/1‑151 (MQ=255)

                 gctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAAc                                                                                                                                      >  6:572640/1‑151 (MQ=255)

                  ctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACt                                                                                                                                     >  7:41180/1‑151 (MQ=255)

                     tgTTTAAAAACTGGGGGCCCCGGCAGGTGCCATTCATTAAACATCATAATTCCCCTGTTATCAAAATCTGATCACCTTGCCATAGCCTAATCCAGGATCACTCATCATCTGAAGGATTTTTGGGGGAGGCGGCCGGGGGCGGATAACttt                                                                                                                                   <  2:380444/150‑1 (MQ=255)

                        ttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCgg                                                                                                                                >  7:1358331/1‑150 (MQ=255)

                        ttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                                >  7:10475/1‑150 (MQ=255)

                          aaCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGCGGGCGGATAACTTTCGGAAt                                                                                                                             >  6:837645/1‑151 (MQ=255)

                           aCAAATGGGGGCACCGGCAGGTGCCATGAATTTAAAATAATAATTCCCCTGTTATTAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATa                                                                                                                            <  1:1007789/151‑1 (MQ=255)

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

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

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

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

                                tGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCgg                                                                                                                       >  7:1488484/1‑151 (MQ=255)

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

                                        ccGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGACTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGAAGGGGGCGGATAACTTTCGGAATAGACGGATCAGGAt                                                                                                               >  3:213497/1‑151 (MQ=255)

                                         cGGCCGGTGCCATTCATTTAACATCATAATGCCCCTGTGATGAAAATCTGACCACCTTGCAATAGCCTAATGCAGGATCACTCATCATCTGAAGGATTTGTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAtt                                                                                                              <  1:277227/151‑1 (MQ=255)

                                         cGGCAGGTGCCATGCATTTAACATAATAAATCCCCTGTTATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACGCATCATCTGCAGGTTATGGGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAtt                                                                                                              <  5:1197797/151‑1 (MQ=255)

                                            cAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcc                                                                                                           >  2:824611/1‑151 (MQ=255)

                                            cAGGGGCCATGCATTAAACATCATAATTCCCCTGTTATCAAAATCTGATCCCCTTCCCATAGCCTACTGCACGATCACTCATCATCTGAAGGATTTTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcc                                                                                                           <  8:1307575/151‑1 (MQ=255)

                                              ggTGCAATGCATTTAACATAATAATTCCCCTTTGATGAAAATCTGATCACCTTCCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTGCGGAATAGCCGGATCAGGATTtcctc                                                                                                         <  7:482210/151‑1 (MQ=255)

                                                  caaTGAATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGa                                                                                                     <  7:18059/149‑1 (MQ=255)

                                                   cattaaTTTAACATCATAATTCCCCTGTGAATAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCTTCTTAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGaa                                                                                                    <  1:904417/146‑1 (MQ=255)

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

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

                                                      gCATATAACATCATAATTCCCGTGTGATGAAACTCTGATCCCCTAGCCATAGCCTACAGCAGGAACACTCATCATCTGAAGGATATGGGGGGGCGGCGCCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAcc                                                                                                  <  7:1441303/150‑1 (MQ=255)

                                                      ccaataaaCATCATAATTCCCCTCTCATCAAAATCTCACCCCCTTCCCCTACCCTAATCCACCATCACTCATCATCTGAAGCTTATCTGCGGGAGGACGCCGCCCCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAcc                                                                                                  <  4:1104968/144‑1 (MQ=255)

                                                      caataaaacaaaataataCCCCTGTGATGAAAATCTGATCCCCTGGCCATAGCCTAATCCAGGATCAATCATCATCTGAAGTATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCg                                                                                                 <  1:1041717/139‑1 (MQ=255)

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

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

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

                                                         aaaaaCAGCATAATCCCCCTGGCATCAAAATCTGACCCCCTTCCCACACCCTACCGCAGCAGCACCCATCCTCTGAAGGATATTGGGGGGCGGCGGCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                              <  7:425901/148‑1 (MQ=255)

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

                                                           tAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                             >  7:1543247/1‑150 (MQ=255)

                                                           aaaCATCATAATACCCCTGAGAGGAAAATCTGATCACCTGCCCATCGCCTAACGCCGCATCACACATCATCTGAAGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTa                                                                                            <  5:1135138/150‑1 (MQ=255)

                                                             aaatcatAGTTCCCCTCAGTTGAAATTCCGATCACCTGGCACTAGCCAATAGCACGATCACTCATTATCTCACTTATTTTGGGGGGCGGCGCCCGGCGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACa                                                                                          <  3:1307543/149‑1 (MQ=255)

                                                             aCATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                           >  3:1059021/1‑150 (MQ=255)

                                                             aCATAATAATTCCCCTGTCATCAAAAGCTGACCCCCTTCCCATACCCTAAGCCAGCACCACTCATCATCTGAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACa                                                                                          <  1:1574197/151‑1 (MQ=255)

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

                                                                 catccctCCCCTCTGATCAAACGCTGACCACCCTCCCATACCCTACTCCACGATCCCTCATCCTCTGACGGATCTTTGGGGGCGGCGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAAc                                                                                      <  8:291178/145‑1 (MQ=255)

                                                                 catAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAAc                                                                                      >  1:941349/1‑151 (MQ=255)

                                                                 catAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGCGGGAGGAGGAAGGGGGCGGATAACTTTCTGAATAGCAGGATCAGGATATCCTCAAGAAGCGGCTTACAGAAc                                                                                      >  6:256595/1‑151 (MQ=255)

                                                                 caatcttCCCCTCTCATCAAAATCTCATCACCTCCCCATACCCTAATGCACGATCCCTCCTCATCTGAAGGATATCTGGCGGAGGCCGCCGCGCCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAAc                                                                                      <  2:601324/146‑1 (MQ=255)

                                                                      ttCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                   >  2:1424847/1‑149 (MQ=255)

                                                                       tCCCCTGCGATCAAAATCTCACCCCCTTGCCATCCACTACTGCCGGCTCACTCCTCCTCTGAAGGATCTGTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAg                                                                                 <  6:1169959/150‑1 (MQ=255)

                                                                        ccccTGTGATGAAAATCTTAACCCCTTGCCATAGCCTAATTCAGGATCACTCATCATCTTAAGGATTTTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGa                                                                               <  4:26013/151‑1 (MQ=255)

                                                                            tgtgATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCATCCTGAACCGGCTTACAGAACGCAAGGATTa                                                                            >  1:483056/1‑150 (MQ=255)

                                                                              tgatgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTg                                                                         >  8:1345855/1‑151 (MQ=255)

                                                                                atgaTAATCTGAACCACTTGCCATAGCCTAATGCAGGACCACTGATCATGTGATGCATTTTTGGGGGCGGCGGGGGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGt                                                                       <  7:1128331/151‑1 (MQ=255)

                                                                                atgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGt                                                                       >  3:856990/1‑151 (MQ=255)

                                                                                   aaaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                    >  4:1356384/1‑151 (MQ=255)

                                                                                    aaaTCTGATCACCTTGCCATAGCCTAATGCCGGATCACTCATCATCTGAATGATATTGGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                   <  2:1030982/151‑1 (MQ=255)

                                                                                     aaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                   >  4:635871/1‑150 (MQ=255)

                                                                                      aTCTGATCACCTTGCCATAGCCTAAAGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGCATTACTGCTTTTc                                                                   >  2:250750/1‑149 (MQ=255)

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

                                                                                          cagccccTCGCCATAGCCTAATGCAGGAGCCCTCAGCATCTGATGGATCTTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                               <  7:611967/144‑1 (MQ=255)

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

                                                                                                  tGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGAt                                                     >  2:942162/1‑151 (MQ=255)

                                                                                                   gCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGCCGATc                                                    >  2:356540/1‑151 (MQ=255)

                                                                                                    ccATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCg                                                   <  1:1249432/151‑1 (MQ=255)

                                                                                                      aTAGCCTAATGCAGGATCACTCATCATCTGAAGGTTATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                                 <  5:1494102/151‑1 (MQ=255)

                                                                                                      aTAGCCTAATCCAGGATCACTCAGCATCTGAAGGATATTGGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                                 <  5:277998/151‑1 (MQ=255)

                                                                                                       tAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTgg                                                >  1:580877/1‑151 (MQ=255)

                                                                                                       aaGCCTAATCCAGGATCACTCATCATCTGAAGTATATTTGGGGCAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTgg                                                <  2:459889/150‑1 (MQ=255)

                                                                                                        aGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggc                                               >  7:284639/1‑151 (MQ=255)

                                                                                                             aaTGCAGCATCACTCATCATCTGATGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGgcggcg                                           <  8:815962/150‑1 (MQ=255)

                                                                                                              aTGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGGACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGtt                                         >  5:1681737/1‑151 (MQ=255)

                                                                                                               cgCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGttt                                        <  5:647271/150‑1 (MQ=255)

                                                                                                                gCAGGATCCCTCATCATATGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAAATTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGATGACGATCTTGGCGGCGttt                                        >  4:1249550/1‑150 (MQ=255)

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

                                                                                                                      tCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATc                                 >  6:911338/1‑151 (MQ=255)

                                                                                                                      tCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAACGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATc                                 >  1:32817/1‑151 (MQ=255)

                                                                                                                      gcACTCATCATCTGCAGGATTTTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATc                                 <  8:101830/150‑1 (MQ=255)

                                                                                                                        aCTCATCATCTGAAGGATATTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                               <  7:1261091/151‑1 (MQ=255)

                                                                                                                            atcatcTGAATGATATTTGGGGGCGCCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCa                           <  6:966970/151‑1 (MQ=255)

                                                                                                                            atcatcTGAACGATATTGGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGGTTAAACATCCAGCCa                           <  3:327107/151‑1 (MQ=255)

                                                                                                                             acatcTGAAGGAGATTGGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCaa                          <  5:260574/150‑1 (MQ=255)

                                                                                                                              catcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCATCCAGCTGACGATCCTGGCGGCGTTTAAACATCCAGCCAAt                         >  5:1597111/1‑151 (MQ=255)

                                                                                                                              catcTGAAGGATATTTGGGGGAGGAGGACGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCaa                          <  4:40507/150‑1 (MQ=255)

                                                                                                                                   gtagtatatTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAAc                     <  2:737604/145‑1 (MQ=255)

                                                                                                                                   gcaggctcgttGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCa                       <  4:1620998/139‑1 (MQ=255)

                                                                                                                                   gAAGGATATTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCaa                      <  8:418942/149‑1 (MQ=255)

                                                                                                                                    aaGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAAc                     >  8:207749/1‑149 (MQ=255)

                                                                                                                                      ggATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGc                  >  7:703956/1‑150 (MQ=255)

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

                                                                                                                                             ggggggCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAAt          <  3:359625/150‑1 (MQ=255)

                                                                                                                                              gggggAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGAACGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATc         <  4:855710/151‑1 (MQ=255)

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

                                                                                                                                               ggggAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATcc        <  2:1274423/151‑1 (MQ=255)

                                                                                                                                               ggggAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATcc        <  3:418452/151‑1 (MQ=255)

                                                                                                                                                  gaggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcg     >  3:561348/1‑151 (MQ=255)

ATGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGA  >  DH1/3448373‑3448672

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