BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,190,316	T→C	6.0%	R82R (CGT→CGC) ‡	PP_3684 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,190,316	0	T	C	6.0%	417.5 / 8.9	132	R82R (CGT→CGC) ‡	PP_3684	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (59/65); new base C (4/4); total (63/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

GTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATA  >  NC_002947/4190167‑4190461

gTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGt                                                                                                                                                   >  1:549508/1‑150 (MQ=255)

gTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGt                                                                                                                                                   <  2:1258815/150‑1 (MQ=255)

gTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGt                                                                                                                                                   <  1:4261450/150‑1 (MQ=255)

  cTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGa                                                                                                                                                 <  1:1105616/150‑1 (MQ=255)

       cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctc                                                                                                                                            >  1:3549016/1‑150 (MQ=255)

       cAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctc                                                                                                                                            >  2:1766932/1‑150 (MQ=255)

        aGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGctct                                                                                                                                           >  2:1573660/1‑150 (MQ=255)

         gACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCtctc                                                                                                                                          >  1:1921134/1‑150 (MQ=255)

         gACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCtctc                                                                                                                                          <  2:1921134/150‑1 (MQ=255)

           ccGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCaa                                                                                                                                        <  1:1453831/150‑1 (MQ=255)

           ccGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCaa                                                                                                                                        <  2:2884453/150‑1 (MQ=255)

           ccGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCaa                                                                                                                                        <  1:4606810/150‑1 (MQ=255)

           ccGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCaa                                                                                                                                        <  1:1805206/150‑1 (MQ=255)

            cGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAg                                                                                                                                       <  2:5240666/150‑1 (MQ=255)

            cGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAg                                                                                                                                       <  1:2239353/150‑1 (MQ=255)

            cGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAg                                                                                                                                       >  1:3929139/1‑150 (MQ=255)

             gAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCTCACGGCTCTCAAgg                                                                                                                                      <  2:3354849/150‑2 (MQ=255)

               agtccAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTTACCGAGGCGGGATGTCGATCGACCGTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCAc                                                                                                                                    >  2:217359/4‑150 (MQ=255)

                  ccAGACGCAGAGATACGGCATGAGCGCCACTTCCTCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCCCTGACCGCGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTg                                                                                                                                 <  1:3780600/150‑1 (MQ=255)

                  ccAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTg                                                                                                                                 <  2:3544765/150‑1 (MQ=255)

                   cAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTgg                                                                                                                                >  1:2263814/1‑150 (MQ=255)

                        gCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCGTGAGGGTCCCTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGAc                                                                                                                           >  2:689400/1‑150 (MQ=255)

                          agagATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGAccc                                                                                                                         >  2:554538/1‑150 (MQ=255)

                          agagATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGAccc                                                                                                                         <  1:1573660/150‑1 (MQ=255)

                                aaGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCa                                                                                                                   <  1:3848622/150‑1 (MQ=255)

                                  ggCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATc                                                                                                                 <  1:3104592/150‑1 (MQ=255)

                                  ggCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATc                                                                                                                 >  1:2068312/1‑150 (MQ=255)

                                  ggCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATc                                                                                                                 <  2:2973730/150‑1 (MQ=255)

                                   gCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCg                                                                                                                >  1:2678908/1‑150 (MQ=255)

                                       gAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGaaa                                                                                                            >  2:46702/1‑150 (MQ=255)

                                       gAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGaaa                                                                                                            <  2:5230468/150‑1 (MQ=255)

                                       gAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGaaa                                                                                                            <  1:1008253/150‑1 (MQ=255)

                                        aGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAg                                                                                                           <  2:1936282/150‑1 (MQ=255)

                                        aGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAg                                                                                                           >  1:5181740/1‑150 (MQ=255)

                                            ccAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGaa                                                                                                       <  2:2450593/150‑1 (MQ=255)

                                            ccAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGaa                                                                                                       >  2:632176/1‑150 (MQ=255)

                                            ccAAGTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGaa                                                                                                       <  2:2930816/150‑1 (MQ=255)

                                             cAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAAc                                                                                                      <  1:2482952/150‑1 (MQ=255)

                                              aaTTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACa                                                                                                     >  2:79488/1‑150 (MQ=255)

                                               aTTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACaa                                                                                                    >  2:146480/1‑150 (MQ=255)

                                                      gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGa                                                                                             <  2:1908070/150‑1 (MQ=255)

                                                      gAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGa                                                                                             >  1:4640092/1‑150 (MQ=255)

                                                      gAGGCGGTGAACAGCCTCAAAACCTGGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGa                                                                                             <  2:2028015/150‑1 (MQ=255)

                                                            gTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCgctg                                                                                       >  2:2152937/1‑150 (MQ=255)

                                                             tGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCgctgc                                                                                      <  1:1226034/150‑1 (MQ=255)

                                                             tGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCgctgc                                                                                      <  2:4189477/150‑1 (MQ=255)

                                                               aaCAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCtt                                                                                    >  2:504392/1‑150 (MQ=255)

                                                                aCAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTtg                                                                                   >  1:2051935/1‑150 (MQ=255)

                                                                     cTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGc                                                                                                          >  1:3629073/1‑122 (MQ=255)

                                                                     cTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGc                                                                                                          <  2:3629073/122‑1 (MQ=255)

                                                                     cTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGtt                                                                              >  1:2113505/1‑150 (MQ=255)

                                                                      tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCg                                                                                                                >  2:1941800/1‑115 (MQ=255)

                                                                      tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGttt                                                                             >  1:1940342/1‑150 (MQ=255)

                                                                      tCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGttt                                                                             >  1:4668264/1‑150 (MQ=255)

                                                                       cAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTc                                                                            >  1:982744/1‑150 (MQ=255)

                                                                        aaaaCCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGaaa                                                                                                            >  2:3282524/1‑117 (MQ=255)

                                                                        aaaaCCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCa                                                                           >  2:3039484/1‑150 (MQ=255)

                                                                          aaCCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGa                                                                         <  2:3688925/150‑1 (MQ=255)

                                                                           aCCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGAt                                                                        <  1:3036308/150‑1 (MQ=255)

                                                                             cTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGAtct                                                                      >  1:2896470/1‑150 (MQ=255)

                                                                              ttGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGAtctc                                                                     >  2:2374883/1‑150 (MQ=255)

                                                                                gAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCaa                                                                   <  2:4595262/150‑1 (MQ=255)

                                                                                gAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCaa                                                                   <  2:5253155/150‑1 (MQ=255)

                                                                                        cTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCaa                                                                   >  2:2919624/1‑142 (MQ=255)

                                                                                          gACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTaaa                                                         <  1:2003814/150‑1 (MQ=255)

                                                                                             cGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTg                                                      >  2:1830646/1‑150 (MQ=255)

                                                                                                 gCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGc                                                  <  1:3913478/150‑1 (MQ=255)

                                                                                                   gggATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAgc                                                                                           >  2:756144/1‑107 (MQ=255)

                                                                                                     gATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACt                                              >  2:1774474/1‑150 (MQ=255)

                                                                                                        gTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCa                                           >  1:3786656/1‑150 (MQ=255)

                                                                                                          cgatcgaCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAAc                                         <  2:448914/150‑1 (MQ=255)

                                                                                                                aCCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTa                                   >  1:3670668/1‑150 (MQ=255)

                                                                                                                 cccTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGAGCTGCTTGTGGTTTCAGATCTAAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTAt                                  >  1:4398359/1‑150 (MQ=255)

                                                                                                                    tGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGc                                      <  2:2343527/143‑1 (MQ=255)

                                                                                                                    tGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTAtttt                               >  1:5322283/1‑150 (MQ=255)

                                                                                                                        gagaTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGc                           <  2:1906390/150‑1 (MQ=255)

                                                                                                                        gagaTACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGc                           <  2:491370/150‑1 (MQ=255)

                                                                                                                             aCGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCttt                      <  1:46702/150‑1 (MQ=255)

                                                                                                                              cGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTg                     <  1:1092354/150‑1 (MQ=255)

                                                                                                                               gtgatgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTgc                    <  1:217359/150‑1 (MQ=255)

                                                                                                                                  atgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCt                 >  1:2278594/1‑150 (MQ=255)

                                                                                                                                   tgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTc                <  2:3809088/150‑1 (MQ=255)

                                                                                                                                   tgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTc                <  1:2722729/150‑1 (MQ=255)

                                                                                                                                   tgtgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTc                <  2:4367356/150‑1 (MQ=255)

                                                                                                                                     tgatCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGACAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGa              >  2:126227/1‑150 (MQ=255)

                                                                                                                                     aagtgcTTGAGAGCCTCACGGCTCTCAAGGACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGa              >  1:4603273/6‑150 (MQ=255)

                                                                                                                                          cTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACt         >  1:4531373/1‑150 (MQ=255)

                                                                                                                                          cTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACt         <  2:3929139/150‑1 (MQ=255)

                                                                                                                                           ttGAGAGCCTCACGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTc        >  1:4279192/1‑150 (MQ=255)

                                                                                                                                           ttGAGAGCCTCACGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTc        >  1:2003637/1‑150 (MQ=255)

                                                                                                                                           ttGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTc        >  2:1397280/1‑150 (MQ=255)

                                                                                                                                           ttGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTc        >  1:322432/1‑150 (MQ=255)

                                                                                                                                            tGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTcc       <  1:554538/150‑1 (MQ=255)

                                                                                                                                              agagCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGt     <  1:689400/150‑1 (MQ=255)

                                                                                                                                               gagaCGTGAGGCTCTCTAGGACTTGGTCGACCCGGCGCATCGGATAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCATCTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGta    <  1:4589989/150‑1 (MQ=255)

                                                                                                                                               gagCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGta    <  1:4937993/150‑1 (MQ=255)

                                                                                                                                               gagCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGta    >  1:2503526/1‑150 (MQ=255)

                                                                                                                                                agCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACt         >  1:3950894/1‑144 (MQ=255)

                                                                                                                                                agCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACt         <  2:3950894/144‑1 (MQ=255)

                                                                                                                                                 gCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGtata  <  1:5138838/150‑1 (MQ=255)

GTCTTCCCAGACCGAAGACCAGAAGCAGAGATAAGGCATGAGCGCCAATTCATCGAGGCGGTGAACAGCCTCAAAACCTTGAGGGTCACTGACCGAGGCGGGATGTCGATCGACCCTGAAGAGATACGTGATGTGATCCTTGAGAGCCGTGAGGCTCTCAAGCACTTGGTCGACCCGGCGCATCGGAAAGCGAACAAAACCTGAGCGCTGCTTGTGGTTTCAGATCTCAAGCCCTCTAAACTGCAGCGACTTCAACTGCGTATTTTGTGCCCTTTGCGCTCGACAACTCCGTATA  >  NC_002947/4190167‑4190461

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