BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,238,518	A→T	5.3%	L74Q (CTG→CAG)	catA‑I ←	catechol 1,2‑dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,238,518	0	A	T	5.3%	576.5 / 14.4	170	L74Q (CTG→CAG)	catA‑I	catechol 1,2‑dioxygenase
Reads supporting (aligned to +/- strand): ref base A (71/90); new base T (4/5); total (75/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 188 total aligned reads displayed.

ATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGTTTCCTTCG  >  NC_002947/4238371‑4238667

  gcgcgcTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGccagcagccagcag                                                                                                                                                   <  1:2262124/150‑1 (MQ=255)

     cgcTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGgcc                                                                                                                                                >  1:324418/1‑150 (MQ=255)

      gcTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGgcct                                                                                                                                               >  2:250838/1‑150 (MQ=255)

        ttCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGgcctgc                                                                                                                                             <  2:1502884/150‑1 (MQ=255)

         tCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGgcctgcc                                                                                                                                            <  2:3868700/150‑1 (MQ=255)

            cccTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCg                                                                                                                                         <  2:1467399/150‑1 (MQ=255)

                 ggCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACg                                                                                                                                    <  1:2297022/150‑1 (MQ=255)

                  gCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACgg                                                                                                                                   >  1:660474/1‑150 (MQ=255)

                  gCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACgg                                                                                                                                   <  2:2824964/150‑1 (MQ=255)

                    cAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGgcc                                                                                                                                 <  1:101161/150‑1 (MQ=255)

                       cggcgCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGgccgcc                                                                                                                              <  2:3988613/150‑1 (MQ=255)

                           gcCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAgg                                                                                                                          >  1:1486565/1‑150 (MQ=255)

                            ccccGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGc                                                                                                                         <  1:5094455/150‑1 (MQ=255)

                             cccGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCg                                                                                                                        <  1:615385/150‑1 (MQ=255)

                                ggCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGtt                                                                                                                     <  2:4818372/150‑1 (MQ=255)

                                 gCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTg                                                                                                                    <  2:3580441/150‑1 (MQ=255)

                                  cAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGa                                                                                                                   <  2:2961552/150‑1 (MQ=255)

                                   aaCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                  <  2:91806/150‑1 (MQ=255)

                                   aaCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                  <  1:4746980/150‑1 (MQ=255)

                                   aaCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                  <  2:1681681/150‑1 (MQ=255)

                                   aaCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                  <  2:3035731/150‑1 (MQ=255)

                                   aaCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                  <  2:3845824/150‑1 (MQ=255)

                                   aaCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                  <  2:617319/150‑1 (MQ=255)

                                       tACAACGGGCCGTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAg                                                                                                              <  2:4359662/150‑1 (MQ=255)

                                            cGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGAc                                                                                                         >  1:3498290/1‑150 (MQ=255)

                                                 cTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGTCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGt                                                                                                    <  1:3171332/150‑1 (MQ=255)

                                                 cTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGt                                                                                                    >  1:587923/1‑150 (MQ=255)

                                                  ttCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAgg                                                                                                                 <  2:5320622/136‑1 (MQ=255)

                                                  ttCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAgg                                                                                                                 >  1:5320772/1‑136 (MQ=255)

                                                  ttCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTg                                                                                                   >  2:3510033/1‑150 (MQ=255)

                                                     gATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTTGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCa                                                                                                <  1:624782/150‑1 (MQ=255)

                                                      aTGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAg                                                                                               <  2:2790714/150‑1 (MQ=255)

                                                       tGGTGCGCGGGGTGCCGCCGCCAAGGCCTGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGa                                                                                              >  2:525109/1‑150 (MQ=255)

                                                        ggTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGaa                                                                                             <  2:673160/150‑1 (MQ=255)

                                                          tGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCa                                                                                                <  2:457771/145‑1 (MQ=255)

                                                          tGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCa                                                                                                >  1:457771/1‑145 (MQ=255)

                                                           gcgcGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTc                                                                                          <  2:4324330/150‑1 (MQ=255)

                                                           gcgcGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTc                                                                                          <  1:91181/150‑1 (MQ=255)

                                                             gcgGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGt                                                                                        >  2:462060/1‑150 (MQ=255)

                                                              cgGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTc                                                                                       <  1:321783/150‑1 (MQ=255)

                                                              aaggatGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                  <  2:2162876/118‑1 (MQ=255)

                                                              aaggatGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                  >  1:2162876/6‑123 (MQ=255)

                                                                  gTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCg                                                                                   >  2:1349241/1‑150 (MQ=255)

                                                                  gTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCg                                                                                   <  1:373975/150‑1 (MQ=255)

                                                                  gTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCg                                                                                   <  1:5419093/150‑1 (MQ=255)

                                                                     ccgccgccAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGc                                                                                                  >  2:5121919/1‑132 (MQ=255)

                                                                     ccgccgccAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGc                                                                                                  <  1:5121919/132‑1 (MQ=255)

                                                                     ccgccgccAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTa                                                                                <  2:5196668/150‑1 (MQ=255)

                                                                       gccgccAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAAt                                                                              <  1:4659792/150‑1 (MQ=255)

                                                                        ccgccAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAAtc                                                                             >  2:1970344/1‑150 (MQ=255)

                                                                          gccAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGATGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGt                                                                                                    <  2:35365/125‑1 (MQ=255)

                                                                          gccAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGATGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGt                                                                                                    >  1:35365/1‑125 (MQ=255)

                                                                            cAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCa                                                                         <  1:1908723/150‑1 (MQ=255)

                                                                            cAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCa                                                                         <  1:2193122/150‑1 (MQ=255)

                                                                             aaGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAg                                                                        >  1:1591565/1‑150 (MQ=255)

                                                                             aaGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAg                                                                        <  1:4711835/150‑1 (MQ=255)

                                                                               ggCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGt                                                                      >  2:4950654/1‑150 (MQ=255)

                                                                                gCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTc                                                                     >  1:5282646/1‑150 (MQ=255)

                                                                                 ccGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCt                                                                    <  2:3208008/150‑1 (MQ=255)

                                                                                  cGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCtt                                                                   <  1:1131044/150‑1 (MQ=255)

                                                                                  cGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCtt                                                                   <  1:4437331/150‑1 (MQ=255)

                                                                                    gCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCg                                                                 <  1:5055683/150‑1 (MQ=255)

                                                                                     cTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGa                                                                <  2:1486565/150‑1 (MQ=255)

                                                                                     cTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGa                                                                <  2:2685636/150‑1 (MQ=255)

                                                                                        cGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCa                                                             >  2:3213162/1‑150 (MQ=255)

                                                                                        cGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCa                                                             <  2:2370706/150‑1 (MQ=255)

                                                                                           cATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAggc                                                          <  1:2894478/150‑1 (MQ=255)

                                                                                           cATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAggc                                                          <  1:2400724/150‑1 (MQ=255)

                                                                                           cATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAggc                                                          <  1:2976986/150‑1 (MQ=255)

                                                                                                ttGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAggcgggcg                                                     >  1:3689789/1‑150 (MQ=255)

                                                                                                ttGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAggcgggcg                                                     >  1:4681106/1‑150 (MQ=255)

                                                                                                 tGGCATCCTGTAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGgcgggcgg                                                    <  1:3882339/150‑1 (MQ=255)

                                                                                                  ggCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAggcgtggtga                                                   <  2:1150072/150‑7 (MQ=255)

                                                                                                  ggCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGgt                                                   <  1:4343798/150‑1 (MQ=255)

                                                                                                     aTCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTc                                                >  1:3585322/1‑150 (MQ=255)

                                                                                                      tCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCt                                               >  1:2455951/1‑150 (MQ=255)

                                                                                                       ccTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGTTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCtt                                              <  1:5155380/150‑1 (MQ=255)

                                                                                                       ccTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCtt                                              >  1:1451296/1‑150 (MQ=255)

                                                                                                        cTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTg                                             >  1:3384841/1‑150 (MQ=255)

                                                                                                                 gTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCa                                                             >  1:2693563/1‑125 (MQ=255)

                                                                                                                 gTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCAcgcg                                    >  2:1619628/1‑150 (MQ=255)

                                                                                                                 gTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCAcgcg                                    >  1:1697584/1‑150 (MQ=255)

                                                                                                                       gAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAAt                              <  2:324418/150‑1 (MQ=255)

                                                                                                                       gAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAAt                              >  1:3323641/1‑150 (MQ=255)

                                                                                                                            gCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTCGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCt                         >  1:4185340/1‑150 (MQ=255)

                                                                                                                                     ccAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTgca                                           <  1:2818527/123‑1 (MQ=255)

                                                                                                                                     ccAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAgc                >  2:2249270/1‑150 (MQ=255)

                                                                                                                                     ccAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAgc                >  1:1993707/1‑150 (MQ=255)

                                                                                                                                      cAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAgcg               >  2:3732891/1‑150 (MQ=255)

                                                                                                                                      cAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAgcg               >  2:4252825/1‑150 (MQ=255)

                                                                                                                                       aGGCCTGCTGGCTGCTGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGaa                               >  2:3235546/1‑133 (MQ=25)

                                                                                                                                       aGGCCTGCTGGCTGCTGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGaa                               <  1:3235546/133‑1 (MQ=25)

                                                                                                                                       aGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAgcgc              <  1:4173118/150‑1 (MQ=255)

                                                                                                                                       aGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAgcgc              <  2:1952286/150‑1 (MQ=255)

                                                                                                                                         gccagcagccagcaGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCgg            >  1:2028016/1‑150 (MQ=255)

                                                                                                                                           cagcagccagcagGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGt          >  1:1073658/1‑150 (MQ=255)

                                                                                                                                            agcagccagcagGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGtt         <  1:2688629/150‑1 (MQ=255)

                                                                                                                                               agccagcagGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGTTTcc      >  2:3092085/1‑150 (MQ=255)

                                                                                                                                               agccagcagGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGTTTcc      <  1:48785/150‑1 (MQ=255)

                                                                                                                                                   agcagGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGTTTCCTTCg  <  2:2414930/150‑1 (MQ=255)

ATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGTTTCCTTCG  >  NC_002947/4238371‑4238667

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