BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,754,857	T→G	9.7%	S283A (TCG→GCG) ‡	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,754,857	0	T	G	9.7%	377.9 / 24.2	123	S283A (TCG→GCG) ‡	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base T (54/57); new base G (6/6); total (60/63)
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.93e-01

Only 100 of 128 total aligned reads displayed.

CGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTG  >  NC_002947/5754709‑5755006

cGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTtctc                                                                                                                                                      <  2:2092762/150‑1 (MQ=255)

  aaGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCgg                                                                                                                                                    <  1:1862490/150‑1 (MQ=255)

    gTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGt                                                                                                                                                  <  1:3428653/150‑1 (MQ=255)

    gTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGt                                                                                                                                                  <  1:121691/150‑1 (MQ=255)

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

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

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

         ggCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAgg                                                                                                                                             <  2:4102836/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:2701070/150‑1 (MQ=255)

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

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

                           gCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGt                                                                                                                                                  >  2:142070/1‑127 (MQ=255)

                           gCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGt                                                                                                                                                  <  1:142070/127‑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                                                                                                               <  1:1229345/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                                                                                                         >  1:2610940/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)

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

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

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

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

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

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

                                                                 tCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACg                                                                                     <  1:4538881/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                      ccGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAg                                                                                <  1:3901893/150‑1 (MQ=255)

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

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

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

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

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

                                                                                  tACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGc                                                                    <  2:520396/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)

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

                                                                                                           aCGTTCACTTCTCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTcct                                           >  2:4225376/1‑148 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                                   ttctCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGt                                   >  1:748965/4‑150 (MQ=255)

                                                                                                                    tGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTc                                  >  1:1887770/1‑150 (MQ=255)

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

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

                                                                                                                     gTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCt                                 >  2:1140819/1‑150 (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)

                                                                                                                            tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                          >  2:3177198/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:902172/150‑1 (MQ=255)

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

                                                                                                                                      gCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAg                >  2:2472470/1‑150 (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)

CGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTG  >  NC_002947/5754709‑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: `‑`