BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,886,958	T→G	6.1%	A53A (GCA→GCC) ‡	PP_5159 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,886,958	0	T	G	6.1%	583.0 / 20.7	181	A53A (GCA→GCC) ‡	PP_5159	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (87/83); new base G (5/6); total (92/89)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

Only 100 of 187 total aligned reads displayed.

AGGAAAAGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAGGGGA  >  NC_002947/5886809‑5887107

aGGAAAAGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCt                                                                                                                                                       >  1:3977660/1‑150 (MQ=255)

    aaaGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCgg                                                                                                                                                   <  1:1049399/150‑1 (MQ=255)

     aaGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCggg                                                                                                                                                  <  1:3829761/150‑1 (MQ=255)

     aaGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCggg                                                                                                                                                  <  2:4107548/150‑1 (MQ=255)

       ggCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGt                                                                                                                                                <  1:2427381/150‑1 (MQ=255)

        gCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCACCCCCTGCGGGGtt                                                                                                                                               >  2:3003214/1‑150 (MQ=255)

        gCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGtt                                                                                                                                               <  2:2037359/150‑1 (MQ=255)

          cTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGtatg                                                                                                                                             >  1:4353792/1‑147 (MQ=255)

          cTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGa                                                                                                                                             <  1:4919494/150‑1 (MQ=255)

          cTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGa                                                                                                                                             >  2:246200/1‑150 (MQ=255)

           ttGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAg                                                                                                                                            >  2:2943887/1‑150 (MQ=255)

             ggTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCGCAGGGGTTgggcg                                                                                                                                          <  2:3417806/150‑5 (MQ=255)

               ttGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCGCAGGGGTTgggcggg                                                                                                                                        <  2:3836581/150‑7 (MQ=255)

                tGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGc                                                                                                                                       <  1:1323833/150‑1 (MQ=255)

                    tGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTc                                                                                                                                   <  1:1281075/150‑1 (MQ=255)

                    tGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTc                                                                                                                                   <  2:2044545/150‑1 (MQ=255)

                    tGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTc                                                                                                                                   <  2:245152/150‑1 (MQ=255)

                     gCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCg                                                                                                                                  <  1:388853/150‑1 (MQ=255)

                         ggCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCAt                                                                                                                              >  2:5381067/1‑150 (MQ=255)

                              aGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACgg                                                                                                                         >  2:2424509/1‑150 (MQ=255)

                                  gATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACAGCGcc                                                                                                                     >  2:2849413/1‑150 (MQ=255)

                                   aTCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGc                                                                                                                                       <  2:4249598/131‑1 (MQ=255)

                                   aTCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGc                                                                                                                                       >  1:4249598/1‑131 (MQ=255)

                                    tCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCgg                                                                                                                   >  1:2798146/1‑150 (MQ=255)

                                            cgccAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGc                                                                                                           >  2:4875556/1‑150 (MQ=255)

                                             gccAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCt                                                                                                          >  1:5138502/1‑150 (MQ=255)

                                             gccAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCt                                                                                                          <  1:5384401/150‑1 (MQ=255)

                                                 ggTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCaa                                                                                                      <  2:37912/150‑1 (MQ=255)

                                                 ggTACCGGTCGCAGGCCTCGGCCTGCCCCGGGACGATCCACTGGATCTGGCCACCGCGAATGGTGTCGCGCACCCCCGGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGCACGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCaa                                                                                                      <  1:294143/150‑1 (MQ=255)

                                                  gTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAg                                                                                                     >  1:2024861/1‑150 (MQ=255)

                                                     cTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgc                                                                                                  <  1:1023992/150‑1 (MQ=255)

                                                      tGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgcg                                                                                                 >  1:4413030/1‑150 (MQ=255)

                                                       ggTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGGTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCCAGTgcgc                                                                                                >  1:2979212/1‑150 (MQ=255)

                                                        gTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACCGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgcgcg                                                                                               >  2:12271/1‑150 (MQ=255)

                                                         tCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgcgcgc                                                                                              <  1:4561982/150‑1 (MQ=255)

                                                          cGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCTCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCa                                                                                             >  1:3174405/1‑150 (MQ=255)

                                                           gCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCAc                                                                                            <  2:4950107/150‑1 (MQ=255)

                                                            cAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCAcc                                                                                           >  2:2385009/1‑150 (MQ=255)

                                                             aGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCa                                                                                          >  2:2502447/1‑150 (MQ=255)

                                                             aGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCa                                                                                          >  1:3633684/1‑150 (MQ=255)

                                                             aGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCa                                                                                          >  2:5229690/1‑150 (MQ=255)

                                                               gCCTCGGCCTGCCCCGGGTCGATCAACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAgg                                                                                        >  2:526293/1‑150 (MQ=255)

                                                                 cTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTc                                                                                      <  1:2453651/150‑1 (MQ=255)

                                                                  tCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCg                                                                                     <  1:2360598/150‑1 (MQ=255)

                                                                  tCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCg                                                                                     >  2:4350152/1‑150 (MQ=255)

                                                                  tCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCg                                                                                     <  1:1221700/150‑1 (MQ=255)

                                                                     gCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggca                                                                                  <  1:4662586/150‑1 (MQ=255)

                                                                          ccccGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggc                                                                                   >  2:625010/1‑144 (MQ=255)

                                                                          ccccGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggc                                                                                   <  1:625010/144‑1 (MQ=255)

                                                                             cGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGaaa                                                                          >  2:212848/1‑150 (MQ=255)

                                                                             cGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGaaa                                                                          <  1:1575390/150‑1 (MQ=255)

                                                                                 tCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTgc                                                                      <  1:2574273/150‑1 (MQ=255)

                                                                                   gATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTgcgc                                                                    <  1:3982235/150‑1 (MQ=255)

                                                                                    aTCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCAGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGcc                                                                   >  1:3724337/1‑150 (MQ=255)

                                                                                     tCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGTTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCct                                                                  >  2:3252099/1‑150 (MQ=255)

                                                                                       cACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggcagg                                                                                <  1:997804/134‑1 (MQ=255)

                                                                                        aCTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCctgct                                                               >  2:988723/1‑150 (MQ=255)

                                                                                        aCTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCctgct                                                               >  1:4914802/1‑150 (MQ=255)

                                                                                           ggATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTggg                                                            >  1:3270329/1‑150 (MQ=255)

                                                                                               cTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCa                                                        <  1:1911239/150‑1 (MQ=255)

                                                                                                     aCCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATg                                                  >  1:4169182/1‑150 (MQ=255)

                                                                                                     aCCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATg                                                  >  2:4611020/1‑150 (MQ=255)

                                                                                                     aCCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATg                                                  >  2:1715620/1‑150 (MQ=255)

                                                                                                      ccGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                                 >  1:3381197/1‑150 (MQ=255)

                                                                                                      ccGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                                 >  1:1376780/1‑150 (MQ=255)

                                                                                                       cgcgAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATggg                                                >  1:461474/1‑150 (MQ=255)

                                                                                                           aaTGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGc                                            <  1:4875556/150‑1 (MQ=255)

                                                                                                             tGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCa                                          >  2:4516000/1‑150 (MQ=255)

                                                                                                              ggTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAg                                         <  2:461474/150‑1 (MQ=255)

                                                                                                                    gcgcACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGAtcgtc                                   >  1:3527281/1‑150 (MQ=255)

                                                                                                                    gcgcACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGAtcgtc                                   <  2:5257655/150‑1 (MQ=255)

                                                                                                                             ctGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGc                          >  1:1677941/1‑150 (MQ=255)

                                                                                                                              tGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGcc                         <  2:3616849/150‑1 (MQ=255)

                                                                                                                               gggTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCg                        >  1:4035619/1‑150 (MQ=255)

                                                                                                                                  tAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCaa                     <  1:5223526/150‑1 (MQ=255)

                                                                                                                                  tAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCaa                     <  1:756987/150‑1 (MQ=255)

                                                                                                                                  tAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCaa                     <  2:4035619/150‑1 (MQ=255)

                                                                                                                                   aGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAg                    >  2:1594948/1‑150 (MQ=255)

                                                                                                                                      tCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCAt                 >  2:3672072/1‑150 (MQ=255)

                                                                                                                                      tCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCAt                 >  2:887747/1‑150 (MQ=255)

                                                                                                                                      accccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCAt                 >  2:2589666/2‑150 (MQ=255)

                                                                                                                                       ccccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATc                >  1:1429852/1‑150 (MQ=255)

                                                                                                                                       ccccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATc                <  2:5172228/150‑1 (MQ=255)

                                                                                                                                         gaactcAACCCCGCAGGGGTTGGGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg              >  2:2756904/6‑150 (MQ=255)

                                                                                                                                         gaactcAACCCCGCAGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg              >  1:4735192/6‑150 (MQ=255)

                                                                                                                                         gaactcAACCCCGCAGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg              >  1:597252/6‑150 (MQ=255)

                                                                                                                                         ccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg              <  2:3776443/150‑1 (MQ=255)

                                                                                                                                         ccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg              <  1:3297909/150‑1 (MQ=255)

                                                                                                                                         ccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg              <  1:210681/150‑1 (MQ=255)

                                                                                                                                           gCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGgt            <  2:3816229/150‑1 (MQ=255)

                                                                                                                                             ccAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGgtgt          >  2:1182969/1‑150 (MQ=255)

                                                                                                                                             ccAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGgtgt          >  2:969617/1‑150 (MQ=255)

                                                                                                                                              cAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGtt         <  1:4051268/150‑1 (MQ=255)

                                                                                                                                               aaCCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTc        <  1:3482060/150‑1 (MQ=255)

                                                                                                                                               aaCCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTc        <  1:747214/150‑1 (MQ=255)

                                                                                                                                                 ccccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAg      <  1:4882375/150‑1 (MQ=255)

                                                                                                                                                  cccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAgg     >  1:1252371/1‑150 (MQ=255)

                                                                                                                                                   ccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAggg    <  2:1376780/150‑1 (MQ=255)

                                                                                                                                                     tGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAGGGGa  >  1:1091888/1‑150 (MQ=255)

                                                                                                                                                     tGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAGGGGa  <  1:2589666/150‑1 (MQ=255)

AGGAAAAGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAGGGGA  >  NC_002947/5886809‑5887107

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