BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,992,317	T→G	8.4%	I110L (ATC→CTC)	ycaC‑III ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,992,317	0	T	G	8.4%	1131.5 / 66.2	341	I110L (ATC→CTC)	ycaC‑III	putative hydrolase
Reads supporting (aligned to +/- strand): ref base T (195/117); new base G (19/10); total (214/127)
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.93e-01

Only 100 of 343 total aligned reads displayed.

CGCCTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGC  >  NC_002947/5992168‑5992466

cGCCTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGAt                                                                                                                                                       <  1:8656620/150‑1 (MQ=255)

  ccTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCa                                                                                                                                                     <  1:2555650/150‑1 (MQ=255)

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

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

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

            gCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGcttctt                                                                                                                                           <  1:5078916/150‑1 (MQ=255)

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

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

                tCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCg                                                                                                                                       <  1:2650133/150‑1 (MQ=255)

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

                   tGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTc                                                                                                                                    >  1:4539442/1‑150 (MQ=255)

                     cgcgGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGc                                                                                                                                  <  1:2218068/150‑1 (MQ=255)

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

                      gcgGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGcc                                                                                                                                 >  2:6130871/1‑150 (MQ=255)

                      gcgGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGcc                                                                                                                                 >  2:3411524/1‑150 (MQ=255)

                       cgGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCt                                                                                                                                >  2:5159536/1‑150 (MQ=255)

                           aTCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCAc                                                                                                                            >  2:3678219/1‑150 (MQ=255)

                               cGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGcc                                                                                                                        >  2:1680270/1‑150 (MQ=255)

                                gggTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCt                                                                                                                       <  2:3122809/150‑1 (MQ=255)

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

                                  gTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGACCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTg                                                                                                                     <  2:2254427/150‑1 (MQ=255)

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

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

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

                                       ttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGaa                                                                                                                <  1:5553576/150‑1 (MQ=255)

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

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

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

                                             gAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACTCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTc                                                                                                          <  2:7306873/150‑1 (MQ=255)

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

                                                      ggAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGACACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCATCGTTGTCCCa                                                                                                 >  1:3010092/1‑150 (MQ=255)

                                                       gAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAg                                                                                                >  1:570682/1‑150 (MQ=255)

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

                                                           cATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGt                                                                                            >  2:1590083/1‑150 (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:4122773/1‑150 (MQ=255)

                                                                     cacGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGcc                                                                                  >  2:4301915/1‑150 (MQ=255)

                                                                       cGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAg                                                                                <  1:4008846/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                    ttcttcCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTc                                                                                                          >  1:7345971/1‑111 (MQ=255)

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

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

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

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

                                                                                           tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCaa                                                                       <  2:6518863/139‑1 (MQ=255)

                                                                                           tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCaa                                                                       >  1:6518863/1‑139 (MQ=255)

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

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

                                                                                                 gCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCCCTGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTGGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCggg                                                      >  2:285247/1‑150 (MQ=255)

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

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

                                                                                                       aGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGGCGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAgg                                                <  2:4833423/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                                                   gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTcc                                    <  1:7701029/150‑1 (MQ=255)

                                                                                                                    ccACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCg                                   >  1:751069/1‑150 (MQ=255)

                                                                                                                     cACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCgg                                  <  1:6800065/150‑1 (MQ=255)

                                                                                                                       cGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCa                                <  1:6173304/150‑1 (MQ=255)

                                                                                                                        gCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCAc                               >  2:4611421/1‑150 (MQ=9)

                                                                                                                        gCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTATGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCAc                               >  1:54893/1‑150 (MQ=255)

                                                                                                                        gCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCAc                               <  2:3764711/150‑1 (MQ=255)

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

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

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

                                                                                                                               tCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCggg                        >  1:3901745/1‑150 (MQ=255)

                                                                                                                               tCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCggg                        >  2:2711738/1‑150 (MQ=255)

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

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

                                                                                                                                  gTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGAGCTGGCCAGTGCGGGCAATGTACGGCGCATCCGGGACCAGGGCTTTCAGTTCCGGCACCAGCGGGCCg                     >  2:7510511/1‑150 (MQ=255)

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

                                                                                                                                     acaacaCCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCAttc                  >  1:830551/1‑149 (MQ=2)

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

                                                                                                                                       aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg                <  1:7073403/150‑1 (MQ=255)

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

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

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

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

                                                                                                                                          acCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCCCCAACGGCCCGCTGGTGcc             <  1:5278437/150‑1 (MQ=18)

                                                                                                                                            cTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGccgt           >  2:1780673/1‑148 (MQ=2)

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

                                                                                                                                              gCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGc         <  1:4477822/150‑1 (MQ=255)

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

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

                                                                                                                                                   gATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGaa    >  1:4644539/1‑150 (MQ=255)

CGCCTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGC  >  NC_002947/5992168‑5992466

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