BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,064,968	T→A	6.0%	E446V (GAG→GTG) ‡	PP_3581 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,064,968	0	T	A	6.0%	608.7 / 20.2	182	E446V (GAG→GTG) ‡	PP_3581	membrane protein
Reads supporting (aligned to +/- strand): ref base T (81/90); new base A (5/6); total (86/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 185 total aligned reads displayed.

CCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGTCC  >  NC_002947/4064821‑4065117

ccGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTcc                                                                                                                                                     <  1:4916237/150‑1 (MQ=255)

       tttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCa                                                                                                                                              >  2:1987982/1‑150 (MQ=255)

       tttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCa                                                                                                                                              <  2:548951/150‑1 (MQ=255)

         tCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCAgt                                                                                                                                            <  1:1641922/150‑2 (MQ=255)

         tCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAgg                                                                                                                                            <  1:934680/150‑1 (MQ=255)

         tCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAgg                                                                                                                                            >  1:345312/1‑150 (MQ=255)

          ccGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgc                                                                                                                                           >  2:615710/1‑150 (MQ=255)

          ccGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgc                                                                                                                                           >  1:5179661/1‑150 (MQ=255)

          ccGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgc                                                                                                                                           >  2:2093332/1‑150 (MQ=255)

              cAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCaa                                                                                                                                       <  2:4945632/150‑1 (MQ=255)

              cAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCaa                                                                                                                                       <  2:4055622/150‑1 (MQ=255)

              cAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCaa                                                                                                                                       <  1:2221301/150‑1 (MQ=255)

              cAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCaa                                                                                                                                       <  1:2373580/150‑1 (MQ=255)

              cAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCaa                                                                                                                                       >  2:2218988/1‑150 (MQ=255)

               aGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAg                                                                                                                                      >  1:652855/1‑150 (MQ=255)

                 gCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTc                                                                                                                                    <  1:3959444/150‑1 (MQ=255)

                  ccGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTcc                                                                                                                                   >  2:2846674/1‑150 (MQ=255)

                   cGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCt                                                                                                                                  >  1:3294630/1‑150 (MQ=255)

                        caGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgcagc                                                                                                                             <  2:1101204/149‑1 (MQ=255)

                          gCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACCTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGcagcag                                                                                                                           <  2:4667151/150‑1 (MQ=255)

                           cGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCagcagc                                                                                                                          >  2:773910/1‑150 (MQ=255)

                             gTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCa                                                                                                                                              >  2:2913903/1‑128 (MQ=255)

                             gTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTg                                                                                                                        <  2:1222856/150‑1 (MQ=255)

                                  gcagcagcaGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAgg                                                                                                                                            >  1:2648152/1‑125 (MQ=255)

                                    agcagcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActgc                                                                                                                 >  1:3428641/1‑150 (MQ=255)

                                    agcagcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActgc                                                                                                                 <  2:319374/150‑1 (MQ=255)

                                    agcagcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActgc                                                                                                                 <  1:5189824/150‑1 (MQ=255)

                                         cagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGcc                                                                                                                                               >  1:1778744/1‑115 (MQ=255)

                                             cccTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTTGCGCGCa                                                                                                        <  1:1008008/150‑1 (MQ=255)

                                             cccTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCa                                                                                                        <  1:1008009/150‑1 (MQ=255)

                                             cccTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCa                                                                                                        <  1:3684687/150‑1 (MQ=255)

                                                tgctgcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTg                                                                                                     >  2:5161564/1‑150 (MQ=255)

                                                tgctgcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTg                                                                                                     >  1:4084197/1‑150 (MQ=255)

                                                  ctgcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCt                                                                                                   <  2:798892/150‑1 (MQ=255)

                                                  ctgcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCt                                                                                                   <  2:3505817/150‑1 (MQ=255)

                                                  ctgcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCt                                                                                                   <  2:2663072/150‑1 (MQ=255)

                                                    gcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGt                                                                                                 <  1:4608522/150‑1 (MQ=255)

                                                        gTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCAt                                                                                             >  1:2147770/1‑150 (MQ=255)

                                                        gTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCAt                                                                                             >  1:750914/1‑150 (MQ=255)

                                                         tGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATc                                                                                            <  2:952851/150‑1 (MQ=255)

                                                           cTCCCAGCGCAGCAACGCTTCGGCAACGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGa                                                                                          >  1:4750028/1‑150 (MQ=255)

                                                            tCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGac                                                                                         >  1:2876926/1‑150 (MQ=255)

                                                            tCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGac                                                                                         <  2:4155996/150‑1 (MQ=255)

                                                             cccAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGaca                                                                                        >  1:4686413/1‑150 (MQ=255)

                                                              ccAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGacac                                                                                       >  1:1602790/1‑150 (MQ=255)

                                                              ccAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGacac                                                                                       >  2:5029557/1‑150 (MQ=255)

                                                                aGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACAt                                                                                     <  1:2727436/150‑1 (MQ=255)

                                                                    cagcaACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaa                                                                                 <  2:1602790/150‑1 (MQ=255)

                                                                     agcaACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaag                                                                                >  1:3738560/1‑150 (MQ=255)

                                                                       caACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaagaa                                                                              >  1:1039454/1‑150 (MQ=255)

                                                                        aaCGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACATATCGaagaag                                                                             >  2:3699320/1‑150 (MQ=255)

                                                                        aaCGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                             >  2:3370754/1‑150 (MQ=255)

                                                                        aaCGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                             <  2:2490290/150‑1 (MQ=255)

                                                                            cTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                             >  1:145649/1‑146 (MQ=255)

                                                                             ttCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTa                                                                        <  2:728272/150‑1 (MQ=255)

                                                                               cGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAAc                                                                      <  2:154392/150‑1 (MQ=255)

                                                                                   aCCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGtt                                                                  <  1:4260802/150‑1 (MQ=255)

                                                                                   aCCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGtt                                                                  >  2:2078304/1‑150 (MQ=255)

                                                                                    ccGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTc                                                                 >  1:309866/1‑150 (MQ=255)

                                                                                    ccGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTc                                                                 <  1:4952538/150‑1 (MQ=255)

                                                                                        tCGGTTGGCAGGCGTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAGTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGc                                                             <  2:5076846/150‑1 (MQ=255)

                                                                                        tCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGc                                                             <  2:60983/150‑1 (MQ=255)

                                                                                           gTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGt                                                                         >  1:4050899/1‑135 (MQ=255)

                                                                                           gTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCgg                                                          >  1:3896207/1‑150 (MQ=255)

                                                                                               gCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTc                                                                                                                                    <  1:4117803/72‑1 (MQ=255)

                                                                                                 aGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTTTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCttt                                                    <  1:4020326/150‑1 (MQ=255)

                                                                                                 aGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCttt                                                    <  2:4365636/150‑1 (MQ=255)

                                                                                                 aGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCttt                                                    >  2:3087243/1‑150 (MQ=255)

                                                                                                  ggCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCtttt                                                   <  1:1621394/150‑1 (MQ=255)

                                                                                                  ggCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCtttt                                                   <  1:4366550/150‑1 (MQ=255)

                                                                                                   gtcTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTg                                                  <  1:4511084/148‑1 (MQ=255)

                                                                                                   gCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGa                                                                                          <  1:4787185/110‑1 (MQ=255)

                                                                                                   gCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGa                                                                                          >  2:4787185/1‑110 (MQ=255)

                                                                                                      tGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCa                                               <  2:1352834/150‑1 (MQ=255)

                                                                                                      tGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCa                                               <  2:1366845/150‑1 (MQ=255)

                                                                                                       gggCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCAt                                              <  1:5029557/150‑1 (MQ=255)

                                                                                                        ggCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATg                                             <  1:1364630/150‑1 (MQ=255)

                                                                                                         gCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACATATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATgg                                            <  1:3699320/150‑1 (MQ=255)

                                                                                                            tCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTAc                                         >  2:5357605/1‑150 (MQ=255)

                                                                                                              gAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACAt                                       >  2:4815990/1‑150 (MQ=255)

                                                                                                                  ttGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCa                                   >  2:1829717/1‑150 (MQ=255)

                                                                                                                    gggCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGc                                 <  1:626132/150‑1 (MQ=255)

                                                                                                                       cTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTc                              <  2:1512540/150‑1 (MQ=255)

                                                                                                                           tAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCg                          >  2:5371187/1‑150 (MQ=255)

                                                                                                                                cAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACg                     >  2:3895128/1‑150 (MQ=255)

                                                                                                                                  aGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACgca                   <  1:785333/150‑1 (MQ=255)

                                                                                                                                   gCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCAGTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACgcag                  >  2:1653673/1‑150 (MQ=255)

                                                                                                                                      gAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGcagcag               >  2:2897689/1‑150 (MQ=255)

                                                                                                                                       gccTGGCGCTGGAGCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACAt                                       >  2:4033949/3‑125 (MQ=255)

                                                                                                                                        aCTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCt             >  2:4097797/1‑150 (MQ=255)

                                                                                                                                         cTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGTCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACACAGCAGCTc            <  2:3791828/150‑1 (MQ=255)

                                                                                                                                         cTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTc            >  1:888604/1‑150 (MQ=255)

                                                                                                                                         cTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTc            >  1:888388/1‑150 (MQ=255)

                                                                                                                                         cTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTc            <  2:888604/150‑1 (MQ=255)

                                                                                                                                           ggCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGt          <  1:1348841/150‑1 (MQ=255)

                                                                                                                                            gcgcTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCgg                            <  2:5109024/131‑1 (MQ=255)

                                                                                                                                            gcgcTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCgg                            >  1:5109024/1‑131 (MQ=255)

                                                                                                                                            gcgcTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTg         <  2:3600557/150‑1 (MQ=255)

                                                                                                                                            gcgcTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATAGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTg         >  2:3699259/1‑150 (MQ=255)

                                                                                                                                                tgctCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTgg     >  2:2094720/1‑150 (MQ=255)

CCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGTCC  >  NC_002947/4064821‑4065117

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