BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,876,600	G→C	5.7%	H258Q (CAC→CAG)	PP_3420 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,876,600	0	G	C	5.7%	826.8 / 29.7	244	H258Q (CAC→CAG)	PP_3420	sensor histidine kinase
Reads supporting (aligned to +/- strand): ref base G (120/110); new base C (6/8); total (126/118)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 256 total aligned reads displayed.

TTCTGCGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTCCACCATGGC  >  NC_002947/3876451‑3876749

ttctGCGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAgg                                                                                                                                                       >  1:3051768/1‑150 (MQ=255)

ttctGCGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAgg                                                                                                                                                       >  1:3909981/1‑150 (MQ=255)

  ctGCGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTg                                                                                                                                                     >  2:1921859/1‑150 (MQ=255)

   tGCGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTgc                                                                                                                                                    <  2:3601292/150‑1 (MQ=255)

   tGCGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGGAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTgc                                                                                                                                                    >  2:5135871/1‑150 (MQ=255)

    gcgcAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTgca                                                                                                                                                   >  1:1018750/1‑150 (MQ=255)

    gcgcAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTgca                                                                                                                                                   >  2:5518739/1‑150 (MQ=255)

    gcgcAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTgca                                                                                                                                                   <  1:5323619/150‑1 (MQ=255)

        agcagcGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCat                                                                                                                                               >  1:2393761/1‑150 (MQ=255)

        agcagcGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCACCTgcagcag                                                                                                                                               <  2:1281676/150‑2 (MQ=255)

         gcagcGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCACCTgcagcagg                                                                                                                                              <  2:4082651/150‑3 (MQ=255)

             cGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATAtgtg                                                                                                                                          <  1:5097812/150‑1 (MQ=255)

              gTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTgc                                                                                                                                         >  1:3757749/1‑150 (MQ=255)

               ttGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTgcc                                                                                                                                        <  2:2784401/150‑1 (MQ=255)

                     gTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCgg                                                                                                                                  >  1:825350/1‑150 (MQ=255)

                       ttCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGgcg                                                                                                                                <  1:763131/150‑1 (MQ=255)

                        tCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGgcgc                                                                                                                               >  1:3375817/1‑150 (MQ=255)

                          cTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCAc                                                                                                                             >  1:2059026/1‑150 (MQ=255)

                           tGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCata                                                                                                                            <  1:2808255/150‑3 (MQ=255)

                           tGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACg                                                                                                                            >  1:2792992/1‑150 (MQ=255)

                           tGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACg                                                                                                                            >  2:4135309/1‑150 (MQ=255)

                            gggCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGa                                                                                                                           >  1:14680/1‑150 (MQ=255)

                              gCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATc                                                                                                                         <  2:891254/150‑1 (MQ=255)

                              gCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATc                                                                                                                         <  2:3721442/150‑1 (MQ=255)

                                   ttCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGaa                                                                                                                    >  1:5199650/1‑150 (MQ=255)

                                   ttCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGaa                                                                                                                    <  2:4382806/150‑1 (MQ=255)

                                    tCGGCGCTGGCGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAAc                                                                                                                   >  1:1777770/1‑150 (MQ=255)

                                     cGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACa                                                                                                                  <  1:188912/150‑1 (MQ=255)

                                       gcgcTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAgg                                                                                                                <  2:5627553/150‑1 (MQ=255)

                                        cgcTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGc                                                                                                               >  2:4935216/1‑150 (MQ=255)

                                            ggAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGaa                                                                                                           >  1:2254299/1‑150 (MQ=255)

                                             gAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAg                                                                                                          >  2:4803355/1‑150 (MQ=255)

                                             gAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAg                                                                                                          <  2:2263687/150‑1 (MQ=255)

                                               gCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTa                                                                                                        <  1:5472274/150‑1 (MQ=255)

                                               gCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTa                                                                                                        <  2:2613523/150‑1 (MQ=255)

                                                 gTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGa                                                                                                      <  1:4651988/150‑1 (MQ=255)

                                                  tATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGCGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGAt                                                                                                     <  1:5547918/150‑1 (MQ=255)

                                                       tGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACg                                                                                                >  2:5348372/1‑150 (MQ=255)

                                                        ggCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGc                                                                                               <  2:2732149/150‑1 (MQ=255)

                                                        ggCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGc                                                                                               >  1:3111572/1‑150 (MQ=255)

                                                         gCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGcc                                                                                              >  2:3462519/1‑150 (MQ=255)

                                                          ccAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGccc                                                                                             <  1:1687365/150‑1 (MQ=255)

                                                            aGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCccca                                                                                           >  1:2065482/1‑150 (MQ=255)

                                                            aGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCccca                                                                                           >  2:3345109/1‑150 (MQ=255)

                                                                 gTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCccacccac                                                                                      <  1:1500925/150‑1 (MQ=255)

                                                                 gTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCccacccac                                                                                      >  1:2180073/1‑150 (MQ=255)

                                                                      gCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAgtg                                                                                 >  2:4017649/1‑150 (MQ=255)

                                                                           cTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAg                                                                            <  1:1373760/150‑1 (MQ=255)

                                                                              tgAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCa                                                                                                                              <  1:810124/97‑1 (MQ=255)

                                                                              tgAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCa                                                                                                                              >  2:811941/1‑97 (MQ=255)

                                                                               gAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGccc                                                                        >  1:3380423/1‑150 (MQ=255)

                                                                                 gCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAg                                                                      <  2:2929685/150‑1 (MQ=255)

                                                                                  ccACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGa                                                                     >  2:2838971/1‑150 (MQ=255)

                                                                                  ccACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGa                                                                     <  1:3127832/150‑1 (MQ=255)

                                                                                   cACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGAt                                                                    >  1:5230850/1‑150 (MQ=255)

                                                                                       gCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGc                                                                <  2:3351175/150‑1 (MQ=255)

                                                                                           acgacATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACAt                                                            >  1:4944134/1‑150 (MQ=255)

                                                                                            cgacATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATa                                                           <  2:271723/150‑1 (MQ=255)

                                                                                                aTTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCcg                                                       >  2:5607730/1‑150 (MQ=255)

                                                                                                 tttCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCcgc                                                      >  1:1287902/1‑150 (MQ=255)

                                                                                                      cccAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCCGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTgg                                                 >  1:5204579/1‑150 (MQ=255)

                                                                                                      cccAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTgg                                                 >  1:2516725/1‑150 (MQ=255)

                                                                                                      cccAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTgg                                                 >  2:3457920/1‑150 (MQ=255)

                                                                                                       ccAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGGAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTggg                                                >  2:778688/1‑149 (MQ=255)

                                                                                                        cAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGAt                                               <  2:2997050/150‑1 (MQ=255)

                                                                                                        cAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGAt                                               <  1:1987698/150‑1 (MQ=255)

                                                                                                        aaGCGAAACCTACGTTGCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCAGCGCTGGAt                                               >  1:590830/2‑150 (MQ=255)

                                                                                                         aGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATg                                              >  2:3800114/1‑150 (MQ=255)

                                                                                                          gcgcAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATgc                                             <  2:2212730/150‑1 (MQ=255)

                                                                                                                cGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATgcgcagaa                                       <  2:3051768/150‑5 (MQ=255)

                                                                                                                   aGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGa                                    <  1:2892118/150‑1 (MQ=255)

                                                                                                                     gTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCccca                                                                                           <  2:1100130/93‑1 (MQ=255)

                                                                                                                     gTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCccca                                                                                           >  1:1100130/1‑93 (MQ=255)

                                                                                                                      tttCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGt                                 >  2:1597754/1‑150 (MQ=255)

                                                                                                                          gCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTg                             <  1:4290410/150‑1 (MQ=255)

                                                                                                                           cTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGc                            <  1:4018381/150‑1 (MQ=255)

                                                                                                                            ttGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCg                           >  1:5161264/1‑150 (MQ=255)

                                                                                                                            ttGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCg                           <  1:503593/150‑1 (MQ=255)

                                                                                                                             tGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCgg                          <  2:4253899/150‑1 (MQ=255)

                                                                                                                                caccaGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTc                       >  1:5502742/1‑150 (MQ=255)

                                                                                                                                     gCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGt                  <  1:3310727/150‑1 (MQ=255)

                                                                                                                                     gCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGt                  <  1:3506285/150‑1 (MQ=255)

                                                                                                                                     gCCTGTCCTGCTGCAGGTGCAGCATATGGGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACAGAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGt                  <  2:3142387/150‑1 (MQ=255)

                                                                                                                                      ccTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTa                 <  2:2463624/150‑1 (MQ=255)

                                                                                                                                      ccTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTa                 <  1:3940121/150‑1 (MQ=255)

                                                                                                                                        tGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGa                                    >  1:5259422/1‑129 (MQ=255)

                                                                                                                                        tGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGa                                    <  2:5259422/129‑1 (MQ=255)

                                                                                                                                        tGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGa               >  1:4803737/1‑150 (MQ=255)

                                                                                                                                        tGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGACGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGa               >  1:3621326/1‑150 (MQ=255)

                                                                                                                                        acatatgctgcACCTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGa               >  2:3889182/6‑150 (MQ=255)

                                                                                                                                          tCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGAtt             >  1:2076459/1‑150 (MQ=255)

                                                                                                                                          tCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGAtt             >  2:1681756/1‑150 (MQ=255)

                                                                                                                                          atatgctgcACCTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCtt                    >  2:5068532/4‑143 (MQ=255)

                                                                                                                                            ctgctgCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTc           <  2:2973033/150‑1 (MQ=255)

                                                                                                                                             tgctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTcc          <  2:2065482/150‑1 (MQ=255)

                                                                                                                                             tgctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTcc          >  1:650809/1‑150 (MQ=255)

                                                                                                                                              gctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTcca         >  1:981416/1‑150 (MQ=255)

                                                                                                                                                   aGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTCCACCATg    >  1:2026703/1‑150 (MQ=255)

                                                                                                                                                     gTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTCCACCATGGc  <  2:2059026/150‑1 (MQ=255)

                                                                                                                                                     gTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTCCACCATGGc  <  2:3375817/150‑1 (MQ=255)

TTCTGCGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTCCACCATGGC  >  NC_002947/3876451‑3876749

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