BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	939,928:1	+C	5.2%	coding (5697/18933 nt)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	939,928	1	.	C	5.2%	389.8 / 20.1	115	coding (5697/18933 nt)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand): ref base . (46/63); new base C (4/2); total (50/65)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 116 total aligned reads displayed.

GGAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAG  >  NC_002947/939781‑940077

ggAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑Gt                                                                                                                                                     <  1:1237917/150‑1 (MQ=255)

 gAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑Gtg                                                                                                                                                    <  1:3845185/150‑1 (MQ=255)

 gAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑Gtg                                                                                                                                                    <  1:3845539/150‑1 (MQ=255)

   ggCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑Gtgct                                                                                                                                                  <  1:5082395/150‑1 (MQ=255)

   ggCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑Gtgct                                                                                                                                                  <  2:60004/150‑1 (MQ=255)

     cGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑Gggctgc                                                                                                                                                >  2:4697343/1‑150 (MQ=255)

      gTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGGCCGCCCCCCGACG‑GGGCTGcc                                                                                                                                               >  1:2387209/1‑150 (MQ=14)

      gTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGcc                                                                                                                                               >  2:2006969/1‑150 (MQ=255)

      gTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGcc                                                                                                                                               <  2:541633/150‑1 (MQ=255)

      gTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGcc                                                                                                                                               >  2:2006794/1‑150 (MQ=255)

         gcctgccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGGCAGCAC‑CGTACGCGTGCTGGCCgc                                                                                                                                            <  2:587181/150‑2 (MQ=14)

          cctgccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGt                                                                                                                                           <  1:3437291/150‑1 (MQ=255)

           ctgccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGtt                                                                                                                                          <  2:5112676/150‑1 (MQ=255)

             gccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCTGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCg                                                                                                                                        >  2:3805114/1‑150 (MQ=255)

               cGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGgc                                                                                                                                      <  1:3666500/150‑1 (MQ=255)

                 cATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGGAAG‑GTGCTGCCCCGTTCGgcgc                                                                                                                                    >  1:4204736/1‑150 (MQ=255)

                      gTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACCCGGACG‑GTGCGGCCCGGGTCGGCGCGGCCg                                                                                                                               >  2:1947386/1‑150 (MQ=255)

                        ggCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGt                                                                                                                             <  2:2170665/150‑1 (MQ=255)

                        ggCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGt                                                                                                                             <  2:2170380/150‑1 (MQ=255)

                        ggCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGt                                                                                                                             <  2:4088174/150‑1 (MQ=255)

                          ccGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGa                                                                                                                           <  2:2883392/150‑1 (MQ=255)

                           cGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGAc                                                                                                                          >  2:3555481/1‑150 (MQ=255)

                           cGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGAc                                                                                                                          >  2:1728194/1‑150 (MQ=255)

                            gTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGGAGCGCCAACGAGCACGCTCCCCAGCACAGAGGTGCCGTCGGCGCCCCGCACGCGTCCC‑GTGCTGCCCGGGGCGGCGCGGCCGGGGAcc                                                                                                                         >  1:4235285/1‑150 (MQ=14)

                              gACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTGCG‑GTGCTGCCCGGTTGGGCGCGGCCGGTGACCAc                                                                                                                       >  1:3367886/1‑150 (MQ=255)

                                 ttGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGGGGCCACCCCGCGGAAG‑GTGCGGGCCCGGTGGGCGCGGCCGGTGGCCAcccg                                                                                                                    >  1:187465/1‑147 (MQ=14)

                                 ttGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCt                                                                                                                    >  2:3102224/1‑150 (MQ=255)

                                 ttGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCACCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCt                                                                                                                    >  2:1449245/1‑150 (MQ=255)

                                 ttGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCACCACCCCGACC‑GGGCGGCCCCGGTAGGCGCGGCCCGGGACCAccca                                                                                                                    >  1:612189/1‑147 (MQ=14)

                                    aaGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGcc                                                                                                                 <  2:36810/150‑1 (MQ=40)

                                      gCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGGGGCCCGCCCGCGTACG‑GTGCGGCCCCGGTGGGCGCGGCCGGGGACCACCCAGCCgg                                                                                                               >  1:3643712/1‑149 (MQ=14)

                                      gCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGGCGGTGTCCCCGCTGCCgg                                                                                                               >  2:1348978/1‑149 (MQ=255)

                                      gCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGt                                                                                                               <  2:4203101/150‑1 (MQ=39)

                                      gCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGt                                                                                                               <  2:4203103/150‑1 (MQ=39)

                                          cTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCCCGCGTACC‑GTGCTGCCCGGTGCGGCGCGGCCGGGGGCCACGCGGCCggcggc                                                                                                           >  2:468017/1‑150 (MQ=11)

                                          cTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGc                                                                                                           <  1:4138757/150‑1 (MQ=38)

                                          cTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGc                                                                                                           <  2:3183947/150‑1 (MQ=38)

                                           tCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGGGGCCCGCACGCGTACG‑GGGCTGCCCGGTTGGGCGCGGCCGGTGGCCACCCGGCCggcggcg                                                                                                          >  1:3127172/1‑149 (MQ=11)

                                               cGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGa                                                                                                      >  2:728133/1‑150 (MQ=34)

                                                gTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGAt                                                                                                     <  2:828224/150‑1 (MQ=34)

                                                   tATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCAc                                                                                                  <  1:219684/150‑1 (MQ=34)

                                                   tATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGAACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCAc                                                                                                  <  2:1675469/150‑1 (MQ=21)

                                                   tATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCCTACG‑GTGCTTCCCGGTTCGGCGCGGCCGGTGCCCACGCAGCCGTCGGCAATTATCAc                                                                                                  >  2:169877/1‑150 (MQ=11)

                                                     tCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCa                                                                                                <  1:3840346/150‑1 (MQ=34)

                                                     tCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCa                                                                                                <  1:3839860/150‑1 (MQ=34)

                                                     tCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCa                                                                                                <  1:3839953/150‑1 (MQ=34)

                                                     tCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGGACG‑GTGCTGCCCGGGTCGGCGCGGCCGGTGACCAAGCTGCCGTCGGCCGTGATAACCa                                                                                                >  1:2074257/1‑150 (MQ=11)

                                                     tCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCACCACGCCTAAC‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCAATGATCAccc                                                                                                >  2:635117/1‑149 (MQ=11)

                                                            cGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTtcggt                                                                                         <  2:1344554/150‑1 (MQ=34)

                                                             gggTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTtcggtc                                                                                        <  1:510321/150‑1 (MQ=34)

                                                                        gATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAgg                                                                             <  2:292415/150‑1 (MQ=34)

                                                                          tGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCCGTTCGGCGCGGCCGGTGACCACGCAGCCGTCGGCATTGATCACCAGTTCGGTCGGGGCGGCAgggg                                                                           >  1:2541630/1‑150 (MQ=255)

                                                                            ccgAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCTAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCg                                                                         <  2:1146208/148‑1 (MQ=21)

                                                                            cTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCg                                                                         <  1:4662399/150‑1 (MQ=34)

                                                                              gAAGCTGCCGGTTGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                       <  1:1419102/150‑1 (MQ=21)

                                                                              gAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                       <  1:2094629/150‑1 (MQ=34)

                                                                                 gCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCACGCGTACG‑GTGCTGCCCGGGTCGGCGCGGCCGGGGACCACGCTGCCGGCGGCATTGATAACCAGTTCGGTCGGGGCGGCAGGGGCgggggt                                                                    >  2:1382071/1‑146 (MQ=11)

                                                                                    gCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATc                                                                 <  1:4292351/150‑1 (MQ=34)

                                                                                        gTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCACGCGGACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCGGCCGTCGGCAGTGATCACCAGTTCGGTCGGGGCGGCCGGGGCGGTGATCTCcggg                                                             >  2:4622568/1‑149 (MQ=11)

                                                                                              ccAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCCGTGACCACGCTGCCGTCGGCATTGATCACCAGTTAGGTCGGTGCGGCCGGGGCGGTGATATCcggtgcggtg                                                       >  2:4999458/1‑150 (MQ=255)

                                                                                               cAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCCCGCTGCCGGCGGCATTGCTCACCAGTTCGGTCGGTGCGGCCGGGGCGGTGATATCcggtgcggtgg                                                      >  2:1606459/1‑149 (MQ=255)

                                                                                               cAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGa                                                      <  2:2195316/150‑1 (MQ=34)

                                                                                               cAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGa                                                      >  2:4372504/1‑150 (MQ=34)

                                                                                                aaCTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGAt                                                     <  1:1728194/150‑1 (MQ=34)

                                                                                                aaCTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGGTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCCTTGATCACCACGTCGGTCGGGGCGGCAGGGGCGGTGATATCcggtgcggtgct                                                     >  2:4094279/1‑148 (MQ=11)

                                                                                                   tATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtg                                                  >  2:2966283/1‑150 (MQ=34)

                                                                                                   tATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGAACACGCTGCCGTCGGCATTGATAAACAGTTCGGTCGGGGCGGCAGGGGCGGTGATATACGGTGCGGTGATTtg                                                  >  2:4562964/1‑150 (MQ=11)

                                                                                                    aTCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCGCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgg                                                 <  2:1850901/150‑1 (MQ=18)

                                                                                                    aTCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGTCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgg                                                 <  2:4987474/150‑1 (MQ=18)

                                                                                                     tCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtggg                                                <  1:4372504/150‑1 (MQ=32)

                                                                                                      cACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGCCCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATCTCCGGTGCGGTGATTtgggt                                               <  2:3530026/150‑1 (MQ=17)

                                                                                                       aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtg                                              <  2:4602750/150‑1 (MQ=32)

                                                                                                       aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtg                                              <  2:4296899/150‑1 (MQ=32)

                                                                                                       aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtg                                              <  2:1318780/150‑1 (MQ=32)

                                                                                                         gCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATACGGTGCGGTGCTTtgggtggg                                            >  2:1630538/1‑150 (MQ=17)

                                                                                                             cccAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGc                                        <  1:3555481/150‑1 (MQ=32)

                                                                                                             cccAGCACAGTGGTGCCGTCGGCCGCCCGCACCCCGACG‑GGGCTGCCCGGGTCGGCGCGGCCGGGGCCCACCCGGCCGGCGGCAGTGATCACCACGTCGGGCGGTGCGGCACGGGCGGGGGTATCCCGTGCGGGGAGTTGGGTGGGCCGc                                        >  1:1174529/1‑150 (MQ=2)

                                                                                                              ccAGCACAGTGGTGCCGTCGGCGGCCAGCACGCCTACG‑GGGCTGCCCCGGTCGGCGCGGCCGGTGTCCACCCGGCCGGCGGCATAGATCACCAGGTCGGTCGGTGCGGCAGGGGCGGGGATATCCGGTGCGGTGATTTGGGGGGGCcgcg                                       >  1:2091245/1‑150 (MQ=11)

                                                                                                               cAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGa                                      >  1:3805796/1‑150 (MQ=32)

                                                                                                                aGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGTCCACGCTGCCGTCGGCCTTGATAACCACTTCGGTCGGTGCGGCACGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAg                                     >  2:4419763/1‑150 (MQ=11)

                                                                                                                aGCACAGTGGTGCCGTCGGCGGCCACCCCGCGTACG‑GTGCTTCCCGGTTGGGCGCGGCCCGGGAACACCCTGCCGTCGGCCATTATGACCACGTGGGGCGGTGCGGCAGGGGCGGTGATATCCGGGGCGGTGATTTGGGTGGGCAGCGAg                                     >  2:4419540/1‑150 (MQ=11)

                                                                                                                 gCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAgg                                    <  1:370466/150‑1 (MQ=32)

                                                                                                                 gCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAgg                                    <  1:2928648/150‑1 (MQ=32)

                                                                                                                 gCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAgg                                    <  1:2928637/150‑1 (MQ=32)

                                                                                                                   acaGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGaa                                  <  2:3870618/150‑1 (MQ=32)

                                                                                                                    caGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGaaa                                 <  2:3643712/150‑1 (MQ=32)

                                                                                                                     aGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAac                                <  1:4287513/150‑1 (MQ=32)

                                                                                                                     aGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAac                                <  1:4287838/150‑1 (MQ=32)

                                                                                                                              gTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGtc                       <  2:187465/150‑1 (MQ=32)

                                                                                                                               tCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGtca                      <  1:2998892/150‑1 (MQ=32)

                                                                                                                                cggcggCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGtcat                     <  1:2006969/150‑1 (MQ=32)

                                                                                                                                cggcggCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGtcat                     <  1:2006794/150‑1 (MQ=32)

                                                                                                                                aaccgggCAGCAC‑CGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGGCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCCGCGAGGAAACACCGCCGtcat                     >  2:366330/4‑150 (MQ=11)

                                                                                                                                  ccgggCAGCAC‑CGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGtcatca                   >  1:2772276/2‑150 (MQ=11)

                                                                                                                                     ggcAGCAC‑CGTACGCGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                       >  1:3015061/3‑95 (MQ=11)

                                                                                                                                     ggcAGCAC‑CGTACGCGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                       <  2:3015061/93‑1 (MQ=11)

                                                                                                                                       cAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGc              >  2:1903809/1‑150 (MQ=32)

                                                                                                                                         gCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGcca            >  1:2895502/1‑150 (MQ=32)

                                                                                                                                                 aCG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTcc    <  1:4697343/150‑1 (MQ=32)

                                                                                                                                                   g‑gTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCcag  >  2:848454/1‑150 (MQ=32)

GGAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAG  >  NC_002947/939781‑940077

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