BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,817,796	C→A	6.8%	A97S (GCA→TCA) ‡	PP_3371 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,817,796	0	C	A	6.8%	582.6 / 23.1	177	A97S (GCA→TCA) ‡	PP_3371	sensor histidine kinase
Reads supporting (aligned to +/- strand): ref base C (67/98); new base A (7/5); total (74/103)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 179 total aligned reads displayed.

CAGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTGCAGCT  >  NC_002947/3817647‑3817945

  ggTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCt                                                                                                                                                     <  2:3870636/150‑1 (MQ=255)

      aaaGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACTGTCGAGCAGTGCCTGCTc                                                                                                                                                 <  2:5083126/150‑1 (MQ=255)

       aaGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCg                                                                                                                                                <  1:660545/150‑1 (MQ=255)

         gCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGaa                                                                                                                                              <  2:1131065/150‑1 (MQ=255)

         gCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGaa                                                                                                                                              <  1:2314252/150‑1 (MQ=255)

         gCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGaa                                                                                                                                              <  2:4748896/150‑1 (MQ=255)

         gCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGaa                                                                                                                                              <  1:3297096/150‑1 (MQ=255)

         gCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGaa                                                                                                                                              <  2:2155373/150‑1 (MQ=255)

         gCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGGCACTGCTCGac                                                                                                                                              <  1:2140904/150‑2 (MQ=255)

          cTGCACAACCGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGGCACTGCTCGacc                                                                                                                                             <  2:5007616/150‑3 (MQ=255)

          cTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGGCACTGCTCGacc                                                                                                                                             <  2:761547/150‑3 (MQ=255)

           tGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAgc                                                                                                                                            >  1:3773912/1‑150 (MQ=255)

           tGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAgc                                                                                                                                            >  1:2892739/1‑150 (MQ=255)

            gCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAgcg                                                                                                                                           >  1:5091523/1‑150 (MQ=255)

             cacaACAGCGCGAACATCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAgcgc                                                                                                                                          >  1:3058324/1‑150 (MQ=255)

             cacaACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAgcgc                                                                                                                                          <  2:139139/150‑1 (MQ=255)

             cacaACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAgcgc                                                                                                                                          <  1:1215810/150‑1 (MQ=255)

              acaacaGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAgcgcg                                                                                                                                         >  2:570226/1‑150 (MQ=255)

                aacagcgcgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGcc                                                                                                                                       <  1:1027452/150‑1 (MQ=255)

                aacagcgcgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGcc                                                                                                                                       >  1:4102273/1‑150 (MQ=255)

                aacagcgcgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGcc                                                                                                                                       >  1:1135386/1‑150 (MQ=255)

                     cgcgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTg                                                                                                                                  <  2:546614/150‑1 (MQ=255)

                     cgcgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTACCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTg                                                                                                                                  <  2:3713786/150‑1 (MQ=255)

                      gcgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGc                                                                                                                                 >  2:2788298/1‑150 (MQ=255)

                           cagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGt                                                                                                                            >  1:3142641/1‑150 (MQ=255)

                            agcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGtt                                                                                                                           >  1:4202694/1‑150 (MQ=255)

                              cgcgcTCATCACCGCCCCAGTCACGAACGGGCCACGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTg                                                                                                                         <  1:3668001/150‑1 (MQ=255)

                              cgcgcTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTg                                                                                                                         <  2:242009/150‑1 (MQ=255)

                               gcgcTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTgg                                                                                                                        <  1:2485053/150‑1 (MQ=255)

                               gcgcTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTgg                                                                                                                        <  2:4708958/150‑1 (MQ=255)

                               gcgcTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTgg                                                                                                                        <  2:4401441/150‑1 (MQ=255)

                                  cTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgcca                                                                                                                     <  1:811142/150‑1 (MQ=255)

                                  cTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgcca                                                                                                                     <  1:811037/150‑1 (MQ=255)

                                  cTCATCACCGCCCCAGTCACGAACGGGCCAAGCTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgcca                                                                                                                     <  2:5094586/150‑1 (MQ=255)

                                   tcatcaCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgccag                                                                                                                    >  1:4642172/1‑150 (MQ=255)

                                    catcaCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgccagc                                                                                                                   >  1:4511081/1‑150 (MQ=255)

                                     atcaCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgccagcc                                                                                                                  >  1:692233/1‑150 (MQ=255)

                                     atcaCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgccagcc                                                                                                                  >  1:1224364/1‑150 (MQ=255)

                                      tcaCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgccagcca                                                                                                                 <  1:4685315/150‑1 (MQ=255)

                                        aCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCcagccagc                                                                                                               >  1:3479634/1‑150 (MQ=255)

                                        aCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCcagccagc                                                                                                               <  2:171723/150‑1 (MQ=255)

                                          cGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCaa                                                                                                             <  1:1540940/150‑1 (MQ=255)

                                             cccAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCg                                                                                                          >  1:1095764/1‑150 (MQ=255)

                                              ccAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCgg                                                                                                         <  1:2775265/150‑1 (MQ=255)

                                                aGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTc                                                                                                       <  1:1969547/150‑1 (MQ=255)

                                                  tCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCtt                                                                                                     <  2:196832/150‑1 (MQ=255)

                                                    aCGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGa                                                                                                   <  1:810815/150‑1 (MQ=255)

                                                         cGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGAtcgtc                                                                                              <  2:4211408/150‑1 (MQ=255)

                                                             ccAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTaa                                                                                          >  1:4260349/1‑150 (MQ=255)

                                                              cAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTaaa                                                                                         >  2:1617671/1‑150 (MQ=255)

                                                                  gTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGcc                                                                                     >  2:581269/1‑150 (MQ=255)

                                                                   ttGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCt                                                                                    >  1:1896316/1‑150 (MQ=255)

                                                                         gggCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGc                                                                              <  2:3200259/150‑1 (MQ=255)

                                                                         gggCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGc                                                                              <  2:2458407/150‑1 (MQ=255)

                                                                         gggCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGc                                                                              <  1:2106285/150‑1 (MQ=255)

                                                                         gggCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGc                                                                              <  2:2073500/150‑1 (MQ=255)

                                                                         gggCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGc                                                                              <  1:4011551/150‑1 (MQ=255)

                                                                           gcgcGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGcac                                                                            >  2:155901/1‑150 (MQ=255)

                                                                                 tGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTg                                                                      <  2:2197806/150‑1 (MQ=255)

                                                                                    cATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCgg                                                                   <  1:3613433/150‑1 (MQ=255)

                                                                                     aTACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGc                                                                  >  1:1833151/1‑150 (MQ=255)

                                                                                     aTACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGc                                                                  >  1:160787/1‑150 (MQ=255)

                                                                                      tACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCt                                                                 >  1:4122447/1‑150 (MQ=255)

                                                                                         cGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCtgtg                                                              >  1:1373404/1‑150 (MQ=255)

                                                                                          gACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTgg                                                             <  1:2471205/150‑1 (MQ=255)

                                                                                          gACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTgg                                                             <  1:3032350/150‑1 (MQ=255)

                                                                                           acacGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGc                                                            <  1:536990/150‑1 (MQ=255)

                                                                                           acacGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGc                                                            <  2:3908135/150‑1 (MQ=255)

                                                                                            cacGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCg                                                           >  1:2409884/1‑150 (MQ=255)

                                                                                             acGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGa                                                          <  2:850604/150‑1 (MQ=255)

                                                                                               gggTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGaca                                                        <  1:1827278/150‑1 (MQ=255)

                                                                                               gggTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGaca                                                        >  2:2551167/1‑150 (MQ=255)

                                                                                               gggTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGaca                                                        >  2:4406542/1‑150 (MQ=255)

                                                                                               gggTGGCCACCGAAGCGCGCAGCACCAGTGGGAGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGACTACTTGCACATCCTGCGGCTGTGGCGaca                                                        >  1:4539303/1‑150 (MQ=255)

                                                                                                 gTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACg                                                      >  2:4086602/1‑150 (MQ=255)

                                                                                                     ccACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCaa                                                  >  1:2015309/1‑150 (MQ=255)

                                                                                                     ccACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCaa                                                  <  1:1278451/150‑1 (MQ=255)

                                                                                                     ccACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCaa                                                  <  2:2775328/150‑1 (MQ=255)

                                                                                                     ccACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCaa                                                  <  1:3569704/150‑1 (MQ=255)

                                                                                                      cACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAg                                                 <  1:708558/150‑1 (MQ=255)

                                                                                                            aGCGCGCAGCACCAGTGGGCGGCCGTCTCGGTCGAGCATTGCCTGCTCGAAGCGCGCCACGTTCACCGTTTGGCCAGCCAGCAAGCGGTATTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACtt                                           <  1:3408035/150‑2 (MQ=255)

                                                                                                             gcgcgcAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTgg                                          >  1:3085506/1‑150 (MQ=255)

                                                                                                               gcgcAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTggcg                                        <  2:3462514/150‑1 (MQ=255)

                                                                                                               gcgcAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTggcg                                        >  2:2494280/1‑150 (MQ=255)

                                                                                                               gcgcAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTggcg                                        >  2:5037431/1‑150 (MQ=255)

                                                                                                                 gcagcaCCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTTGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTggcggc                                      >  2:3488705/1‑150 (MQ=255)

                                                                                                                 gcagcaCCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTggcggc                                      >  2:4044816/1‑150 (MQ=255)

                                                                                                                  cagcaCCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGcggcc                                     >  2:1561257/1‑150 (MQ=255)

                                                                                                                     caCCAGTGGGCGGCCATCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGcggccggc                                  <  2:2068601/150‑1 (MQ=255)

                                                                                                                           tGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTc                            <  1:3170592/150‑1 (MQ=255)

                                                                                                                                 gcCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACGTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATg                      <  2:105742/150‑1 (MQ=255)

                                                                                                                                       aCGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTc                <  2:3350505/150‑1 (MQ=255)

                                                                                                                                       aCGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTc                <  1:1339222/150‑1 (MQ=255)

                                                                                                                                         ggTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGa              >  2:5021509/1‑150 (MQ=255)

                                                                                                                                         cttCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGa              >  1:3696282/3‑150 (MQ=255)

                                                                                                                                         cttCGAGCAGGCACTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGa              >  2:26391/3‑150 (MQ=255)

                                                                                                                                          ttCGAGCAGGCACTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAg             >  1:4439051/2‑150 (MQ=255)

                                                                                                                                          ttCGAGCAGGCACTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCGGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAg             >  2:904590/2‑150 (MQ=255)

                                                                                                                                             gAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCAt          <  2:4230809/150‑1 (MQ=255)

                                                                                                                                                cAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGGCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGTCCGGCCCGCTCGCCATGCCGTTCGAGCATTTg       <  2:3176620/150‑1 (MQ=255)

CAGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTGCAGCT  >  NC_002947/3817647‑3817945

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