BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,113,955	T→A	6.3%	H118Q (CAT→CAA)	PP_5362 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,113,955	0	T	A	6.3%	565.1 / 20.3	175	H118Q (CAT→CAA)	PP_5362	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (76/88); new base A (5/6); total (81/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 176 total aligned reads displayed.

GGGCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAACC  >  NC_002947/6113807‑6114103

gggCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAtt                                                                                                                                                     >  1:1351675/1‑150 (MQ=255)

 ggCACAAACCACCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTg                                                                                                                                                    >  1:4645594/1‑150 (MQ=255)

 ggCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTg                                                                                                                                                    >  1:975794/1‑150 (MQ=255)

   cacaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCa                                                                                                                                                  <  2:3306171/150‑1 (MQ=255)

      aaaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCt                                                                                                                                               >  1:533922/1‑150 (MQ=255)

       aaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCtt                                                                                                                                              <  2:110663/150‑1 (MQ=255)

       aaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCtt                                                                                                                                              >  1:4854120/1‑150 (MQ=255)

        aCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTg                                                                                                                                             >  2:500441/1‑150 (MQ=255)

        aCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTg                                                                                                                                             >  1:1673390/1‑150 (MQ=255)

        aCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTg                                                                                                                                             <  2:2386957/150‑1 (MQ=255)

          cAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGaa                                                                                                                                           <  2:3126916/150‑2 (MQ=255)

          cAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGAt                                                                                                                                           >  2:4874013/1‑150 (MQ=255)

          cAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGAAGCTTGAt                                                                                                                                           >  2:456681/1‑150 (MQ=255)

          cAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaa                                                                                                                                           <  1:681075/150‑2 (MQ=255)

           aaCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATc                                                                                                                                          >  1:2096549/1‑150 (MQ=255)

           aaCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCt                                                                                                                                               <  2:593707/145‑1 (MQ=255)

              cgCCGCTATGGCAACTTGACTTTTTTTAGCGGCCTGGCCCGCGTTGGTGTGGCCTGGGTGGGCTACAGCCGCCGGGTAGGCTGGACCACCTTGAACCGCGTATCCTGTTCCAGCGGGGAGCGCTTCAAGCTGGATTGCAGCTTGATCGAg                                                                                                                                       <  1:4760697/149‑1 (MQ=255)

              aGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAg                                                                                                                                       <  1:4135392/150‑1 (MQ=255)

              aGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAg                                                                                                                                       <  1:3563971/150‑1 (MQ=255)

               gCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAgg                                                                                                                                      >  1:3944330/1‑150 (MQ=255)

               gCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAgg                                                                                                                                      <  1:804583/150‑1 (MQ=255)

               gCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAgg                                                                                                                                      >  1:4523032/1‑150 (MQ=255)

               gCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAgg                                                                                                                                      <  2:3680834/150‑1 (MQ=255)

                  gCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGt                                                                                                                                   <  2:27057/150‑1 (MQ=255)

                    tATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGc                                                                                                                                 >  2:2039/1‑150 (MQ=255)

                          acaTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGc                                                                                                                                                >  2:2093496/1‑129 (MQ=255)

                          acaTTGACTTTTTGTCCCGGCCGGGCCAGCGAGGGTGAGGCCGGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACGTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATGGCAGCTTGATCGAGGGGTGCCGTTGa                                                                                                                           <  1:4695058/150‑1 (MQ=255)

                            aTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaag                                                                                                                                          >  2:2661904/1‑131 (MQ=255)

                            aTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaag                                                                                                                                          <  1:2661904/133‑3 (MQ=255)

                                aCTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTTAACCGCGTATCCAGTTCAAGCGGGGAGCGCTTCAAGCTGCATTGAAGCTTGATCGAGGGGTGCCGTTGACGcctc                                                                                                                     <  1:1486324/150‑1 (MQ=255)

                                  tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctcct                                                                                                                   >  2:933251/1‑150 (MQ=255)

                                  tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctcct                                                                                                                   >  2:3332315/1‑150 (MQ=255)

                                  tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctcct                                                                                                                   >  2:4115885/1‑150 (MQ=255)

                                      tttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCa                                                                                                               <  2:5312367/150‑1 (MQ=255)

                                       ttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAt                                                                                                              >  1:3138455/1‑150 (MQ=255)

                                             gCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGaa                                                                                                        >  2:3308294/1‑150 (MQ=255)

                                             gCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGaa                                                                                                        >  2:5140803/1‑150 (MQ=255)

                                              ccTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAAc                                                                                                       <  1:500441/150‑1 (MQ=255)

                                                  gCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGttt                                                                                                   <  1:2806040/150‑1 (MQ=255)

                                                    cAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGttttt                                                                                                 >  1:2263332/1‑150 (MQ=255)

                                                     aGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGtttttt                                                                                                >  2:3082667/1‑150 (MQ=255)

                                                     aGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGtttttt                                                                                                <  1:291460/150‑1 (MQ=255)

                                                       cGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTtct                                                                                              >  2:1388014/1‑150 (MQ=255)

                                                         aTGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACAACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTTCCGTTGAAGCATCCTTTCATTTTGATCGTTTTTTctct                                                                                            <  2:3115501/150‑1 (MQ=255)

                                                            gTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCAAGCGGGGAGCGCTTCAAGCTGCATTGCAGATTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctgct                                                                                         <  1:3738333/150‑1 (MQ=255)

                                                              gAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCtgctgc                                                                                       <  1:4905660/150‑1 (MQ=255)

                                                                 gCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAgg                                                                                                                                      <  2:1618324/100‑1 (MQ=255)

                                                                 gCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAgg                                                                                                                                      >  1:1618324/1‑100 (MQ=255)

                                                                 gCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTa                                                                                    >  2:688264/1‑150 (MQ=255)

                                                                    tgggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAttt                                                                                 >  2:2436079/1‑150 (MQ=255)

                                                                    tgggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAttt                                                                                 <  2:747618/150‑1 (MQ=255)

                                                                      ggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAtttt                                                                                >  1:1745055/1‑149 (MQ=255)

                                                                      ggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTg                                                                               >  1:544193/1‑150 (MQ=255)

                                                                       gtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGt                                                                              <  1:633239/150‑1 (MQ=255)

                                                                       gtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGt                                                                              <  1:3359816/150‑1 (MQ=255)

                                                                         gggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGTGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTaa                                                                            <  1:1097016/150‑1 (MQ=255)

                                                                         gggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTaa                                                                            <  1:601995/150‑1 (MQ=255)

                                                                             tACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctg                                                                                           >  1:2064234/1‑131 (MQ=255)

                                                                             tACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctg                                                                                           <  2:2064234/131‑1 (MQ=255)

                                                                             tACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCa                                                                        >  1:3915650/1‑150 (MQ=255)

                                                                                 gccgccGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCtt                                                                    <  2:3757551/150‑1 (MQ=255)

                                                                                 gccgccGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGAAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCtt                                                                    >  1:751643/1‑150 (MQ=255)

                                                                                    gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                                 <  2:2967623/150‑1 (MQ=255)

                                                                                    gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                                 <  1:4572748/150‑1 (MQ=255)

                                                                                     ccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAg                                                                <  1:4582861/150‑1 (MQ=255)

                                                                                       gcgTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCa                                                              >  2:3194024/1‑150 (MQ=255)

                                                                                            ggCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACgcgc                                                         <  2:3656912/150‑1 (MQ=255)

                                                                                             gCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCt                                                        <  2:1514466/150‑1 (MQ=255)

                                                                                               tGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGa                                                      >  1:2643175/1‑150 (MQ=255)

                                                                                                  accaccTTGAACCGCGTAGCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATCGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaaa                                                   <  1:4676891/150‑1 (MQ=255)

                                                                                                    caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAaca                                                 >  1:1347639/1‑150 (MQ=255)

                                                                                                     accTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGc                                                               >  1:3565108/1‑135 (MQ=255)

                                                                                                      ccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaaaggga                                               >  1:4930623/1‑146 (MQ=255)

                                                                                                      ccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACAcc                                               >  1:3208393/1‑150 (MQ=255)

                                                                                                       cTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCa                                              >  1:2998077/1‑150 (MQ=255)

                                                                                                          gAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGcc                                           <  1:2558590/150‑1 (MQ=255)

                                                                                                            aCCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCtt                                         <  1:1054929/150‑1 (MQ=255)

                                                                                                            aCCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCtt                                         <  1:1054401/150‑1 (MQ=255)

                                                                                                              cgcgTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGa                                       >  2:4678141/1‑150 (MQ=255)

                                                                                                                          tCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGt                           <  2:655067/150‑1 (MQ=255)

                                                                                                                          tCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGt                           <  1:262717/150‑1 (MQ=255)

                                                                                                                            cAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCa                         >  1:4535666/1‑150 (MQ=255)

                                                                                                                               cGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGAc                      <  1:3099442/150‑1 (MQ=255)

                                                                                                                                     gcgcTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGaaaa                <  2:2186204/150‑1 (MQ=255)

                                                                                                                                     gcgcTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGaaaa                <  2:2580482/150‑1 (MQ=255)

                                                                                                                                     gcgcTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGaaaa                <  2:2186203/150‑1 (MQ=255)

                                                                                                                                     gcgcTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGaaaa                <  1:4190851/150‑1 (MQ=255)

                                                                                                                                         ttCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAttt                                                                                 >  1:5190320/1‑81 (MQ=255)

                                                                                                                                         ttCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTa            >  2:3543852/1‑150 (MQ=255)

                                                                                                                                         ttCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTa            <  1:2424539/150‑1 (MQ=255)

                                                                                                                                          tCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTAt           <  2:2634529/150‑1 (MQ=255)

                                                                                                                                             aGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGcc        <  1:4316139/150‑1 (MQ=255)

                                                                                                                                             aGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGcc        <  2:2025385/150‑1 (MQ=255)

                                                                                                                                              gCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCg       >  2:4825097/1‑150 (MQ=255)

                                                                                                                                               cTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGt      >  1:4295335/1‑150 (MQ=255)

                                                                                                                                                 gCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTaa                                                                            <  2:1766349/78‑1 (MQ=255)

                                                                                                                                                 gCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTaa                                                                            >  1:1766349/1‑78 (MQ=255)

                                                                                                                                                 gCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTaa    >  2:1528507/1‑150 (MQ=255)

                                                                                                                                                  cATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAAc   >  2:4581776/1‑150 (MQ=255)

                                                                                                                                                   aTTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAAcc  <  1:3679331/150‑1 (MQ=255)

GGGCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAACC  >  NC_002947/6113807‑6114103

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