BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,764,033	G→C	5.6%	D41E (GAC→GAG)	PP_2412 ←	ParA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,764,033	0	G	C	5.6%	777.6 / 24.4	229	D41E (GAC→GAG)	PP_2412	ParA family protein
Reads supporting (aligned to +/- strand): ref base G (91/125); new base C (8/5); total (99/130)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

TAGCGGGTTTCGGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCGGGGTCAA  >  NC_002947/2763884‑2764182

tAGCGGGTTTCGGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGggtcgagg                                                                                                                                                       >  2:2047832/1‑150 (MQ=255)

  gCGGGTTTCGGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGggtcgaggtc                                                                                                                                                     >  2:3069523/1‑150 (MQ=255)

  gCGGGTTTCGGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGggtcgaggtc                                                                                                                                                     <  1:578109/150‑1 (MQ=255)

  gCGGGTTTCGGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGggtcgaggtc                                                                                                                                                     >  1:2956261/1‑150 (MQ=255)

  gCGGGTTTCGGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGggtcgaggtc                                                                                                                                                     >  2:1372705/1‑150 (MQ=255)

   cGGGTTTCGGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGggtcgaggtcg                                                                                                                                                    >  1:4655413/1‑150 (MQ=255)

    gggTTTCGGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGggtcgaggtcga                                                                                                                                                   <  2:1501772/150‑1 (MQ=255)

        ttCGGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAg                                                                                                                                               <  1:1372705/150‑1 (MQ=255)

           ggTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCaa                                                                                                                                            >  2:5290426/1‑150 (MQ=255)

           ggTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCaa                                                                                                                                            >  2:1873931/1‑150 (MQ=255)

           ggTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCaa                                                                                                                                            >  2:2720699/1‑150 (MQ=255)

            gTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAg                                                                                                                                           >  2:3604436/1‑150 (MQ=255)

               aTTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGtt                                                                                                                                        >  2:4678227/1‑150 (MQ=255)

                ttGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTc                                                                                                                                       <  2:1868395/150‑1 (MQ=255)

                 tGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCg                                                                                                                                      <  1:4946239/150‑1 (MQ=255)

                   ccACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGt                                                                                                                                    >  1:1795941/1‑150 (MQ=255)

                    cACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACTAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTa                                                                                                                                   <  2:3322626/150‑1 (MQ=255)

                    cACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTa                                                                                                                                   >  1:2177242/1‑150 (MQ=255)

                       gcgGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGcc                                                                                                                                >  2:2884249/1‑150 (MQ=255)

                       gcgGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGcc                                                                                                                                >  2:620412/1‑150 (MQ=255)

                         ggTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGCGGTCGACCAGCCAGGTGCGGTAGCCCt                                                                                                                              >  2:15510/1‑150 (MQ=255)

                         ggTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAAGGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCt                                                                                                                              >  2:3298339/1‑150 (MQ=255)

                              ttcttGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTcggcg                                                                                                                         <  2:461855/150‑1 (MQ=255)

                              ttcttGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTcggcg                                                                                                                         <  1:3475715/150‑1 (MQ=255)

                              ttcttGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTcggcg                                                                                                                         <  1:2293108/150‑1 (MQ=255)

                              ttcttGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCATcggcg                                                                                                                         >  2:682550/1‑150 (MQ=255)

                               tcttGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCTGGTCGACCTCGACCCGCag                                                                                                                                            <  1:884514/130‑2 (MQ=255)

                               tcttGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCTGGTCGACCTCGACCCGCag                                                                                                                                            >  2:3713625/1‑129 (MQ=255)

                                   gCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCt                                                                                                                    <  1:807773/150‑1 (MQ=255)

                                   gCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCt                                                                                                                    <  1:3635975/150‑1 (MQ=255)

                                    ccGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGt                                                                                                                                    <  1:3759729/133‑1 (MQ=255)

                                    ccGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGt                                                                                                                                    >  2:3759729/1‑133 (MQ=255)

                                    ccGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTg                                                                                                                   <  2:3945422/150‑1 (MQ=255)

                                     cggcggcggTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTgg                                                                                                                  >  2:5333223/1‑150 (MQ=255)

                                       gcggcggTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCa                                                                                                                >  1:4713305/1‑150 (MQ=255)

                                         ggcggTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAgc                                                                                                              <  1:5129873/150‑1 (MQ=255)

                                         ggcggTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAgc                                                                                                              <  1:4563832/150‑1 (MQ=255)

                                           cggTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAgcgg                                                                                                            <  1:5066899/150‑1 (MQ=255)

                                              tAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGa                                                                                                         <  1:2837934/150‑1 (MQ=255)

                                                      gACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAg                                                                                                 <  1:2313350/150‑1 (MQ=255)

                                                       aCAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAgg                                                                                                >  2:3695175/1‑150 (MQ=255)

                                                       aCAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAgg                                                                                                <  1:2911870/150‑1 (MQ=255)

                                                       aCAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAgg                                                                                                <  1:389842/150‑1 (MQ=255)

                                                          aGGTCTGGCGGAAGAAGTCGGCGCTGCCGGGCGGACTGGTGTCGTGGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGa                                                                                             <  2:4782656/150‑1 (MQ=255)

                                                            gTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATg                                                                                           >  1:1845516/1‑150 (MQ=255)

                                                             tCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGc                                                                                          >  2:5094806/1‑150 (MQ=255)

                                                             tCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGc                                                                                          >  1:5372108/1‑150 (MQ=255)

                                                              cTGGCGGAAGCAGTCGGCGATGCCGGGTGGAATGGTGTCGGTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCt                                                                                         <  2:2127090/150‑1 (MQ=255)

                                                               tGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTg                                                                                        <  1:473745/150‑1 (MQ=255)

                                                                    gaagaaGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGAttt                                                                                   <  2:2243932/150‑1 (MQ=255)

                                                                    gaagaaGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGAttt                                                                                   >  1:4167680/1‑150 (MQ=255)

                                                                      agaagTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGc                                                                                 <  2:4002830/150‑1 (MQ=255)

                                                                      agaagTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGc                                                                                 <  1:255339/150‑1 (MQ=255)

                                                                        aagTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCg                                                                               <  1:1054207/150‑1 (MQ=255)

                                                                             ggCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGcc                                                                          >  2:3864396/1‑150 (MQ=255)

                                                                             ggCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGcc                                                                          >  1:2598338/1‑150 (MQ=255)

                                                                              gCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCa                                                                         >  1:3024491/1‑150 (MQ=255)

                                                                                      gggTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTg                                                                 >  1:3909700/1‑150 (MQ=255)

                                                                                       ggTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGa                                                                <  2:4083039/150‑1 (MQ=255)

                                                                                        gTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGAt                                                               <  1:3101686/150‑1 (MQ=255)

                                                                                         tGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGAtt                                                              <  2:5568845/150‑1 (MQ=255)

                                                                                              tGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGaaaa                                                         <  2:2811026/150‑1 (MQ=255)

                                                                                              tGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGaaaa                                                         <  1:4451182/150‑1 (MQ=255)

                                                                                              tGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGaaaa                                                         <  1:4493795/150‑1 (MQ=255)

                                                                                               ggTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAAc                                                        <  1:5173162/150‑1 (MQ=255)

                                                                                                 tgtCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACca                                                      <  2:4413811/150‑1 (MQ=255)

                                                                                                   tCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACcaca                                                    >  1:5106030/1‑150 (MQ=255)

                                                                                                     gTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACacg                                                  >  2:3122615/1‑150 (MQ=255)

                                                                                                     gTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACacg                                                  >  1:3536835/1‑150 (MQ=255)

                                                                                                     gTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACAACacg                                                  >  1:5540005/1‑150 (MQ=255)

                                                                                                      ttGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACacga                                                 <  2:4214862/150‑1 (MQ=255)

                                                                                                       tGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACacgac                                                <  2:5056745/150‑1 (MQ=255)

                                                                                                         aCCCGCCCGGTGAGGAAGCAGGGGGCATTGGCCTGGGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAACCCACACGACGc                                              <  2:2956320/150‑1 (MQ=255)

                                                                                                           cAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCAt                                            <  2:4816705/150‑1 (MQ=255)

                                                                                                            aGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATg                                           <  2:1143781/150‑1 (MQ=255)

                                                                                                              cccGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATgag                                         >  1:1715100/1‑150 (MQ=255)

                                                                                                                       gTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATc                                >  1:2462974/1‑150 (MQ=255)

                                                                                                                       gTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATc                                <  2:3909605/150‑1 (MQ=255)

                                                                                                                          gAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCAt                             <  1:1364164/150‑1 (MQ=255)

                                                                                                                          gAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCAt                             >  1:4812820/1‑150 (MQ=255)

                                                                                                                           aaGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGaaaa                                                         >  2:1743522/1‑121 (MQ=255)

                                                                                                                           aaGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATg                            >  2:4468531/1‑150 (MQ=255)

                                                                                                                             ggTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGaa                          <  1:1835357/150‑1 (MQ=255)

                                                                                                                              gTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGAGTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAg                         >  2:4886205/1‑150 (MQ=255)

                                                                                                                               tGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGt                        >  1:3233628/1‑150 (MQ=255)

                                                                                                                                ggCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTa                       <  2:2324761/150‑1 (MQ=255)

                                                                                                                                 gCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTaa                      >  1:5003988/1‑150 (MQ=255)

                                                                                                                                   aTTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTaaaa                    >  1:4059212/1‑150 (MQ=255)

                                                                                                                                   aTTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTaaaa                    >  1:4451342/1‑150 (MQ=255)

                                                                                                                                    ttGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTaaaaa                   <  2:419912/150‑1 (MQ=255)

                                                                                                                                    gaaccatGCTGGTCGACCACGACCCGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGACGTAGCAGGCGATGCTGGATATGCCGACCCCACCTTACAGATTGAAAACCACACGACGCATGAGAACCTCATCCATGAAGAataaa                   >  2:1200746/7‑150 (MQ=255)

                                                                                                                                       ccttGCTGGTCGACCTCGACCCGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGt                >  1:3225090/4‑150 (MQ=38)

                                                                                                                                       ccttGCTGGTCGACCTCGACCCGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGt                >  1:4888336/4‑150 (MQ=38)

                                                                                                                                        ccTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTc               <  1:2720699/150‑1 (MQ=255)

                                                                                                                                         ttGCTGGTCGACCTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATTTct              >  1:4263216/2‑149 (MQ=255)

                                                                                                                                             gggTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTc          >  1:3500829/1‑150 (MQ=255)

                                                                                                                                               gtcgaggtcgaCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCgg        >  2:3224100/1‑150 (MQ=255)

                                                                                                                                                tcgaggtcgaCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCggg       <  1:4092636/150‑1 (MQ=255)

                                                                                                                                                tcgaggtcgaCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCggg       <  1:1166187/150‑1 (MQ=255)

                                                                                                                                                 cgaggtcgaCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCgggg      <  1:2030715/150‑1 (MQ=255)

TAGCGGGTTTCGGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCGGGGTCAA  >  NC_002947/2763884‑2764182

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