BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,886,960	C→A	5.5%	A53S (GCA→TCA) ‡	PP_5159 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,886,960	0	C	A	5.5%	618.0 / 13.5	185	A53S (GCA→TCA) ‡	PP_5159	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (90/84); new base A (5/6); total (95/90)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 186 total aligned reads displayed.

GAAAAGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAGGGGA  >  NC_002947/5886811‑5887107

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

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

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

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

     ggCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGt                                                                                                                                                >  2:2280009/1‑150 (MQ=255)

     ggCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGt                                                                                                                                                >  1:5172228/1‑150 (MQ=255)

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

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

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

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

        cTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGa                                                                                                                                             <  1:3691889/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)

            gTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTc                                                                                                                                         <  1:3227385/150‑1 (MQ=255)

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

              tGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCTCAACCCCGCAGGGGTTgggcgggg                                                                                                                                       <  2:3718118/150‑8 (MQ=255)

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

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

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

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

                    ccAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCgg                                                                                                                                 >  1:4001028/1‑150 (MQ=255)

                    ccAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCgg                                                                                                                                 >  1:4442648/1‑150 (MQ=255)

                      aGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCa                                                                                                                               >  1:3005562/1‑150 (MQ=255)

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

                            aGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACgg                                                                                                                         >  2:2424509/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)

                                    cATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCa                                                                                                                 <  1:2424509/150‑1 (MQ=255)

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

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

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

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

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

                                                 tACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGt                                                                                                    <  2:1114281/150‑1 (MQ=255)

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

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

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

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

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

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

                                                             gCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCggg                                                                                                                                                  <  1:4050380/92‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                              gTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTg                                                                       >  1:3776443/1‑150 (MQ=255)

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

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

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

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

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

                                                                                         ggATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTggg                                                            <  1:5398862/150‑1 (MQ=255)

                                                                                           aTCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGAc                                                          <  2:277491/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                    ccGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                                 >  1:3381197/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                                   <  2:5257655/150‑1 (MQ=255)

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

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

                                                                                                                    gcACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTcga                                 >  2:1268412/1‑150 (MQ=255)

                                                                                                                         tcctGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGacgacg                            >  2:4446278/1‑150 (MQ=255)

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

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

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

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

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

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

                                                                                                                                    tCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGGCGACGACGGCCGCAAGCAt                 >  2:1211049/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)

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

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

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

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

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

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

                                                                                                                                         actcAACCCCGCAGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCaa                                                                                                      >  1:4594174/4‑60 (MQ=37)

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

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

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

                                                                                                                                              aCCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCa       <  2:2875528/150‑1 (MQ=255)

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

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

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

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

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

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

GAAAAGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAGGGGA  >  NC_002947/5886811‑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: `‑`