BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,629,798	G→C	5.2%	G34A (GGC→GCC)	PP_3198 →	2Fe‑2S ferredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,629,798	0	G	C	5.2%	986.5 / 30.4	290	G34A (GGC→GCC)	PP_3198	2Fe‑2S ferredoxin
Reads supporting (aligned to +/- strand): ref base G (134/141); new base C (6/9); total (140/150)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 306 total aligned reads displayed.

ATGGGCACCCAGCGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCAC  >  NC_002947/3629649‑3629946

 tGGGCACCCAGCGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGc                                                                                                                                                     >  2:7903833/1‑150 (MQ=255)

    gCACCCAGCGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAg                                                                                                                                                  >  2:8503092/1‑150 (MQ=255)

      aCCCAGCGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGt                                                                                                                                                <  2:426210/150‑1 (MQ=255)

        ccAGCGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTAc                                                                                                                                              <  2:1331557/150‑1 (MQ=255)

          aGCGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTAcgc                                                                                                                                            <  1:6772415/150‑1 (MQ=255)

          aGCGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCCGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTAcga                                                                                                                                            <  2:7266803/150‑2 (MQ=255)

           gCGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTAcgaa                                                                                                                                           <  2:3681187/150‑3 (MQ=255)

            cGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGCCCAGGTAccgaa                                                                                                                                          <  2:3754028/150‑5 (MQ=255)

            cGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGCCCAGGTAccgaa                                                                                                                                          <  1:7714660/150‑5 (MQ=255)

             ggCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTc                                                                                                                                         <  1:2559544/150‑1 (MQ=255)

             ggCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTc                                                                                                                                         <  1:2046032/150‑1 (MQ=255)

             ggCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTc                                                                                                                                         >  2:8296099/1‑150 (MQ=255)

              gCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGCCCAGGTAccgaagg                                                                                                                                        <  1:8337132/150‑7 (MQ=255)

                ccACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTAc                                                                                                                                      <  2:4827189/150‑1 (MQ=255)

                    cTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCa                                                                                                                                  >  1:3828331/1‑150 (MQ=255)

                    cTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCa                                                                                                                                  >  1:4497427/1‑150 (MQ=255)

                      gCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAAc                                                                                                                                >  2:7514774/1‑150 (MQ=255)

                      gCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAAc                                                                                                                                <  2:186745/150‑1 (MQ=255)

                      gCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAAc                                                                                                                                <  2:1961613/150‑1 (MQ=255)

                      gCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAAc                                                                                                                                >  2:2241695/1‑150 (MQ=255)

                            cAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGc                                                                                                                          <  2:2646724/150‑1 (MQ=255)

                              gggTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGccc                                                                                                                        <  2:8755351/150‑1 (MQ=255)

                               ggTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCg                                                                                                                       >  1:2678738/1‑150 (MQ=255)

                                gTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGc                                                                                                                      >  1:4755237/1‑150 (MQ=255)

                                    aaaaaTGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCAttt                                                                                                                  >  2:8712445/1‑150 (MQ=255)

                                    aaaaaTGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCAttt                                                                                                                  >  2:1826999/1‑150 (MQ=255)

                                      aaaTTGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGa                                                                                                                <  2:4856588/150‑1 (MQ=255)

                                      aaaTGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGa                                                                                                                >  1:3156224/1‑150 (MQ=255)

                                      aaaTGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGa                                                                                                                <  1:7018564/150‑1 (MQ=255)

                                        aTGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATg                                                                                                              <  1:1797860/150‑1 (MQ=255)

                                        aTGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATg                                                                                                              <  1:3076304/150‑1 (MQ=255)

                                        aTGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATg                                                                                                              >  1:8518737/1‑150 (MQ=255)

                                         tGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATgg                                                                                                             >  2:589918/1‑150 (MQ=255)

                                             gtgggtggAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGcc                                                                                                         >  1:5171329/1‑150 (MQ=255)

                                              tgggtggAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGccc                                                                                                        >  1:3569945/1‑150 (MQ=255)

                                                ggtggAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCt                                                                                                      >  2:3046135/1‑150 (MQ=255)

                                                ggtggAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCt                                                                                                      >  1:4344413/1‑150 (MQ=255)

                                                  tggAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTgg                                                                                                    <  1:5098891/150‑1 (MQ=255)

                                                  tggAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGATTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTgg                                                                                                    <  2:3235266/150‑1 (MQ=255)

                                                   ggAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGa                                                                                                   <  1:589918/150‑1 (MQ=255)

                                                           gCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCa                                                                                           <  1:8604786/150‑1 (MQ=255)

                                                            ccGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCaa                                                                                          >  1:669434/1‑150 (MQ=255)

                                                            ccGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCaa                                                                                          <  2:1611766/150‑1 (MQ=255)

                                                             cGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAg                                                                                         >  1:5806765/1‑150 (MQ=255)

                                                              ggCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGa                                                                                        >  2:7681865/1‑150 (MQ=255)

                                                               gCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGaa                                                                                       <  1:7083210/150‑1 (MQ=255)

                                                               gCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGaa                                                                                       >  2:6575394/1‑150 (MQ=255)

                                                                cTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAAc                                                                                      <  2:6426430/150‑1 (MQ=255)

                                                                   gATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGc                                                                                   >  2:4084608/1‑150 (MQ=255)

                                                                     tGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCtt                                                                                 <  2:7448483/150‑1 (MQ=255)

                                                                       aGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCtttt                                                                               <  1:3918265/150‑1 (MQ=255)

                                                                        gCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCttttt                                                                              <  2:6874238/150‑1 (MQ=255)

                                                                        gCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCAGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCttttt                                                                              <  2:2145470/150‑1 (MQ=255)

                                                                          tGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGa                                                                            <  2:4053788/150‑1 (MQ=255)

                                                                           gggCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAg                                                                           <  1:5632798/150‑1 (MQ=255)

                                                                            ggCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAgc                                                                          <  2:7058218/150‑1 (MQ=255)

                                                                              cGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAgcgc                                                                        >  1:7838429/1‑150 (MQ=255)

                                                                               gAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGcc                                                                       >  2:3135455/1‑150 (MQ=255)

                                                                                  ggTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGAt                                                                    >  2:5687855/1‑150 (MQ=255)

                                                                                  ggTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGAt                                                                    >  1:7545198/1‑150 (MQ=255)

                                                                                    tCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGAggg                                                                  >  2:7618937/1‑150 (MQ=255)

                                                                                    tCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATgg                                                                  <  2:5437989/150‑1 (MQ=255)

                                                                                    tCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATgg                                                                  >  2:3299784/1‑150 (MQ=255)

                                                                                     cAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGc                                                                 <  2:2294816/150‑1 (MQ=255)

                                                                                     cAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGc                                                                 <  1:797923/150‑1 (MQ=255)

                                                                                        gcgcgcGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAg                                                              >  1:6504519/1‑150 (MQ=255)

                                                                                        gcgcgcGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAg                                                              >  2:2950498/1‑150 (MQ=255)

                                                                                                tttCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTTCCGCAAGCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGt                                                      <  2:6211191/150‑1 (MQ=255)

                                                                                                 ttCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTc                                                     <  1:589685/150‑1 (MQ=255)

                                                                                                 ttCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTc                                                     <  2:1467882/150‑1 (MQ=255)

                                                                                                  tCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCt                                                    >  1:5383570/1‑150 (MQ=255)

                                                                                                    gACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGc                                                  >  2:6959663/1‑150 (MQ=255)

                                                                                                    gACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGc                                                  >  2:4637457/1‑150 (MQ=255)

                                                                                                       ccACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTAc                                               >  1:3765414/1‑150 (MQ=255)

                                                                                                        cacaGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACg                                              >  2:5001748/1‑150 (MQ=255)

                                                                                                         acaGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGc                                             <  2:4336215/150‑1 (MQ=255)

                                                                                                         acaGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGc                                             >  1:1191577/1‑150 (MQ=255)

                                                                                                             ggcagggGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCg                                                                                             >  2:4491258/1‑98 (MQ=255)

                                                                                                                    gcgcGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGgcgcgc                                  <  2:2376208/150‑1 (MQ=255)

                                                                                                                    gcgcGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGgcgcgc                                  >  1:7347908/1‑150 (MQ=255)

                                                                                                                      gcgGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACATGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGCTGGCCAGTGGATCGTCTg                                                   >  2:2603543/1‑131 (MQ=255)

                                                                                                                               cTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGccc                       <  1:2858327/150‑1 (MQ=255)

                                                                                                                                     gCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGAc                 >  2:3359361/1‑150 (MQ=255)

                                                                                                                                       ccTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACgg               >  1:975888/1‑150 (MQ=255)

                                                                                                                                         ttCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCa             >  1:8418035/1‑150 (MQ=255)

                                                                                                                                           gcgTACCTGGCCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACt           >  2:8623399/3‑150 (MQ=255)

                                                                                                                                           cGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACt           >  1:5288265/1‑150 (MQ=255)

                                                                                                                                            ggTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTGCGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACtt          <  2:8142213/150‑1 (MQ=255)

                                                                                                                                            ggTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACtt          <  1:6167735/150‑1 (MQ=255)

                                                                                                                                            ggTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACtt          <  2:7203323/150‑1 (MQ=255)

                                                                                                                                              tACCTGGGCCAGGCACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGa                  >  1:4085098/1‑140 (MQ=255)

                                                                                                                                               aCCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGCAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACATGCg       >  2:7868129/1‑150 (MQ=255)

                                                                                                                                                ccTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGt      >  1:8898917/1‑150 (MQ=255)

                                                                                                                                                ccTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGt      <  2:7737807/150‑1 (MQ=255)

                                                                                                                                                  tGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTc    >  2:2329699/1‑150 (MQ=255)

                                                                                                                                                   gggCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCa   >  1:1671773/1‑150 (MQ=255)

                                                                                                                                                   gggCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCa   <  1:140275/150‑1 (MQ=255)

                                                                                                                                                   gggCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCa   >  2:672380/1‑150 (MQ=255)

                                                                                                                                                   gggCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCa   <  1:6271376/150‑1 (MQ=255)

                                                                                                                                                    ggCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCAc  <  2:355881/150‑1 (MQ=255)

ATGGGCACCCAGCGGCCCACCTGCGAAGCAGGGTTGAAAAATGGCGTGGGTGGAGCTGTGCCGGCTCGATGAGCTGGGCGAGGGTCAAGCGCGCGGTTTCGACCCACAGGGCAGGGGCGCGGACAGCCTGTTCGCCCTTCGGTACCTGGGCCAGGTACGCGTCTACCGCAACCTGTGCCCGCATTTGATGGTGCCCCTGGAGTACCGCAAGAACCGCTTTTTGAGCGCCGATGGCCAGTGGATCGTCTGCTACGCCCACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCAC  >  NC_002947/3629649‑3629946

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