BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,319,178	C→G	5.5%	P608P (CCG→CCC)	treSA ←	trehalose synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,319,178	0	C	G	5.5%	502.4 / 15.0	146	P608P (CCG→CCC)	treSA	trehalose synthase A
Reads supporting (aligned to +/- strand): ref base C (62/76); new base G (5/3); total (67/79)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 147 total aligned reads displayed.

CGACCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCACGCGG  >  NC_002947/3319030‑3319327

  aCCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCggg                                                                                                                                                    <  1:3929967/150‑1 (MQ=255)

  aCCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCggg                                                                                                                                                    <  2:2742138/150‑1 (MQ=255)

   ccGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGCCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgc                                                                                                                                                   >  1:964976/1‑150 (MQ=255)

    cGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgct                                                                                                                                                  <  1:1495928/150‑1 (MQ=255)

     gTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctg                                                                                                                                                 <  1:1673104/150‑1 (MQ=255)

         cGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgctgg                                                                                                                                             <  1:2220751/150‑3 (MQ=255)

         cGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgctgg                                                                                                                                             <  1:992562/150‑3 (MQ=255)

         cGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgctgg                                                                                                                                             <  1:1483520/150‑3 (MQ=255)

          gTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCAc                                                                                                                                            >  1:337419/1‑150 (MQ=255)

          gTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCAc                                                                                                                                            >  1:1468041/1‑150 (MQ=255)

               tCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                                       >  2:3108378/1‑150 (MQ=255)

               tCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                                       <  2:928121/150‑1 (MQ=255)

               tCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                                       <  1:2278281/150‑1 (MQ=255)

               tCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                                       >  2:463127/1‑150 (MQ=255)

                 aTGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGt                                                                                                                                     >  2:323103/1‑150 (MQ=255)

                    tCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTAcc                                                                                                                                  >  1:3204400/1‑150 (MQ=255)

                    tCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTAcc                                                                                                                                  <  2:862278/150‑1 (MQ=255)

                     cGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCa                                                                                                                                 <  2:3927411/150‑1 (MQ=255)

                     cGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCa                                                                                                                                 <  2:3927107/150‑1 (MQ=255)

                     cGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCa                                                                                                                                 <  2:4264334/150‑1 (MQ=255)

                      gACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAg                                                                                                                                >  1:690351/1‑150 (MQ=255)

                      gACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAg                                                                                                                                >  1:157189/1‑150 (MQ=255)

                        ccaccGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCa                                                                                                                              <  2:831614/150‑1 (MQ=255)

                         caccGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCAc                                                                                                                             <  1:1759460/150‑1 (MQ=255)

                         caccGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCAc                                                                                                                             >  2:2310209/1‑150 (MQ=255)

                         caccGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCAc                                                                                                                             <  2:3902024/150‑1 (MQ=255)

                             ggCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGACACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGa                                                                                                                         >  1:173649/1‑150 (MQ=255)

                              gCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGAc                                                                                                                        <  2:2658005/150‑1 (MQ=255)

                                 cAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCCTCACCAGCAGCCCCGGGCTGCTCACCTCAGGGACCAGCACCTGACGGc                                                                                                                     <  2:3693935/150‑1 (MQ=255)

                                  aGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCt                                                                                                                    >  1:859681/1‑150 (MQ=255)

                                   ggTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCtg                                                                                                                   <  2:1823201/150‑1 (MQ=255)

                                   ggTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCtg                                                                                                                   >  1:229385/1‑150 (MQ=255)

                                    gtgtGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctg                                                                                                                                                 >  2:489776/1‑119 (MQ=255)

                                    gtgtGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctg                                                                                                                                                 <  1:489884/119‑1 (MQ=255)

                                    gtgtGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctg                                                                                                                                                 <  1:489776/119‑1 (MQ=255)

                                     tgtgAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCtggt                                                                                                                 <  1:424918/150‑1 (MQ=255)

                                       tgAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCtggtgg                                                                                                               <  2:2604654/150‑1 (MQ=255)

                                         aaCCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCg                                                                                                             >  1:1213849/1‑150 (MQ=255)

                                           cccGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGAt                                                                                                           >  1:867682/1‑150 (MQ=255)

                                           cccGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGAt                                                                                                           >  2:2984122/1‑150 (MQ=255)

                                           cccGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGAt                                                                                                           >  2:2984197/1‑150 (MQ=255)

                                               gTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCg                                                                                                       >  1:3784496/1‑150 (MQ=255)

                                                tCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGt                                                                                                      <  1:1583034/150‑1 (MQ=255)

                                                 cAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTa                                                                                                     <  2:399608/150‑1 (MQ=255)

                                                  aaCAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAg                                                                                                    >  2:2923974/1‑150 (MQ=255)

                                                   aCAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAgg                                                                                                   >  2:501550/1‑150 (MQ=255)

                                                     agcagTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGcc                                                                                                 >  1:1956964/1‑150 (MQ=255)

                                                      gcagTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCt                                                                                                >  2:2348342/1‑150 (MQ=255)

                                                         gTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGgc                                                                                             <  2:1754991/150‑1 (MQ=255)

                                                         gTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGgc                                                                                             <  1:2195548/150‑1 (MQ=255)

                                                         gTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGgc                                                                                             <  1:4016763/150‑1 (MQ=255)

                                                            ccTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGca                                                                                          <  1:2802815/150‑1 (MQ=255)

                                                                ggCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCAcagc                                                                                      <  1:1267588/150‑1 (MQ=255)

                                                                ggCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCAcagc                                                                                      <  1:1235229/150‑1 (MQ=255)

                                                                ggCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCAcagc                                                                                      <  1:1188252/150‑1 (MQ=255)

                                                                 gCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCAcagcc                                                                                     <  1:1741673/150‑1 (MQ=255)

                                                                     tCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCcagca                                                                                 <  2:838094/150‑1 (MQ=255)

                                                                          tCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttc                                                                            <  1:2386312/150‑1 (MQ=255)

                                                                            cTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttctt                                                                          <  2:586703/150‑1 (MQ=255)

                                                                            cTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttctt                                                                          <  1:2310421/150‑1 (MQ=255)

                                                                               gCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCAc                                                                       >  2:4206743/1‑150 (MQ=255)

                                                                                ccGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCAcc                                                                      <  2:1511339/150‑1 (MQ=255)

                                                                                ccGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCAcc                                                                      <  2:1511684/150‑1 (MQ=255)

                                                                                 cGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCt                                                                     >  2:1351544/1‑150 (MQ=255)

                                                                                 cGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCt                                                                     >  1:191960/1‑150 (MQ=255)

                                                                                  gAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCtt                                                                    <  2:2466924/150‑1 (MQ=255)

                                                                                        cAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGc                                                              <  1:4112888/150‑1 (MQ=255)

                                                                                            gcAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAAt                                                          <  2:1014406/150‑1 (MQ=255)

                                                                                             cAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATg                                                         <  1:3410684/150‑1 (MQ=255)

                                                                                             cAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATg                                                         <  1:3410594/150‑1 (MQ=255)

                                                                                                tGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAAt                                                      <  2:1977052/150‑1 (MQ=255)

                                                                                                    cTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCAct                                                  <  1:2610503/150‑1 (MQ=255)

                                                                                                     tGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCCGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCActc                                                 <  2:2546885/150‑1 (MQ=255)

                                                                                                      ggATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCActct                                                >  2:4272715/1‑150 (MQ=255)

                                                                                                          aCCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTc                                            >  1:214072/1‑150 (MQ=255)

                                                                                                          aCCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTc                                            <  1:3904408/150‑1 (MQ=255)

                                                                                                            cgcgCCCGGCTGGCAGCTCATGCACCATCACCCGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCa                                          >  2:68883/1‑150 (MQ=255)

                                                                                                            cgcgCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCAc                                                                                                                                            >  1:3130822/1‑52 (MQ=255)

                                                                                                            cgcgCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCACATGAATCACTCTCGTCCa                                          >  2:3490978/1‑150 (MQ=255)

                                                                                                            cgcgCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCa                                          >  2:291529/1‑150 (MQ=255)

                                                                                                              cgcCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGc                                        >  1:2678994/1‑150 (MQ=255)

                                                                                                               gcCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCt                                       <  2:3029419/150‑1 (MQ=255)

                                                                                                                 ccGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGc                                     <  1:4272715/150‑1 (MQ=255)

                                                                                                                        gCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCa                              >  1:2926210/1‑150 (MQ=255)

                                                                                                                          aGCTCATGTACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCCGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCCGCTGCCTGTCCAgg                            >  2:3741110/1‑150 (MQ=255)

                                                                                                                          aGCTCATGTACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAgg                            >  2:664554/1‑150 (MQ=255)

                                                                                                                            cTCATGCACCATAACCAGCAGCACCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGAGCACAGACAGCAGTTTCTTCACATTGAAGGCAAATGAATCACTCTCTTCAAGCTGCCTGTCCAGTCt                          >  2:1358279/1‑150 (MQ=255)

                                                                                                                             tCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTg                         >  1:3927137/1‑150 (MQ=255)

                                                                                                                               aTGCACCCTCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGcc                       <  2:3385777/150‑1 (MQ=255)

                                                                                                                                tGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCg                      >  2:2019674/1‑150 (MQ=255)

                                                                                                                                    ccATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGtata                  >  1:1397778/1‑150 (MQ=255)

                                                                                                                                    ccATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGtata                  <  1:4185869/150‑1 (MQ=255)

                                                                                                                                    ccATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGtata                  <  1:2310209/150‑1 (MQ=255)

                                                                                                                                    ccATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGtata                  <  1:4185824/150‑1 (MQ=255)

                                                                                                                                     cATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGACGTATAg                 >  2:130898/1‑150 (MQ=255)

                                                                                                                                         gtgagcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGcc             >  1:2499493/4‑150 (MQ=255)

                                                                                                                                           cagcagCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCg           >  1:676384/1‑150 (MQ=255)

                                                                                                                                            agcagcCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCgg          >  1:1517873/1‑150 (MQ=255)

                                                                                                                                             gcagcCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCggg         >  2:1195838/1‑150 (MQ=255)

                                                                                                                                                    cGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCACGCgg  <  2:4225997/150‑1 (MQ=255)

CGACCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCACGCGG  >  NC_002947/3319030‑3319327

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