BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,649,375	T→G	5.2%	E576A (GAG→GCG) ‡	PP_1446 ←	TonB‑dependent receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,649,375	0	T	G	5.2%	517.1 / 13.2	153	E576A (GAG→GCG) ‡	PP_1446	TonB‑dependent receptor
Reads supporting (aligned to +/- strand): ref base T (74/70); new base G (6/2); total (80/73)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 157 total aligned reads displayed.

CCAAGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAG  >  NC_002947/1649226‑1649524

ccAAGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCt                                                                                                                                                       <  1:3401366/150‑1 (MQ=255)

 cAAGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTc                                                                                                                                                      <  1:3377729/150‑1 (MQ=255)

   aGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGc                                                                                                                                                    <  2:2962251/150‑1 (MQ=255)

   aGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGc                                                                                                                                                    >  2:232836/1‑150 (MQ=255)

    ggtggtGTAGCCATCGGTGGTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGGTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGcc                                                                                                                                                   >  2:2680521/1‑150 (MQ=255)

     gtggtgTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCa                                                                                                                                                  <  1:2765360/150‑1 (MQ=255)

           tAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCg                                                                                                                                            >  2:593065/1‑150 (MQ=255)

            aGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGa                                                                                                                                           <  1:110781/150‑1 (MQ=255)

            aGCCAACGGTGCTGGTGTCGGTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGGCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGACAGTGGCGAGCCACTgccac                                                                                                                                           <  2:4719550/150‑5 (MQ=255)

             gCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAg                                                                                                                                          <  1:3347115/150‑1 (MQ=255)

              ccATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAgg                                                                                                                                         >  2:2575859/1‑150 (MQ=255)

                aTCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGtt                                                                                                                                       >  1:1790406/1‑150 (MQ=255)

                aTCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGtt                                                                                                                                       >  2:2594836/1‑150 (MQ=255)

                  cGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCt                                                                                                                                     >  1:1542064/1‑150 (MQ=255)

                   ggTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCtt                                                                                                                                    <  1:983131/150‑1 (MQ=255)

                    gTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTg                                                                                                                                   <  2:2746612/150‑1 (MQ=255)

                      gCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGc                                                                                                                                 >  1:4279055/1‑150 (MQ=255)

                        tGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTc                                                                                                                               <  1:165152/150‑1 (MQ=255)

                        tGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTc                                                                                                                               >  2:3093223/1‑150 (MQ=255)

                           tttCGTTGGCGGGTTTGTGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCggg                                                                                                                            <  1:4113416/150‑1 (MQ=255)

                           tttCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCggg                                                                                                                            >  2:957240/1‑150 (MQ=255)

                           gttCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCggg                                                                                                                            <  2:3342297/149‑1 (MQ=255)

                               gTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                        <  1:4749639/150‑1 (MQ=255)

                                  ggCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTc                                                                                                                     <  1:4728765/150‑1 (MQ=255)

                                   gCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCTGGCGCAGCGGGGCGATGCGGGGCAGTGGGTCGCCACTGTCGTGGTTCTTGGCTGGGGttttgtga                                                                                                                    >  1:4341156/1‑143 (MQ=255)

                                   gCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCa                                                                                                                    <  1:1080693/150‑1 (MQ=255)

                                   gCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCa                                                                                                                    >  2:1846286/1‑150 (MQ=255)

                                      ggTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCg                                                                                                                 <  1:1509648/150‑1 (MQ=255)

                                       gTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGGATAGTCACCgg                                                                                                                >  2:2843045/1‑150 (MQ=255)

                                            cGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAAc                                                                                                           <  2:3195831/150‑1 (MQ=255)

                                             ggTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACt                                                                                                          >  1:1882609/1‑150 (MQ=255)

                                             ggTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACt                                                                                                          >  2:4501099/1‑150 (MQ=255)

                                                  tgcGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAgc                                                                                                     <  1:4474806/150‑1 (MQ=255)

                                                      aCGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGa                                                                                                 >  1:3596491/1‑150 (MQ=255)

                                                       cGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGaa                                                                                                <  2:1775318/150‑1 (MQ=255)

                                                            cGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGCCCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGc                                                                                           <  1:1187623/150‑1 (MQ=255)

                                                                 tgcACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTaa                                                                                      >  2:3533775/1‑150 (MQ=255)

                                                                 tgcACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTaa                                                                                      >  2:2143635/1‑150 (MQ=255)

                                                                 tgcACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTaa                                                                                      >  1:5062133/1‑150 (MQ=255)

                                                                  gcACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAAc                                                                                     >  2:4258624/1‑150 (MQ=255)

                                                                   cACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCg                                                                                         <  2:3499296/145‑1 (MQ=255)

                                                                    aCGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACgg                                                                                   >  1:2922392/1‑150 (MQ=255)

                                                                          aCCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCg                                                                             <  1:3225270/150‑1 (MQ=255)

                                                                           cccGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCgg                                                                            <  2:2501155/150‑1 (MQ=255)

                                                                             cgcgCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGcc                                                                          <  1:3135950/150‑1 (MQ=255)

                                                                              gcgcCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGGGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCa                                                                         >  1:4682416/1‑150 (MQ=255)

                                                                                    gCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGcc                                                                   >  1:809362/1‑150 (MQ=255)

                                                                                      cAGCGTTCCAGCTCCCAGACCAGGCCGCTGTTCACGCGCAGCGGGGCGATGCGCGGCAGTGGCACGACAAAGTCTAGCTTCTTGTCACTGGTATAGTCACCTGACAACTCAAGCCAGAAGCTCCCCTAACGCTTCTCGGCCAAAAGCCAg                                                                 <  2:45819/150‑1 (MQ=255)

                                                                                      cAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAg                                                                 >  2:4508705/1‑150 (MQ=255)

                                                                                      cAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAg                                                                 >  2:3948821/1‑150 (MQ=255)

                                                                                      cAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAg                                                                 >  1:1998204/1‑150 (MQ=255)

                                                                                       aGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGc                                                                >  1:82161/1‑150 (MQ=255)

                                                                                        gCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCg                                                               >  2:187590/1‑150 (MQ=255)

                                                                                          ggTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAg                                                                 <  2:797659/146‑1 (MQ=255)

                                                                                             ccAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCt                                                          <  2:809362/150‑1 (MQ=255)

                                                                                                   cccAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCt                                                    >  1:3494655/1‑150 (MQ=255)

                                                                                                     cAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCg                                                  >  2:101786/1‑150 (MQ=255)

                                                                                                      aGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGa                                                 >  2:4605030/1‑150 (MQ=255)

                                                                                                      aGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGa                                                 >  1:4935564/1‑150 (MQ=255)

                                                                                                      aGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGa                                                 >  1:1355993/1‑150 (MQ=255)

                                                                                                       gACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGAt                                                <  2:4971958/150‑1 (MQ=255)

                                                                                                        aCCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATg                                               <  1:4208955/150‑1 (MQ=255)

                                                                                                         ccAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGc                                              <  2:4016886/150‑1 (MQ=255)

                                                                                                          cAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGcc                                             <  2:1784112/150‑1 (MQ=255)

                                                                                                            ggCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGt                                           <  1:3645120/150‑1 (MQ=255)

                                                                                                            ggCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGt                                           <  1:3645204/150‑1 (MQ=255)

                                                                                                             gCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTa                                          >  1:764116/1‑150 (MQ=255)

                                                                                                             gCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTa                                          <  1:1435158/150‑1 (MQ=255)

                                                                                                               cGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGa                                        <  2:727951/150‑1 (MQ=255)

                                                                                                               cGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGa                                        >  1:1998539/1‑150 (MQ=255)

                                                                                                                gCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGaa                                       >  1:3376686/1‑150 (MQ=255)

                                                                                                                gCTGTTCAGGCGCAGCGGGGCGATGAGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGaa                                       >  1:3102039/1‑150 (MQ=255)

                                                                                                                  tGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAgc                                     >  1:4991952/1‑150 (MQ=255)

                                                                                                                    ttCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAgcgc                                   <  1:2222712/150‑1 (MQ=255)

                                                                                                                    ttCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAgcgc                                   <  2:4274402/150‑1 (MQ=255)

                                                                                                                       aGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGcc                                >  1:4578940/1‑150 (MQ=255)

                                                                                                                       aGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGcc                                >  2:3686980/1‑150 (MQ=255)

                                                                                                                                 gggCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCtt                      <  1:2845629/150‑1 (MQ=255)

                                                                                                                                  ggCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTg                     <  1:2097567/150‑1 (MQ=255)

                                                                                                                                  ggCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTg                     <  1:1316550/150‑1 (MQ=255)

                                                                                                                                       tGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACt                <  2:38383/150‑1 (MQ=255)

                                                                                                                                        gCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTg               <  2:2880273/150‑1 (MQ=255)

                                                                                                                                        gCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTg               <  2:2880268/150‑1 (MQ=255)

                                                                                                                                        gCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTg               <  1:593065/150‑1 (MQ=255)

                                                                                                                                         ctcgacAGTGGCGAGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGCTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTgg              >  2:2148471/6‑150 (MQ=255)

                                                                                                                                         ctcgacAGTGGCGAGCCACTGCCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTgg              >  2:2459607/6‑150 (MQ=255)

                                                                                                                                          tcgacAGTGGCGAGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGt             >  1:3796533/5‑150 (MQ=255)

                                                                                                                                           cgacAGTGGCGAGCCACTGCCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTa            >  1:4133232/4‑150 (MQ=255)

                                                                                                                                            gacAGTGGCGAGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTAt           >  1:4940165/3‑150 (MQ=255)

                                                                                                                                             caagaaccTCGACAGTGGCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTAtt          >  2:2944163/8‑150 (MQ=255)

                                                                                                                                               aGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTcc        <  2:1569974/150‑1 (MQ=255)

                                                                                                                                                gTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCg       <  1:961835/150‑1 (MQ=255)

                                                                                                                                                gTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCg       <  1:3359100/150‑1 (MQ=255)

                                                                                                                                                gTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCg       <  1:4360963/150‑1 (MQ=255)

                                                                                                                                                 tGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCgg      <  1:2078967/150‑1 (MQ=255)

                                                                                                                                                 tGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCgg      <  1:683761/150‑1 (MQ=255)

                                                                                                                                                   gCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGa    >  2:1937681/1‑150 (MQ=255)

                                                                                                                                                    cTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGaa   <  2:4854022/150‑1 (MQ=255)

                                                                                                                                                    cTCGACAGTGGCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGaa   <  1:2944163/150‑1 (MQ=255)

                                                                                                                                                     tCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAg  >  2:3932196/1‑150 (MQ=255)

CCAAGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAG  >  NC_002947/1649226‑1649524

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