BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,323,653	A→G	5.0%	A301A (GCA→GCG)	PP_2039 →	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,323,653	0	A	G	5.0%	1163.1 / 27.7	355	A301A (GCA→GCG)	PP_2039	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (199/138); new base G (14/4); total (213/142)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 360 total aligned reads displayed.

CATTGCCCCACCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCG  >  NC_002947/2323504‑2323802

 aTTGCCCCACCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCAc                                                                                                                                                      >  1:2233851/1‑150 (MQ=255)

    gCCCCACCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACgca                                                                                                                                                   <  2:5596518/150‑1 (MQ=255)

       ccACCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACgcagca                                                                                                                                                >  1:6206277/1‑150 (MQ=255)

       ccACCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTAGGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACgcagca                                                                                                                                                <  1:7998943/150‑1 (MQ=255)

         aCCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATg                                                                                                                                              <  2:7643227/150‑1 (MQ=255)

         aCCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTCCATGCTGCGTGCAGCATg                                                                                                                                              <  1:766198/150‑1 (MQ=255)

         aCCTTGCAGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATg                                                                                                                                              >  2:5921196/1‑150 (MQ=255)

            ttGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTc                                                                                                                                           <  1:4883707/150‑1 (MQ=255)

              gCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTa                                                                                                                                         <  2:3333513/150‑1 (MQ=255)

                tGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACg                                                                                                                                       >  2:5443674/1‑150 (MQ=255)

                    gCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCAc                                                                                                                                   >  1:5224378/1‑150 (MQ=255)

                    gCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCAc                                                                                                                                   >  1:4963545/1‑150 (MQ=255)

                      tGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACgg                                                                                                                                 <  2:4270394/150‑1 (MQ=255)

                           ggccTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCg                                                                                                                            <  2:4887925/150‑1 (MQ=255)

                             ccTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGt                                                                                                                          >  2:20861/1‑150 (MQ=255)

                               tGGCACTGAACCTCTCGATCTACTATGCCCGCAAGCGCAAGCAGTTCGGCCAGCCTATCGCCAAGGTCCAGTCGATCAACCAAATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGgaga                                                                                                                        >  1:7351487/1‑150 (MQ=255)

                               tGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGa                                                                                                                        <  1:5668447/150‑1 (MQ=255)

                               tGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGa                                                                                                                        <  1:1444498/150‑1 (MQ=255)

                                ggCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGaa                                                                                                                       <  2:7437801/150‑1 (MQ=255)

                                   aCTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGa                                                                                                                    >  2:5061001/1‑150 (MQ=255)

                                           tctcGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCcggcg                                                                                                            >  2:5696801/1‑150 (MQ=255)

                                           tctcGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCcggcg                                                                                                            >  1:6936918/1‑150 (MQ=255)

                                           tctcGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCcggcg                                                                                                            >  2:8265044/1‑150 (MQ=255)

                                            ctcGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCcggcgg                                                                                                           >  2:481967/1‑150 (MQ=255)

                                                atcgaCTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCaag                                                                                                       >  2:7046918/1‑149 (MQ=255)

                                                atcgaCTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAat                                                                                                       >  1:6793692/1‑150 (MQ=255)

                                                 tcgaCTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAata                                                                                                      >  2:2242074/1‑150 (MQ=255)

                                                   gaCTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATc                                                                                                    >  1:1939525/1‑150 (MQ=255)

                                                    aCTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATcc                                                                                                   <  2:7550516/150‑1 (MQ=255)

                                                      tATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGGCGCCGGCGGCAATATCCGc                                                                                                 >  2:7532813/1‑150 (MQ=255)

                                                      tATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGc                                                                                                 >  1:760739/1‑150 (MQ=255)

                                                           ccGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCa                                                                                            >  2:1838158/1‑150 (MQ=255)

                                                           ccGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCa                                                                                            >  2:6594912/1‑150 (MQ=255)

                                                                   gcAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGa                                                                                    >  1:950597/1‑150 (MQ=255)

                                                                   gcAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGa                                                                                    >  1:5095434/1‑150 (MQ=255)

                                                                     aaGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGAtg                                                                                  >  2:3464857/1‑150 (MQ=255)

                                                                     aaGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGAtg                                                                                  >  1:8044671/1‑150 (MQ=255)

                                                                      aGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGAtgt                                                                                 <  2:801819/150‑1 (MQ=255)

                                                                        cAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGc                                                                               <  2:6856104/150‑1 (MQ=255)

                                                                            tCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGc                                                                           <  1:6475248/150‑1 (MQ=255)

                                                                              ggCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCACGCTg                                                                         >  1:2713856/1‑150 (MQ=255)

                                                                              ggCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTg                                                                         >  2:1443408/1‑150 (MQ=255)

                                                                               gCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGt                                                                        <  1:2858002/150‑1 (MQ=255)

                                                                                ccAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGtt                                                                       <  2:2651275/150‑1 (MQ=255)

                                                                                 cAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTc                                                                      <  2:3557108/150‑1 (MQ=255)

                                                                                 cAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTc                                                                      >  1:8125606/1‑150 (MQ=255)

                                                                                   gCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGccc                                                                                       <  1:378275/131‑1 (MQ=255)

                                                                                   gCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGccc                                                                                       >  2:378275/1‑131 (MQ=255)

                                                                                    ccTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGtt                                                                   >  1:3551645/1‑150 (MQ=255)

                                                                                     cTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGttt                                                                  <  1:1279246/150‑1 (MQ=255)

                                                                                                ttCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGcc                                                       <  2:439061/150‑1 (MQ=255)

                                                                                                ttCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGcc                                                       <  1:792463/150‑1 (MQ=255)

                                                                                                 tCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCt                                                      >  1:534310/1‑150 (MQ=255)

                                                                                                  ccAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCtt                                                     >  1:7738117/1‑150 (MQ=255)

                                                                                                   cAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTc                                                    >  1:3025773/1‑150 (MQ=255)

                                                                                                    aGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACACAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTcc                                                   >  1:3401983/1‑150 (MQ=255)

                                                                                                     gTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCg                                                  >  2:6226174/1‑150 (MQ=255)

                                                                                                      tCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGc                                                 >  2:2498717/1‑150 (MQ=255)

                                                                                                        gATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCAt                                               >  1:4618926/1‑150 (MQ=255)

                                                                                                            aaCCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGcc                                           >  1:2151596/1‑150 (MQ=255)

                                                                                                             cccATATGATCGCCGACGTGGCCACCGAGCTGCATGCTGCACGCAGCTTGGTCGACGCCACGGCGGCGGTGCATGACGCCGGCGGCACTATCCGCACCCAGGCCCCGAGGTGCGCGCTGTTCGTGTCGGAGACGGCCTTCCGCAGCGCCg                                          >  2:4020313/2‑150 (MQ=255)

                                                                                                             aCCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCg                                          >  2:7708411/1‑150 (MQ=255)

                                                                                                              ccATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGa                                         <  2:1652156/150‑1 (MQ=255)

                                                                                                               cATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGAt                                        >  1:5762959/1‑150 (MQ=255)

                                                                                                                atatGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATc                                       <  1:4157/150‑1 (MQ=255)

                                                                                                                 tatGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATcg                                      >  1:1194141/1‑150 (MQ=255)

                                                                                                                  atGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATcgc                                     >  2:3193157/1‑150 (MQ=255)

                                                                                                                   tGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATcgcg                                    >  2:1151570/1‑150 (MQ=255)

                                                                                                                    gATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCgcgc                                   >  1:3985074/1‑150 (MQ=255)

                                                                                                                     attgCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGcc                                  >  2:3601287/4‑150 (MQ=255)

                                                                                                                     aTCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGcc                                  <  1:1929109/150‑1 (MQ=255)

                                                                                                                       cGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGt                                <  2:8651985/150‑1 (MQ=255)

                                                                                                                       cGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGt                                >  2:8181569/1‑150 (MQ=255)

                                                                                                                       cGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCAACGCCGATCGCGCCgg                                <  1:669285/150‑2 (MQ=255)

                                                                                                                         ccGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGc                              <  1:441862/150‑1 (MQ=255)

                                                                                                                         ccGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGc                              >  1:83136/1‑150 (MQ=255)

                                                                                                                          cGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCa                             <  2:4616994/150‑1 (MQ=255)

                                                                                                                           gACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAg                            <  1:6124382/150‑1 (MQ=255)

                                                                                                                              aTGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATa                         >  2:1248837/1‑150 (MQ=255)

                                                                                                                              aTGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATa                         >  1:1648397/1‑150 (MQ=255)

                                                                                                                                ggCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGAGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGAGCCGTGCAGATaca                       >  2:5802502/1‑150 (MQ=255)

                                                                                                                                    aCCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACAcggc                   <  1:5060591/150‑1 (MQ=255)

                                                                                                                                      cGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACAcggcgg                 >  1:1020879/1‑150 (MQ=255)

                                                                                                                                      cGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACAcggcgg                 >  1:7596040/1‑150 (MQ=255)

                                                                                                                                        tagaccATGCTGCGTGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTg               >  1:5136898/6‑150 (MQ=255)

                                                                                                                                          cTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGcc             <  2:6739547/150‑1 (MQ=255)

                                                                                                                                          cTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGcc             <  2:5198398/150‑1 (MQ=255)

                                                                                                                                          cTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGcc             >  2:1047396/1‑150 (MQ=255)

                                                                                                                                           tGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCg            >  2:8404349/1‑150 (MQ=255)

                                                                                                                                           accATGCTGCGTGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCg            >  1:8508595/3‑150 (MQ=255)

                                                                                                                                           accATGCTGCGTGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCg            >  1:5597162/3‑150 (MQ=255)

                                                                                                                                            gCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCgg           <  1:3909723/150‑1 (MQ=255)

                                                                                                                                             cATGCTGCGTGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCggg          >  1:8602050/1‑150 (MQ=255)

                                                                                                                                             cATGCTGCGTGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCggg          >  1:2514416/1‑150 (MQ=255)

                                                                                                                                             cATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCggg          >  2:7328396/1‑150 (MQ=255)

                                                                                                                                                  tgcACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGgctgct     <  1:3292758/150‑1 (MQ=255)

                                                                                                                                                   gcACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGctgctg    >  1:709362/1‑150 (MQ=255)

                                                                                                                                                   gcACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGctgctg    >  1:8685633/1‑150 (MQ=255)

                                                                                                                                                     aCGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCg  >  1:8069381/1‑150 (MQ=255)

                                                                                                                                                     aCGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCg  >  2:2905660/1‑150 (MQ=255)

CATTGCCCCACCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCG  >  NC_002947/2323504‑2323802

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