BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,650	A→C	6.7%	Y196S (TAC→TCC) ‡	rpsC →	30S ribosomal protein S3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	554,650	0	A	C	6.7%	546.7 / 20.6	173	Y196S (TAC→TCC) ‡	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand): ref base A (59/102); new base C (4/8); total (63/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 175 total aligned reads displayed.

GGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGCCAAATGTTGCAAC  >  NC_002947/554501‑554799

ggTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTa                                                                                                                                                       >  1:4863801/1‑150 (MQ=255)

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

 gTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTAc                                                                                                                                                      >  1:2426357/1‑150 (MQ=255)

       aaCGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTCTCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCt                                                                                                                                                >  2:3571398/1‑150 (MQ=255)

       aaCGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCt                                                                                                                                                >  2:3867025/1‑150 (MQ=255)

         cGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCTTGCCGATATCGACTACAACACCTACGAAGCTca                                                                                                                                              <  2:3336814/150‑1 (MQ=255)

         cGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTca                                                                                                                                              <  2:4151082/150‑1 (MQ=255)

          gCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTcac                                                                                                                                             >  2:1719370/1‑150 (MQ=255)

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

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

                  tATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTa                                                                                                                                     >  2:2412299/1‑150 (MQ=255)

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

                            aaGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATc                                                                                                                           <  1:3030023/150‑1 (MQ=255)

                             aGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCg                                                                                                                          <  2:99448/150‑1 (MQ=255)

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

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

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

                                gCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCggcg                                                                                                                       <  2:3194735/150‑3 (MQ=255)

                                  aTCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGtgtg                                                                                                                     >  2:1117442/1‑150 (MQ=255)

                                  aTCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGtgtg                                                                                                                     >  2:1916494/1‑150 (MQ=255)

                                   tCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGa                                                                                                                    <  2:4504061/150‑1 (MQ=255)

                                   tCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGa                                                                                                                    <  2:3902990/150‑1 (MQ=255)

                                   tCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGa                                                                                                                    <  1:351426/150‑1 (MQ=255)

                                      aGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAgg                                                                                                                 >  2:188173/1‑150 (MQ=255)

                                         tCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGTTCGTGTGCCTCTGCACACCCTGCGTGCCGTTATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGttt                                                                                                              <  2:3678635/150‑1 (MQ=255)

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

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

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

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

                                          ccAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTg                                                                                                             <  2:3557085/150‑1 (MQ=255)

                                            aGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGa                                                                                                           <  2:193830/150‑1 (MQ=255)

                                                  gCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCa                                                                                                     >  2:3089063/1‑150 (MQ=255)

                                                   cGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCaa                                                                                                    <  2:1998485/150‑1 (MQ=255)

                                                    ggTCGTCTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCaag                                                                                                   <  2:2287628/150‑2 (MQ=255)

                                                    ggTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCaaa                                                                                                   <  1:2701612/150‑1 (MQ=255)

                                                      tcgtcTTGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCAAggg                                                                                                 <  1:3311039/150‑1 (MQ=255)

                                                      tcgtcTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAgg                                                                                                 >  2:4305277/1‑150 (MQ=255)

                                                      tcgtcTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAgg                                                                                                 >  2:3951318/1‑150 (MQ=255)

                                                         tcttGGCGGCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCAAGGGCGa                                                                                              >  1:4038895/1‑150 (MQ=255)

                                                         tctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGa                                                                                              <  2:2267106/150‑1 (MQ=255)

                                                         tctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGa                                                                                              <  2:2267132/150‑1 (MQ=255)

                                                          ctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTCCGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGaa                                                                                             <  2:836752/150‑1 (MQ=255)

                                                          ctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGaa                                                                                             <  2:1274372/150‑1 (MQ=255)

                                                          ctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGaa                                                                                             <  1:2318747/150‑1 (MQ=255)

                                                          ctcGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGaa                                                                                             <  2:3166057/150‑1 (MQ=255)

                                                              gcggcgCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCAAGGGCGagatcc                                                                                         <  2:24053/150‑6 (MQ=255)

                                                               cggCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCAAGGGCGagatcct                                                                                        <  2:4466956/150‑7 (MQ=255)

                                                               cggCGCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTGGGTATTCAAGGGCGagatcct                                                                                        <  2:3305162/150‑7 (MQ=255)

                                                                ggTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTAtt                                                                                       <  2:1722938/150‑1 (MQ=255)

                                                                ggTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTAtt                                                                                       <  2:2496037/150‑1 (MQ=255)

                                                                 gTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtg                                                                                      <  2:564557/150‑1 (MQ=255)

                                                                 gTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtg                                                                                      <  1:2630505/150‑1 (MQ=255)

                                                                 gTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCGACGAAGCTCACACCACTTCCGGTGTGATCGGTGTGAAGGTTTGGATCTTCAACGGCGAAGTTATtg                                                                                      <  2:564618/150‑1 (MQ=255)

                                                                   gCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTgg                                                                                                            <  2:1306714/126‑1 (MQ=255)

                                                                   gCTGAAATCGCACGTATGGAATGGTATCGCGAAGGTCGCGTTCCGCTTCACACGCTGCGTGCCGATATCGACTACGGTACGGCAGAAGCGCAGACCGCTTACGGCATCTGCGGCGTGAAGGTTTgg                                                                                                            >  1:1306714/1‑126 (MQ=255)

                                                                     tGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATtggtgg                                                                                  >  2:789830/1‑150 (MQ=255)

                                                                         aTTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGc                                                                              <  2:2880940/150‑1 (MQ=255)

                                                                         aTTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGc                                                                              <  1:2864161/150‑1 (MQ=255)

                                                                              tCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCaaga                                                                         >  2:2431395/1‑150 (MQ=255)

                                                                                  aCCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAgaagaa                                                                     <  1:1111423/150‑1 (MQ=255)

                                                                                   ccGAGTGGTCGCGCGAAGGTCGTGTGCCTCTGCACCCCCTGCGTGCCGAGATCGACTACAACACCTACGAAGCTCACACCACTTACGGGGTGATCGGTGTGAAGGGTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCACGAAGAAc                                                                    >  2:4576269/1‑150 (MQ=255)

                                                                                    cGAGTTGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGTCTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACt                                                                   <  1:1912871/150‑1 (MQ=255)

                                                                                     gAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTg                                                                  >  2:1829985/1‑150 (MQ=255)

                                                                                         ggTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAAc                                                              <  1:100444/150‑1 (MQ=255)

                                                                                             tCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACcaca                                                          >  1:2789411/1‑150 (MQ=255)

                                                                                              cgcgAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGTAGTTATTGGTGGTCGCCAAGAAGAACTGTAACCACaa                                                         <  2:4827371/150‑1 (MQ=255)

                                                                                               gcgAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAg                                                        >  2:2082618/1‑150 (MQ=255)

                                                                                               gcgAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAg                                                        >  1:4212197/1‑150 (MQ=255)

                                                                                                 gAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCa                                                      >  1:3799763/1‑150 (MQ=255)

                                                                                                 gAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTCACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCa                                                      >  1:1216855/1‑150 (MQ=255)

                                                                                                    ggTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCa                                                   >  1:5063600/1‑150 (MQ=255)

                                                                                                         tgtgCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGcc                                              >  2:2541856/1‑150 (MQ=255)

                                                                                                         tgtgCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGcc                                              <  2:4879029/150‑1 (MQ=255)

                                                                                                         tgtgCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGcc                                              >  1:3021150/1‑150 (MQ=255)

                                                                                                          gtgCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCt                                             <  1:3446098/150‑1 (MQ=255)

                                                                                                            gCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCg                                           <  2:753645/150‑1 (MQ=255)

                                                                                                             ccTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGt                                          <  1:3951139/150‑1 (MQ=255)

                                                                                                              ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAAcc                                                             <  2:932734/130‑1 (MQ=255)

                                                                                                              ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAAcc                                                             >  1:932734/1‑130 (MQ=255)

                                                                                                              ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTa                                         >  1:4450452/1‑150 (MQ=255)

                                                                                                              ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTa                                         >  2:5031231/1‑150 (MQ=255)

                                                                                                              ctctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTa                                         >  1:548960/1‑150 (MQ=255)

                                                                                                               tctGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTaa                                        <  2:1803863/150‑1 (MQ=255)

                                                                                                                  gCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTaaaaa                                     <  1:727764/150‑1 (MQ=255)

                                                                                                                  gCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTaaaaa                                     >  2:653568/1‑150 (MQ=255)

                                                                                                                  gCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTaaaaa                                     >  2:53703/1‑150 (MQ=255)

                                                                                                                      acCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAgct                                 <  1:1928389/150‑1 (MQ=255)

                                                                                                                        ccTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAgctgc                               <  1:1320908/150‑1 (MQ=255)

                                                                                                                        ccTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAgctgc                               <  1:4585211/150‑1 (MQ=255)

                                                                                                                                 cGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAgg                      <  2:2571904/150‑1 (MQ=255)

                                                                                                                                  gATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAggg                     <  1:1719370/150‑1 (MQ=255)

                                                                                                                                        gACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCTGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGc               >  1:2411082/1‑150 (MQ=255)

                                                                                                                                        gACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGc               <  2:4863801/150‑1 (MQ=255)

                                                                                                                                        gACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGc               >  1:2887378/1‑150 (MQ=255)

                                                                                                                                         aCTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGcc              >  1:1327770/1‑150 (MQ=255)

                                                                                                                                               acacCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGCCAAATGt        <  1:868543/150‑1 (MQ=255)

                                                                                                                                               acacCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGCCAAATGt        >  1:370056/1‑150 (MQ=255)

                                                                                                                                                cacCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGCCAAATGtt       <  1:3143963/150‑1 (MQ=255)

                                                                                                                                                    tACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGCCAAATGTTGCaa   >  1:4968967/1‑150 (MQ=255)

                                                                                                                                                     aCGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGCCAAATGTTGCAAc  >  1:5063647/1‑150 (MQ=255)

GGTACAGAACGCCATGCGTATTGGTGCCAAGGGCATCAAGATCCAGGTGAGCGGTCGTCTCGGCGGTGCTGAGATTGCTCGTACCGAGTGGTATCGCGAAGGTCGTGTGCCTCTGCACACCCTGCGTGCCGATATCGACTACAACACCTACGAAGCTCACACCACTTACGGTGTGATCGGTGTGAAGGTTTGGATCTTCAAAGGCGAAGTTATTGGTGGTCGCCAAGAAGAACTGAAACCACAAGCACCAGCGCCTCGTAAAAAAGCTGCTAAGTAAGGGGTACGCCAAATGTTGCAAC  >  NC_002947/554501‑554799

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