BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,992,320	G→A	8.3%	L109L (CTG→TTG)	ycaC‑III ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,992,320	0	G	A	8.3%	1074.4 / 56.5	339	L109L (CTG→TTG)	ycaC‑III	putative hydrolase
Reads supporting (aligned to +/- strand): ref base G (191/119); new base A (18/10); total (210/129)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 347 total aligned reads displayed.

CTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGG  >  NC_002947/5992175‑5992469

cTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAgctgc                                                                                                                                                   >  1:1083953/1‑150 (MQ=255)

 tCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAgctgct                                                                                                                                                  >  2:9089552/1‑150 (MQ=255)

  cATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGctt                                                                                                                                                 <  2:3892677/150‑1 (MQ=255)

   aTGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGcttc                                                                                                                                                <  2:6930867/150‑1 (MQ=255)

    tGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGcttct                                                                                                                                               >  1:3684844/1‑150 (MQ=255)

       ggTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCACACCTCGGTCACAACACCTGCCATGATCAGCTGCTTCTTGc                                                                                                                                            >  1:9059442/1‑150 (MQ=255)

        gTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGcc                                                                                                                                           >  1:246915/1‑150 (MQ=255)

           aTGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGt                                                                                                                                        >  2:1135372/1‑150 (MQ=255)

              cgcgGCATCACGGGTAATTTCATTGAAGGTGCCGGAAGCATCGTTCACCACGAACACCTCAAATTCTTACTCCAAGGACGACAGCGCAGGGAACGGCACGCAAAACTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGAAGGTCGc                                                                                                                                     >  2:2356479/1‑150 (MQ=255)

                        cGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCACCTTCACCGcc                                                                                                                           >  1:4687434/1‑150 (MQ=255)

                          ggTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCtt                                                                                                                         <  2:6218541/150‑1 (MQ=255)

                           gTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTg                                                                                                                        >  1:5700511/1‑150 (MQ=255)

                           gTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTg                                                                                                                        <  1:2684793/150‑1 (MQ=255)

                            tcatttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGa                                                                                                                       >  1:6309798/1‑150 (MQ=255)

                              atttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACg                                                                                                                     >  1:5897324/1‑150 (MQ=255)

                               tttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGa                                                                                                                    >  1:554019/1‑150 (MQ=255)

                                ttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGaa                                                                                                                   >  2:79438/1‑150 (MQ=255)

                                 tcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCACACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGaaa                                                                                                                  >  2:1002944/1‑150 (MQ=255)

                                   attGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATc                                                                                                                >  2:610450/1‑150 (MQ=255)

                                   attGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATc                                                                                                                >  1:6113810/1‑150 (MQ=255)

                                       aaGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCg                                                                                                            >  1:2946662/1‑150 (MQ=255)

                                        aGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGt                                                                                                           <  1:1812758/150‑1 (MQ=255)

                                                gAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAg                                                                                                   <  2:5113590/150‑1 (MQ=255)

                                                gAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAg                                                                                                   <  2:8904280/150‑1 (MQ=255)

                                                    cATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGt                                                                                               >  2:1590083/1‑150 (MQ=255)

                                                    cATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGt                                                                                               >  1:8330152/1‑150 (MQ=255)

                                                     aacGGTCACCACGAACACCACAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGtt                                                                                              <  1:7926273/148‑1 (MQ=255)

                                                      tCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTg                                                                                             >  2:8316988/1‑150 (MQ=255)

                                                       cGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGa                                                                                            <  1:1050486/150‑1 (MQ=255)

                                                         gTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATc                                                                                          >  1:1130253/1‑150 (MQ=255)

                                                          tCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCt                                                                                         >  2:6125922/1‑150 (MQ=255)

                                                           caccacGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTg                                                                                        <  2:5461813/150‑1 (MQ=255)

                                                              cacGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGt                                                                                                           >  1:752774/1‑128 (MQ=255)

                                                                 gAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgg                                                                                  >  1:8275009/1‑150 (MQ=255)

                                                                       cTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggcgggc                                                                            >  2:6740586/1‑150 (MQ=255)

                                                                        tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                           >  2:3347290/1‑150 (MQ=255)

                                                                        tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                           >  2:6401398/1‑150 (MQ=255)

                                                                         cAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCaa                                                                          >  1:5595771/1‑150 (MQ=255)

                                                                         cAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCaa                                                                          >  2:2988416/1‑150 (MQ=255)

                                                                          aaaTTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTAACAACACCTGCAATGATCAGCTGATTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCATGCGATGATCTGGCCAGGGCGGGCAAt                                                                         >  1:2235041/1‑150 (MQ=255)

                                                                                ttcCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACgg                                                                   >  1:653866/1‑150 (MQ=255)

                                                                                 tcctccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGgc                                                                  >  1:7114525/1‑150 (MQ=255)

                                                                                   ctccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGgcgc                                                                >  2:3275994/1‑150 (MQ=255)

                                                                                    tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCa                                                               >  2:4115535/1‑150 (MQ=255)

                                                                                    tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCa                                                               >  1:7029676/1‑150 (MQ=255)

                                                                                    tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCa                                                               >  2:4709461/1‑150 (MQ=255)

                                                                                    tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCa                                                               >  1:3371522/1‑150 (MQ=255)

                                                                                     ccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCAt                                                              >  2:7468252/1‑150 (MQ=255)

                                                                                       aaGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATcc                                                            <  1:1995893/150‑1 (MQ=255)

                                                                                            cGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACAGCGCATCCGGGaa                                                       <  2:2482433/150‑1 (MQ=255)

                                                                                             gACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAAc                                                      >  1:7387101/1‑150 (MQ=255)

                                                                                              aCAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACa                                                     <  1:6644448/150‑1 (MQ=255)

                                                                                              aCAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACa                                                     <  1:151082/150‑1 (MQ=255)

                                                                                               cAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAg                                                    <  1:6394607/150‑1 (MQ=255)

                                                                                                  cgcCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGc                                                 <  1:1323111/150‑1 (MQ=255)

                                                                                                  cgcCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGc                                                 <  1:7457703/150‑1 (MQ=255)

                                                                                                   gcCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGaatgtttct                                                <  1:8838747/150‑8 (MQ=14)

                                                                                                     cGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCttt                                              <  2:4494784/150‑1 (MQ=255)

                                                                                                      gggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTc                                             >  1:2101246/1‑150 (MQ=255)

                                                                                                      gggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTc                                             <  1:4525845/150‑1 (MQ=255)

                                                                                                       ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                            <  2:356101/150‑1 (MQ=255)

                                                                                                       ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                            <  2:8649672/150‑1 (MQ=255)

                                                                                                       ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                            <  2:7759557/150‑1 (MQ=255)

                                                                                                       ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                            <  2:2822733/150‑1 (MQ=255)

                                                                                                       ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                            <  1:6861779/150‑1 (MQ=255)

                                                                                                            gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTcc                                       >  2:6648070/1‑150 (MQ=255)

                                                                                                            gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTcc                                       >  1:758725/1‑150 (MQ=255)

                                                                                                            gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTcc                                       >  2:131231/1‑150 (MQ=255)

                                                                                                              cACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCgg                                     >  2:2667585/1‑150 (MQ=255)

                                                                                                              cACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCgg                                     >  2:37018/1‑150 (MQ=255)

                                                                                                                 gCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCAc                                  <  1:8591899/150‑1 (MQ=9)

                                                                                                                   aaaCCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCa                                >  2:7227516/1‑150 (MQ=255)

                                                                                                                   aaaCCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCa                                >  1:3968973/1‑150 (MQ=255)

                                                                                                                     aCCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGc                              >  1:263800/1‑150 (MQ=255)

                                                                                                                      cttCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCaatg                             <  2:8137616/148‑4 (MQ=9)

                                                                                                                       cTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCgg                            <  2:2057225/150‑1 (MQ=255)

                                                                                                                           gTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGcca                        >  2:7036451/1‑149 (MQ=9)

                                                                                                                           gTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCg                        >  2:7582085/1‑150 (MQ=255)

                                                                                                                            tCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGccat                       >  1:1608516/1‑148 (MQ=9)

                                                                                                                              acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                     >  2:4901561/1‑150 (MQ=255)

                                                                                                                              acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                     <  2:428245/150‑1 (MQ=255)

                                                                                                                               caacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTgg                    >  1:8093054/1‑150 (MQ=255)

                                                                                                                               caacaCCAGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTgg                    <  1:5090160/150‑1 (MQ=255)

                                                                                                                                aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg                   >  1:922333/1‑150 (MQ=255)

                                                                                                                                aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg                   >  2:1449256/1‑150 (MQ=255)

                                                                                                                                 acacCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTcggg                  <  1:387560/150‑1 (MQ=2)

                                                                                                                                     cTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCtt              >  2:1748584/1‑150 (MQ=255)

                                                                                                                                     cTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCtt              <  1:5393684/150‑1 (MQ=255)

                                                                                                                                      tGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTg             <  2:1968011/150‑1 (MQ=255)

                                                                                                                                       gcgaTGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGccgttt            >  2:5650462/4‑146 (MQ=2)

                                                                                                                                       gcgaTGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGccgttt            <  1:4895631/147‑5 (MQ=2)

                                                                                                                                        cAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCt           >  2:8788004/1‑150 (MQ=255)

                                                                                                                                         aaTGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCACTGTACGGCGCATCCGGGAACAGGGCTGTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTc          >  2:9012743/1‑150 (MQ=255)

                                                                                                                                         aaTGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTc          >  1:4640926/1‑150 (MQ=255)

                                                                                                                                              tCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGc     >  1:4000746/1‑150 (MQ=255)

                                                                                                                                               cAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCt    >  2:7310755/1‑150 (MQ=255)

                                                                                                                                                aGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtg   <  2:1611044/150‑1 (MQ=255)

                                                                                                                                                 gctgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtgg  >  1:4382307/1‑150 (MQ=255)

                                                                                                                                                 gctgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtgg  >  1:3060364/1‑150 (MQ=255)

                                                                                                                                                 gctgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtgg  >  2:2448074/1‑150 (MQ=255)

CTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGG  >  NC_002947/5992175‑5992469

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