BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,529	A→G	8.9%	K156E (AAG→GAG) ‡	rpsC →	30S ribosomal protein S3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	554,529	0	A	G	8.9%	1037.7 / 61.8	349	K156E (AAG→GAG) ‡	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand): ref base A (161/155); new base G (19/12); total (180/169)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 358 total aligned reads displayed.

GCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGAT  >  NC_002947/554381‑554677

gCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCaa                                                                                                                                                     <  1:9096164/150‑1 (MQ=255)

gCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCaa                                                                                                                                                     <  1:4470718/150‑1 (MQ=255)

 ccTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCGGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAg                                                                                                                                                    <  1:3814068/150‑1 (MQ=255)

  cTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAgg                                                                                                                                                   >  2:1849656/1‑150 (MQ=255)

  cTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAgg                                                                                                                                                   >  2:6567725/1‑150 (MQ=255)

  cTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAgg                                                                                                                                                   <  2:3645939/150‑1 (MQ=255)

  cTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAgg                                                                                                                                                   >  2:5996812/1‑150 (MQ=255)

  cTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCACAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAgg                                                                                                                                                   <  2:7888946/150‑1 (MQ=255)

        acaTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCa                                                                                                                                             >  2:4172279/1‑150 (MQ=255)

        acaTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCa                                                                                                                                             >  1:4963507/1‑150 (MQ=255)

        acaTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCa                                                                                                                                             >  2:3117918/1‑150 (MQ=255)

         caTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCaa                                                                                                                                            <  1:5057088/150‑1 (MQ=255)

          aTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAg                                                                                                                                           >  2:6185808/1‑150 (MQ=255)

          aTCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAg                                                                                                                                           >  2:8238926/1‑150 (MQ=255)

           tCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGa                                                                                                                                          <  2:3078509/150‑1 (MQ=255)

           tCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGa                                                                                                                                          <  1:2276059/150‑1 (MQ=255)

                aTCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAg                                                                                                                                     <  1:100890/150‑1 (MQ=255)

                aTCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAg                                                                                                                                     >  1:7727520/1‑150 (MQ=255)

                aTCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAg                                                                                                                                     <  2:6926149/150‑1 (MQ=255)

                aTCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAg                                                                                                                                     <  2:213998/150‑1 (MQ=255)

                 tCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGAAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAgg                                                                                                                                    <  2:6090528/150‑1 (MQ=255)

                   gAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTg                                                                                                                                  >  2:2698129/1‑150 (MQ=255)

                   gAAGAGATCCGCAAGCCAGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATAAATATCAAGGTg                                                                                                                                  >  1:5299767/1‑150 (MQ=255)

                    aaGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGa                                                                                                                                 <  2:3272731/150‑1 (MQ=255)

                    aaGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGa                                                                                                                                 >  1:1930446/1‑150 (MQ=255)

                    aaGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGa                                                                                                                                 <  1:2654963/150‑1 (MQ=255)

                       agaTCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCg                                                                                                                              <  2:7634486/150‑1 (MQ=255)

                         aTCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGgt                                                                                                                            >  1:3580784/1‑150 (MQ=255)

                            cGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCaa                                                                                                                                                     <  1:4307944/122‑1 (MQ=255)

                            cGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCaa                                                                                                                                                     >  2:4307944/1‑122 (MQ=255)

                            cGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGgtcgt                                                                                                                         >  1:7566345/1‑150 (MQ=255)

                             gCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGgtcgtc                                                                                                                        >  1:1670833/1‑150 (MQ=255)

                              cAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGtct                                                                                                                       >  1:2828048/1‑150 (MQ=255)

                                   cGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTcggcg                                                                                                                  >  2:6799596/1‑150 (MQ=255)

                                    ggAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTcggcgg                                                                                                                 <  1:8287635/150‑1 (MQ=255)

                                       aCTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGc                                                                                                              >  2:5151385/1‑150 (MQ=255)

                                         tCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTg                                                                                                            >  1:3272111/1‑150 (MQ=255)

                                          cGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTga                                                                                                           <  2:392/150‑1 (MQ=255)

                                           gACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTgag                                                                                                          >  1:1271972/1‑150 (MQ=255)

                                           gACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTgag                                                                                                          >  1:556603/1‑150 (MQ=255)

                                           gACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTgag                                                                                                          >  2:4460200/1‑150 (MQ=255)

                                           gACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTgag                                                                                                          >  2:6398428/1‑150 (MQ=255)

                                           gACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTgag                                                                                                          >  2:5240185/1‑150 (MQ=255)

                                           gACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGCCTCGGCGGTGCTgag                                                                                                          >  1:6570490/1‑150 (MQ=255)

                                            aCGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTgaga                                                                                                         <  1:6162773/150‑1 (MQ=255)

                                            aCGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTgaga                                                                                                         <  2:7727881/150‑1 (MQ=255)

                                                tATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGc                                                                                                     >  1:8993774/1‑150 (MQ=255)

                                                tATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGc                                                                                                     >  2:8746069/1‑150 (MQ=255)

                                                tATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGc                                                                                                     <  2:1754844/150‑1 (MQ=255)

                                                  tGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTc                                                                                                   >  1:7531268/1‑150 (MQ=255)

                                                    cTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGt                                                                                                 >  1:7396684/1‑150 (MQ=255)

                                                      ggTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTAc                                                                                               <  1:5959908/150‑1 (MQ=255)

                                                         tGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCACCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGa                                                                                            <  1:1726768/150‑1 (MQ=255)

                                                            gcAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTg                                                                                         <  2:4542657/150‑1 (MQ=255)

                                                            gcAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTg                                                                                         >  1:4807694/1‑150 (MQ=255)

                                                            gcAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTg                                                                                         >  1:7155757/1‑150 (MQ=255)

                                                                tcgatcgCTCAGCAGCTCGAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTAt                                                                                     >  2:8392482/7‑150 (MQ=255)

                                                                      gCTCAGGAGGTGGAACGCCGTTTAATGTTCAGTCTCGCCATGAAGCGTGCGGTACATATCGCCTTGCGTATTGGTTCATAGGGAATCAAGATCGAGGTGAGCGGTCGTCTCGGCTGTGCTGAGTTTGCTAGTAGCTAGTGGTTTCGCGaa                                                                               <  2:1304361/150‑1 (MQ=255)

                                                                        tCAGCAGCTCGAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAgg                                                                             >  1:4382578/1‑150 (MQ=255)

                                                                           gcagcTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCg                                                                          >  2:7803381/1‑150 (MQ=255)

                                                                               ctcgAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCgcgtt                                                                      >  2:7209842/4‑149 (MQ=255)

                                                                                tGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGc                                                                     >  2:1918104/1‑150 (MQ=255)

                                                                                 ggAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTcc                                                                    >  1:1446202/1‑150 (MQ=255)

                                                                                 ggAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGcc                                                                    >  2:4437608/1‑150 (MQ=255)

                                                                                  gAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTAc                                                                                               <  1:6436829/122‑1 (MQ=255)

                                                                                  gAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCct                                                                   <  2:1451301/150‑1 (MQ=255)

                                                                                   aaCGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTccgc                                                                  <  2:1003817/150‑3 (MQ=255)

                                                                                    aCGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCcgct                                                                 >  2:4238227/1‑150 (MQ=255)

                                                                                    aCGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGCCCGCCATGCGTCTTGGTGCCAAGGGCATCACGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCctct                                                                 <  2:7608383/150‑1 (MQ=255)

                                                                                    aCGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTTGTGTGcagag                                                                 >  2:5388717/1‑146 (MQ=255)

                                                                                    aCGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCctct                                                                 <  1:2503165/150‑1 (MQ=255)

                                                                                        cGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGCTTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGcac                                                             >  2:8156513/1‑150 (MQ=255)

                                                                                             ggcattCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCt                                                        <  2:668120/146‑1 (MQ=255)

                                                                                               cattCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGc                                                      <  1:5827995/148‑1 (MQ=255)

                                                                                                   ccGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGc                                                  <  2:5262654/150‑1 (MQ=255)

                                                                                                    cGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGcc                                                 <  1:5706758/150‑1 (MQ=255)

                                                                                                      tCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGa                                               >  2:3435799/1‑150 (MQ=255)

                                                                                                        gcgcCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGata                                             >  1:1328603/1‑150 (MQ=255)

                                                                                                         cgcCCTGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCACGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGatat                                            <  2:6448277/150‑1 (MQ=255)

                                                                                                            cATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGa                                         <  1:2946988/150‑1 (MQ=255)

                                                                                                            cATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGa                                         <  1:2550704/150‑1 (MQ=255)

                                                                                                            cATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGa                                         <  1:5382181/150‑1 (MQ=255)

                                                                                                            cATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGa                                         <  2:7112708/150‑1 (MQ=255)

                                                                                                             aTGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGAc                                        >  1:6793038/1‑150 (MQ=255)

                                                                                                              tGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACt                                       >  2:2508113/1‑150 (MQ=255)

                                                                                                              tGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACt                                       <  1:8746069/150‑1 (MQ=255)

                                                                                                                   cGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTacaac                                  <  2:4951818/150‑1 (MQ=255)

                                                                                                                   cGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGATCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTacaac                                  <  2:27832/150‑1 (MQ=255)

                                                                                                                           aCAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAGGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGa                          >  2:7335066/1‑150 (MQ=255)

                                                                                                                           aCAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGa                          >  1:1112558/1‑150 (MQ=255)

                                                                                                                              gAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGc                       >  1:3992278/1‑150 (MQ=255)

                                                                                                                                   ccATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTcaca                  <  1:449239/150‑1 (MQ=255)

                                                                                                                                   ccATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTcaca                  <  1:41317/150‑1 (MQ=255)

                                                                                                                                    cATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCacac                 <  1:5178598/150‑1 (MQ=255)

                                                                                                                                      tGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCAcacca               <  2:5700847/150‑1 (MQ=255)

                                                                                                                                      tGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCAcacca               >  1:6654929/1‑150 (MQ=255)

                                                                                                                                            ttGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCAcaccgtaagtg         <  1:6536938/150‑8 (MQ=255)

                                                                                                                                               gTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGgtg      >  2:5600519/1‑150 (MQ=255)

                                                                                                                                               gTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGgtg      <  1:7350100/150‑1 (MQ=255)

                                                                                                                                                  ccAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGa   >  2:8945253/1‑150 (MQ=255)

GCCTGTGCACATCAACATCGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGAT  >  NC_002947/554381‑554677

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