BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,889,985	T→G	5.8%	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,889,985	0	T	G	5.8%	1395.7 / 50.8	413	A53A (GCA→GCC) ‡	PP_5159	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (188/201); new base G (16/8); total (204/209)
Fisher's exact test for biased strand distribution p-value = 9.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 415 total aligned reads displayed.

GAAAAGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAGGGGA  >  NC_002947/5889838‑5890134

  aaaGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCgg                                                                                                                                                   >  1:2618800/1‑150 (MQ=255)

    aGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCgggg                                                                                                                                                 <  2:4850365/150‑1 (MQ=255)

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

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

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

      gCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGtt                                                                                                                                               >  1:7291523/1‑150 (MQ=255)

       ccTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTg                                                                                                                                              >  2:3497297/1‑150 (MQ=255)

       ccTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTg                                                                                                                                              <  1:2533946/150‑1 (MQ=255)

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

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

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

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

         ttGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCGCAGGGGTTggg                                                                                                                                            <  2:2096733/150‑3 (MQ=255)

            gTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAgcgg                                                                                                                                         <  1:8520503/150‑4 (MQ=255)

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

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

            gTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTc                                                                                                                                         >  2:3712329/1‑150 (MQ=255)

            gTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCGCAGGGGTTgggcgg                                                                                                                                         <  2:8242664/150‑6 (MQ=255)

             ttGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTTGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCgggg                                                                                                                                                 <  1:2431131/141‑1 (MQ=255)

             ttGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTTGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCgggg                                                                                                                                                 >  2:2431131/1‑141 (MQ=255)

             ttGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCg                                                                                                                                        >  1:5022538/1‑150 (MQ=255)

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

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

                   gCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCg                                                                                                                                  >  2:188347/1‑150 (MQ=255)

                     cAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGc                                                                                                                                >  1:8774054/1‑150 (MQ=255)

                        gcgcAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATc                                                                                                                             <  1:2024617/150‑1 (MQ=255)

                          gcAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCAc                                                                                                                           <  1:6039995/150‑1 (MQ=255)

                                     aTGGCCGCCAGGAACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCAc                                                                                                                >  2:784981/1‑150 (MQ=255)

                                         ccgccAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCgg                                                                                                            <  2:3844453/150‑1 (MQ=255)

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

                                                 tACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGt                                                                                                    >  2:1316385/1‑150 (MQ=255)

                                                  aCTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTg                                                                                                   >  2:5155373/1‑150 (MQ=255)

                                                  aCTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTg                                                                                                   >  1:3266941/1‑150 (MQ=255)

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

                                                      gTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgcgcg                                                                                               >  1:5329682/1‑150 (MQ=255)

                                                      gTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgcgcg                                                                                               >  2:5878944/1‑150 (MQ=255)

                                                       tCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgcgcgc                                                                                              >  2:6916759/1‑150 (MQ=255)

                                                             gCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAgg                                                                                        >  1:3603339/1‑150 (MQ=255)

                                                             gCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAgg                                                                                        >  2:3008178/1‑150 (MQ=255)

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

                                                                  ggCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggc                                                                                   <  2:6743292/150‑1 (MQ=255)

                                                                   gCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggca                                                                                  >  1:6190851/1‑150 (MQ=255)

                                                                    ccTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggcag                                                                                 <  2:1120964/150‑1 (MQ=255)

                                                                    ccTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCgg                                                                                    <  1:4789242/147‑1 (MQ=255)

                                                                      tGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggcaggc                                                                               <  2:1549654/150‑1 (MQ=255)

                                                                       gCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggcaggca                                                                              >  2:2841180/1‑150 (MQ=255)

                                                                         cccGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGa                                                                            >  1:7395574/1‑150 (MQ=255)

                                                                         cccGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGa                                                                            <  1:4597769/150‑1 (MQ=255)

                                                                              gTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgcg                                                                                                 >  1:4174715/1‑124 (MQ=255)

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

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

                                                                                        tGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTgg                                                             >  1:2728817/1‑150 (MQ=255)

                                                                                          gATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGTGCCTGCTGGGa                                                           <  1:6717070/150‑1 (MQ=255)

                                                                                          gATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGa                                                           >  2:4036989/1‑150 (MQ=255)

                                                                                          gATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGa                                                           <  2:7483522/150‑1 (MQ=255)

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

                                                                                             cTGGTCACCGCGACTGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCa                                                        <  2:4681718/150‑1 (MQ=255)

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

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

                                                                                               ggTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGACAGTGCGCCTGCTGGGCCCAGc                                                      >  2:5076791/1‑150 (MQ=255)

                                                                                                 tCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCa                                                    >  1:1680575/1‑150 (MQ=255)

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

                                                                                                     cgcgAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCACGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATggg                                                >  1:3339292/1‑150 (MQ=255)

                                                                                                       cgAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTg                                              <  1:2277358/150‑1 (MQ=255)

                                                                                                        gAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTgg                                             >  1:7157053/1‑150 (MQ=255)

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

                                                                                                          aTGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGcc                                           <  2:5825946/150‑1 (MQ=255)

                                                                                                              tttCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGAt                                       >  1:5104077/1‑150 (MQ=255)

                                                                                                               ttCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGAtc                                      <  1:8103524/150‑1 (MQ=255)

                                                                                                                 cgcgCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGAtcgt                                    <  2:559313/150‑1 (MQ=255)

                                                                                                                  gcgcACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTgc                                                                      <  1:8717258/115‑1 (MQ=255)

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

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

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

                                                                                                                    gcACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTcga                                 <  2:8191353/150‑1 (MQ=255)

                                                                                                                     cACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTcgac                                <  2:6402710/150‑1 (MQ=255)

                                                                                                                       cctcctGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTcgacga                              <  1:1486990/150‑1 (MQ=255)

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

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

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

                                                                                                                               gTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCa                      <  1:3497297/150‑1 (MQ=255)

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

                                                                                                                                      cccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCg               >  2:5000936/1‑150 (MQ=255)

                                                                                                                                      cccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCg               >  1:7673325/1‑150 (MQ=255)

                                                                                                                                      cccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCg               >  1:4897091/1‑150 (MQ=255)

                                                                                                                                       gaactcAACCCCGCAGGGGTTGGGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg              >  1:3831397/6‑150 (MQ=255)

                                                                                                                                       gaactcAACCCCGCAGGGGTTGGGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg              >  1:5477399/6‑150 (MQ=255)

                                                                                                                                       gaactcAACCCCGCAGGGGTTGGGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg              >  1:5260926/6‑150 (MQ=255)

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

                                                                                                                                       ccGCCCACCCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGACTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCgg              >  2:1593772/1‑150 (MQ=255)

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

                                                                                                                                          ctcAACCCCGCAGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGACAGTGCGCCTGCTGGGACCAGCCAGGGGTGGCCCGATCGTCGACGACGGCCGCACGCATCGGgtg           >  2:7607569/3‑150 (MQ=255)

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

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

                                                                                                                                                cccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggca                                                                                  >  2:675164/1‑73 (MQ=255)

                                                                                                                                                cccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAgg     <  2:8726273/150‑1 (MQ=255)

                                                                                                                                                cccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAgg     <  2:5329682/150‑1 (MQ=255)

                                                                                                                                                cccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAgg     <  1:15731/150‑1 (MQ=255)

                                                                                                                                                cccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAgg     <  2:3389085/150‑1 (MQ=255)

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

GAAAAGGCCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAGGGGA  >  NC_002947/5889838‑5890134

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