BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,551	G→C	9.3%	S163T (AGC→ACC) ‡	rpsC →	30S ribosomal protein S3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	554,551	0	G	C	9.3%	1089.3 / 77.2	334	S163T (AGC→ACC) ‡	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand): ref base G (145/158); new base C (19/12); total (164/170)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 337 total aligned reads displayed.

AGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCA  >  NC_002947/554402‑554700

agagATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAg                                                                                                                                                       <  2:1978999/150‑1 (MQ=255)

     tCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGgtc                                                                                                                                                  <  1:5920035/150‑1 (MQ=255)

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

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

            gCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTcgg                                                                                                                                           <  2:7890944/150‑1 (MQ=255)

                 aaCTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTg                                                                                                                                      <  2:3580784/150‑1 (MQ=255)

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

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

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

                    tCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGAGGTCGTCTCGGCGGTGCTg                                                                                                                                   >  1:3310108/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                         gCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGAtt                                                                                                                              >  2:4251888/1‑150 (MQ=255)

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

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

                            aTGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCt                                                                                                                           <  1:3146511/150‑1 (MQ=255)

                                    tGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGa                                                                                                                   <  1:1153504/150‑1 (MQ=255)

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

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

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

                                                 gCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGaa                                                                                                      <  2:8889900/150‑1 (MQ=255)

                                                 gCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGaa                                                                                                      <  2:1363752/150‑1 (MQ=255)

                                                 gCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGaa                                                                                                      >  2:5227613/1‑150 (MQ=255)

                                                  cTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCt                                                                                                                           <  2:1822878/128‑1 (MQ=255)

                                                  cTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAg                                                                                                     >  2:1896246/1‑150 (MQ=255)

                                                   tCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAgg                                                                                                    >  2:7685022/1‑150 (MQ=255)

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

                                                      gcagcTCGAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCg                                                                                                 >  1:7483754/1‑150 (MQ=255)

                                                        agcTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCgtg                                                                                               <  1:8288523/150‑1 (MQ=255)

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

                                                            cgAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTcc                                                                                           <  1:6604352/149‑1 (MQ=255)

                                                             gAACGCCGTGTGGCATTCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTccg                                                                                          >  2:8313519/1‑149 (MQ=255)

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

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

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

                                                                       tggcattCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCacacgc                                                                                >  2:2871986/6‑148 (MQ=255)

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

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

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

                                                                          tGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGc                                                                             <  1:5151385/150‑1 (MQ=255)

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

                                                                           gTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCg                                                                            >  2:7247465/1‑150 (MQ=255)

                                                                           gTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCg                                                                            >  2:4231423/1‑150 (MQ=255)

                                                                           gTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCg                                                                            >  1:444472/1‑150 (MQ=255)

                                                                           attCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCg                                                                            >  1:611220/2‑150 (MQ=255)

                                                                            ttCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGt                                                                           >  1:7580722/1‑150 (MQ=255)

                                                                            ttCCGTCGCGCCATGAAGCGAGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGt                                                                           <  2:8569376/150‑1 (MQ=255)

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

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

                                                                                 tCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGa                                                                      <  2:2243454/150‑1 (MQ=255)

                                                                                  cgcgCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGat                                                                     <  2:1604712/150‑1 (MQ=255)

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

                                                                                     gcCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATGTc                                                                  >  2:4391600/1‑150 (MQ=255)

                                                                                     gcCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATc                                                                  >  1:9105871/1‑150 (MQ=255)

                                                                                       cATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGa                                                                >  1:1102822/1‑150 (MQ=255)

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

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

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

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

                                                                                                 gCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACAcc                                                      <  2:2779292/150‑1 (MQ=255)

                                                                                                  cGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCt                                                     <  1:2356378/150‑1 (MQ=255)

                                                                                                     tACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACg                                                  >  1:6831599/1‑150 (MQ=255)

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

                                                                                                          aaCGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGAGGGGCTGAGATTGCTCTTACCGAGTCGTATCGCGAAGGTCGTGTGCCTCTGCACACACTGCGAGACGATATCGACTACAACACCTACGAAGca                                             >  1:7372041/1‑149 (MQ=255)

                                                                                                              ccATGCGTATTGGTGCCTAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTcaca                                         <  2:6118005/150‑1 (MQ=255)

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

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

                                                                                                                    gTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACtt                                   >  1:4328776/1‑150 (MQ=255)

                                                                                                                    gTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACtt                                   >  1:2877079/1‑150 (MQ=255)

                                                                                                                       ttGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACg                                <  1:3759094/150‑1 (MQ=255)

                                                                                                                         ggTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTAGTACCGAGTTGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGgt                              <  2:6578930/150‑1 (MQ=255)

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

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

                                                                                                                           tGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGgtgt                            <  1:5340089/150‑1 (MQ=255)

                                                                                                                                aGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCg                       >  2:50615/1‑150 (MQ=255)

                                                                                                                                 gggCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCgg                      <  2:4781078/150‑1 (MQ=255)

                                                                                                                                  ggCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgt                     <  2:3132413/150‑1 (MQ=255)

                                                                                                                                   gCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgtg                    <  1:7041785/150‑1 (MQ=255)

                                                                                                                                    cATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgtgt                   <  2:2982202/150‑1 (MQ=255)

                                                                                                                                     aTCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGtgtg                  <  2:4879948/150‑1 (MQ=255)

                                                                                                                                        aaGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAg               <  1:6486564/150‑1 (MQ=255)

                                                                                                                                          gATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGt             >  2:6882994/1‑150 (MQ=255)

                                                                                                                                            tCCAGGTGAGCGGTCGTCTCGGCGGTGCTTAGATTGCTAGTACCGAGTGGTATCGCGAATGTAGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACTAAGCTCACACCACTTACTGTGTGATCGGTGTGAAGGttt           <  1:5636007/150‑1 (MQ=255)

                                                                                                                                            tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGACGGttt           >  2:5707651/1‑150 (MQ=255)

                                                                                                                                            tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt           <  1:5996812/150‑1 (MQ=255)

                                                                                                                                            tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt           <  1:6567725/150‑1 (MQ=255)

                                                                                                                                             ccAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTATAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTg          <  1:4252420/150‑1 (MQ=255)

                                                                                                                                             ccAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTg          >  1:6189610/1‑150 (MQ=255)

                                                                                                                                              cAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCGGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTgg         >  1:2622237/1‑150 (MQ=255)

                                                                                                                                               aGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGa        >  1:7947945/1‑150 (MQ=255)

                                                                                                                                                   gAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCtt    >  1:4977043/1‑150 (MQ=255)

                                                                                                                                                    aGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTc   <  1:6774021/150‑1 (MQ=255)

                                                                                                                                                    aGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTc   <  2:1271972/150‑1 (MQ=255)

                                                                                                                                                    aGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTc   <  1:1849656/150‑1 (MQ=255)

AGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCA  >  NC_002947/554402‑554700

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