BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,810,509	A→C	9.1%	noncoding (416/434 nt)	PP_mr35 ←	cspA regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,810,509	0	A	C	9.1%	662.1 / 39.3	228	noncoding (416/434 nt)	PP_mr35	cspA regulatory region
Reads supporting (aligned to +/- strand): ref base A (95/111); new base C (11/11); total (106/122)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

TTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTA  >  NC_002947/2810360‑2810658

ttGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCca                                                                                                                                                       >  2:5164025/1‑150 (MQ=255)

   cgcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACTCGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACaa                                                                                                                                                    <  2:4684632/150‑1 (MQ=255)

   cgcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACaa                                                                                                                                                    <  2:4684295/150‑1 (MQ=255)

     cAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAgg                                                                                                                                                  <  2:703140/150‑1 (MQ=255)

       aaaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGc                                                                                                                                                <  2:8103305/150‑1 (MQ=255)

         aTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCgg                                                                                                                                              <  2:5389456/150‑1 (MQ=255)

            tgAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAACGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCACTCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCCAAGGGCGGGtt                                                                                                                                           >  2:6487128/1‑150 (MQ=255)

             gAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttt                                                                                                                                          >  2:1711208/1‑150 (MQ=255)

              aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGGtttt                                                                                                                                         >  2:2366936/1‑150 (MQ=39)

              aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAAACCGCCACAAGGGCGGGtttt                                                                                                                                         >  1:5632943/1‑150 (MQ=255)

               aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGGttttt                                                                                                                                        <  2:3915594/150‑1 (MQ=39)

               aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  2:1960129/150‑1 (MQ=255)

               aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  1:4002522/150‑1 (MQ=255)

               aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  1:4612709/150‑1 (MQ=255)

               aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  2:7662745/150‑1 (MQ=255)

               aaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        <  2:693400/150‑1 (MQ=255)

                aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGGtttttt                                                                                                                                       <  1:932701/150‑1 (MQ=39)

                aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttttt                                                                                                                                       <  1:2886013/150‑1 (MQ=255)

                aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttttt                                                                                                                                       <  2:3103989/150‑1 (MQ=255)

                  gCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTa                                                                                                                                     <  1:3405060/150‑1 (MQ=255)

                  gCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTa                                                                                                                                     >  2:9122298/1‑150 (MQ=255)

                    cTGCGCGGTGGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACa                                                                                                                                   <  1:1853876/150‑1 (MQ=255)

                    cTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACa                                                                                                                                   <  1:5761263/150‑1 (MQ=255)

                     tGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAg                                                                                                                                  >  1:5971853/1‑150 (MQ=255)

                          ggTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCACTCCCGCCCCTTCAGATGTTGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATc                                                                                                                             >  1:2641552/1‑150 (MQ=255)

                                   gTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACt                                                                                                                    <  1:5993303/150‑1 (MQ=255)

                                           gggTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGa                                                                                                            <  1:6489807/150‑1 (MQ=255)

                                           gggTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGa                                                                                                            <  1:5431417/150‑1 (MQ=255)

                                           gggTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGa                                                                                                            <  1:3012672/150‑1 (MQ=255)

                                            ggTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGaa                                                                                                           <  2:331373/150‑1 (MQ=255)

                                               tCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCt                                                                                                        <  2:3835119/150‑1 (MQ=255)

                                                   tgaatggcacGGGCTTTGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCGTCAGATGATGCATAACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGc                                                                                                    >  1:2500899/8‑150 (MQ=255)

                                                   ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATTCATCACCCAAAAAAAAAACCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGc                                                                                                    >  2:2018895/1‑150 (MQ=255)

                                                   ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                        >  1:7351000/1‑114 (MQ=255)

                                                   ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGAAGAACCTGTGc                                                                                                    >  1:8865460/1‑150 (MQ=255)

                                                   ccGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACACGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGc                                                                                                    >  2:5603620/1‑150 (MQ=255)

                                                       tggcacGGGCTTTGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAAc                                                                                                          <  2:6604462/137‑1 (MQ=255)

                                                       tggcacGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGcc                                                                                                >  2:8251125/4‑150 (MQ=255)

                                                         acacGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAATCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTg                                                                                              >  1:3744922/1‑150 (MQ=255)

                                                            cGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCAt                                                                                           <  2:6090428/150‑1 (MQ=255)

                                                             gggCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCgg                                                                                                                                              <  2:635044/98‑1 (MQ=255)

                                                             gggCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGACTTACTCCGGACGAACCTGTGCAGCCTGCATg                                                                                          >  2:6557528/1‑150 (MQ=255)

                                                             gggCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATg                                                                                          <  2:2259121/150‑1 (MQ=255)

                                                               gCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGcc                                                                                        >  2:7588216/1‑150 (MQ=255)

                                                                 aaaGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCt                                                                                      >  1:6965322/1‑150 (MQ=255)

                                                                   aGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCGGCATGCCCTTc                                                                                    <  1:4139859/150‑1 (MQ=255)

                                                                    gCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCt                                                                                   <  1:3844149/150‑1 (MQ=255)

                                                                          gCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATACAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCcgc                                                                             <  2:5839691/150‑1 (MQ=255)

                                                                          gCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCcgc                                                                             >  2:5536597/1‑150 (MQ=255)

                                                                            cATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCgcgc                                                                           >  2:678703/1‑150 (MQ=255)

                                                                                     tACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCAACCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCAGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcgacg                                                                  >  2:4524934/1‑150 (MQ=255)

                                                                                        aGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAg                                                               <  2:1066872/150‑1 (MQ=255)

                                                                                        aGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCAGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAg                                                               >  2:3021804/1‑150 (MQ=255)

                                                                                          tCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGa                                                             >  1:3280074/1‑150 (MQ=255)

                                                                                          tCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGa                                                             >  1:7896294/1‑150 (MQ=255)

                                                                                           cGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGaa                                                            <  1:4850376/150‑1 (MQ=255)

                                                                                            ggTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAACAAACCCGCCACAAGGGCGGGTTTTTTTACAGCGATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGaaa                                                           <  1:8334970/150‑1 (MQ=255)

                                                                                             gTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAAc                                                          >  2:5712553/1‑150 (MQ=255)

                                                                                             gTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCGTCTGGCCGCGCTCGGCGACGACGGacac                                                          >  1:2759016/1‑150 (MQ=255)

                                                                                                 cGCAGGCAATCCCGCCCCTTCAGAGGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTc                                                      <  2:6220706/150‑1 (MQ=255)

                                                                                                  gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACCGTCt                                                     >  1:6485432/1‑150 (MQ=255)

                                                                                                  gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCt                                                     >  1:1038673/1‑150 (MQ=255)

                                                                                                  gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCt                                                     >  2:3469128/1‑150 (MQ=255)

                                                                                                  gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCt                                                     >  2:6758936/1‑150 (MQ=255)

                                                                                                  gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCGGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCt                                                     >  1:4475468/1‑150 (MQ=255)

                                                                                                   cAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                    >  1:1927956/1‑150 (MQ=255)

                                                                                                   cAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                    <  2:2944586/150‑1 (MQ=255)

                                                                                                   cAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                    <  2:1439248/150‑1 (MQ=255)

                                                                                                    aGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTgg                                                   <  2:4758361/150‑1 (MQ=255)

                                                                                                     gtcAATCCCGCCCCTTCAGATGATGCATCAACCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACCGCTCTCAGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGc                                                  <  1:8381454/148‑1 (MQ=255)

                                                                                                       cAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACCGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCct                                                <  2:1264106/150‑1 (MQ=255)

                                                                                                        aaTCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCctt                                               <  1:5164025/150‑1 (MQ=255)

                                                                                                        aaTCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAAGTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGACACAGTCTGGCctt                                               >  2:7293942/1‑150 (MQ=255)

                                                                                                           cccGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCcttctt                                            >  1:2101160/1‑150 (MQ=255)

                                                                                                              gCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                         >  1:2244062/1‑150 (MQ=255)

                                                                                                              gCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                         <  1:2047640/150‑1 (MQ=255)

                                                                                                              gCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                         <  1:6119738/150‑1 (MQ=255)

                                                                                                              gCCCCGTCAGATGATGCATCACCCACACCACACCCCGCCCCAAGGGCGGGTTTTTTTACAGCTCTCCGAATTACGCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                         <  1:8494332/150‑1 (MQ=255)

                                                                                                               ccccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAg                                        <  1:1711208/150‑1 (MQ=255)

                                                                                                               ccccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTCTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAg                                        <  1:8512920/150‑1 (MQ=255)

                                                                                                                   ttCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGctc                                    >  1:2834070/1‑150 (MQ=255)

                                                                                                                    tCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGctct                                   <  1:6387897/150‑1 (MQ=255)

                                                                                                                        atgatgCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGaa                               <  1:962366/150‑1 (MQ=255)

                                                                                                                        atgatgCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGaa                               <  1:5498176/150‑1 (MQ=255)

                                                                                                                               aTCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCa                        <  2:1498675/150‑1 (MQ=255)

                                                                                                                                     tgtAAAAAAACCCGCCCTTGTGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCg                  >  1:2540593/3‑150 (MQ=39)

                                                                                                                                         aaaaaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTACCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATgg              >  2:1381162/1‑150 (MQ=255)

                                                                                                                                          aaaaaCCCGCCCTTGTGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGcc                                                                               >  1:6323625/1‑84 (MQ=37)

                                                                                                                                          aaaaaCCCGCCCTTGTGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGc             >  1:1405536/1‑150 (MQ=39)

                                                                                                                                          aaaaaCCCGCCCTTGTGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTCTCAGGCTCTTGAAGCCGTCAGATTCGATGGc             >  2:2790461/1‑150 (MQ=255)

                                                                                                                                          aaaaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGc             >  1:5534055/1‑150 (MQ=255)

                                                                                                                                            aatcccGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCtt           >  2:3779371/4‑150 (MQ=255)

                                                                                                                                            aaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCtt           >  1:4103516/1‑150 (MQ=255)

                                                                                                                                            aaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCtt           <  2:8123590/150‑1 (MQ=255)

                                                                                                                                             aaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCttt          >  2:2163850/1‑150 (MQ=255)

                                                                                                                                               cccGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGa        <  1:3165392/150‑1 (MQ=255)

                                                                                                                                                  gCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAgt     <  2:2536702/150‑1 (MQ=255)

                                                                                                                                                     acaAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCTGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTa  >  2:8737281/1‑150 (MQ=255)

                                                                                                                                                     acaAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTa  <  1:3656921/150‑1 (MQ=255)

                                                                                                                                                     acaAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTa  <  1:1470696/150‑1 (MQ=255)

TTGCGCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTA  >  NC_002947/2810360‑2810658

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`