BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,754,850	C→A	9.7%	N280K (AAC→AAA) ‡	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,754,850	0	C	A	9.7%	396.2 / 25.2	124	N280K (AAC→AAA) ‡	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base C (55/57); new base A (6/6); total (61/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.98e-01

Only 100 of 130 total aligned reads displayed.

AAGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCG  >  NC_002947/5754701‑5754997

 agGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAaca                                                                                                                                                    <  1:3861447/150‑1 (MQ=255)

 agGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAaca                                                                                                                                                    <  1:3861344/150‑1 (MQ=255)

 agGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAaca                                                                                                                                                    <  2:1368677/150‑1 (MQ=255)

   gCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACAcc                                                                                                                                                  >  2:223332/1‑150 (MQ=255)

      gACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTtc                                                                                                                                               >  1:2385661/1‑150 (MQ=255)

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

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

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

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

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

                  gcgcTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGc                                                                                                                                   >  2:1615779/1‑150 (MQ=255)

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

                    gcTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTa                                                                                                                                 >  2:2845382/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)

                                cgtGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTc                                                                                                                     <  1:3198246/150‑1 (MQ=255)

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

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

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

                                      caACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGacac                                                                                                               >  2:3194737/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:2328743/1‑150 (MQ=255)

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

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

                                                     ttACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGACCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCAt                                                                                                >  2:3810816/1‑150 (MQ=255)

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

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

                                                          gCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcggcg                                                                                           <  2:1979172/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)

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

                                                                    cACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCaa                                                                                 <  2:2867769/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)

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

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

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

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

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

                                                                           aTGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGt                                                                          >  2:4331578/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:520476/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                                 gCACGTTCACTTCTCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGAGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTc                                    <  1:2886203/150‑1 (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:1281151/143‑1 (MQ=255)

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

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

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

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

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

                                                                                                                           aTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGt                          >  1:2055028/1‑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                        >  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)

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

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

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

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

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

                                                                                                                                       aaaGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGACCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGc              >  1:1727053/1‑150 (MQ=255)

                                                                                                                                        aaGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCt             >  2:1101259/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)

                                                                                                                                             cgcaatcacacCTTCTCGGGTGTAGGCTATTCCGGCCTTTCCGACACCGGCCTTTTATTAATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa        <  2:199086/143‑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)

AAGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCG  >  NC_002947/5754701‑5754997

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