BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	DH1	2,688,963	A→C	K404T (AAG→ACG)	purB →	adenylosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	2,688,963	0	A	C	99.0%	1590.4 / 0.8	412	K404T (AAG→ACG)	purB	adenylosuccinate lyase
Reads supporting (aligned to +/- strand): ref base A (2/2); new base C (213/194); total (215/197)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 412 total aligned reads displayed.

ATTGCATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGAAGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATGGTTGA  >  DH1/2688815‑2689110

aTTGCATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGc                                                                                                                                                   <  8:284206/151‑1 (MQ=255)

 ttGCATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGACCCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGcg                                                                                                                                                  >  8:1076870/1‑150 (MQ=255)

  tGCATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTg                                                                                                                                                 >  3:1534625/1‑151 (MQ=255)

    cATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGCACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGa                                                                                                                                                >  7:1056591/1‑150 (MQ=255)

    cATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGa                                                                                                                                                >  8:402927/1‑150 (MQ=255)

    cATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGa                                                                                                                                                >  3:1390788/1‑150 (MQ=255)

       atCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAaga                                                                                                                                            <  4:521618/151‑1 (MQ=255)

         cAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGc                                                                                                                                          >  6:1103414/1‑151 (MQ=255)

          aaTCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCt                                                                                                                                         >  6:204703/1‑151 (MQ=255)

                cccTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTcg                                                                                                                                   >  8:7805/1‑151 (MQ=255)

                cccTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTcg                                                                                                                                   >  5:354040/1‑151 (MQ=255)

                cccTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTcg                                                                                                                                   <  6:1244359/151‑1 (MQ=255)

                 ccTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTcgc                                                                                                                                  >  5:1613339/1‑151 (MQ=255)

                      aaGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTaa                                                                                                                             <  4:1481892/151‑1 (MQ=255)

                          cGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGCACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAgcgc                                                                                                                         >  3:713028/1‑151 (MQ=255)

                           gTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGcgcg                                                                                                                        >  5:1681987/1‑151 (MQ=255)

                             gAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGtt                                                                                                                      <  5:960358/151‑1 (MQ=255)

                             gAGCAAACTGGAAGGGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGt                                                                                                                       <  1:1012973/150‑1 (MQ=255)

                              aGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTg                                                                                                                     <  4:4165/151‑1 (MQ=255)

                              aGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTg                                                                                                                     <  7:846036/151‑1 (MQ=255)

                               gCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGAAGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTg                                                                                                                     >  7:750259/1‑150 (MQ=255)

                                cAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGAc                                                                                                                   <  3:1105819/151‑1 (MQ=255)

                                 aaaCTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACg                                                                                                                  >  1:379528/1‑151 (MQ=255)

                                   aCTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGcc                                                                                                                <  6:330425/151‑1 (MQ=255)

                                      ggAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCg                                                                                                               <  1:1464998/149‑1 (MQ=255)

                                      ggAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGa                                                                                                              <  3:1420811/150‑1 (MQ=255)

                                       gAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGATGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGCCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAg                                                                                                            <  4:629046/151‑1 (MQ=255)

                                       gAAGTGAACCGTGACCATCTGCTGGAGGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAg                                                                                                            <  8:655457/151‑1 (MQ=255)

                                         aGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGc                                                                                                          <  7:205335/151‑1 (MQ=255)

                                         aGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGc                                                                                                          >  2:645953/1‑151 (MQ=255)

                                         aGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGc                                                                                                          >  2:41013/1‑151 (MQ=255)

                                         aGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGc                                                                                                          >  8:1288597/1‑151 (MQ=255)

                                                cGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGCAGc                                                                                                   >  5:259019/1‑151 (MQ=255)

                                                cGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGc                                                                                                   <  6:946518/151‑1 (MQ=255)

                                                      cATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTa                                                                                             <  5:791709/151‑1 (MQ=255)

                                                          tgctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCg                                                                                          >  2:1268008/1‑150 (MQ=255)

                                                           gctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCg                                                                                          <  5:642265/149‑1 (MQ=255)

                                                           gctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCg                                                                                          <  7:27230/149‑1 (MQ=255)

                                                           gctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGa                                                                                         <  3:828820/150‑1 (MQ=255)

                                                           gctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGAt                                                                                        >  2:1414101/1‑151 (MQ=255)

                                                            ctgGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATg                                                                                       >  7:1480898/1‑151 (MQ=255)

                                                              ggATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATg                                                                                       <  7:1621584/149‑1 (MQ=255)

                                                              ggATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATg                                                                                       <  1:652243/149‑1 (MQ=255)

                                                               gATGAACTGGATCACAACTGGGATGTGCTGTCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTc                                                                                    <  2:851221/151‑1 (MQ=255)

                                                                aTGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCt                                                                                   >  6:1352770/1‑151 (MQ=255)

                                                                 tGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTg                                                                                  >  6:65359/1‑151 (MQ=255)

                                                                  gAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTgg                                                                                 >  4:858926/1‑151 (MQ=255)

                                                                  gAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTgg                                                                                 <  4:1132668/151‑1 (MQ=255)

                                                                      tGGATCACAAATGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCTAAAAACCTAACGAGACGCTGAAAGAGCTGACTCGCGGTAAGAGCGTTGACGCAGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGtt                                                                             >  3:626981/1‑151 (MQ=255)

                                                                       ggATCACACCTGGGAACTGCTGGCTGAACCAAACCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAACCGCGTTGCCGCCGAAGGCATGAAGCAGTTTATCGTTGGTCTGGAGtgc                                                                            >  6:57519/1‑149 (MQ=255)

                                                                       ggATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGtt                                                                             <  8:537110/150‑1 (MQ=255)

                                                                        gATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGc                                                                           >  6:996540/1‑151 (MQ=255)

                                                                           cacaACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGcc                                                                          >  1:798625/1‑149 (MQ=255)

                                                                             caACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCaga                                                                       >  7:1294173/1‑150 (MQ=255)

                                                                              aaCTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCagaag                                                                     <  1:1154028/151‑1 (MQ=255)

                                                                              aaCTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCagaag                                                                     <  5:836656/151‑1 (MQ=255)

                                                                               aCTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCaga                                                                       >  5:1128255/1‑148 (MQ=255)

                                                                                cTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGttt                                                                                              >  7:457911/1‑124 (MQ=255)

                                                                                  ggtttgTGCTGGCTGAACCAATCCAGAGTGTTATGCGTCGCTATTGCATCGAAAAACCGTTCGAGACGCTGAAAGAGCTGACTCTCTGTAACCGCGTTGTCGCCGAATGCATGAATCAGTTTATCGATGGTCTGGCGTTGCCagaagaaga                                                                 <  8:576825/146‑1 (MQ=255)

                                                                                     aaGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAgaga                                                               >  4:1200833/1‑150 (MQ=255)

                                                                                     aaGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGaa                                                              <  5:886742/151‑1 (MQ=255)

                                                                                         gctggctgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGcc                                                          >  7:968180/1‑151 (MQ=255)

                                                                                          ctggctgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGccc                                                         >  2:373640/1‑151 (MQ=255)

                                                                                           tggctgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGcccc                                                        >  3:1003424/1‑150 (MQ=255)

                                                                                           tggctgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGCTCTGGCGTTGCCAGAAGAAGAGAAAGGCCg                                                        >  4:1088893/1‑151 (MQ=255)

                                                                                            ggcttAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGc                                                       <  1:978711/151‑1 (MQ=255)

                                                                                            ggctgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGccc                                                         <  7:361675/149‑1 (MQ=255)

                                                                                              ctgAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCt                                                     >  8:1195021/1‑151 (MQ=255)

                                                                                                 aaCCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGa                                                   >  3:267251/1‑150 (MQ=255)

                                                                                                  aCCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGaa                                                  >  1:1326267/1‑150 (MQ=255)

                                                                                                     aaTCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCg                                              <  5:1457187/151‑1 (MQ=255)

                                                                                                            aCAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGc                                       >  4:1009912/1‑151 (MQ=255)

                                                                                                            aCAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGc                                       >  2:819033/1‑151 (MQ=255)

                                                                                                             cAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGcc                                      >  3:1654496/1‑151 (MQ=255)

                                                                                                             cAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGcc                                      >  6:1514583/1‑151 (MQ=255)

                                                                                                              aGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCg                                     >  5:1372352/1‑151 (MQ=255)

                                                                                                                 tATGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCt                                  >  1:1390212/1‑151 (MQ=255)

                                                                                                                  aTGCGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACCCCGGCTa                                 <  2:311501/151‑1 (MQ=255)

                                                                                                                     cGTCGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACt                              <  2:1516561/151‑1 (MQ=255)

                                                                                                                       gcGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACtat                            <  7:664282/150‑1 (MQ=255)

                                                                                                                        cGCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACtat                            >  8:443535/1‑150 (MQ=255)

                                                                                                                         gCTATGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTatat                          <  1:45870/151‑1 (MQ=255)

                                                                                                                            aTGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTgg                       <  4:1323375/151‑1 (MQ=255)

                                                                                                                            aTGGCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTgg                       >  8:310535/1‑151 (MQ=255)

                                                                                                                              ggCATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTc                     <  3:470638/151‑1 (MQ=255)

                                                                                                                                cATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCgc                   >  8:202027/1‑150 (MQ=255)

                                                                                                                                cATCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGa                   <  6:145286/151‑1 (MQ=255)

                                                                                                                                 aTCGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAg                  <  6:76080/151‑1 (MQ=255)

                                                                                                                                   cGAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCt                <  1:1660574/151‑1 (MQ=255)

                                                                                                                                    gAAAAACCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCt                >  2:1149171/1‑150 (MQ=255)

                                                                                                                                      aaaaCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCCTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATc             >  8:1670913/1‑151 (MQ=255)

                                                                                                                                       aaaCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATc             >  5:924324/1‑150 (MQ=255)

                                                                                                                                       aaaCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCa            <  4:209655/151‑1 (MQ=255)

                                                                                                                                        aaCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCAc           <  1:222630/151‑1 (MQ=255)

                                                                                                                                        aaCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCAc           >  7:560110/1‑151 (MQ=255)

                                                                                                                                        aaCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCAc           <  6:1613000/151‑1 (MQ=255)

                                                                                                                                         aCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACg          <  8:1441031/151‑1 (MQ=255)

                                                                                                                                         aCCGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACg          >  5:579846/1‑151 (MQ=255)

                                                                                                                                           cGTACGAGACGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGAt        >  2:179326/1‑151 (MQ=255)

                                                                                                                                                 agaCGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGCAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATGGtgg   >  8:1507346/1‑148 (MQ=255)

ATTGCATATCAATCCACCCTGAAAGGCGTGAGCAAACTGGAAGTGAACCGTGACCATCTGCTGGATGAACTGGATCACAACTGGGAAGTGCTGGCTGAACCAATCCAGACAGTTATGCGTCGCTATGGCATCGAAAAACCGTACGAGAAGCTGAAAGAGCTGACTCGCGGTAAGCGCGTTGACGCCGAAGGCATGAAGCAGTTTATCGATGGTCTGGCGTTGCCAGAAGAAGAGAAAGCCCGCCTGAAAGCGATGACGCCGGCTAACTATATTGGTCGAGCTATCACGATGGTTGA  >  DH1/2688815‑2689110

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