BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,718,787	A→C	5.0%	Q675P (CAG→CCG) ‡	PP_4171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,718,787	0	A	C	5.0%	793.3 / 19.3	229	Q675P (CAG→CCG) ‡	PP_4171	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (113/104); new base C (8/4); total (121/108)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

GATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTTGTCCGC  >  NC_002947/4718638‑4718936

gATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgca                                                                                                                                                       <  2:2891754/150‑1 (MQ=255)

   gcgcGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgcagca                                                                                                                                                    <  2:867001/150‑1 (MQ=255)

   gcgcGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgcagca                                                                                                                                                    <  2:2505741/150‑1 (MQ=255)

     gcgAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCagcagc                                                                                                                                                  <  2:5044669/150‑1 (MQ=255)

     gcgAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCagcagc                                                                                                                                                  >  2:5874862/1‑150 (MQ=255)

       gAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAgcgc                                                                                                                                                >  2:4453186/1‑150 (MQ=255)

        aaTGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCGACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAGCGcc                                                                                                                                               <  2:5608614/150‑1 (MQ=255)

        aaTGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGGACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAGCGcc                                                                                                                                               <  2:2762050/150‑1 (MQ=255)

          tGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTg                                                                                                                                             <  1:1541274/150‑1 (MQ=255)

           ggggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTgg                                                                                                                                            <  1:4862094/150‑1 (MQ=255)

           ggggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTgg                                                                                                                                            <  2:6324598/150‑1 (MQ=255)

           ggggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAGCGccgcc                                                                                                                                            <  2:3355978/150‑4 (MQ=255)

            gggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGa                                                                                                                                           <  2:5372096/150‑1 (MQ=255)

             ggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGAc                                                                                                                                          >  2:2411799/1‑150 (MQ=255)

               tGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGa                                                                                                                                        <  1:5218283/150‑1 (MQ=255)

               tGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGa                                                                                                                                        <  1:6477616/150‑1 (MQ=255)

                   ggCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGaa                                                                                                                                    <  1:2601334/150‑1 (MQ=255)

                   ggCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGaa                                                                                                                                    <  2:6722585/150‑1 (MQ=255)

                     cTGGTCGAGCGCGCTGTTCTATTCCGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAAcc                                                                                                                                  <  2:6277387/150‑1 (MQ=255)

                       ggTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCACGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCa                                                                                                                                <  2:4390484/150‑1 (MQ=255)

                       ggTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCa                                                                                                                                <  1:3190578/150‑1 (MQ=255)

                        gTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCAc                                                                                                                               >  1:6764172/1‑150 (MQ=255)

                        gTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGGTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCAc                                                                                                                               >  2:6230457/1‑150 (MQ=255)

                           gAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTc                                                                                                                            >  1:1397117/1‑150 (MQ=255)

                                  cTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTg                                                                                                                     <  2:2252343/150‑1 (MQ=255)

                                   tGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGc                                                                                                                    >  2:3872519/1‑150 (MQ=255)

                                   tGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGc                                                                                                                    <  1:4643463/150‑1 (MQ=255)

                                   tGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGc                                                                                                                    >  2:3232475/1‑150 (MQ=255)

                                   tGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGc                                                                                                                    >  2:4795030/1‑150 (MQ=255)

                                    gTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCa                                                                                                                   <  2:4237162/150‑1 (MQ=255)

                                    gTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCa                                                                                                                   <  1:982265/150‑1 (MQ=255)

                                      tCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGa                                                                                                                 <  1:3692858/150‑1 (MQ=255)

                                       cTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGaa                                                                                                                >  1:6621693/1‑150 (MQ=255)

                                        tATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAAc                                                                                                               >  2:6837741/1‑150 (MQ=255)

                                        tATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAAc                                                                                                               >  1:62209/1‑150 (MQ=255)

                                          ttACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCg                                                                                                             >  1:4195302/1‑150 (MQ=255)

                                              ggTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGACCTGCAGAACCGCTAc                                                                                                         <  2:1739295/150‑1 (MQ=255)

                                                    gATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGc                                                                                                   >  1:3852736/1‑150 (MQ=255)

                                                      tGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCg                                                                                                 >  1:3267593/1‑150 (MQ=255)

                                                        cccTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTc                                                                                               <  1:5919605/150‑1 (MQ=255)

                                                         ccTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAgcgc                                                                                                                                                >  2:2964592/1‑100 (MQ=255)

                                                          cTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGTACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAc                                                                                             >  1:4969664/1‑150 (MQ=255)

                                                          cTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAc                                                                                             <  1:1232830/150‑1 (MQ=255)

                                                           tCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAcc                                                                                            >  2:2674133/1‑150 (MQ=255)

                                                           tCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAcc                                                                                            >  1:5513763/1‑150 (MQ=255)

                                                           tCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAcc                                                                                            >  1:577019/1‑150 (MQ=255)

                                                           tCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAcc                                                                                            >  2:1622051/1‑150 (MQ=255)

                                                             aaCCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCACTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGc                                                                                          <  2:3189219/150‑1 (MQ=255)

                                                                 aGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtga                                                                                      >  2:4895949/1‑150 (MQ=255)

                                                                     ccGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgt                                                                                  >  2:1163343/1‑150 (MQ=255)

                                                                      cGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtg                                                                                 >  1:186551/1‑150 (MQ=255)

                                                                       gCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtga                                                                                <  1:6080811/150‑1 (MQ=255)

                                                                       gCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtga                                                                                <  2:4313222/150‑1 (MQ=255)

                                                                        cTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtgag                                                                               >  2:6515473/1‑150 (MQ=255)

                                                                                 tCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCt                                                                      >  2:6053742/1‑150 (MQ=255)

                                                                                 tCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCt                                                                      <  1:3764334/150‑1 (MQ=255)

                                                                                    ggACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTg                                                                     <  1:2992783/148‑1 (MQ=255)

                                                                                    ggACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTg                                                                     >  2:2992783/1‑148 (MQ=255)

                                                                                     gACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAg                                                                  <  1:5332913/150‑1 (MQ=255)

                                                                                     gACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAg                                                                  <  1:450890/150‑1 (MQ=255)

                                                                                       ccTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGtt                                                                <  1:267647/150‑1 (MQ=255)

                                                                                        cTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTc                                                               >  1:4445715/1‑150 (MQ=255)

                                                                                         tGCGCCTGGCCAAGCGTTTTCGGGTGGATGGCAGCAAACTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCt                                                              >  2:563936/1‑150 (MQ=255)

                                                                                            gcCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGAt                                                           >  1:6206028/1‑150 (MQ=255)

                                                                                            gcCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGAt                                                           >  2:352906/1‑150 (MQ=255)

                                                                                             ccTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATg                                                          >  1:1273852/1‑150 (MQ=255)

                                                                                                ggCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGcc                                                       <  1:1470694/150‑1 (MQ=255)

                                                                                                 gCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGccc                                                      >  1:1046442/1‑150 (MQ=255)

                                                                                                   cAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctg                                                    >  2:1134839/1‑150 (MQ=255)

                                                                                                   cAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctg                                                    >  1:2841429/1‑150 (MQ=255)

                                                                                                    aaGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgc                                                   >  1:4007689/1‑150 (MQ=255)

                                                                                                     aGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgct                                                  >  2:3074415/1‑150 (MQ=255)

                                                                                                     aGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgct                                                  >  1:483553/1‑150 (MQ=255)

                                                                                                        gTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgctgct                                               >  2:2046436/1‑150 (MQ=255)

                                                                                                         ttttCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgctgctg                                              >  1:4743327/1‑150 (MQ=255)

                                                                                                         ttttCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgctgctg                                              >  2:2591050/1‑150 (MQ=255)

                                                                                                             cGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGcc                                          >  1:2272859/1‑150 (MQ=255)

                                                                                                              gggTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGc                                                                                          <  1:3314554/101‑1 (MQ=255)

                                                                                                               ggTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtgagc                                                                              <  2:1250193/112‑1 (MQ=255)

                                                                                                                gtgtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGt                                       <  2:1205105/150‑1 (MQ=255)

                                                                                                                  gtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtg                                     >  2:4642921/1‑150 (MQ=255)

                                                                                                                            aaCCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGcc                           >  1:5881642/1‑150 (MQ=255)

                                                                                                                                tGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCtgtg                       >  2:1581882/1‑150 (MQ=255)

                                                                                                                                   aGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGcc                    <  2:6764172/150‑1 (MQ=255)

                                                                                                                                     gtccaggcgCTGCTGCCACAGAGCCTGGAAGATGAACCCACCTCGGCCCTGCAGAACCGCGACGACAGCCGTCAACGCCTGAGCGTGAGCGCGGACCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGGTGCTTCTCTGCCAGTGGCCGt                  >  1:251543/6‑150 (MQ=255)

                                                                                                                                      tGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCACGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGtt                 >  2:3681629/1‑150 (MQ=255)

                                                                                                                                      tGGCGGCGCTGTGGCAGCAGCGCCTGGACGAAGAACCAACCTCGGCCCTTCAGAACCGCTACGAAAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGATTTCTTATGGACCTGCTGCTGCGCCTGTTGATGCTCTGCCTGTTGCCGtt                 >  2:3245267/1‑150 (MQ=255)

                                                                                                                                        gcggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTgg               <  2:1518558/150‑1 (MQ=255)

                                                                                                                                         aggcgCTGCTGCCACAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTggg              >  1:1343060/2‑150 (MQ=255)

                                                                                                                                          ggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGcc                                                       >  1:2090651/1‑108 (MQ=255)

                                                                                                                                          ggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtgctgc                                 >  1:5784730/1‑130 (MQ=255)

                                                                                                                                          ggcgCTGCTGCCACAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTgggg             >  2:2331359/1‑150 (MQ=255)

                                                                                                                                             gcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCt          <  1:5088195/150‑1 (MQ=255)

                                                                                                                                             gcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCt          <  1:1321178/150‑1 (MQ=255)

                                                                                                                                             gcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCt          >  2:6466874/1‑150 (MQ=255)

                                                                                                                                               tgtgGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCttt        <  1:1664220/150‑1 (MQ=255)

                                                                                                                                               tgtgGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCttt        <  2:5161702/150‑1 (MQ=255)

                                                                                                                                                  ggCAGCAGCGCCTGGACGCTGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTTGTc     <  2:4045695/150‑1 (MQ=255)

                                                                                                                                                  ggCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTTGTc     <  2:186551/150‑1 (MQ=255)

                                                                                                                                                   gcagcaGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTTGTcc    <  2:4696681/150‑1 (MQ=255)

GATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTTGTCCGC  >  NC_002947/4718638‑4718936

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