BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,471,047	G→C	7.3%	noncoding (105/119 nt)	PP_mr39 ←	P15 ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,471,047	0	G	C	7.3%	1086.9 / 55.1	324	noncoding (105/119 nt)	PP_mr39	P15 ncRNA
Reads supporting (aligned to +/- strand): ref base G (165/135); new base C (9/15); total (174/150)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

CGATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAG  >  NC_002947/3470899‑3471195

  aTCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAg                                                                                                                                                   <  2:469808/150‑1 (MQ=255)

   tCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGc                                                                                                                                                  <  1:4430591/150‑1 (MQ=255)

      gAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGa                                                                                                                                               <  2:3595131/150‑1 (MQ=255)

       aGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGAc                                                                                                                                              >  2:213276/1‑150 (MQ=255)

        ggcggcAGGCTGGCCCGATCGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATACGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAACGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGAcc                                                                                                                                             <  2:7308956/150‑1 (MQ=255)

         gcggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAATTCCCCCGGTCGGGAAGCCGACCg                                                                                                                                            >  2:1651601/1‑150 (MQ=255)

           ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCggg                                                                                                                                          <  2:4848282/150‑1 (MQ=255)

           ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCggg                                                                                                                                          >  1:6943352/1‑150 (MQ=255)

           ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCggg                                                                                                                                          <  1:4822167/150‑1 (MQ=255)

           ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCggg                                                                                                                                          <  1:6324934/150‑1 (MQ=255)

           ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCggg                                                                                                                                          <  2:5855506/150‑1 (MQ=255)

                gCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGttt                                                                                                                                     <  2:1588841/150‑1 (MQ=255)

                 cTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTg                                                                                                                                    >  1:541690/1‑150 (MQ=255)

                 cTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTg                                                                                                                                    >  1:466088/1‑150 (MQ=255)

                 cTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCGGGGGtttt                                                                                                                                    <  2:33108/150‑2 (MQ=255)

                 cTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCGGGGGtttt                                                                                                                                    <  2:2498436/150‑2 (MQ=255)

                  tGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGa                                                                                                                                   <  1:6507264/150‑1 (MQ=255)

                    gCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGt                                                                                                                                 >  2:4394218/1‑150 (MQ=255)

                    gCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGt                                                                                                                                 >  2:1504812/1‑150 (MQ=255)

                      ccGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTAt                                                                                                                               >  1:1000353/1‑150 (MQ=255)

                            gggATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCAt                                                                                                                         >  2:283668/1‑150 (MQ=255)

                             ggATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATg                                                                                                                        >  1:1145108/1‑150 (MQ=255)

                             ggATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATg                                                                                                                        <  2:792285/150‑1 (MQ=255)

                              gATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATCCTCaaac                                                                                                                       <  1:4938183/150‑4 (MQ=255)

                                taggCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATCAAACTGTCCTCAAAGTGAATCGGTGTTGCTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTc                                                                                                                     >  2:3803524/3‑150 (MQ=255)

                                   gCTAAATCGCCAGGCGTAGGTGTTTGCGTTGTAATAAGTGGCTTATTTCATGAATGTGTCCTCATAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggc                                                                                                                  <  1:307181/150‑1 (MQ=255)

                                   gCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggc                                                                                                                  >  1:6662781/1‑150 (MQ=255)

                                     tAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggcgg                                                                                                                >  2:5291665/1‑150 (MQ=255)

                                     tAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggcgg                                                                                                                >  1:6129326/1‑150 (MQ=255)

                                       aaTCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGcc                                                                                                              >  2:2796329/1‑150 (MQ=255)

                                         tCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTCCGGCGTccac                                                                                                            <  1:3776706/150‑3 (MQ=255)

                                           gCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTc                                                                                                          >  1:2372333/1‑150 (MQ=255)

                                            ccAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCg                                                                                                         >  2:5585330/1‑150 (MQ=255)

                                                gCGTAGGTGTTTGAGTAGTAATAAGTGGCGAAGTTCATGAATCTGTCCTCAAAGTGAATCGTTGTTGGTCAAAGGCCTCAAAAACAAAACCACCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTACGGCGGCCCCCAGAAgg                                                                                                     >  2:6679620/1‑150 (MQ=255)

                                                 cGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGt                                                                                                    <  2:4069634/150‑1 (MQ=255)

                                                   tAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTgg                                                                                                  >  2:3530299/1‑150 (MQ=255)

                                                    aGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGa                                                                                                 <  1:6915780/150‑1 (MQ=255)

                                                        gtTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGcc                                                                                             >  2:6256750/1‑150 (MQ=255)

                                                             cGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtg                                                                                        <  1:3452335/150‑1 (MQ=255)

                                                               tagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtggg                                                                                      >  1:1431577/1‑150 (MQ=255)

                                                                   aataaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGcc                                                                                                                                                 >  2:2256035/1‑87 (MQ=255)

                                                                   aataaGGGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATc                                                                                  <  2:3896716/150‑1 (MQ=255)

                                                                    ataaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCg                                                                                 <  2:6943352/150‑1 (MQ=255)

                                                                    ataaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCg                                                                                 >  1:7206813/1‑150 (MQ=255)

                                                                        gTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgcgc                                                                             >  1:7349705/1‑150 (MQ=255)

                                                                              tAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTg                                                                       <  1:6845629/150‑1 (MQ=255)

                                                                              tAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTg                                                                       <  2:2386973/150‑1 (MQ=255)

                                                                               aGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTgg                                                                      <  1:2218720/150‑1 (MQ=255)

                                                                                   cATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaa                                                                  >  2:4626432/1‑150 (MQ=255)

                                                                                   cATGAATCTGTCCTCAAAGTGAATCGGTGGTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTCTTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaa                                                                  >  2:5208151/1‑150 (MQ=255)

                                                                                    aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaaa                                                                 <  2:5735022/150‑1 (MQ=255)

                                                                                         tCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAAt                                                            <  2:3082398/150‑1 (MQ=255)

                                                                                             tCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGt                                                                                                                                 >  1:4054872/1‑77 (MQ=255)

                                                                                               cTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaa                                                      >  1:2265740/1‑150 (MQ=255)

                                                                                                tCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaaa                                                     >  2:2628177/1‑150 (MQ=255)

                                                                                                tCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaaa                                                     >  2:5806304/1‑150 (MQ=255)

                                                                                                tCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaaa                                                     >  1:805290/1‑150 (MQ=255)

                                                                                                tCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaaa                                                     >  1:4243130/1‑150 (MQ=255)

                                                                                                tCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaaa                                                     >  1:6358623/1‑150 (MQ=255)

                                                                                                tCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaaa                                                     >  2:1163164/1‑150 (MQ=255)

                                                                                                     gTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAAt                                                >  1:3215275/1‑150 (MQ=255)

                                                                                                      tGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATa                                               <  2:96331/150‑1 (MQ=255)

                                                                                                         aTCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCt                                            >  1:4600453/1‑150 (MQ=255)

                                                                                                          tCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTg                                           <  2:1564026/150‑1 (MQ=255)

                                                                                                           cGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTgg                                          <  1:3755250/150‑1 (MQ=255)

                                                                                                            ggTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCAt                                                                                                                         <  2:4971407/70‑1 (MQ=255)

                                                                                                             gtgtTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGa                                                                                                        <  2:6706164/86‑1 (MQ=255)

                                                                                                             gtgtTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCa                                        >  2:7215178/1‑150 (MQ=255)

                                                                                                             gtgtTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCa                                        >  1:2467159/1‑150 (MQ=255)

                                                                                                             gtgtTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCa                                        >  1:4197442/1‑150 (MQ=255)

                                                                                                                ttGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGt                                     >  1:6738302/1‑150 (MQ=255)

                                                                                                                ttGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGt                                     >  1:3501825/1‑150 (MQ=255)

                                                                                                                ttGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGt                                     >  1:4162323/1‑150 (MQ=255)

                                                                                                                    tCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTgg                                                                      >  2:1372865/1‑113 (MQ=255)

                                                                                                                     cAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCAc                                <  1:1966998/150‑1 (MQ=255)

                                                                                                                      aaaGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATa                                               >  1:1563721/1‑134 (MQ=255)

                                                                                                                      aaaGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACg                               >  2:6367797/1‑150 (MQ=255)

                                                                                                                      aaaGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACg                               >  2:5536165/1‑150 (MQ=255)

                                                                                                                      aaaGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTAGTCTCAGGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACg                               >  2:6096626/1‑150 (MQ=255)

                                                                                                                           ccTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGcc                          <  1:6000674/150‑1 (MQ=255)

                                                                                                                            cTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCg                         >  1:5160228/1‑150 (MQ=255)

                                                                                                                             tGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCgg                        >  1:5748672/1‑150 (MQ=255)

                                                                                                                              gAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGc                       >  1:2293318/1‑150 (MQ=255)

                                                                                                                                 aaacaaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGACACGCAGCCGGCAAc                    >  2:1044706/1‑150 (MQ=255)

                                                                                                                                  aacaaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACg                   >  2:184970/1‑150 (MQ=255)

                                                                                                                                    caaatcCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGc                 >  1:4859190/1‑150 (MQ=255)

                                                                                                                                     aaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGcc                >  2:1748713/1‑150 (MQ=255)

                                                                                                                                      aaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCa               <  1:5789501/150‑1 (MQ=255)

                                                                                                                                        acCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGc             >  1:3897242/1‑150 (MQ=255)

                                                                                                                                         cccccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGca            <  2:5969098/150‑1 (MQ=255)

                                                                                                                                          ccccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATa                                                           <  1:555669/102‑1 (MQ=255)

                                                                                                                                          ccccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcac           >  1:5825386/1‑150 (MQ=255)

                                                                                                                                          ccccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcac           >  1:5296093/1‑150 (MQ=255)

                                                                                                                                          ccccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcac           >  1:2583874/1‑150 (MQ=255)

                                                                                                                                          ccccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcac           >  1:6258108/1‑150 (MQ=255)

                                                                                                                                          ccccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcac           >  2:2112839/1‑150 (MQ=255)

                                                                                                                                            ccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAg         <  1:2799638/150‑1 (MQ=255)

                                                                                                                                             cGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTgg                                                                      <  2:2945594/88‑1 (MQ=255)

                                                                                                                                             cGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAgt        <  1:5959324/150‑2 (MQ=255)

                                                                                                                                             cGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAgg        <  2:4285650/150‑1 (MQ=255)

CGATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAG  >  NC_002947/3470899‑3471195

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