BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,612,408	A→T	5.4%	K29M (AAG→ATG) ‡	PP_4079 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,612,408	0	A	T	5.4%	1193.6 / 33.8	346	K29M (AAG→ATG) ‡	PP_4079	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (170/157); new base T (11/8); total (181/165)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 353 total aligned reads displayed.

GGACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGC  >  NC_002947/4612259‑4612556

ggACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCaa                                                                                                                                                      <  1:4231638/150‑1 (MQ=255)

 gACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAg                                                                                                                                                     <  1:2470164/150‑1 (MQ=255)

   ccAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGt                                                                                                                                                   >  2:1184/1‑150 (MQ=255)

   ccAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGt                                                                                                                                                   <  1:158294/150‑1 (MQ=255)

    cAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTg                                                                                                                                                  >  2:6072829/1‑150 (MQ=255)

    cAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTg                                                                                                                                                  >  2:4119598/1‑150 (MQ=255)

    cAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTg                                                                                                                                                  >  1:4510357/1‑150 (MQ=255)

       cTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGgcg                                                                                                                                               >  1:1986305/1‑150 (MQ=255)

        ttCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGgctc                                                                                                                                              <  2:1249487/150‑3 (MQ=255)

        ttCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGgcgc                                                                                                                                              >  1:2087826/1‑150 (MQ=255)

         tcgtcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCt                                                                                                                                             <  2:5607338/150‑1 (MQ=255)

          cgtcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTg                                                                                                                                            <  2:587713/150‑1 (MQ=255)

           gtcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTgg                                                                                                                                           <  2:6383630/150‑1 (MQ=255)

            tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                          <  1:3439952/150‑1 (MQ=255)

            tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                          >  2:6915717/1‑150 (MQ=255)

            tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                          >  1:4028659/1‑150 (MQ=255)

            tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                          <  2:714094/150‑1 (MQ=255)

                  aCCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTc                                                                                                                                    <  1:6832310/150‑1 (MQ=255)

                    cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTc                                                                                                                                    >  2:6616425/1‑148 (MQ=255)

                    cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTc                                                                                                                                    <  1:6616425/148‑1 (MQ=255)

                    aaCGATGGCCAAATCAGTCGCGTTGTCCGCTTCCGGGGGCGGTCATGAGGCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGt                                                                                                                                  >  2:2561582/2‑150 (MQ=255)

                       gATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTAc                                                                                                                               <  1:5859810/150‑1 (MQ=255)

                         tGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCTTGGTGGCTCTgc                                                                                                                                           >  2:1143090/1‑135 (MQ=255)

                              aaaTCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCgg                                                                                                                        <  1:2236695/150‑1 (MQ=255)

                              aaaTCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCgg                                                                                                                        <  1:4112986/150‑1 (MQ=255)

                                aTCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGgc                                                                                                                      >  2:2887331/1‑150 (MQ=255)

                                   aGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGgccgc                                                                                                                   >  2:885876/1‑150 (MQ=255)

                                     tCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTg                                                                                                                                            >  2:3797607/1‑123 (MQ=255)

                                       ccGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCaa                                                                                                               >  1:7193594/1‑150 (MQ=255)

                                       ccGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCaa                                                                                                               <  2:5914658/150‑1 (MQ=255)

                                       ccGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCaa                                                                                                               <  2:7257392/150‑1 (MQ=255)

                                       ccGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCaa                                                                                                               >  1:2659356/1‑150 (MQ=255)

                                        cGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAg                                                                                                              <  1:4857211/150‑1 (MQ=255)

                                          ttGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGa                                                                                                            >  1:4333532/1‑150 (MQ=255)

                                          ttGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGa                                                                                                            >  1:5059719/1‑150 (MQ=255)

                                             tCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                          >  1:5975490/1‑117 (MQ=255)

                                             tCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaa                                                                                                         >  2:5783763/1‑150 (MQ=255)

                                              ccGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaac                                                                                                        <  1:6540924/150‑1 (MQ=255)

                                                gCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaaccc                                                                                                      >  2:4887919/1‑150 (MQ=255)

                                                  ttCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAc                                                                                                           <  2:2372679/143‑1 (MQ=255)

                                                  ttCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAc                                                                                                           >  1:2372679/1‑143 (MQ=255)

                                                    ccGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTAc                                                                                                                               >  1:2865572/1‑121 (MQ=255)

                                                              tCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGt                                                                                        <  2:1529428/150‑1 (MQ=255)

                                                                 tGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCa                                                                                     <  1:3500697/150‑1 (MQ=255)

                                                                 tGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCa                                                                                     <  1:3219179/150‑1 (MQ=255)

                                                                 tGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCa                                                                                     <  1:935045/150‑1 (MQ=255)

                                                                    gccgccTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCa                                                                                  <  1:6134746/150‑1 (MQ=255)

                                                                         cTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAAc                                                                             <  2:6261081/150‑1 (MQ=255)

                                                                               gTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGc                                                                       <  2:752278/150‑1 (MQ=255)

                                                                                tCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGcc                                                                      >  1:6436350/1‑150 (MQ=255)

                                                                                       tGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCt                                                               <  1:7389511/150‑1 (MQ=255)

                                                                                           ccTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAAc                                                           >  2:6143987/1‑150 (MQ=255)

                                                                                            ctctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACa                                                          <  1:1693584/150‑1 (MQ=255)

                                                                                            ctctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACa                                                          <  2:3366751/150‑1 (MQ=255)

                                                                                            ctctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACa                                                          <  1:1162489/150‑1 (MQ=255)

                                                                                            ctctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACa                                                          >  1:6327062/1‑150 (MQ=255)

                                                                                               tgctggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCa                                                       >  1:518852/1‑150 (MQ=255)

                                                                                                 ctggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAAc                                                     <  2:3656658/150‑1 (MQ=255)

                                                                                                  tggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAAcc                                                    >  1:7246758/1‑150 (MQ=255)

                                                                                                      cgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGaa                                                <  1:1230035/150‑1 (MQ=255)

                                                                                                          ctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAgt                                            <  2:248559/150‑1 (MQ=255)

                                                                                                            ggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAgtgt                                          <  2:735931/150‑1 (MQ=255)

                                                                                                             gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTc                                         >  1:1038701/1‑150 (MQ=255)

                                                                                                                ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGAc                                      >  1:4973444/1‑150 (MQ=255)

                                                                                                                    tGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCaa                                                      >  1:3751009/1‑130 (MQ=255)

                                                                                                                    tGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTgg                                  >  2:4793881/1‑150 (MQ=255)

                                                                                                                      ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcg                                >  2:4347949/1‑150 (MQ=255)

                                                                                                                      ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcg                                >  2:3893383/1‑150 (MQ=255)

                                                                                                                       cAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTc                                                              <  2:2597208/119‑1 (MQ=255)

                                                                                                                       cAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcgg                               >  2:3254228/1‑150 (MQ=255)

                                                                                                                             gCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATcc                         <  2:2940711/150‑1 (MQ=255)

                                                                                                                             gCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATcc                         <  1:2186714/150‑1 (MQ=255)

                                                                                                                             gCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATcc                         <  1:5478647/150‑1 (MQ=255)

                                                                                                                              cACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATccc                        >  1:1008403/1‑150 (MQ=255)

                                                                                                                               aCGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCt                       >  2:5239305/1‑150 (MQ=255)

                                                                                                                               aCGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCt                       >  2:4747525/1‑150 (MQ=255)

                                                                                                                                 gTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga                     <  2:1805751/150‑1 (MQ=255)

                                                                                                                                 gTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga                     <  1:5086006/150‑1 (MQ=255)

                                                                                                                                  tCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgag                    >  1:3022134/1‑150 (MQ=255)

                                                                                                                                   ccTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga                   <  1:677928/150‑1 (MQ=255)

                                                                                                                                   ccTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga                   <  2:1877850/150‑1 (MQ=255)

                                                                                                                                   ccTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga                   <  2:5069516/150‑1 (MQ=255)

                                                                                                                                         gCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATcc             >  2:351785/1‑150 (MQ=255)

                                                                                                                                         ccAGCGCCACCTTGGTGGCTCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATcc             >  1:3812066/2‑150 (MQ=255)

                                                                                                                                          cAGCGCCACCTTGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCCACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCt            >  1:4841813/1‑150 (MQ=255)

                                                                                                                                           agagCCACCAAGGTGGCGCTGGACCTCTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTg           >  2:1695769/1‑150 (MQ=255)

                                                                                                                                             agCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGa         >  1:2591373/1‑150 (MQ=255)

                                                                                                                                              gCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGaa        <  2:1395888/150‑1 (MQ=255)

                                                                                                                                              gCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGaa        <  1:4694364/150‑1 (MQ=255)

                                                                                                                                               ccaccaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAc       <  2:1642930/150‑1 (MQ=255)

                                                                                                                                                caccaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAcc      <  1:7086325/150‑1 (MQ=255)

                                                                                                                                                 accaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAccc     <  2:3071016/150‑1 (MQ=255)

                                                                                                                                                 accaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAccc     >  2:2261048/1‑150 (MQ=255)

                                                                                                                                                  ccaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTc                                         >  1:166702/1‑113 (MQ=255)

                                                                                                                                                  ccaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAcccc    >  1:1439799/1‑149 (MQ=255)

                                                                                                                                                  ccaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCa    <  1:4147077/150‑1 (MQ=255)

                                                                                                                                                  ccaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGGTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCa    >  2:298195/1‑150 (MQ=255)

                                                                                                                                                   caAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAg   >  1:4716932/1‑150 (MQ=255)

                                                                                                                                                   caAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAg   >  1:7442842/1‑150 (MQ=255)

                                                                                                                                                   caAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAg   >  1:6221581/1‑150 (MQ=255)

GGACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGC  >  NC_002947/4612259‑4612556

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