BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,238,521	A→T	5.3%	L73Q (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,521	0	A	T	5.3%	579.6 / 14.8	170	L73Q (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 = 8.95e-01

Only 100 of 188 total aligned reads displayed.

TGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGTTTCCTTCGG  >  NC_002947/4238372‑4238668

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

    cgcTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGgcc                                                                                                                                                 >  1:324418/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)

                   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)

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

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

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

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

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

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

                                  aaCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                   <  1:4746980/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                                                                                                                  >  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)

                                                    gATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTTGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCa                                                                                                 <  1:624782/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                                                                                                                   >  1:2162876/6‑123 (MQ=255)

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

                                                                         gccAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGATGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGt                                                                                                     <  2:35365/125‑1 (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:1591565/1‑150 (MQ=255)

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

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

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

                                                                                ccGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCt                                                                     >  2:2527074/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                               ttGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAggcgggcg                                                      >  1:3689789/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)

                                                                                                               ggTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCAcgc                                      >  2:3104302/1‑150 (MQ=255)

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

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

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

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

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

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

                                                                                                                            cTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTg                         <  2:2946165/150‑1 (MQ=255)

                                                                                                                            cTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTg                         <  2:215629/150‑1 (MQ=255)

                                                                                                                                    ccAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                                   <  2:3228576/52‑1 (MQ=255)

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

                                                                                                                                    ccAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTgca                                            >  2:2818527/1‑123 (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:3235622/1‑133 (MQ=25)

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

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

                                                                                                                                      aGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAgcgc               <  1:4173118/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       <  1:48785/150‑1 (MQ=255)

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

TGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGGGTTTCCTTCGG  >  NC_002947/4238372‑4238668

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