BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,548	T→G	9.1%	V162G (GTG→GGG) ‡	rpsC →	30S ribosomal protein S3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	554,548	0	T	G	9.1%	1083.2 / 76.8	342	V162G (GTG→GGG) ‡	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand): ref base T (152/159); new base G (19/12); total (171/171)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 345 total aligned reads displayed.

CGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCT  >  NC_002947/554399‑554697

cGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGt                                                                                                                                                       <  1:6528739/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

         ccGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGgtcg                                                                                                                                              >  2:3500781/1‑150 (MQ=255)

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

              aGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTcg                                                                                                                                         >  1:455508/1‑150 (MQ=255)

              aGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTcg                                                                                                                                         <  1:3837909/150‑1 (MQ=255)

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

                ccGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTcggc                                                                                                                                       >  2:2679006/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                      ttGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCg                                                                                                                 >  2:2299772/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                     cTCAGCAGCTAGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCt                                                                                                                                 >  2:6446061/1‑119 (MQ=255)

                                                     cTCAGCAGCTAGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCt                                                                                                                                 <  1:6446061/119‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                        tgtAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCacac                                                                               <  1:735199/150‑1 (MQ=255)

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

                                                                          tAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACAccc                                                                             <  2:1626993/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                               ttCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGc                                                                                         <  1:517454/133‑1 (MQ=255)

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

                                                                               ttCCGTCGCGCCATGAAGCGTGCCGTTCAGTCGGCAATGCGTCTCGGTGCTGAAGGTATCCGTATCAACTGCTCCGGTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGt                                                                        >  1:6728003/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                                                                      ggTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTa                                                 <  2:7727520/150‑1 (MQ=255)

                                                                                                       gTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTAc                                                <  2:6925354/150‑1 (MQ=255)

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

                                                                                                         aCAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACCCCCTGCGTGCCGATCTCGACTACAACACCTACGa                                              >  1:5012552/1‑150 (MQ=255)

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

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

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

                                                                                                                          ttGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACg                             <  1:3759094/150‑1 (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:4781050/150‑1 (MQ=255)

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

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

                                                                                                                                     ggCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGAGTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgt                  <  2:3291531/150‑1 (MQ=255)

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

                                                                                                                                       cATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGgtgt                >  2:2226513/1‑150 (MQ=255)

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

                                                                                                                                               tCCAGGTGAGCGGTTGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt        <  2:5299767/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                                                     tGAGCGGTCGTCTCGGCGGTTCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCt  <  2:6676876/150‑1 (MQ=255)

CGAAGAGATCCGCAAGCCGGAACTCGACGCTATGCTGGTTGCGCAGAGCGTAGCTCAGCAGCTGGAACGCCGTGTAATGTTCCGTCGCGCCATGAAGCGTGCGGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCT  >  NC_002947/554399‑554697

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