BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,609,200	C→G	5.0%	P12R (CCG→CGG)	PP_4076 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,609,200	0	C	G	5.0%	946.0 / 28.8	278	P12R (CCG→CGG)	PP_4076	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (130/134); new base G (8/6); total (138/140)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 283 total aligned reads displayed.

TCTCGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGAC  >  NC_002947/4609051‑4609348

  tcGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGt                                                                                                                                                    >  2:5324606/1‑150 (MQ=255)

    gggcGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTcc                                                                                                                                                  >  1:1529299/4‑150 (MQ=255)

     gACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCg                                                                                                                                                 <  2:3332998/150‑1 (MQ=255)

       cGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGTAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGcc                                                                                                                                               <  1:2259941/150‑1 (MQ=255)

          ggCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATc                                                                                                                                            <  2:4077456/150‑1 (MQ=255)

          ggCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATc                                                                                                                                            >  1:4386505/1‑150 (MQ=255)

          ggCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATc                                                                                                                                            >  2:3577161/1‑150 (MQ=255)

           gCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATcc                                                                                                                                           >  2:1756143/1‑150 (MQ=255)

              ggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCggg                                                                                                                                        >  2:5583598/1‑150 (MQ=255)

              ggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCggg                                                                                                                                        >  1:4086245/1‑150 (MQ=255)

              ggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCggg                                                                                                                                        <  2:1953512/150‑1 (MQ=255)

               gTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGa                                                                                                                                       >  2:630234/1‑150 (MQ=255)

               gTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGa                                                                                                                                       <  1:2992389/150‑1 (MQ=255)

               gTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGa                                                                                                                                       <  2:3738850/150‑1 (MQ=255)

                tGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAg                                                                                                                                      <  1:5188587/150‑1 (MQ=255)

                tGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAg                                                                                                                                      >  2:1939509/1‑150 (MQ=255)

                 gATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAgg                                                                                                                                     <  2:3406281/150‑1 (MQ=255)

                   tCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCg                                                                                                                                   >  1:4385176/1‑150 (MQ=255)

                     cGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtc                                                                                                                                 <  1:1374139/150‑1 (MQ=255)

                     cGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtc                                                                                                                                 <  1:738122/150‑1 (MQ=255)

                                cAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGaaa                                                                                                                      >  2:224137/1‑150 (MQ=255)

                                cAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCCGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGaaa                                                                                                                      <  1:6272/150‑1 (MQ=255)

                                    tctcGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAgcgc                                                                                                                  >  2:4385263/1‑150 (MQ=255)

                                       cgctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTc                                                                                                               <  2:1081964/150‑1 (MQ=255)

                                        gctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCg                                                                                                              >  1:5313243/1‑150 (MQ=255)

                                         ctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGc                                                                                                             <  1:1839966/150‑1 (MQ=255)

                                          tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGca                                                                                                            <  1:2101174/150‑1 (MQ=255)

                                          tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGca                                                                                                            >  2:3029728/1‑150 (MQ=255)

                                            gcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGcaca                                                                                                          >  2:4012948/1‑150 (MQ=255)

                                             ccGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACaa                                                                                                         <  2:4992589/150‑1 (MQ=255)

                                             ccGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACaa                                                                                                         <  1:2326877/150‑1 (MQ=255)

                                                    cccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTcgcg                                                                                                  <  1:3858696/150‑1 (MQ=255)

                                                     ccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCgg                                                                                                 >  2:475254/1‑150 (MQ=255)

                                                     ccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGGTCGCgg                                                                                                 >  1:1266946/1‑150 (MQ=255)

                                                       gCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGcc                                                                                               >  2:3268808/1‑150 (MQ=255)

                                                        cACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCa                                                                                              <  2:2707233/150‑1 (MQ=255)

                                                        cACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCa                                                                                              <  2:5180986/150‑1 (MQ=255)

                                                         aCTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGCACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCAt                                                                                             >  2:769639/1‑150 (MQ=255)

                                                           ttCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCa                                                                                           >  2:1552459/1‑150 (MQ=255)

                                                             cAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAg                                                                                         <  2:4428672/150‑1 (MQ=255)

                                                              aaGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                        <  1:5594748/150‑1 (MQ=255)

                                                              aaGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                        <  2:259224/150‑1 (MQ=255)

                                                              aaGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                        <  2:2077948/150‑1 (MQ=255)

                                                              aaGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                        <  2:1868588/150‑1 (MQ=255)

                                                                gTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGc                                                                                      <  2:2465993/150‑1 (MQ=255)

                                                                gTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGc                                                                                      <  2:1022425/150‑1 (MQ=255)

                                                                   cATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCAc                                                                                   <  1:50794/150‑1 (MQ=255)

                                                                   cATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCAc                                                                                   <  2:551108/150‑1 (MQ=255)

                                                                   cATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCAc                                                                                   >  1:1147858/1‑150 (MQ=255)

                                                                      gggCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTc                                                                                <  2:5415937/150‑1 (MQ=255)

                                                                         cgcgAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGaa                                                                             <  2:4149226/150‑1 (MQ=255)

                                                                            gAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACg                                                                          >  1:2469212/1‑150 (MQ=255)

                                                                            gAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACg                                                                          >  2:2945/1‑150 (MQ=255)

                                                                            gAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACg                                                                          >  1:1589835/1‑150 (MQ=255)

                                                                              aCTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGtt                                                                        <  1:1693008/150‑1 (MQ=255)

                                                                              aCTTCGAAATGTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGtt                                                                        <  2:5272182/150‑1 (MQ=255)

                                                                               cTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                       <  2:1301704/150‑1 (MQ=255)

                                                                               cTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                       <  1:1977645/150‑1 (MQ=255)

                                                                               cTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                       <  2:2851499/150‑1 (MQ=255)

                                                                                ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      <  1:1910759/150‑1 (MQ=255)

                                                                                ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      >  1:93206/1‑150 (MQ=255)

                                                                                ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      >  2:5609695/1‑150 (MQ=255)

                                                                                ttCGAAATTTCCCTGGTAGGCCGACTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      >  2:683923/1‑150 (MQ=255)

                                                                                 tCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCt                                                                     <  1:4767756/150‑1 (MQ=255)

                                                                                    aaaTTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGaa                                                                  <  1:3547266/150‑1 (MQ=255)

                                                                                          cccTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCg                                                            <  2:247645/150‑1 (MQ=255)

                                                                                              ggTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCaa                                                        >  2:4632590/1‑150 (MQ=255)

                                                                                                   gCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTg                                                   <  2:3675759/150‑1 (MQ=255)

                                                                                                      ggCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTc                                                >  1:2673615/1‑150 (MQ=255)

                                                                                                      ggCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTc                                                >  1:4484504/1‑150 (MQ=255)

                                                                                                        cTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGACATCAGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACATCGAAACGTTGCTGAACGCTCGCCAAGGTAGCTccc                                              >  2:29831/1‑150 (MQ=255)

                                                                                                         tCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCa                                             <  1:2818988/150‑1 (MQ=255)

                                                                                                              acGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGagca                                        >  1:3203610/1‑150 (MQ=255)

                                                                                                                       gTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCTGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGc                               >  1:510308/1‑150 (MQ=255)

                                                                                                                         ggTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTg                             <  2:89560/150‑1 (MQ=255)

                                                                                                                         ggTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTg                             <  2:3097212/150‑1 (MQ=255)

                                                                                                                         ggTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTg                             >  1:3006570/1‑150 (MQ=255)

                                                                                                                          gtgtcttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACACTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                            <  1:837653/147‑1 (MQ=255)

                                                                                                                          gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                            <  1:4632590/150‑1 (MQ=255)

                                                                                                                            ttcttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGc                          <  2:2523286/150‑1 (MQ=255)

                                                                                                                            tacttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGc                          >  2:5216160/3‑150 (MQ=255)

                                                                                                                               ttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTg                       >  2:3490401/1‑150 (MQ=255)

                                                                                                                               ttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTg                       >  2:2180490/1‑150 (MQ=255)

                                                                                                                                tcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTgc                      <  1:3577161/150‑1 (MQ=255)

                                                                                                                                 cGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTgcg                     >  1:3642068/1‑150 (MQ=255)

                                                                                                                                  gATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTgcgc                    <  2:4828921/150‑1 (MQ=255)

                                                                                                                                   aTCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGcgcg                   >  2:3917581/1‑150 (MQ=255)

                                                                                                                                    tCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGa                  <  1:2977618/150‑1 (MQ=255)

                                                                                                                                      gCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACt                >  1:3452166/1‑150 (MQ=255)

                                                                                                                                       ccggATGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTg                                                   >  2:4063209/4‑114 (MQ=255)

                                                                                                                                       ccggATGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTg                                                   <  1:4063209/111‑1 (MQ=255)

                                                                                                                                       cTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACtt               >  2:2328329/1‑150 (MQ=255)

                                                                                                                                        ttGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTc              >  2:1423897/1‑150 (MQ=255)

                                                                                                                                          gATGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCaa            >  1:312911/1‑150 (MQ=255)

                                                                                                                                          gATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCaa            >  2:3831897/1‑150 (MQ=255)

                                                                                                                                          gATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCaa            >  2:1403422/1‑150 (MQ=255)

                                                                                                                                          gATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCaa            >  2:4268592/1‑150 (MQ=255)

                                                                                                                                             ggCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACgg         <  1:2160173/150‑1 (MQ=255)

                                                                                                                                              gCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACggg        >  1:4539931/1‑150 (MQ=255)

                                                                                                                                                    ccGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGAc  >  1:2143574/1‑150 (MQ=255)

TCTCGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGAC  >  NC_002947/4609051‑4609348

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