BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	DH1	3,668,924	(T)_5→6	pseudogene (353/455 nt)	lpxA ←	hydroxymyristoyl‑ACP dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	3,668,919	1	.	T	98.1%	1127.0 / 8.9	315	pseudogene (358/455 nt)	lpxA	hydroxymyristoyl‑ACP dehydratase
Reads supporting (aligned to +/- strand): ref base . (4/2); new base T (146/162); total (151/164)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 346 total aligned reads displayed.

CCTCTTCCACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑TTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGAATACCTGTTGCCTGT  >  DH1/3668770‑3669064

ccTCTTCCACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑t                                                                                                                                                   >  2:172253/1‑151 (MQ=255)

     tCCACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑tttttt                                                                                                                                              >  8:1268681/1‑150 (MQ=255)

     tCCACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑ttttt                                                                                                                                               <  8:547479/150‑1 (MQ=255)

     tCCACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑ttttt                                                                                                                                               <  8:367749/150‑1 (MQ=255)

     tCCACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑tttt                                                                                                                                                <  8:1044665/149‑1 (MQ=255)

      ccACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑ttttttc                                                                                                                                             <  1:474424/151‑3 (MQ=255)

      ccACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑ttttttc                                                                                                                                             >  4:592798/1‑149 (MQ=255)

      ccACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑ttttttc                                                                                                                                             <  3:212247/151‑3 (MQ=255)

      ccACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑ttttttc                                                                                                                                             >  2:591334/1‑149 (MQ=255)

        acaATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑ttttttcga                                                                                                                                           <  1:591225/151‑5 (MQ=255)

        acaATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑ttttttcga                                                                                                                                           <  5:685698/151‑5 (MQ=255)

        acaATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑ttttttcg                                                                                                                                            >  4:296304/1‑147 (MQ=255)

             ggCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGCGG‑TTTTT‑CCAAAGTg                                                                                                                                      >  7:290837/1‑151 (MQ=255)

             ggCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTA‑CATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGa                                                                                                                                     <  8:1121502/151‑1 (MQ=255)

                    ggATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGGTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTcc                                                                                                                                <  8:859346/150‑1 (MQ=255)

                    ggATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTcc                                                                                                                                <  3:876525/150‑1 (MQ=255)

                    ggATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCa                                                                                                                               <  4:407432/151‑1 (MQ=255)

                      aTGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATg                                                                                                                             <  1:763570/151‑1 (MQ=255)

                       tGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATg                                                                                                                             >  3:990247/1‑150 (MQ=255)

                        gCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGa                                                                                                                            <  1:1017332/150‑1 (MQ=255)

                         cacaAAGGCGGATTTATCAATCACTTATCAGGCCTCCCGGGTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGAAACCCCTTTAAAACGGGTCAGGCCGCCGCGC‑G‑TTTTTACGAAAGTGACATCCa                                                                                                                               >  3:116151/1‑146 (MQ=255)

                           caAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCAt                                                                                                                        >  7:325994/1‑151 (MQ=255)

                           aaaaGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCAt                                                                                                                        <  3:587374/150‑1 (MQ=255)

                                 cGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATc                                                                                                                  >  6:483418/1‑151 (MQ=255)

                                 cGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATc                                                                                                                  >  6:690093/1‑151 (MQ=255)

                                   gATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGc                                                                                                                <  4:447368/151‑1 (MQ=255)

                                        atcaatcaCGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGc                                                                                                           <  2:904242/151‑1 (MQ=255)

                                              caCGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGAcc                                                                                                     <  4:849900/151‑1 (MQ=255)

                                               aCGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGCCCg                                                                                                    >  7:1233142/1‑151 (MQ=255)

                                                cGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCgg                                                                                                   >  6:1167394/1‑151 (MQ=255)

                                                 gTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCggg                                                                                                  <  4:915517/151‑1 (MQ=255)

                                                      aGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCtt                                                                                             >  5:1087152/1‑151 (MQ=255)

                                                       ggCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTg                                                                                            <  5:1165793/151‑1 (MQ=255)

                                                        gCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTg                                                                                            >  4:831571/1‑150 (MQ=255)

                                                           tCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAgc                                                                                        >  4:817163/1‑151 (MQ=255)

                                                            cccGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAgc                                                                                        <  8:947776/150‑1 (MQ=255)

                                                            cccGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAg                                                                                         <  3:426186/149‑1 (MQ=255)

                                                              cGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGCCCGGGCGCTTGAAgcgcg                                                                                     >  7:963463/1‑151 (MQ=255)

                                                              cGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAgcgcg                                                                                     <  8:501764/151‑1 (MQ=255)

                                                                gCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCt                                                                                   <  1:856147/151‑1 (MQ=255)

                                                                   aCGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCtt                                                                                  <  2:357382/149‑1 (MQ=255)

                                                                   aCGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTc                                                                                 <  8:143879/150‑1 (MQ=255)

                                                                   aCGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTc                                                                                 <  3:1292013/150‑1 (MQ=255)

                                                                           acaTCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATAcc                                                                        <  4:1170613/151‑1 (MQ=255)

                                                                           acaTCATCGTTGCTTCCCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATAcc                                                                        <  5:1222746/151‑1 (MQ=255)

                                                                             atcatcGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAg                                                                      <  3:808992/151‑1 (MQ=255)

                                                                               catcGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGc                                                                     <  6:1332754/150‑1 (MQ=255)

                                                                                  cGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAag                                                                 <  8:422155/151‑1 (MQ=255)

                                                                                    ttGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAagt                                                                >  1:116191/1‑150 (MQ=255)

                                                                                      gCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCCCCACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAagtagt                                                             >  7:488403/1‑151 (MQ=255)

                                                                                      gCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAagtag                                                              <  2:755359/150‑1 (MQ=255)

                                                                                       cTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAagtagt                                                             >  3:372619/1‑150 (MQ=255)

                                                                                         tcccAAACTACTTTACCACCGACCAGAGCAACCCCTTAAAAACGGGTCCGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAgtagta                                                            <  6:300442/146‑1 (MQ=255)

                                                                                          cGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACa                                                          <  8:1066829/150‑1 (MQ=255)

                                                                                          cGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACa                                                          <  6:925660/150‑1 (MQ=255)

                                                                                          cGCAAACTACTTTACCATCGACCAGACCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGATTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACa                                                          <  5:849666/150‑1 (MQ=255)

                                                                                            cAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCt                                                       >  5:917795/1‑151 (MQ=255)

                                                                                             aaaCTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTc                                                      <  4:941659/151‑1 (MQ=255)

                                                                                             aaaCTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTc                                                      >  1:961629/1‑151 (MQ=255)

                                                                                                ctactTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑GGTAAGTGCCTTCCATGTTGAATTTATCGCCACGCCCGCCCGCGCGCATGGAGCCCGCGTCGGCAAAACCCGCGAAGTACGTCCGCTCAcc                                                   >  7:609288/1‑151 (MQ=255)

                                                                                                  actTTACCCTCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCAccc                                                  <  2:977201/150‑1 (MQ=255)

                                                                                                  actTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCAccc                                                  <  5:342918/150‑1 (MQ=255)

                                                                                                  actTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCAccc                                                  <  4:615999/150‑1 (MQ=255)

                                                                                                  actTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCAcc                                                   <  8:1205985/149‑1 (MQ=255)

                                                                                                  actTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCg                                                 >  3:24756/1‑151 (MQ=255)

                                                                                                  actTTACCATCGACCAGAGCAACCCCTATAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGAGAAGTAGTACAGCTCACCCg                                                 >  6:863215/1‑151 (MQ=255)

                                                                                                    tttACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCggg                                               >  7:694398/1‑150 (MQ=255)

                                                                                                         cATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTcc                                            >  1:297120/1‑149 (MQ=255)

                                                                                                          aTCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGt                                         <  8:793301/151‑1 (MQ=255)

                                                                                                          aTCGACCAGAGCAACCCCTTTAAAACGGGCCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGt                                         <  7:547336/151‑1 (MQ=255)

                                                                                                             gACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGCCCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGtttt                                      >  2:1140193/1‑151 (MQ=255)

                                                                                                             gACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGCAGCGCGCGTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGtttt                                      >  3:1077106/1‑151 (MQ=255)

                                                                                                             gACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGtttt                                      <  4:166999/151‑1 (MQ=255)

                                                                                                              aCCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGTAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTc                                     >  3:717729/1‑151 (MQ=255)

                                                                                                              aCCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTc                                     >  5:255380/1‑151 (MQ=255)

                                                                                                               ccAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTcc                                    <  7:878896/151‑1 (MQ=255)

                                                                                                               ccAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTcc                                    <  7:769728/151‑1 (MQ=255)

                                                                                                                     cAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTACTACGCt                              >  2:261021/1‑151 (MQ=255)

                                                                                                                       aCCCCTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCttt                            >  2:521605/1‑151 (MQ=255)

                                                                                                                        ccccTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCtttt                           <  2:1063692/151‑1 (MQ=255)

                                                                                                                          ccTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCCTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCACTTTTCCTACGCTTTTa                          >  5:806022/1‑150 (MQ=255)

                                                                                                                           cTTTAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTaa                         >  8:1327919/1‑150 (MQ=255)

                                                                                                                            tttAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAAc                       >  2:833667/1‑151 (MQ=255)

                                                                                                                             ttAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTaaa                        <  6:528074/149‑1 (MQ=255)

                                                                                                                             ttAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAAc                       <  1:172015/150‑1 (MQ=255)

                                                                                                                              tAAAACGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGc                     >  2:682475/1‑151 (MQ=255)

                                                                                                                               aaaaCGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGc                     <  5:899857/150‑1 (MQ=255)

                                                                                                                                aaaCGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCa                   <  8:676822/151‑1 (MQ=255)

                                                                                                                                 aaCGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTGCCCGGTTTCCTACGCTTTTAAACGcc                    >  6:1144341/1‑149 (MQ=255)

                                                                                                                                   cGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGaa                <  7:1112078/151‑1 (MQ=255)

                                                                                                                                   cGGGTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGa                 >  7:1028338/1‑150 (MQ=255)

                                                                                                                                    gggTCAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGAAt               <  1:1203920/151‑1 (MQ=255)

                                                                                                                                        cAGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTGCTACGCTTTTAAACGCCAGAATACCt           >  8:268407/1‑151 (MQ=255)

                                                                                                                                         aGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTGACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGAATACCTg          >  3:620235/1‑151 (MQ=255)

                                                                                                                                         aGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGAATACCTg          <  6:812308/151‑1 (MQ=255)

                                                                                                                                         aGGCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGAATACCTg          <  4:710840/151‑1 (MQ=255)

                                                                                                                                          ggCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGAATACCTGt         <  1:526759/151‑1 (MQ=255)

                                                                                                                                          ggCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCCCGCTTCGTCAATACCAGCGAAGTAGTACACCTCACCCGGATCCAGTTTTACTACGCTTTTAAACCCCAGAATAACTGt         >  8:609277/1‑151 (MQ=255)

                                                                                                                                           gCCGCGGCGC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGAATACCTGtt        <  1:425551/151‑1 (MQ=255)

                                                                                                                                                ggcgC‑GTTTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGAATACCTGTTGCCtg   <  5:622004/151‑1 (MQ=255)

CCTCTTCCACAATGGCGGTTGGATGCACAAAGGCGGATTTATCAATCACGTATCAGGCCTCCCGGCTACGAGCACACATCATCGTTGCTTCGCAAACTACTTTACCATCGACCAGAGCAACCCCTTTAAAACGGGTCAGGCCGCGGCGC‑G‑TTTTT‑CGAAAGTGACTTCCATGATCATTTGATCGCCAGGCACGACCGGGCGCTTGAAGCGCGCTTCGTCAATACCAGCGAAGTAGTACAGCTCACCCGGTTCCAGTTTTCCTACGCTTTTAAACGCCAGAATACCTGTTGCCTGT  >  DH1/3668770‑3669064

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