BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,238,516	C→G	5.3%	A75P (GCT→CCT) ‡	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,516	0	C	G	5.3%	569.9 / 15.9	170	A75P (GCT→CCT) ‡	catA‑I	catechol 1,2‑dioxygenase
Reads supporting (aligned to +/- strand): ref base C (71/90); new base G (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.84e-01

Only 100 of 186 total aligned reads displayed.

GTCCATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGTTTCCTT  >  NC_002947/4238367‑4238665

gtcCATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGgccagcagc                                                                                                                                                       <  1:3826793/150‑1 (MQ=255)

gtcCATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGgccagcagc                                                                                                                                                       >  2:4670524/1‑150 (MQ=255)

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

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

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

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

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

            ttCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGgcctgc                                                                                                                                           >  1:3208008/1‑150 (MQ=255)

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

                cccTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCTGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCg                                                                                                                                       <  2:1467408/150‑1 (MQ=255)

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

                 ccTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGt                                                                                                                                      >  2:3531213/1‑150 (MQ=255)

                  cTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCTGCTGGCTGCTGgcctgggtatcg                                                                                                                                     <  2:3778030/150‑8 (MQ=25)

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

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

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

                       ccAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGgc                                                                                                                                >  1:2790714/1‑150 (MQ=255)

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

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

                             gcgcCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGATTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCa                                                                                                                          <  2:3540000/150‑1 (MQ=255)

                               gcCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCCGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAgg                                                                                                                        >  2:5112300/1‑150 (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)

                                    ggCAACGTACAACGGGCATTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGtt                                                                                                                   <  2:3053503/150‑1 (MQ=255)

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

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

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

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

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

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

                                       aaCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCATCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                <  1:416036/150‑1 (MQ=255)

                                                    ccTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCg                                                                                                   <  1:2845598/150‑1 (MQ=255)

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

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

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

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

                                                       tCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGc                                                                                                <  2:2341089/150‑1 (MQ=255)

                                                        cGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGcc                                                                                               <  1:347705/150‑1 (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                                                                                              >  1:457771/1‑145 (MQ=255)

                                                              tGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACt                                                                                         >  1:350588/1‑150 (MQ=255)

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

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

                                                                  aaggatGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                <  2:2162883/118‑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:5419093/150‑1 (MQ=255)

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

                                                                       tGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACg                                                                                                      <  1:1953438/128‑1 (MQ=255)

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

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

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

                                                                          cgccgccAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTaa                                                                             <  1:5008556/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                  aGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAgg                                                                     >  2:3017928/1‑150 (MQ=255)

                                                                                    gCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTc                                                                   >  1:5282646/1‑150 (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:2685636/150‑1 (MQ=255)

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

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

                                                                                                 tCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAggcgg                                                      <  1:2897559/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                agcagGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTgcagc                                       >  1:3575882/1‑150 (MQ=255)

                                                                                                                    ggTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCAcgc                                   >  2:3104302/1‑150 (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)

                                                                                                                               tGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATc                        <  1:4420992/150‑1 (MQ=255)

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

                                                                                                                                    gATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCtt                   >  2:4567404/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                           aGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAgcgc            >  1:2421598/1‑150 (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)

GTCCATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGTTTCCTT  >  NC_002947/4238367‑4238665

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