BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,357,033	A→G	6.1%	R48R (CGT→CGC)	PP_2950 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,357,033	0	A	G	6.1%	991.4 / 31.0	310	R48R (CGT→CGC)	PP_2950	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (159/132); new base G (16/3); total (175/135)
Fisher's exact test for biased strand distribution p-value = 1.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 315 total aligned reads displayed.

TCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGT  >  NC_002947/3356884‑3357182

tcgtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTa                                                                                                                                                       <  1:5998934/150‑1 (MQ=255)

      cACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCAt                                                                                                                                                 <  1:8883925/150‑1 (MQ=255)

       aCTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCAtt                                                                                                                                                >  1:1906464/1‑150 (MQ=255)

        cTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTg                                                                                                                                               >  1:3071629/1‑150 (MQ=255)

        cTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTg                                                                                                                                               <  2:8475216/150‑1 (MQ=255)

         tCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTgg                                                                                                                                              <  2:3369344/150‑1 (MQ=255)

         tCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTgg                                                                                                                                              <  1:8963699/150‑1 (MQ=255)

           gTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGgt                                                                                                                                            <  1:475754/150‑1 (MQ=255)

             aGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGgtgt                                                                                                                                          >  1:2750473/1‑150 (MQ=255)

                      ggACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTAGCGAGTTAGGAAGTACCGAGCTTCTTCAAGGTTGTTTGCAAACTCGGCACTGTTACTCCATTGGGTGTCCACCGtat                                                                                                                                 >  2:3539441/1‑150 (MQ=255)

                          cgcgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTgg                                                                                                                                              >  1:7263176/1‑133 (MQ=255)

                          cgcgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGc                                                                                                                             >  1:2821490/1‑150 (MQ=255)

                            cgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCa                                                                                                                           <  2:2428689/150‑1 (MQ=255)

                            cgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCa                                                                                                                           <  1:5281499/150‑1 (MQ=255)

                             gATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAg                                                                                                                          <  1:7570807/150‑1 (MQ=255)

                                        cGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCATACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCtttt                                                                                                               <  2:2835404/150‑1 (MQ=255)

                                         ggTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTc                                                                                                              <  2:4255422/150‑1 (MQ=255)

                                          gTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCg                                                                                                             >  2:3747280/1‑150 (MQ=255)

                                           tGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCgg                                                                                                            <  1:8414794/150‑1 (MQ=255)

                                            gAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCggg                                                                                                           >  1:1812276/1‑150 (MQ=255)

                                             aaGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTAGGGt                                                                                                          >  1:9165700/1‑150 (MQ=255)

                                                  ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                     >  1:7934174/1‑150 (MQ=255)

                                                  ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                     >  2:420494/1‑150 (MQ=255)

                                                  ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                     >  1:101337/1‑150 (MQ=255)

                                                   cAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCgg                                                                                                    <  2:3448367/150‑1 (MQ=255)

                                                     gagCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGt                                                                                                  >  1:537008/1‑150 (MQ=255)

                                                     gagCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGt                                                                                                  >  1:2434523/1‑150 (MQ=255)

                                                       gCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                                                                                <  2:6900060/150‑1 (MQ=255)

                                                       gCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                                                                                >  1:7413657/1‑150 (MQ=255)

                                                        ctctTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAtt                                                                                               >  2:4000554/1‑150 (MQ=255)

                                                         tctTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTc                                                                                              <  2:7648900/150‑1 (MQ=255)

                                                           tttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCa                                                                                            >  1:8149253/1‑150 (MQ=255)

                                                                tCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCAt                                                                                       >  2:6299684/1‑150 (MQ=255)

                                                                  aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                                                                     >  1:2875830/1‑150 (MQ=255)

                                                                  aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                                                                     >  2:8329264/1‑150 (MQ=255)

                                                                   aCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCt                                                                                    >  2:9142436/1‑150 (MQ=255)

                                                                    ccGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTc                                                                                   <  1:6461282/150‑1 (MQ=255)

                                                                    ccGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTc                                                                                   <  2:3837772/150‑1 (MQ=255)

                                                                     cggtactCCAGTTCCTTATTCCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTcc                                                                                  <  1:7678563/144‑2 (MQ=255)

                                                                     cGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCa                                                                                  >  2:6830836/1‑150 (MQ=255)

                                                                      gATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAg                                                                                 >  2:8904025/1‑150 (MQ=255)

                                                                       gtactCCAGTTCCTTATTCCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGc                                                                                >  1:7229481/5‑150 (MQ=255)

                                                                          gTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAg                                                                             >  1:8265565/1‑150 (MQ=255)

                                                                          gTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAg                                                                             >  2:4383047/1‑150 (MQ=255)

                                                                           tCCAGTTCCTTATTCCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAgg                                                                            >  1:5070848/1‑149 (MQ=255)

                                                                             cAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCg                                                                          >  1:1934175/1‑150 (MQ=255)

                                                                               gTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGa                                                                        <  2:5508275/150‑1 (MQ=255)

                                                                                ttctttctTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGAc                                                                       <  2:1819196/150‑1 (MQ=255)

                                                                                   tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                                                                    <  2:2216847/150‑1 (MQ=255)

                                                                                   cttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtt                                                                    >  2:1619703/5‑150 (MQ=255)

                                                                                   cttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtt                                                                    >  2:1425408/5‑150 (MQ=255)

                                                                                   cttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtt                                                                    >  1:464019/5‑150 (MQ=255)

                                                                                    ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                                                                   <  1:4866419/150‑1 (MQ=255)

                                                                                    ttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCttt                                                                   >  2:7278749/4‑150 (MQ=255)

                                                                                     tattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtttt                                                                  <  2:77904/148‑1 (MQ=255)

                                                                                      cttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTc                                                                 >  2:5994744/1‑150 (MQ=255)

                                                                                       ttcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCTTTTcg                                                                >  2:7205101/1‑149 (MQ=255)

                                                                                       ttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCa                                                                <  2:452263/150‑1 (MQ=255)

                                                                                         ccACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAAt                                                              >  1:3839940/1‑150 (MQ=255)

                                                                                          cACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATc                                                             >  2:2577482/1‑150 (MQ=255)

                                                                                          cACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATc                                                             >  1:5771161/1‑150 (MQ=255)

                                                                                           aCTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCg                                                            >  2:4880557/1‑150 (MQ=255)

                                                                                               gCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGa                                                        >  2:1512059/1‑150 (MQ=255)

                                                                                                cGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGAt                                                       <  2:5268225/150‑1 (MQ=255)

                                                                                                 gAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGAtt                                                      >  2:2354668/1‑150 (MQ=255)

                                                                                                 gAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGAtt                                                      >  2:6972451/1‑150 (MQ=255)

                                                                                                    ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc                                                   >  2:5006492/1‑150 (MQ=255)

                                                                                                        gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtttt                                                                  <  2:6636007/131‑1 (MQ=255)

                                                                                                        gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                               <  1:5171175/150‑1 (MQ=255)

                                                                                                        gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                               >  2:5529155/1‑150 (MQ=255)

                                                                                                        gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                               >  1:4595796/1‑150 (MQ=255)

                                                                                                        gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                               >  2:6524965/1‑150 (MQ=255)

                                                                                                            tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                           <  1:1712646/150‑1 (MQ=255)

                                                                                                            tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                           >  2:2047823/1‑150 (MQ=255)

                                                                                                             cccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAag                                          <  2:323400/150‑1 (MQ=255)

                                                                                                             cccGAGCTTCTTCACGGTTGTTCGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAag                                          <  2:7902719/150‑1 (MQ=255)

                                                                                                              ccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAggg                                                                                                                        <  1:4198699/71‑1 (MQ=255)

                                                                                                               cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                                        <  1:6885362/150‑1 (MQ=255)

                                                                                                               cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                                        <  2:1270700/150‑1 (MQ=255)

                                                                                                                gAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAaggag                                       <  2:877221/150‑1 (MQ=255)

                                                                                                                gAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAaggag                                       <  1:8635875/150‑1 (MQ=255)

                                                                                                                  gCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAggagga                                     <  1:3226551/150‑1 (MQ=255)

                                                                                                                    ttcttcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAg                                   >  2:4608910/1‑150 (MQ=255)

                                                                                                                      cttcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTg                                 <  1:8837307/150‑1 (MQ=255)

                                                                                                                        tcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGc                               <  2:2821487/150‑1 (MQ=255)

                                                                                                                        tcACGGTTGTTAGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATACCAGCCAACCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGc                               >  2:7572308/1‑150 (MQ=255)

                                                                                                                          aCGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCgg                             >  2:8671766/1‑150 (MQ=255)

                                                                                                                            ggTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGACGTGGCGGAc                           >  1:4433971/1‑150 (MQ=255)

                                                                                                                            ggTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGAc                           >  1:9035653/1‑150 (MQ=255)

                                                                                                                              ttgttTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGa                         <  2:3082489/150‑1 (MQ=255)

                                                                                                                                gttTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgt                       <  1:1234662/150‑1 (MQ=255)

                                                                                                                                gttTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgt                       <  1:7740941/150‑1 (MQ=255)

                                                                                                                                  ttGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgtgt                     <  2:9079699/150‑1 (MQ=255)

                                                                                                                                      aaaCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTc                 >  1:2926160/1‑150 (MQ=255)

                                                                                                                                      aaaCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTc                 <  1:8012039/150‑1 (MQ=255)

                                                                                                                                        aCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGt               >  1:6268789/1‑150 (MQ=255)

                                                                                                                                            ggCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGa           <  1:5495868/150‑1 (MQ=255)

                                                                                                                                             gCACGGTTGCGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGAc          >  1:6602155/1‑150 (MQ=255)

                                                                                                                                                   ttACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACAtggtg    <  2:2659993/150‑1 (MQ=255)

                                                                                                                                                   ttACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACAtggtg    <  2:1354169/150‑1 (MQ=255)

TCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGT  >  NC_002947/3356884‑3357182

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