BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,718,783	T→C	5.1%	W674R (TGG→CGG) ‡	PP_4171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,718,783	0	T	C	5.1%	755.7 / 17.7	233	W674R (TGG→CGG) ‡	PP_4171	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (115/106); new base C (8/4); total (123/110)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 244 total aligned reads displayed.

GGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTTGT  >  NC_002947/4718634‑4718932

ggTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCtgt                                                                                                                                                       <  2:223436/150‑1 (MQ=255)

ggTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCtgt                                                                                                                                                       <  2:5885549/150‑1 (MQ=255)

 gTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGAGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCtgtg                                                                                                                                                      >  2:5103016/1‑150 (MQ=255)

   ggATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgc                                                                                                                                                    <  1:553771/150‑1 (MQ=255)

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

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

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

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

             aTGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCt                                                                                                                                          <  1:5013032/150‑1 (MQ=255)

             aTGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCt                                                                                                                                          <  2:3114197/150‑1 (MQ=255)

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

                 ggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGAc                                                                                                                                      <  2:4286018/150‑1 (MQ=255)

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

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

                        gCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAAc                                                                                                                               <  2:2832019/150‑1 (MQ=255)

                        gCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAAc                                                                                                                               <  2:2976907/150‑1 (MQ=255)

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

                          tGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAAccc                                                                                                                             >  2:3196212/1‑150 (MQ=255)

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

                             tCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCAcc                                                                                                                          >  2:5637193/1‑150 (MQ=255)

                              cGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCt                                                                                                                         <  1:2172380/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                             aTTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAAcc                                                                                                          <  2:63927/150‑1 (MQ=255)

                                              ttACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCg                                                                                                         <  1:6401864/150‑1 (MQ=255)

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

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

                                                 cGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTa                                                                                                      <  2:432572/150‑1 (MQ=255)

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

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

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

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

                                                                    cAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtg                                                                                   >  2:3767777/1‑150 (MQ=255)

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

                                                                     aGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtga                                                                                  <  1:6606974/150‑1 (MQ=255)

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

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

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

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

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

                                                                                        ggACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGa                                                               <  2:2281795/150‑1 (MQ=255)

                                                                                        ggACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGa                                                               <  1:514303/150‑1 (MQ=255)

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

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

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

                                                                                                gcCGGGCCCAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCGGAGTGTGAGCTCGGAGCTGGAGTTCTGAt                                                       <  2:5213071/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                                        aaGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAgcgc                                                                                                                                            <  1:4537624/57‑1 (MQ=37)

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

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

                                                                                                           cGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgctgc                                            >  2:4498484/1‑150 (MQ=255)

                                                                                                              tttCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCtgctgctgc                                         <  1:2417512/150‑1 (MQ=255)

                                                                                                              tttCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCtgctgctgc                                         <  2:2709334/150‑1 (MQ=255)

                                                                                                                 cGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGcc                                      <  1:3720373/150‑1 (MQ=255)

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

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

                                                                                                                   ggTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTg                                    <  1:3842774/150‑1 (MQ=255)

                                                                                                                    gtgtATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGt                                   >  2:3464518/1‑150 (MQ=255)

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

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

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

                                                                                                                           gcagcaACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTgctgct                            >  1:3145049/1‑150 (MQ=255)

                                                                                                                               caACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGc                        <  2:5379935/150‑1 (MQ=255)

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

                                                                                                                                 aCCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCt                      <  1:2724412/150‑1 (MQ=255)

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

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

                                                                                                                                    tGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCtgtg                   <  1:4453186/150‑1 (MQ=255)

                                                                                                                                       aGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGcc                >  1:6028947/1‑150 (MQ=255)

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

                                                                                                                                         cTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGt              <  1:1831306/150‑1 (MQ=255)

                                                                                                                                          tccaggcgCTGCTGCCACAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGtt             >  1:3618995/5‑150 (MQ=255)

                                                                                                                                           ggcggcGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTg            <  2:4445715/150‑1 (MQ=255)

                                                                                                                                            gcggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTgg           >  1:4777206/1‑150 (MQ=255)

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

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

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

                                                                                                                                              ggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtgctgc                             <  2:5784700/130‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

GGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTTGT  >  NC_002947/4718634‑4718932

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