BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,471,050	G→C	7.5%	noncoding (102/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,050	0	G	C	7.5%	1057.9 / 56.2	318	noncoding (102/119 nt)	PP_mr39	P15 ncRNA
Reads supporting (aligned to +/- strand): ref base G (164/130); new base C (9/15); total (173/145)
Fisher's exact test for biased strand distribution p-value = 9.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 323 total aligned reads displayed.

ATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTG  >  NC_002947/3470901‑3471199

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

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

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

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

    gAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTCATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGa                                                                                                                                                   >  2:2218720/1‑150 (MQ=255)

     aGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGAc                                                                                                                                                  <  1:5778136/150‑1 (MQ=255)

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

        cggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCgg                                                                                                                                               <  2:5929646/150‑1 (MQ=255)

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

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

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

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

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

           cAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTTATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCggggg                                                                                                                                            <  1:464233/150‑1 (MQ=255)

           cAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCggggg                                                                                                                                            <  2:3675731/150‑1 (MQ=255)

               cTGGCCCGACCGGGAGTGGCTACAGCGCCAGGCGTAGGTGTTTGCGTAGTACGACGTGGCGTAGTTCCTGAATCGGTCCTCACCGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCGGGGGtttt                                                                                                                                        <  2:5481910/150‑2 (MQ=255)

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

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

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

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

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

                tGGCCCGAACGGGATTGGCTAAATCGCAAGGCGTTGGTGTTTGCGTAGTAATTAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAATGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGa                                                                                                                                       <  2:3553663/150‑1 (MQ=255)

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

                    ccGAACGGGATTGGCTAAATCGCCAGGCGGAGGTGTTTGCGTAGTCATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAAGCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTAt                                                                                                                                   >  2:1747513/1‑150 (MQ=255)

                       aaCGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGATTGAGTATtct                                                                                                                                <  2:1094818/150‑1 (MQ=255)

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

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

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

                            gATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCGGGGGttttg                                                                                                                                       >  2:3498025/1‑136 (MQ=39)

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

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

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

                                    aaaTCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCggcggc                                                                                                                   >  2:2541930/1‑150 (MQ=255)

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

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

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

                                           cAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGa                                                                                                            >  1:2588702/1‑150 (MQ=255)

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

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

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

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

                                                          gCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggt                                                                                             >  2:5044704/1‑150 (MQ=255)

                                                          gCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggt                                                                                             <  2:2925289/150‑1 (MQ=255)

                                                            gtagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtgg                                                                                           >  1:1921390/1‑150 (MQ=255)

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

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

                                                                  ataaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGAc                                                                                                    <  1:444259/135‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                               ttCATGAATCTGTCATCAAAGTGAATCGGTGTTGGTAAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtgg                                                                                           <  1:2494265/131‑1 (MQ=255)

                                                                                  attaaTCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaaa                                                                     <  1:5212658/147‑1 (MQ=255)

                                                                                   tGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAAc                                                                    >  2:7295214/1‑150 (MQ=255)

                                                                                    gAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAAcc                                                                   >  1:7436736/1‑150 (MQ=255)

                                                                                      aTCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCaa                                                                 <  2:4445344/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                   gtgtATCGGTGTTGGTCAAAGGCCTGAAAAACAAAAGCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAAt                                                    <  2:4690461/150‑1 (MQ=255)

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

                                                                                                      aaTCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCGATACCCAAAGTAATAGc                                                 <  2:5794021/150‑1 (MQ=255)

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

                                                                                                           gtgtTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCa                                            <  2:1634591/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                                     aaGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGc                                  <  1:7234458/150‑1 (MQ=255)

                                                                                                                      agtccTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGTAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCATAGTAATAGCTGGCAGGTGCCACGCa                                 <  1:340448/147‑1 (MQ=255)

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

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

                                                                                                                            gAAAAACAAATCCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCACAGTAATAGCTGGCAGGTGCCACGCAGCCGGc                           >  2:1575296/1‑150 (MQ=255)

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

                                                                                                                                  tcaaacCCCCGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgcgc                                                                                 >  1:3657217/3‑90 (MQ=37)

                                                                                                                                  caaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGc                     >  2:5731590/1‑150 (MQ=255)

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

                                                                                                                                   aaaacCCCCGGTCGGGAAGCCGACCGGGGGTGTGAGTATACTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCAAAAGTAATAGCTCGCAGGTGACACGAAGCCGGCAACGGcc                    >  2:6056705/1‑150 (MQ=255)

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

                                                                                                                                    aaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGAAGGTGCCACGCAGCCGGCAACGGCCa                   >  1:422238/1‑150 (MQ=255)

                                                                                                                                       cccccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCACCGGCCAGca                >  2:1397761/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                         cccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcaca              >  2:5229271/1‑150 (MQ=255)

                                                                                                                                         cccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcaca              >  2:4941978/1‑150 (MQ=255)

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

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

                                                                                                                                           cGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAgg            >  2:2964182/1‑150 (MQ=255)

                                                                                                                                                 ggAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAg      >  1:3474134/1‑150 (MQ=255)

                                                                                                                                                     gCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTg  <  1:5291665/150‑1 (MQ=255)

ATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGGCTG  >  NC_002947/3470901‑3471199

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`