BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,189,961	A→T	7.3%	I166F (ATC→TTC)	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,189,961	0	A	T	7.3%	971.9 / 44.5	298	I166F (ATC→TTC)	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (131/145); new base T (10/12); total (141/157)
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 317 total aligned reads displayed.

TCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCACCGC  >  NC_002947/5189812‑5190110

 cGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCAt                                                                                                                                                      <  2:7097152/150‑1 (MQ=255)

   gCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTcc                                                                                                                                                    <  2:1756072/150‑1 (MQ=255)

   gCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATcc                                                                                                                                                    >  1:6860046/1‑150 (MQ=255)

     tGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATcccc                                                                                                                                                  <  2:6000359/150‑1 (MQ=255)

        gcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGaa                                                                                                                                               <  1:966315/150‑1 (MQ=38)

        gcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                               <  2:740908/150‑1 (MQ=255)

        gcCTGTGAAACAGACCAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGAGACCCTGCTCTACTACCTGCTCACCCAGCGCGACGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGaa                                                                                                                                               <  2:305643/150‑1 (MQ=37)

         ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAag                                                                                                                                              >  1:1013119/1‑150 (MQ=255)

               aaaCCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTcc                                                                                                                                        >  2:7966860/1‑150 (MQ=255)

                   cGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTa                                                                                                                                    <  1:5881047/150‑1 (MQ=255)

                          gCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTa                                                                                                                             >  2:3648914/1‑150 (MQ=255)

                          gCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTa                                                                                                                             >  1:7601231/1‑150 (MQ=255)

                              tGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCAGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCt                                                                                                                         >  1:2096048/1‑150 (MQ=255)

                                aaGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTc                                                                                                                       <  2:1255481/150‑1 (MQ=255)

                                aaGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTc                                                                                                                       <  2:2193406/150‑1 (MQ=255)

                                aaGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTc                                                                                                                       <  1:7991034/150‑1 (MQ=255)

                                 aGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGGTACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTct                                                                                                                      <  2:1534414/150‑2 (MQ=255)

                                 aGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCg                                                                                                                      >  2:5435308/1‑150 (MQ=255)

                                   aTGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGc                                                                                                                    <  1:8489883/150‑1 (MQ=255)

                                      ttGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAg                                                                                                                 <  2:7899857/150‑1 (MQ=255)

                                        ggCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGgc                                                                                                               >  1:6667749/1‑150 (MQ=255)

                                             tCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGcgcgcg                                                                                                          >  1:3542306/1‑150 (MQ=255)

                                               gAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcggc                                                                                                        >  2:6871643/1‑150 (MQ=255)

                                               gAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcggc                                                                                                        <  1:1651539/150‑1 (MQ=255)

                                               gAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcggc                                                                                                        >  1:7260941/1‑150 (MQ=255)

                                               gAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcggc                                                                                                        <  2:7702767/150‑1 (MQ=255)

                                                aGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcggcg                                                                                                       <  1:6721704/150‑1 (MQ=255)

                                                 ggTCCCGCCGCCACGCTCGCTCGAGACCCTGCTCTACGCCCTGCTCACCCAGCGCGACGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGcggcgg                                                                                                      <  2:2502113/150‑1 (MQ=255)

                                                 ggTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGcggcgg                                                                                                      <  2:5808367/150‑1 (MQ=255)

                                                  gTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCggg                                                                                                     >  2:324787/1‑150 (MQ=255)

                                                  gTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCggg                                                                                                     >  1:272325/1‑150 (MQ=255)

                                                   tCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGc                                                                                                    >  1:1127195/1‑150 (MQ=255)

                                                   tCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGc                                                                                                    >  1:9064783/1‑150 (MQ=255)

                                                     ccgccgCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCtt                                                                                                  >  2:5578728/1‑150 (MQ=255)

                                                     ccgccgCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCtt                                                                                                  >  1:297275/1‑150 (MQ=255)

                                                     ccgccgCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCtt                                                                                                  >  1:4303840/1‑150 (MQ=255)

                                                         cgcAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGa                                                                                              <  2:714016/150‑1 (MQ=255)

                                                           cAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGaaa                                                                                            <  1:3190847/150‑1 (MQ=255)

                                                              cgctcgctCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAacac                                                                                         <  1:6263219/150‑1 (MQ=255)

                                                                 tcgctcgATACCCTGCTCTCCTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                                      <  1:5623022/150‑1 (MQ=255)

                                                                 tcgctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                                      >  1:6003242/1‑150 (MQ=255)

                                                                 tcgctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                                      <  1:2241465/150‑1 (MQ=255)

                                                                 tcgctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                                      >  1:1727482/1‑150 (MQ=255)

                                                                   gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAca                                                                                    >  2:8667304/1‑150 (MQ=255)

                                                                   gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAca                                                                                    >  2:8759536/1‑150 (MQ=255)

                                                                     tcgATACCCTGCTCTGCTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcagc                                                                                  >  2:6162537/1‑150 (MQ=255)

                                                                       gATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcagcag                                                                                <  2:8478761/150‑1 (MQ=255)

                                                                                   tactacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGAt                                                                    >  1:3297332/1‑150 (MQ=255)

                                                                                   tactacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGAt                                                                    <  2:2248732/150‑1 (MQ=255)

                                                                                      tacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAg                                                                 >  1:3783565/1‑150 (MQ=255)

                                                                                       acCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGc                                                                >  1:8654770/1‑150 (MQ=255)

                                                                                       acCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGc                                                                <  1:1327939/150‑1 (MQ=255)

                                                                                       acCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGc                                                                <  1:5183246/150‑1 (MQ=255)

                                                                                        ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACagcagc                                                                               >  1:4813226/1‑134 (MQ=255)

                                                                                        ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACagcagc                                                                               <  2:4813226/134‑1 (MQ=255)

                                                                                        ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGcc                                                               <  2:8296854/150‑1 (MQ=255)

                                                                                        ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGcc                                                               >  2:3798852/1‑150 (MQ=255)

                                                                                        ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGcc                                                               <  2:1882225/150‑1 (MQ=255)

                                                                                           gCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCa                                                            <  2:4656709/150‑1 (MQ=255)

                                                                                             tCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATg                                                          >  1:3095895/1‑150 (MQ=255)

                                                                                              cACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGc                                                         >  1:1166114/1‑150 (MQ=255)

                                                                                               aCCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                                        >  1:7014718/1‑150 (MQ=255)

                                                                                                 ccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCca                                                      <  2:1695248/150‑1 (MQ=255)

                                                                                                  cAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcag                                                     >  2:4291303/1‑150 (MQ=255)

                                                                                                   aGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcagc                                                    >  1:4928314/1‑150 (MQ=255)

                                                                                                     cgcgAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcagcag                                                  >  2:657787/1‑150 (MQ=255)

                                                                                                              ggCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCgg                                         >  2:626295/1‑150 (MQ=255)

                                                                                                               gCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGt                                        <  1:5498384/150‑1 (MQ=255)

                                                                                                                 ctgctgGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGc                                      <  2:4121779/150‑1 (MQ=255)

                                                                                                                   gctgGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCa                                    <  1:2146970/150‑1 (MQ=255)

                                                                                                                   gctgGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCa                                    <  1:5745008/150‑1 (MQ=255)

                                                                                                                      ggACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACATGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGc                                 <  1:7698266/150‑1 (MQ=255)

                                                                                                                      ggACTTCGAGGAAGGCGTGGTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGc                                 <  1:4887094/150‑1 (MQ=255)

                                                                                                                       gACTTTGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCg                                <  2:2956790/150‑1 (MQ=255)

                                                                                                                       gACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCg                                <  1:7843436/150‑1 (MQ=255)

                                                                                                                       gACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCg                                <  1:4527491/150‑1 (MQ=255)

                                                                                                                          ttCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                             >  2:1045932/1‑150 (MQ=38)

                                                                                                                          ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCCGGCGGtt                             >  2:8484601/1‑150 (MQ=255)

                                                                                                                          ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                             >  1:3586374/1‑150 (MQ=255)

                                                                                                                          ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                             >  1:4075376/1‑150 (MQ=255)

                                                                                                                          ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                             >  1:821944/1‑150 (MQ=255)

                                                                                                                          ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                             >  1:1873468/1‑150 (MQ=255)

                                                                                                                          ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                             >  1:8729131/1‑150 (MQ=255)

                                                                                                                           tCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTa                            >  1:4669324/1‑150 (MQ=255)

                                                                                                                            cGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAc                           <  2:781647/150‑1 (MQ=255)

                                                                                                                              gggATGGCGTGGCGCAACACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCg                         >  1:5357565/2‑150 (MQ=255)

                                                                                                                                gAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGaa                       <  2:3130913/150‑1 (MQ=255)

                                                                                                                                     cGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCa                  >  2:568476/1‑150 (MQ=255)

                                                                                                                                     cGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCa                  <  1:791431/150‑1 (MQ=255)

                                                                                                                                     cGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCa                  <  1:1376604/150‑1 (MQ=255)

                                                                                                                                     cGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCa                  <  2:3566246/150‑1 (MQ=255)

                                                                                                                                      gtgtTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAg                 <  1:8083037/150‑1 (MQ=255)

                                                                                                                                         ttGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCaa              <  1:3013020/150‑1 (MQ=255)

                                                                                                                                         ttGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCaa              <  2:1458099/150‑1 (MQ=255)

                                                                                                                                             gccacgccaTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTcc          <  2:275752/150‑1 (MQ=255)

                                                                                                                                             gccacgccaTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTcc          <  1:7099541/150‑1 (MQ=255)

                                                                                                                                                acgccaTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATc       >  1:2457946/1‑150 (MQ=255)

                                                                                                                                                 cgccaTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGGAGGCTGAGCGGTGCCAGGCGGTTACGGATGATCAGGTAGTCCATCa      <  2:8961072/150‑1 (MQ=255)

                                                                                                                                                  gccaTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCAc     >  2:4890351/1‑150 (MQ=255)

                                                                                                                                                    caTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCACCg   <  2:6813830/150‑1 (MQ=255)

TCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCCATCACCGC  >  NC_002947/5189812‑5190110

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