BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,754,863	G→A	9.8%	E285K (GAA→AAA) ‡	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,754,863	0	G	A	9.8%	345.2 / 21.8	122	E285K (GAA→AAA) ‡	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base G (55/55); new base A (6/6); total (61/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 122 total aligned reads displayed.

CAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTG  >  NC_002947/5754716‑5755006

cAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCTGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGaa                                                                                                                                               <  1:5101607/150‑1 (MQ=255)

cAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGaa                                                                                                                                               <  1:3741233/150‑1 (MQ=255)

  ggCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAgg                                                                                                                                             <  2:4102836/150‑1 (MQ=255)

     gcTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTa                                                                                                                                          <  2:4309598/150‑1 (MQ=255)

     gcTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTa                                                                                                                                          >  2:2845382/1‑150 (MQ=255)

     gcTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTa                                                                                                                                          >  2:1130506/1‑150 (MQ=255)

        gAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGc                                                                                                                                       >  2:137359/1‑150 (MQ=255)

         aaaTACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCCTGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGcc                                                                                                                                      >  2:3692842/1‑150 (MQ=255)

         aaaTACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGcc                                                                                                                                      >  1:2922177/1‑150 (MQ=255)

           aTACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCgg                                                                                                                                    <  1:2201109/150‑1 (MQ=255)

            tACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGc                                                                                                                                   <  2:2701108/150‑1 (MQ=255)

            tACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGc                                                                                                                                   <  2:2701070/150‑1 (MQ=255)

               gtcgtGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTg                                                                                                                                >  2:627922/1‑150 (MQ=255)

                  gtGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTcc                                                                                                                             <  1:3333180/150‑1 (MQ=255)

                    gCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGa                                                                                                                           <  1:4482660/150‑1 (MQ=255)

                      acaACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGaca                                                                                                                         >  1:2922998/1‑150 (MQ=255)

                       caACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGacac                                                                                                                        >  2:3194737/1‑150 (MQ=255)

                        aaCGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACAcc                                                                                                                       >  1:1331829/1‑150 (MQ=255)

                                cGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTa                                                                                                               <  1:184687/150‑1 (MQ=255)

                                cGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTa                                                                                                               <  2:2921260/150‑1 (MQ=255)

                                    gCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTc                                                                                                           >  2:1738925/1‑150 (MQ=255)

                                      ttACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCAt                                                                                                         >  2:4958914/1‑150 (MQ=255)

                                      ttACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCAt                                                                                                         >  1:2610682/1‑150 (MQ=255)

                                      ttACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCAt                                                                                                         >  2:2328743/1‑150 (MQ=255)

                                        aCGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcg                                                                                                       <  2:2443958/150‑1 (MQ=255)

                                           gCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcggcg                                                                                                    <  1:2845782/150‑1 (MQ=255)

                                           gCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcggcg                                                                                                    <  2:1979172/150‑1 (MQ=255)

                                           gCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcggcg                                                                                                    <  1:1816231/150‑1 (MQ=255)

                                           gCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcggcg                                                                                                    <  1:2360895/150‑1 (MQ=255)

                                           gCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcggcg                                                                                                    <  1:734427/150‑1 (MQ=255)

                                            ccGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcggcgg                                                                                                   >  1:1220456/1‑150 (MQ=255)

                                            ccGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcggcgg                                                                                                   >  1:1220317/1‑150 (MQ=255)

                                             cGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGt                                                                                                  <  1:2982788/150‑1 (MQ=255)

                                                     cACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCaa                                                                                          <  2:1932058/150‑1 (MQ=255)

                                                        cTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCa                                                                                       >  1:1876748/1‑150 (MQ=255)

                                                         ttCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCCTGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCAc                                                                                      >  1:1098345/1‑150 (MQ=255)

                                                         ttCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCAc                                                                                      <  1:2800606/150‑1 (MQ=255)

                                                            aTGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGt                                                                                   >  1:1747332/1‑150 (MQ=255)

                                                            aTGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGt                                                                                   >  2:4331578/1‑150 (MQ=255)

                                                            aTGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGt                                                                                   <  1:1130506/150‑1 (MQ=255)

                                                            aTGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGt                                                                                   <  2:1345652/150‑1 (MQ=255)

                                                             tGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTa                                                                                  <  1:1056468/150‑1 (MQ=255)

                                                              gCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTaa                                                                                 <  1:5106901/150‑1 (MQ=255)

                                                                  aaGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGcc                                                                             <  2:3748672/150‑1 (MQ=255)

                                                                    gCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCt                                                                           >  2:1550705/1‑150 (MQ=255)

                                                                       aCTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGaa                                                                        >  1:140013/1‑150 (MQ=255)

                                                                       aCTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGaa                                                                        <  2:2370551/150‑1 (MQ=255)

                                                                       aCTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGaa                                                                        >  1:2570695/1‑150 (MQ=255)

                                                                        cTGTACGGCGACAACGGCTCGGGCATGCACGAGCACATGTCGATCTGGAAAGACGGCAAGAACACCTACTCGGGTGAAGGCTATGCCCGCCTGTCCGACACCCCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAAc                                                                       <  2:2862336/150‑1 (MQ=255)

                                                                           tACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGc                                                                    <  2:520447/150‑1 (MQ=255)

                                                                           tACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGc                                                                    <  1:964063/150‑1 (MQ=255)

                                                                           tACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGc                                                                    <  2:520476/150‑1 (MQ=255)

                                                                              ggCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTc                                                                 >  1:3212508/1‑150 (MQ=255)

                                                                                cGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAATACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCAc                                                               <  2:385564/150‑1 (MQ=255)

                                                                                         cTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTc                                                      <  1:3874127/150‑1 (MQ=255)

                                                                                          tCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCg                                                     <  2:3697175/150‑1 (MQ=255)

                                                                                              gCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCa                                                 >  1:2870125/1‑150 (MQ=255)

                                                                                                aTGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAAc                                               <  1:4569946/150‑1 (MQ=255)

                                                                                                  gCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTc                                             <  1:4604866/150‑1 (MQ=255)

                                                                                                     cGTTCACTTCTCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTa                                          <  1:4940439/150‑1 (MQ=255)

                                                                                                     cGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTa                                          <  2:2936445/150‑1 (MQ=255)

                                                                                                        tcacttctCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATaa                                       <  2:1281151/143‑1 (MQ=255)

                                                                                                        tcacttctCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATaa                                       <  2:177099/143‑1 (MQ=255)

                                                                                                          acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACaatc                                     >  2:416701/1‑148 (MQ=255)

                                                                                                          acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                     >  2:416811/1‑150 (MQ=255)

                                                                                                          acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                     >  2:448607/1‑150 (MQ=255)

                                                                                                          acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                     >  2:448582/1‑150 (MQ=255)

                                                                                                          acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                     >  2:696580/1‑150 (MQ=255)

                                                                                                          acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                     >  2:1256511/1‑150 (MQ=255)

                                                                                                          acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                     <  1:137359/150‑1 (MQ=255)

                                                                                                            aTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGt                                   >  1:2055028/1‑150 (MQ=255)

                                                                                                             tGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTc                                  <  1:2094032/150‑1 (MQ=255)

                                                                                                             tGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGATTCACCAACCCGTCGACCAACTCGTACAAGCGTc                                  <  1:1250506/150‑1 (MQ=255)

                                                                                                              ctCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCt                                 >  1:2028944/2‑150 (MQ=255)

                                                                                                              ctCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCt                                 >  2:4582686/2‑150 (MQ=255)

                                                                                                                 gATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGt                              >  1:4843381/1‑150 (MQ=255)

                                                                                                                 gATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGt                              >  1:482366/1‑150 (MQ=255)

                                                                                                                  aTCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTc                             >  2:4402304/1‑150 (MQ=255)

                                                                                                                  aTCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTc                             >  2:1049903/1‑150 (MQ=255)

                                                                                                                   tCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTcc                            >  1:2351797/1‑150 (MQ=255)

                                                                                                                    cTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTccc                           <  1:4160340/150‑1 (MQ=255)

                                                                                                                    cTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTccc                           <  1:4160247/150‑1 (MQ=255)

                                                                                                                     tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                          >  2:3177198/1‑150 (MQ=255)

                                                                                                                      ggAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAg                         >  1:2611590/1‑150 (MQ=255)

                                                                                                                         aaGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCt                      >  2:1101259/1‑150 (MQ=255)

                                                                                                                          aGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCtt                     >  1:4139359/1‑150 (MQ=255)

                                                                                                                              ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGTGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                 <  2:2306260/150‑1 (MQ=255)

                                                                                                                              ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGcccggca                 <  2:2385661/150‑7 (MQ=255)

                                                                                                                              ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                 <  1:1633489/150‑1 (MQ=255)

                                                                                                                              ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                 <  1:902172/150‑1 (MQ=255)

                                                                                                                                 aaGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGcc              >  1:1999939/1‑150 (MQ=255)

                                                                                                                                  aGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGccc             <  1:3076734/150‑1 (MQ=255)

                                                                                                                                   gAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGcccc            >  2:4290315/1‑150 (MQ=255)

                                                                                                                                   gAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGcccc            >  2:4290363/1‑150 (MQ=255)

                                                                                                                                    aaCACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCg           <  2:3949514/150‑1 (MQ=255)

                                                                                                                                       acCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTa        >  1:4928409/1‑150 (MQ=255)

                                                                                                                                           tctcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGc    >  2:3810017/1‑150 (MQ=255)

                                                                                                                                           tctcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGc    >  2:2964333/1‑150 (MQ=255)

                                                                                                                                             tcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg  <  1:2845382/150‑1 (MQ=255)

                                                                                                                                             tcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg  <  1:3692842/150‑1 (MQ=255)

CAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTG  >  NC_002947/5754716‑5755006

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