BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,357,015	T→A	6.1%	A54A (GCA→GCT) ‡	PP_2950 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,357,015	0	T	A	6.1%	1021.1 / 36.5	304	A54A (GCA→GCT) ‡	PP_2950	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (154/130); new base A (16/3); total (170/134)
Fisher's exact test for biased strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 313 total aligned reads displayed.

AACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGT  >  NC_002947/3356866‑3357163

aaCCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGttt                                                                                                                                                      >  2:5495868/1‑150 (MQ=255)

        gCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTc                                                                                                                                              <  2:6170312/150‑1 (MQ=255)

         cccGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCg                                                                                                                                             >  1:7552385/1‑150 (MQ=255)

           cGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGc                                                                                                                                           >  1:5893351/1‑150 (MQ=255)

                cgtcgtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGt                                                                                                                                      <  1:5999407/150‑1 (MQ=255)

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

                   cgtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTAc                                                                                                                                   <  1:8356757/150‑1 (MQ=255)

                    gtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACg                                                                                                                                  >  1:7316176/1‑150 (MQ=255)

                    gtGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACg                                                                                                                                  <  1:6123976/150‑1 (MQ=255)

                      gCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGcc                                                                                                                                >  2:710702/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                 gCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTcc                                                                                                                     >  2:8814638/1‑150 (MQ=255)

                                      aCGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGt                                                                                                                >  2:7816409/1‑150 (MQ=255)

                                      aCGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGAAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGt                                                                                                                >  1:7977398/1‑150 (MQ=255)

                                       cGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTATGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGta                                                                                                               >  2:1565074/1‑150 (MQ=255)

                                        ggACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGtat                                                                                                              >  1:8218471/1‑150 (MQ=255)

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

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

                                            cgcgATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGc                                                                                                          <  2:7396662/150‑1 (MQ=255)

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

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

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

                                               gATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAg                                                                                                       >  2:1280795/1‑150 (MQ=255)

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

                                                           ggTGAAGGCCCAGAGCTCTTTGATCAACCGAGAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTc                                                                                           <  2:6387510/150‑1 (MQ=255)

                                                           ggTGAAGGCCCAGAGCACATTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTc                                                                                           <  1:4771186/150‑1 (MQ=255)

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

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

                                                                    ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                  <  2:237990/150‑1 (MQ=255)

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

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

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

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

                                                                       gagCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCTGGTTCGCGGGt                                                                               >  1:6982953/1‑150 (MQ=255)

                                                                           tctTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTc                                                                           >  1:7437472/1‑150 (MQ=255)

                                                                             tttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGgt                                                                                                                         >  1:9062882/1‑102 (MQ=255)

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

                                                                             tttGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCa                                                                         <  1:4584052/150‑1 (MQ=255)

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

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

                                                                               tGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGc                                                                       >  2:2270646/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                                                          tactCCAGTTCCTTATTCCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGca                                                            >  1:4254847/4‑149 (MQ=255)

                                                                                           aGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGa                                                           <  2:6587529/150‑1 (MQ=255)

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

                                                                                             tCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGt                                                         <  1:7800667/150‑1 (MQ=255)

                                                                                             tCCAGTTCTTTCTTCAACTTGCGAGTTATGAAGTCACGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCAGTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAga                                                         >  2:86514/1‑149 (MQ=255)

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

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

                                                                                                     tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                                                 >  2:3821212/1‑150 (MQ=255)

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

                                                                                                      ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGTCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                                                <  1:7583285/150‑1 (MQ=255)

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

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

                                                                                                        cttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTc                                              >  2:5994744/1‑150 (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)

                                                                                                                  cGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGag                                    >  1:2437606/1‑149 (MQ=255)

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

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

                                                                                                                      ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc                                <  1:4952041/150‑1 (MQ=255)

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

                                                                                                                         ggAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCATTTCAATCGCCGATTCGCGAg                             >  2:7755081/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                                                ccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                                                             <  1:182504/95‑1 (MQ=255)

                                                                                                                                ccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                                                             >  2:182504/1‑95 (MQ=255)

                                                                                                                                 cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                     >  2:692695/1‑150 (MQ=255)

                                                                                                                                 cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                     >  2:4207656/1‑150 (MQ=255)

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

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

                                                                                                                                         ttcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTgg             >  1:345938/1‑150 (MQ=255)

                                                                                                                                          tcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGc            >  2:7572396/1‑150 (MQ=255)

                                                                                                                                          tcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGc            >  2:2772944/1‑150 (MQ=255)

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

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

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

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

AACCAGTTGCCCGAATCGTCGTGCCACTCGTAGGCGTAACGGACCGCGATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGT  >  NC_002947/3356866‑3357163

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