BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,245,566	A→T	6.0%	D44V (GAC→GTC)	PP_4618 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,245,566	0	A	T	6.0%	1040.9 / 36.5	315	D44V (GAC→GTC)	PP_4618	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (141/155); new base T (4/15); total (145/170)
Fisher's exact test for biased strand distribution p-value = 3.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 318 total aligned reads displayed.

ACGTGTGTGGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACTGC  >  NC_002947/5245417‑5245715

 cGTGTGTGGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGAc                                                                                                                                                      <  2:3235252/150‑1 (MQ=255)

  gtgtgtGGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGAcc                                                                                                                                                     <  2:2134180/150‑1 (MQ=255)

      gtgGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCa                                                                                                                                                 <  1:1550112/150‑1 (MQ=255)

       tgGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAg                                                                                                                                                >  2:8000721/1‑150 (MQ=255)

        ggAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgc                                                                                                                                               >  2:4424865/1‑150 (MQ=255)

          agaTACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcgc                                                                                                                                             >  1:8533139/1‑150 (MQ=255)

           gaTACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAGGTCGTGCAgcgca                                                                                                                                            <  2:2395909/150‑2 (MQ=255)

            aTACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAGGTCGTGCAgcgcag                                                                                                                                           <  2:2709827/150‑3 (MQ=255)

              aCCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCgaga                                                                                                                                         >  1:7253228/1‑150 (MQ=255)

              aCCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCgaga                                                                                                                                         <  1:2745001/150‑1 (MQ=255)

               ccccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGAt                                                                                                                                        >  1:5355426/1‑150 (MQ=255)

               ccccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGAt                                                                                                                                        >  2:4059431/1‑150 (MQ=255)

                 ccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATaa                                                                                                                                      >  1:393844/1‑150 (MQ=255)

                   aTGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAgc                                                                                                                                    <  1:3431010/150‑1 (MQ=255)

                   aTGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAgc                                                                                                                                    <  1:6975496/150‑1 (MQ=255)

                   aTGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGAGTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAgc                                                                                                                                    <  1:1496418/150‑1 (MQ=255)

                        gCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCaa                                                                                                                               <  1:6348258/150‑1 (MQ=255)

                        gCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCaa                                                                                                                               >  2:3709639/1‑150 (MQ=255)

                         ccACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAAc                                                                                                                              >  2:1672589/1‑150 (MQ=255)

                         ccACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAAc                                                                                                                              <  1:5971078/150‑1 (MQ=255)

                         ccACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAAc                                                                                                                              >  1:7293375/1‑150 (MQ=255)

                           aCGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCt                                                                                                                            >  1:1578326/1‑150 (MQ=255)

                           aCGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCt                                                                                                                            >  2:4790566/1‑150 (MQ=255)

                              tCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGcc                                                                                                                         >  2:2818506/1‑150 (MQ=255)

                                   cctgcctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGc                                                                                                                    <  1:4512759/150‑1 (MQ=255)

                                   cctgcctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGc                                                                                                                    >  2:7240128/1‑150 (MQ=255)

                                      gcctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTc                                                                                                                 <  1:3861870/150‑1 (MQ=255)

                                       cctggcctgcctggAACGCGAGCCTGTCGCCCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTcc                                                                                                                <  2:5651704/150‑1 (MQ=255)

                                        ctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCa                                                                                                               <  2:1358696/150‑1 (MQ=255)

                                         tggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCAc                                                                                                              >  1:114485/1‑150 (MQ=255)

                                            cctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTAt                                                                                                           <  2:7110303/150‑1 (MQ=255)

                                              tgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGc                                                                                                         <  1:5013334/150‑1 (MQ=255)

                                              tgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGc                                                                                                         <  2:8177331/150‑1 (MQ=255)

                                              tgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGc                                                                                                         <  2:4525732/150‑1 (MQ=255)

                                                 ctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGaa                                                                                                      <  2:8374943/150‑1 (MQ=255)

                                                  tggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAAc                                                                                                     >  2:3922262/1‑150 (MQ=255)

                                                   ggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACt                                                                                                    >  2:2314229/1‑150 (MQ=255)

                                                    gAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTg                                                                                                   <  2:3308794/150‑1 (MQ=255)

                                                      aCGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGc                                                                                                 <  2:3643919/150‑1 (MQ=255)

                                                        gcgAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGccc                                                                                               <  1:6912493/150‑1 (MQ=255)

                                                         cgAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCa                                                                                              <  2:6952354/150‑1 (MQ=255)

                                                            gCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGcc                                                                                           <  2:7492234/150‑1 (MQ=255)

                                                               tGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCAct                                                                                        >  1:6875683/1‑150 (MQ=255)

                                                                 tCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCActgc                                                                                      <  2:2998775/150‑1 (MQ=255)

                                                                  cGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCActgct                                                                                     <  1:3642865/150‑1 (MQ=255)

                                                                   gACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCActgctg                                                                                    >  2:2536691/1‑150 (MQ=255)

                                                                       tgctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGcc                                                                                >  2:1991119/1‑150 (MQ=255)

                                                                       tgctgGAAGCCGCACTGTGGAGTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGcc                                                                                >  2:8700151/1‑150 (MQ=255)

                                                                        gctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCa                                                                               <  2:6218855/150‑1 (MQ=255)

                                                                        gctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCa                                                                               <  2:6628340/150‑1 (MQ=255)

                                                                          tgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGc                                                                             >  1:4183233/1‑150 (MQ=255)

                                                                          tgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGc                                                                             <  2:114485/150‑1 (MQ=255)

                                                                           ggAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCt                                                                            >  2:5214773/1‑150 (MQ=255)

                                                                             aaGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGa                                                                          <  2:7306858/150‑1 (MQ=255)

                                                                              aGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGaa                                                                         >  2:2605777/1‑150 (MQ=255)

                                                                               gCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAAt                                                                        <  2:3377943/150‑1 (MQ=255)

                                                                               gCCGCACGGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAAt                                                                        >  2:2447126/1‑150 (MQ=255)

                                                                                     cTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTg                                                                  <  2:6025130/150‑1 (MQ=255)

                                                                                      tgtgGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTgg                                                                 <  1:5283703/150‑1 (MQ=255)

                                                                                       gtgGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTggg                                                                <  1:231488/150‑1 (MQ=255)

                                                                                          gATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCtt                                                             >  1:7666137/1‑150 (MQ=255)

                                                                                           aTTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCttt                                                            <  2:4985891/150‑1 (MQ=255)

                                                                                            ttGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTc                                                           <  2:8533139/150‑1 (MQ=255)

                                                                                            ttGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTc                                                           <  1:8685521/150‑1 (MQ=255)

                                                                                             tgctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTca                                                          <  1:2147723/150‑1 (MQ=255)

                                                                                             tgctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTca                                                          <  1:1005356/150‑1 (MQ=255)

                                                                                                tgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagca                                                       >  2:8231072/1‑150 (MQ=255)

                                                                                                tgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagca                                                       <  1:1672589/150‑1 (MQ=255)

                                                                                                tgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagca                                                       >  1:8091845/1‑150 (MQ=255)

                                                                                                tgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagca                                                       <  2:1323097/150‑1 (MQ=255)

                                                                                                 gcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagcag                                                      <  2:3516295/150‑1 (MQ=255)

                                                                                                   cGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGa                                                    >  1:926708/1‑150 (MQ=255)

                                                                                                    gAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGc                                                                                                                          <  2:7715200/79‑1 (MQ=255)

                                                                                                    gAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGAt                                                   >  1:2498866/1‑150 (MQ=255)

                                                                                                     aGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTGTCAGCAGGATg                                                  >  2:5072758/1‑150 (MQ=255)

                                                                                                          gACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAatat                                             >  1:5429587/1‑150 (MQ=255)

                                                                                                           aCCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATc                                            <  1:4924451/150‑1 (MQ=255)

                                                                                                               tGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCAccc                                        >  1:748584/1‑150 (MQ=255)

                                                                                                                      gAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTgcgc                                 >  2:2285892/1‑150 (MQ=255)

                                                                                                                       aGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGcc                                >  1:7728763/1‑150 (MQ=255)

                                                                                                                       aGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGcc                                >  1:3058951/1‑150 (MQ=255)

                                                                                                                        gCCCGCTGCCATCCTTGCTGCGCTGCACGACCTTCAGCGCGAGATAAGCGACAACCTGCCTATGCTTCCACTATGCTAACTGGCCCAGCCTCTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGccc                               <  1:4353239/150‑1 (MQ=255)

                                                                                                                          ccGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGAGGAATATCACCCCTTGCGccccc                             >  2:8590675/1‑150 (MQ=255)

                                                                                                                                 cAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATgcggc                      <  1:4467729/150‑1 (MQ=255)

                                                                                                                                  aGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATgcggcg                     >  1:3108301/1‑150 (MQ=255)

                                                                                                                                   gCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCga                    <  1:4027453/150‑1 (MQ=255)

                                                                                                                                   gCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCga                    <  1:2214311/150‑1 (MQ=255)

                                                                                                                                   gCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCga                    >  1:403308/1‑150 (MQ=255)

                                                                                                                                    ccTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCgat                   <  1:3354043/150‑1 (MQ=255)

                                                                                                                                    ccTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCgat                   >  1:7873506/1‑150 (MQ=255)

                                                                                                                                      ttGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCgatga                 >  1:2948151/1‑150 (MQ=255)

                                                                                                                                        gctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGatgatg               <  1:6817577/150‑1 (MQ=255)

                                                                                                                                        gctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGatgatg               <  2:4466895/150‑1 (MQ=255)

                                                                                                                                        gctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGatgatg               <  1:8397920/150‑1 (MQ=255)

                                                                                                                                        ctcgcgctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGatgatg               >  2:5496646/4‑150 (MQ=255)

                                                                                                                                         tcgcgctgcAGGTCGTGCAGCGCAAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATgg              >  2:1834996/3‑150 (MQ=255)

                                                                                                                                         ctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATgg              <  1:8000721/150‑1 (MQ=255)

                                                                                                                                          tgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGa             >  1:3189086/1‑150 (MQ=255)

                                                                                                                                                    gACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACTg   >  1:5491024/1‑150 (MQ=255)

                                                                                                                                                    gACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACTg   >  2:1562915/1‑150 (MQ=255)

ACGTGTGTGGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACTGC  >  NC_002947/5245417‑5245715

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