BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,067,994	C→G	5.1%	E446D (GAG→GAC) ‡	PP_3581 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,067,994	0	C	G	5.1%	801.0 / 24.0	234	E446D (GAG→GAC) ‡	PP_3581	membrane protein
Reads supporting (aligned to +/- strand): ref base C (124/98); new base G (7/5); total (131/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 242 total aligned reads displayed.

AGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGTC  >  NC_002947/4067845‑4068143

aGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCgctgc                                                                                                                                                       >  1:2728999/1‑150 (MQ=255)

 gCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCgctgct                                                                                                                                                      >  1:4565627/1‑150 (MQ=255)

    cGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCa                                                                                                                                                   >  2:5020490/1‑150 (MQ=255)

    cGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAAATTGGGCTGGTAATGCAAGCTGAACTGGCGCTGCTCCa                                                                                                                                                   <  2:937919/150‑1 (MQ=255)

        cTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAgcgc                                                                                                                                               <  1:3941185/150‑1 (MQ=255)

         ttttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGcc                                                                                                                                              >  2:1589015/1‑150 (MQ=255)

          tttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCa                                                                                                                                             <  2:3816702/150‑1 (MQ=255)

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

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

           ttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGACCTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAg                                                                                                                                            >  1:495892/1‑150 (MQ=255)

               gCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgcgc                                                                                                                                        >  1:1133118/1‑150 (MQ=255)

               gCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgcgc                                                                                                                                        >  1:1512949/1‑150 (MQ=255)

               gCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGAGCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgcgc                                                                                                                                        >  1:109143/1‑150 (MQ=255)

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

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

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

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

                    gCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTc                                                                                                                                   >  1:2353049/1‑150 (MQ=255)

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

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

                      cGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTAGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCt                                                                                                                                 >  1:382296/1‑150 (MQ=255)

                       gATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTg                                                                                                                                >  1:2774647/1‑150 (MQ=255)

                       gATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTg                                                                                                                                >  1:4276085/1‑150 (MQ=255)

                        aTGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgc                                                                                                                               <  1:2856363/150‑1 (MQ=255)

                        aTGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgc                                                                                                                               >  2:921387/1‑150 (MQ=255)

                        aTGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgc                                                                                                                               <  1:3194530/150‑1 (MQ=255)

                        aTGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgc                                                                                                                               >  2:4806552/1‑150 (MQ=255)

                         tGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGTTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgca                                                                                                                              <  2:2016227/150‑1 (MQ=255)

                          gAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgcag                                                                                                                             >  2:5233298/1‑150 (MQ=255)

                             gCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGcagcag                                                                                                                          <  2:4053067/150‑1 (MQ=255)

                                    acccgcagcagGGCCTGCTGCTGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGACCTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAc                                                                                                                   >  2:4537337/5‑150 (MQ=255)

                                     gcagcagcaGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAct                                                                                                                  >  1:972203/1‑150 (MQ=255)

                                     gcagcagcaGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAct                                                                                                                  >  2:3516927/1‑150 (MQ=255)

                                      ccgcagcagGGCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActg                                                                                                                 >  2:5128289/3‑150 (MQ=255)

                                      cagcagcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActg                                                                                                                 <  2:4878431/150‑1 (MQ=255)

                                        gcagcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActgct                                                                                                               >  1:2649321/1‑150 (MQ=255)

                                          agcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTgg                                                                                                             >  2:1001603/1‑150 (MQ=255)

                                            cagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgcg                                                                                                           >  2:3039259/1‑150 (MQ=255)

                                               gCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgcgcgc                                                                                                        >  2:4323028/1‑150 (MQ=255)

                                                 ccTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCAt                                                                                                      >  1:4426556/1‑150 (MQ=255)

                                                 ccTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCAt                                                                                                      >  2:2473785/1‑150 (MQ=255)

                                                   tgctgcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTg                                                                                                    <  1:4293876/150‑1 (MQ=255)

                                                        cGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGtt                                                                                               <  1:1943758/150‑1 (MQ=255)

                                                        cGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGtt                                                                                               <  2:3626054/150‑1 (MQ=255)

                                                         gggTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTc                                                                                              >  2:2384649/1‑150 (MQ=255)

                                                         gggTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTc                                                                                              >  2:4927410/1‑150 (MQ=255)

                                                            tGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTACTGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATc                                                                                           <  2:4569442/150‑1 (MQ=255)

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

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

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

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

                                                                  cAGCGCAGCACCGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGAcaca                                                                                     <  2:5580136/150‑1 (MQ=255)

                                                                  cAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGAcaca                                                                                     >  1:2135095/1‑150 (MQ=255)

                                                                  cAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGAcaca                                                                                     >  2:5343453/1‑150 (MQ=255)

                                                                   aGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACAt                                                                                    >  2:3524664/1‑150 (MQ=255)

                                                                    gcgcAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATc                                                                                   <  1:5240265/150‑1 (MQ=255)

                                                                          caACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaagaa                                                                             <  1:654616/150‑1 (MQ=255)

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

                                                                              gCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTg                                                                         >  1:1640521/1‑150 (MQ=255)

                                                                               cTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGt                                                                        <  2:278343/150‑1 (MQ=255)

                                                                                ttCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCACGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCGGGTCAGCGACCCATCGAAGCCGCTGtc                                                                       >  2:3869545/1‑149 (MQ=255)

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

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

                                                                                        cGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCt                                                               <  2:3195304/150‑1 (MQ=255)

                                                                                         gATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCtt                                                              >  1:2531123/1‑150 (MQ=255)

                                                                                         gATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCtt                                                              >  1:72429/1‑150 (MQ=255)

                                                                                           tCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGc                                                            >  1:1691819/1‑150 (MQ=255)

                                                                                            cGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGcc                                                           <  1:475587/150‑1 (MQ=255)

                                                                                             ggTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCGACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCg                                                          >  1:463290/1‑150 (MQ=255)

                                                                                              gTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATTCAATCGGAACTTGCGCTGCTCCAGCGCCAGGCGCAACTCCTGCAGCAGCTGCAACAGCTGGCGCGCAATGCTGTTCATCGAAACATCGAAGAATCTGTAAACCTTCTTGCCgg                                                         >  1:689698/1‑150 (MQ=255)

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

                                                                                                      gCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTg                                                 <  2:1133118/150‑1 (MQ=255)

                                                                                                      gCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTg                                                 <  1:3482129/150‑1 (MQ=255)

                                                                                                       ccTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTgg                                                <  2:488168/150‑1 (MQ=255)

                                                                                                            gCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATgg                                           >  2:5522369/1‑150 (MQ=255)

                                                                                                              gTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTa                                         >  2:2813995/1‑150 (MQ=255)

                                                                                                                 gAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACAt                                      <  1:1589015/150‑1 (MQ=255)

                                                                                                                   acgtGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATgg                                    <  1:3039945/147‑1 (MQ=255)

                                                                                                                     ttGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCa                                  <  1:809154/150‑1 (MQ=255)

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

                                                                                                                          cTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTc                             >  1:2266385/1‑150 (MQ=255)

                                                                                                                          cTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTc                             >  1:219911/1‑150 (MQ=255)

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

                                                                                                                            ggTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCt                                                                          >  2:2975660/1‑103 (MQ=255)

                                                                                                                                   cAAGCGGAACTGGTGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACg                    <  1:2453338/150‑1 (MQ=255)

                                                                                                                                   cAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACg                    <  1:4419092/150‑1 (MQ=255)

                                                                                                                                    aaGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACgc                   >  2:4555675/1‑150 (MQ=255)

                                                                                                                                    aaGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACgc                   >  2:308404/1‑150 (MQ=255)

                                                                                                                                      gCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTg                                                 >  2:313779/1‑118 (MQ=255)

                                                                                                                                      gCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTg                                                 <  1:313779/118‑1 (MQ=255)

                                                                                                                                      gCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACgcag                 >  1:4998419/1‑150 (MQ=255)

                                                                                                                                      gCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTGTTGGCATGGTACATGGCAGCGTCGGCGTTACgcag                 >  2:1506920/1‑150 (MQ=255)

                                                                                                                                           ccTGGCGCTGGAGCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCt            >  1:1037228/2‑150 (MQ=255)

                                                                                                                                           ccTGGCGCTGGAGCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCt            >  1:924380/2‑150 (MQ=255)

                                                                                                                                            cTGGCGCTGGAGCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTc           >  1:5576222/1‑150 (MQ=255)

                                                                                                                                             tGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGACGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCg          >  1:1953424/1‑150 (MQ=255)

                                                                                                                                               gcgcTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTGTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTg        >  2:5168119/1‑150 (MQ=255)

                                                                                                                                                  ctgctCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTg     <  1:5567408/150‑1 (MQ=255)

                                                                                                                                                    gctCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGt   >  1:1540862/1‑150 (MQ=255)

                                                                                                                                                     ctCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGTc  <  2:1512949/150‑1 (MQ=255)

AGCCCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGTC  >  NC_002947/4067845‑4068143

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