BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,987,828	G→T	6.1%	G118G (GGG→GGT)	PP_3510 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,987,828	0	G	T	6.1%	980.6 / 36.7	294	G118G (GGG→GGT)	PP_3510	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (130/145); new base T (3/15); total (133/161)
Fisher's exact test for biased strand distribution p-value = 1.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 306 total aligned reads displayed.

ACGCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTG  >  NC_002947/3987680‑3987977

aCGCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGc                                                                                                                                                      >  2:4432492/1‑150 (MQ=255)

  gCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCAt                                                                                                                                                    >  1:1028484/1‑150 (MQ=255)

  gCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCAt                                                                                                                                                    >  2:128880/1‑150 (MQ=255)

  gCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCAt                                                                                                                                                    >  1:3072959/1‑150 (MQ=255)

   cAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCAtt                                                                                                                                                   <  2:2065359/150‑1 (MQ=255)

       ggCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGAcc                                                                                                                                               >  2:2033999/1‑150 (MQ=255)

           cGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGc                                                                                                                                           <  2:6683128/150‑1 (MQ=255)

            gTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCCGCAGGGTCAATGCCCCTGca                                                                                                                                          <  2:2697839/150‑2 (MQ=255)

             tATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcag                                                                                                                                         <  2:676459/150‑3 (MQ=255)

             tATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcag                                                                                                                                         <  1:755177/150‑3 (MQ=255)

             tATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGcag                                                                                                                                         <  2:5745646/150‑3 (MQ=255)

              aTGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcagg                                                                                                                                        <  2:5622798/150‑4 (MQ=255)

               tGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGtt                                                                                                                                       <  2:4612665/150‑1 (MQ=255)

               tGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCCGCAGGGTCAATGCCCCTGcaggc                                                                                                                                       <  2:4759766/150‑5 (MQ=255)

                gCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGttt                                                                                                                                      <  2:4239450/150‑1 (MQ=255)

                gCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGttt                                                                                                                                      <  2:6006093/150‑1 (MQ=255)

                  tGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGa                                                                                                                                    >  1:6808564/1‑150 (MQ=255)

                  tGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGa                                                                                                                                    <  1:6705922/150‑1 (MQ=255)

                   ggCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGAc                                                                                                                                   <  2:5903118/150‑1 (MQ=255)

                    gCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACg                                                                                                                                  <  1:879669/150‑1 (MQ=255)

                      aCTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGAt                                                                                                                                <  1:540231/150‑1 (MQ=255)

                       cTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATa                                                                                                                               <  2:615238/150‑1 (MQ=255)

                        tGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAg                                                                                                                              >  1:2676125/1‑150 (MQ=255)

                        tGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAg                                                                                                                              >  2:1024445/1‑150 (MQ=255)

                          aCCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTa                                                                                                                            >  2:5166844/1‑150 (MQ=255)

                          aCCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTa                                                                                                                            >  2:1957591/1‑150 (MQ=255)

                           ccTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTAt                                                                                                                           <  1:5281728/150‑1 (MQ=255)

                            cTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATc                                                                                                                          <  2:6189663/150‑1 (MQ=255)

                             tGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCg                                                                                                                         >  1:1818806/1‑150 (MQ=255)

                                 aGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTc                                                                                                                     >  1:1660748/1‑150 (MQ=255)

                                 aGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTc                                                                                                                     <  2:464200/150‑1 (MQ=255)

                                     aTTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCAc                                                                                                                 >  2:2273537/1‑150 (MQ=255)

                                        ggACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCt                                                                                                              <  2:6276036/150‑1 (MQ=255)

                                         gACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTa                                                                                                             >  2:5280881/1‑150 (MQ=255)

                                         gACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTa                                                                                                             >  1:2420919/1‑150 (MQ=255)

                                          aCCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTAt                                                                                                            <  1:2038420/150‑1 (MQ=255)

                                           ccACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATc                                                                                                           >  2:4382007/1‑150 (MQ=255)

                                             aCGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCaa                                                                                                         >  2:1751938/1‑150 (MQ=255)

                                             aCGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCaa                                                                                                         <  1:4343804/150‑1 (MQ=255)

                                              cGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAAt                                                                                                        <  2:6855825/150‑1 (MQ=255)

                                               gCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGAGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAAtt                                                                                                       >  1:3711134/1‑150 (MQ=255)

                                                ccGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGTGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTc                                                                                                      <  1:5314785/150‑1 (MQ=255)

                                                 cGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCg                                                                                                     >  2:3873652/1‑150 (MQ=255)

                                                 cGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCg                                                                                                     <  2:3057209/150‑1 (MQ=255)

                                                  ggAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGt                                                                                                    <  1:1482615/150‑1 (MQ=255)

                                                        cAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTgg                                                                                              <  1:5186133/150‑1 (MQ=255)

                                                               aCCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTcc                                                                                       <  2:2751518/150‑1 (MQ=255)

                                                               aCCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTcc                                                                                       <  1:1969647/150‑1 (MQ=255)

                                                                 cATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATc                                                                                                                          >  1:2155051/1‑113 (MQ=255)

                                                                  aTCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGa                                                                                    >  2:3999207/1‑150 (MQ=255)

                                                                   tCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGAt                                                                                   <  2:853771/150‑1 (MQ=255)

                                                                           ggCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAg                                                                           <  2:6566991/150‑1 (MQ=255)

                                                                             cATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGTGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGa                                                                         >  1:4226326/1‑150 (MQ=255)

                                                                              aTGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGAt                                                                        >  1:2937867/1‑150 (MQ=255)

                                                                              aTGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGAt                                                                        <  2:4894654/150‑1 (MQ=255)

                                                                              aTGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGAt                                                                        <  2:2460441/150‑1 (MQ=255)

                                                                              aTGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGAt                                                                        >  2:3435979/1‑150 (MQ=255)

                                                                                 tACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTa                                                                     <  1:594079/150‑1 (MQ=255)

                                                                                    gCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTa                                                                  >  2:1499931/1‑150 (MQ=255)

                                                                                     cTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTaa                                                                 <  2:1598918/150‑1 (MQ=255)

                                                                                      tATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTc                                                                                        <  1:987775/126‑1 (MQ=255)

                                                                                        tGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGt                                                              <  1:250336/150‑1 (MQ=255)

                                                                                        tGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGt                                                              <  1:3794418/150‑1 (MQ=255)

                                                                                        tGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGACATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGt                                                              <  2:5499613/150‑1 (MQ=255)

                                                                                         gACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTg                                                             <  2:4353263/150‑1 (MQ=255)

                                                                                         gACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTg                                                             <  1:3280111/150‑1 (MQ=255)

                                                                                          aCAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTgg                                                            <  1:6387856/150‑1 (MQ=255)

                                                                                               gAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATgg                                                       >  1:4132563/1‑150 (MQ=255)

                                                                                                     gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGATGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCa                                                 >  1:403088/1‑150 (MQ=255)

                                                                                                     gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCa                                                 >  1:1382127/1‑150 (MQ=255)

                                                                                                       cTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCata                                               <  1:4838420/150‑1 (MQ=255)

                                                                                                        tgtgCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCatat                                              >  2:4172292/1‑150 (MQ=255)

                                                                                                        tgtgCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCatat                                              >  2:4664394/1‑150 (MQ=255)

                                                                                                          tgCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATaa                                            <  1:6640125/150‑1 (MQ=255)

                                                                                                            cACTGGTGCATGTCGGGGACTTGGGAAAACGCCTGCAGGGGCATTGACCCTGCTGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAca                                          <  2:1881848/150‑1 (MQ=255)

                                                                                                            cACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAca                                          <  1:121165/150‑1 (MQ=255)

                                                                                                            cACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAca                                          >  2:6474186/1‑150 (MQ=255)

                                                                                                             aCGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcag                                         <  1:4785707/150‑1 (MQ=255)

                                                                                                             aCGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcag                                         <  1:5853870/150‑1 (MQ=255)

                                                                                                               gggTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcagca                                       >  2:4935894/1‑150 (MQ=255)

                                                                                                                   gCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGccc                                   <  2:1214032/150‑1 (MQ=255)

                                                                                                                        gCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCg                              >  1:4526825/1‑150 (MQ=255)

                                                                                                                           gggACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATc                           <  2:1066600/150‑1 (MQ=255)

                                                                                                                            ggACTGGCGAAAACGCCTGCAGGGGCATAGACCCGGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCAAAGTATGTATTAAGGTGGCATGTTGCGCATATAACAGCAGCCCGAGCGATCg                          <  1:3633275/150‑1 (MQ=255)

                                                                                                                             gACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCgg                         >  1:6450223/1‑150 (MQ=255)

                                                                                                                              aCTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGATCGCTCggg                        <  1:5127251/150‑1 (MQ=255)

                                                                                                                                tGGGGAAAACGCCTGCAGGGGCATCGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCAATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCa                      >  1:4192181/1‑150 (MQ=255)

                                                                                                                                 ggggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCAc                     <  1:4756098/150‑1 (MQ=255)

                                                                                                                                   ggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACtt                   <  2:4132671/150‑1 (MQ=255)

                                                                                                                                    gAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTgcgc                                                  <  2:3701537/118‑1 (MQ=255)

                                                                                                                                    gAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTg                  >  2:3310516/1‑150 (MQ=255)

                                                                                                                                      aaaCGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGATCGCTCGGGctgctgtt                <  1:1575733/150‑8 (MQ=255)

                                                                                                                                         aaaccgCAGGGTCAATGCCCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAATGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGa             >  1:2268111/6‑150 (MQ=255)

                                                                                                                                          gCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGAc            >  1:570427/1‑150 (MQ=255)

                                                                                                                                             tGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAAt         <  1:547667/150‑1 (MQ=255)

                                                                                                                                                 ggggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAgg     <  1:3519894/150‑1 (MQ=255)

                                                                                                                                                 ggggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAgg     <  2:3948603/150‑1 (MQ=255)

                                                                                                                                                  gggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAggg    >  2:1421608/1‑150 (MQ=255)

                                                                                                                                                  gggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAggg    >  2:4084856/1‑150 (MQ=255)

                                                                                                                                                  gggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAggg    >  2:4971231/1‑150 (MQ=255)

ACGCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAGGGTG  >  NC_002947/3987680‑3987977

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