BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,754,841	C→A	5.2%	D277E (GAC→GAA)	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,754,841	0	C	A	5.2%	911.1 / 27.2	268	D277E (GAC→GAA)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base C (102/152); new base A (8/6); total (110/158)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 269 total aligned reads displayed.

TGGTGAAGAAGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGA  >  NC_002947/5754693‑5754990

tGGTGAAGAAGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACg                                                                                                                                                      <  2:3787966/150‑1 (MQ=255)

  gTGAAGAAGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGc                                                                                                                                                    <  1:8757968/150‑1 (MQ=255)

  gTGAAGAAGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGc                                                                                                                                                    <  1:4553573/150‑1 (MQ=255)

  gTGAAGAAGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGc                                                                                                                                                    <  2:5893685/150‑1 (MQ=255)

    gaagaaGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCCTGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCaa                                                                                                                                                  <  1:781330/150‑1 (MQ=255)

    gaagaaGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCaa                                                                                                                                                  <  1:6216638/150‑1 (MQ=255)

    gaagaaGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCaa                                                                                                                                                  >  1:2849829/1‑150 (MQ=255)

     aagaagGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAg                                                                                                                                                 >  2:4053472/1‑150 (MQ=255)

       gaagGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCCTGCACGTGCACATGTCGATCTGGAAAGACGGCAAGaa                                                                                                                                               >  2:4953638/1‑150 (MQ=255)

       gaagGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGaa                                                                                                                                               <  1:3726779/150‑1 (MQ=255)

       gaagGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGaa                                                                                                                                               <  2:5943466/150‑1 (MQ=255)

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

          ggCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAacac                                                                                                                                            >  1:3673014/1‑150 (MQ=255)

           gCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACAcc                                                                                                                                           <  1:5388932/150‑1 (MQ=255)

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

             tGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCtt                                                                                                                                         >  1:6004696/1‑150 (MQ=255)

                cGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTtctc                                                                                                                                      >  1:3240893/1‑150 (MQ=255)

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

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

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

                  aaGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCgg                                                                                                                                    >  1:2658258/1‑150 (MQ=255)

                      aCAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGa                                                                                                                                <  1:2504621/150‑1 (MQ=255)

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

                        aGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAg                                                                                                                              <  1:4268173/150‑1 (MQ=255)

                        aGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAg                                                                                                                              <  2:3240893/150‑1 (MQ=255)

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

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

                              tGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAg                                                                                                                              <  2:5620340/144‑1 (MQ=255)

                              tGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATg                                                                                                                        <  1:2683266/150‑1 (MQ=255)

                                 aaTACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCg                                                                                                                     <  2:1992771/150‑1 (MQ=255)

                                 aaTACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCg                                                                                                                     <  1:8002711/150‑1 (MQ=255)

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

                                    aCGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGcc                                                                                                                  <  2:2658258/150‑1 (MQ=255)

                                     cgtcgtGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCt                                                                                                                 <  1:1086412/150‑1 (MQ=255)

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

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

                                         gtGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTcc                                                                                                             >  1:7309628/1‑150 (MQ=255)

                                                 cGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGc                                                                                                     >  2:7796225/1‑150 (MQ=255)

                                                   ttGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGccc                                                                                                   >  1:475105/1‑150 (MQ=255)

                                                    tGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCt                                                                                                  <  2:8207076/150‑1 (MQ=255)

                                                      ccGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGt                                                                                                <  2:8070114/150‑1 (MQ=255)

                                                      ccGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGt                                                                                                >  2:8270193/1‑150 (MQ=255)

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

                                                        gACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTAc                                                                                              <  2:8885988/150‑1 (MQ=255)

                                                        gACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTAc                                                                                              >  1:6673800/1‑150 (MQ=255)

                                                        gACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTAc                                                                                              >  2:5105179/1‑150 (MQ=255)

                                                         aCGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACt                                                                                             <  1:3031655/150‑1 (MQ=255)

                                                            cTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCa                                                                                          <  2:4477915/150‑1 (MQ=255)

                                                               aCGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcg                                                                                       >  1:4112765/1‑150 (MQ=255)

                                                                 ggCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcggc                                                                                     >  2:8075651/1‑150 (MQ=255)

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

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

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

                                                                  gCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATcggcg                                                                                    >  1:4400664/1‑150 (MQ=255)

                                                                       aCCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTatc                                                                               <  2:1452978/150‑1 (MQ=255)

                                                                       aCCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTatc                                                                               <  2:7381151/150‑1 (MQ=255)

                                                                           ccACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTAtcatca                                                                           <  1:8794644/150‑1 (MQ=255)

                                                                           ccACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTAtcatca                                                                           <  1:2980471/150‑1 (MQ=255)

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

                                                                             aCCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAg                                                                         >  1:5530853/1‑150 (MQ=255)

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

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

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

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

                                                                                          aGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGccc                                                            <  2:6633019/150‑1 (MQ=255)

                                                                                             cACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTCAGGCCCTGa                                                         >  2:4341728/1‑150 (MQ=255)

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

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

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

                                                                                                   aCGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCt                                                   >  1:9141147/1‑150 (MQ=255)

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

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

                                                                                                      gCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCa                                                <  2:3412995/150‑1 (MQ=255)

                                                                                                       cGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCAc                                               <  2:7090937/150‑1 (MQ=255)

                                                                                                             cGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAAccc                                         <  1:9163473/150‑1 (MQ=255)

                                                                                                             cGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAAccc                                         >  2:4515886/1‑150 (MQ=255)

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

                                                                                                                 tCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCg                                     >  2:8121459/1‑150 (MQ=255)

                                                                                                                   gggCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGAc                                   <  1:5918865/150‑1 (MQ=255)

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

                                                                                                                      taTGCACGTTCACTTCTCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCaa                                <  2:6782531/149‑1 (MQ=255)

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

                                                                                                                         gCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTc                             >  2:6135711/1‑150 (MQ=255)

                                                                                                                          cACGTTCACTTCTCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTcc                            >  1:2240410/1‑149 (MQ=255)

                                                                                                                           aCGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGt                           >  2:485209/1‑150 (MQ=255)

                                                                                                                               gCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTCTCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACac                       >  2:2221201/1‑149 (MQ=255)

                                                                                                                                 atttctCGATCTGGAAAGAAGGCAAGCCAACCTTCGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTTTACTTCATCGGTGGCGTCATCAAGCACGCCAAGGCCGTGAATGCCTTCACCAACCCGTCGACCAACTCTTACAAGc                     >  1:4428474/6‑150 (MQ=255)

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

                                                                                                                                   ttctCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGt                   >  2:4643576/4‑150 (MQ=255)

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

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

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

                                                                                                                                     ggcGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCt                 <  1:8189541/148‑1 (MQ=255)

                                                                                                                                     gTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCt                 >  1:5631167/1‑150 (MQ=255)

                                                                                                                                      tCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTg                <  1:8024644/150‑1 (MQ=255)

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

                                                                                                                                              gAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAgg        >  2:7942569/1‑150 (MQ=255)

                                                                                                                                              gAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAgg        >  2:179831/1‑150 (MQ=255)

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

                                                                                                                                                    cGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCATGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGATTCACCAACCCTTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGa  <  2:6313624/150‑1 (MQ=255)

TGGTGAAGAAGGCTGACGAAGTACAGGCGCTGAAATACGTCGTGCACAACGTTGCCGACGCTTACGGCCGTACCGCCACCTTCATGCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGA  >  NC_002947/5754693‑5754990

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