BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,817,794	T→G	6.4%	A97A (GCA→GCC) ‡	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,794	0	T	G	6.4%	606.6 / 21.2	181	A97A (GCA→GCC) ‡	PP_3371	sensor histidine kinase
Reads supporting (aligned to +/- strand): ref base T (69/100); new base G (7/5); total (76/105)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

GTCAGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTGC  >  NC_002947/3817645‑3817941

gTCAGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGt                                                                                                                                                     <  2:1298208/150‑1 (MQ=255)

 tCAGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTg                                                                                                                                                    >  2:3779861/1‑150 (MQ=255)

       aaaaGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCt                                                                                                                                              >  2:2474667/1‑150 (MQ=255)

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

        aaaGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTggtgctgc                                                                                                                                             <  2:3747515/150‑8 (MQ=255)

           gCTGCACAACCGCGCGAACCGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCCAGGGTGCGTGGGCGCGGTGCCCTACGGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCGTGCTCGaa                                                                                                                                          <  2:3664568/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                             cagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGAGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGt                                                                                                                        >  1:3142650/1‑150 (MQ=255)

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

                               gcgcgcTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGttt                                                                                                                      <  1:686882/150‑1 (MQ=255)

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

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

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

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

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

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

                                  cgcTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgc                                                                                                                   <  2:1112942/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)

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

                                          aCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCcggccagc                                                                                                           >  1:4680700/1‑150 (MQ=255)

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

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

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

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

                                               cccAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCGCCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCg                                                                                                      <  2:2844308/150‑1 (MQ=255)

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

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

                                                   gTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCt                                                                                                  <  1:3615674/150‑1 (MQ=255)

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

                                                      aCGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGa                                                                                               >  2:3198890/1‑150 (MQ=255)

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

                                                      aCGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGACGAGCAGTGCCTGCTCGAAGCGCGCCACCTGCACCGGTTGGCCAGCCAGCAAGCCGTCTAGa                                                                                               >  2:4239341/1‑150 (MQ=255)

                                                           cGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCTGTGCCTGCTCGATGCGCGCCTCGTGCACCTTTTGGCCTGGCAGCAAGCGTTCTTGAtcgtc                                                                                          <  2:4936462/150‑1 (MQ=255)

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

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

                                                                 aaGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAg                                                                                    <  1:4909079/150‑1 (MQ=255)

                                                                 aaGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAg                                                                                    >  1:2141060/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                   gTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACg                                                  <  1:3016820/150‑1 (MQ=255)

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

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

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

                                                                                                         aCCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCTAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGc                                            >  1:2392060/1‑150 (MQ=255)

                                                                                                           cGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCaa                                          >  1:252099/1‑150 (MQ=255)

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

                                                                                                              aGCGCGCAGCACCAGTGGGCGGCCGTCTCGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTg                                       <  2:4729312/150‑1 (MQ=255)

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

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

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

                                                                                                                   gcagcaCCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAGAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTggcggc                                  >  2:4045047/1‑150 (MQ=255)

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

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

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

                                                                                                                         ccAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGccc                            >  1:2466172/1‑150 (MQ=255)

                                                                                                                         ccAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGccc                            >  2:2506063/1‑150 (MQ=255)

                                                                                                                         ccAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGccc                            >  2:4219440/1‑150 (MQ=255)

                                                                                                                                           ggTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGa          <  1:772670/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                                                                  cAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTg   >  2:3204347/1‑150 (MQ=255)

GTCAGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTGC  >  NC_002947/3817645‑3817941

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