BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,655,675	T→A	6.3%	H173L (CAT→CTT) ‡	PP_4959 ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,655,675	0	T	A	6.3%	1361.6 / 53.6	409	H173L (CAT→CTT) ‡	PP_4959	diguanylate cyclase
Reads supporting (aligned to +/- strand): ref base T (212/171); new base A (9/17); total (221/188)
Fisher's exact test for biased strand distribution p-value = 4.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 412 total aligned reads displayed.

GCCCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGATT  >  NC_002947/5655526‑5655824

  cccTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATgg                                                                                                                                                     <  1:4905329/150‑1 (MQ=255)

    cTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGc                                                                                                                                                   >  2:131040/1‑150 (MQ=255)

      ggTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCAt                                                                                                                                                 <  1:7047006/150‑1 (MQ=255)

       gTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATc                                                                                                                                                >  2:6334509/1‑150 (MQ=255)

          ggCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGTTGATGCCCATGGCATCatc                                                                                                                                             <  2:4029678/150‑3 (MQ=255)

          ggCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCCATGGCATCatc                                                                                                                                             <  1:5190624/150‑3 (MQ=255)

          ggCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCCATGGCATCatc                                                                                                                                             <  1:317980/150‑3 (MQ=255)

          ggCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAac                                                                                                                                             <  1:1082119/150‑1 (MQ=255)

             cTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGTTGATGCCCATGGCATCAtccgt                                                                                                                                          <  1:5512324/150‑3 (MQ=255)

             cTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCCATGGCATCAtccgt                                                                                                                                          <  1:6052774/150‑3 (MQ=255)

               tGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAAccacca                                                                                                                                        >  2:1551579/1‑150 (MQ=255)

                  aGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGaa                                                                                                                                     >  2:2150946/1‑150 (MQ=255)

                  aGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGaa                                                                                                                                     <  2:3182864/150‑1 (MQ=255)

                      gTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGcc                                                                                                                                 >  1:6923543/1‑150 (MQ=255)

                         cGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTc                                                                                                                              >  1:972258/1‑150 (MQ=255)

                         cGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTc                                                                                                                              >  2:238969/1‑150 (MQ=255)

                           aCTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTcgg                                                                                                                            <  2:4353211/150‑1 (MQ=255)

                            cTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCACTTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTcggc                                                                                                                           >  2:232364/1‑150 (MQ=255)

                            cTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTcggc                                                                                                                           <  2:6951807/150‑1 (MQ=255)

                             tCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTcggcg                                                                                                                          >  1:5011961/1‑150 (MQ=255)

                               gcgcCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCggcggc                                                                                                                        >  1:2175251/1‑150 (MQ=255)

                                   cacaGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCAttt                                                                                                                    >  1:2366348/1‑150 (MQ=255)

                                   cacaGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCAttt                                                                                                                    >  2:2598428/1‑150 (MQ=255)

                                   cacaGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCAttt                                                                                                                    >  1:493986/1‑150 (MQ=255)

                                       ggCTGGTACTGGTCAGTTGCACAACGTCGAGTAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCa                                                                                                                <  2:7846145/150‑1 (MQ=255)

                                       ggCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCa                                                                                                                <  2:582471/150‑1 (MQ=255)

                                       ggCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCa                                                                                                                <  1:719369/150‑1 (MQ=255)

                                       ggCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCa                                                                                                                <  2:3295426/150‑1 (MQ=255)

                                          tGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGa                                                                                                             <  2:7636912/150‑1 (MQ=255)

                                           ggTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAACCGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGAt                                                                                                            <  2:1691201/150‑1 (MQ=255)

                                            gTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATg                                                                                                           <  2:828080/150‑1 (MQ=255)

                                             tACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATgg                                                                                                          >  1:8128021/1‑150 (MQ=255)

                                              aCTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGa                                                                                                         >  1:8487016/1‑150 (MQ=255)

                                               cTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGAc                                                                                                        >  2:3831621/1‑150 (MQ=255)

                                                tGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACg                                                                                                       >  1:1958097/1‑150 (MQ=255)

                                                  gTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCa                                                                                                     >  2:8650263/1‑150 (MQ=255)

                                                      gTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGa                                                                                                 >  2:6671700/1‑150 (MQ=255)

                                                       ttGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGaa                                                                                                >  1:5637195/1‑150 (MQ=255)

                                                       ttGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGaa                                                                                                <  2:5779277/150‑1 (MQ=255)

                                                       ttGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGaa                                                                                                >  1:4093788/1‑150 (MQ=255)

                                                             aaCGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCt                                                                                          >  2:8529363/1‑150 (MQ=255)

                                                             aaCGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCt                                                                                          >  1:453824/1‑150 (MQ=255)

                                                             aaCGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCt                                                                                          >  2:8139480/1‑150 (MQ=255)

                                                               cGTCGAGCCATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTgg                                                                                        <  2:5137056/150‑1 (MQ=255)

                                                               cGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTgg                                                                                        <  1:5889999/150‑1 (MQ=255)

                                                                gTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGc                                                                                       >  1:5362387/1‑150 (MQ=255)

                                                                gTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGc                                                                                       >  1:7917938/1‑150 (MQ=255)

                                                                gTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGc                                                                                       <  2:6418504/150‑1 (MQ=255)

                                                                   gAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGt                                                                                    <  2:8717981/150‑1 (MQ=255)

                                                                   gAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGt                                                                                    <  2:1164309/150‑1 (MQ=255)

                                                                   gAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGt                                                                                    >  1:922168/1‑150 (MQ=255)

                                                                   gAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGt                                                                                    >  1:7782229/1‑150 (MQ=255)

                                                                       aaTTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTcc                                                                                <  1:5593184/150‑1 (MQ=255)

                                                                        aTTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCa                                                                               <  1:294482/150‑1 (MQ=255)

                                                                            ggTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTc                                                                           >  2:6063948/1‑150 (MQ=255)

                                                                              tGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCAc                                                                         >  1:7613877/1‑150 (MQ=255)

                                                                                 ccTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGAc                                                                      >  1:7416970/1‑150 (MQ=255)

                                                                                  cTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACt                                                                     >  1:3897681/1‑150 (MQ=255)

                                                                                       ggTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGAGGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAgct                                                                >  2:4710868/1‑150 (MQ=255)

                                                                                           ttGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGctgctg                                                            >  2:4197322/1‑150 (MQ=255)

                                                                                             gCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgca                                                          >  1:8685869/1‑150 (MQ=255)

                                                                                              cACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgcag                                                         >  1:6187535/1‑150 (MQ=255)

                                                                                              cACCGGTGCCGCCAACAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgcag                                                         <  1:8576949/150‑1 (MQ=255)

                                                                                                        gccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCgg                                               >  2:1808279/1‑150 (MQ=255)

                                                                                                             cAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTgcg                                          >  1:2772347/1‑150 (MQ=255)

                                                                                                             cAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTgcg                                          <  2:7221397/150‑1 (MQ=255)

                                                                                                               accGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgct                                        <  2:5520369/148‑1 (MQ=255)

                                                                                                                aCGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgctg                                       >  1:2630213/1‑150 (MQ=255)

                                                                                                                 cGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGATTGCgctgc                                      <  1:4698574/150‑1 (MQ=255)

                                                                                                                   ggAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGctgctg                                    <  1:3350020/150‑1 (MQ=255)

                                                                                                                          gCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAgc                                                                 <  1:2146677/114‑1 (MQ=255)

                                                                                                                            cGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGagcag                           >  1:6488008/1‑150 (MQ=255)

                                                                                                                             ggATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTCGACCagcagc                          <  1:823297/150‑1 (MQ=255)

                                                                                                                              gATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTCGACCAGCAgcg                         <  1:147137/150‑1 (MQ=255)

                                                                                                                                 tGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCt                      <  2:6488008/150‑1 (MQ=255)

                                                                                                                                   gCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCTGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTgg                    >  2:5382591/1‑150 (MQ=255)

                                                                                                                                   gCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTgg                    >  2:4626792/1‑150 (MQ=255)

                                                                                                                                   gCAAAACGGATGATGCCATGGGCATCAACAACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCACCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTgg                    >  1:127255/1‑150 (MQ=255)

                                                                                                                                    cAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGa                   >  2:3442467/1‑150 (MQ=255)

                                                                                                                                     aaaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGAc                  >  2:565511/1‑150 (MQ=255)

                                                                                                                                     aaaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGAc                  >  2:1994663/1‑150 (MQ=255)

                                                                                                                                      aaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATAGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGAcc                 >  2:3751998/1‑150 (MQ=255)

                                                                                                                                        gtggttgatgCCCATGGCATCATCCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCtt               >  2:7492386/6‑150 (MQ=255)

                                                                                                                                        gtggttgatgCCCATGGCATCATCCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCtt               >  1:670683/6‑150 (MQ=255)

                                                                                                                                        gtggttgatgCCCATGGCATCATCCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCtt               >  2:6799819/6‑150 (MQ=255)

                                                                                                                                         cGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTg              <  1:1721397/150‑1 (MQ=255)

                                                                                                                                          ggttgatgCCCATGGCATCATCCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTgg             >  2:5314616/4‑150 (MQ=255)

                                                                                                                                           gttgatgCCCATGGCATCATCCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTg                     >  2:6111798/3‑141 (MQ=255)

                                                                                                                                           gttgatgCCCATGGCATCATCCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTggg            >  1:4636754/3‑150 (MQ=255)

                                                                                                                                                tgCCATGGGCATCAACCACCAGACTGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGa       <  2:4859399/150‑1 (MQ=255)

                                                                                                                                                tgCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGa       <  1:3631343/150‑1 (MQ=255)

                                                                                                                                                tgCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGa       <  1:4400510/150‑1 (MQ=255)

                                                                                                                                                tgCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGa       <  1:68376/150‑1 (MQ=255)

                                                                                                                                                 gCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAg      >  1:7676647/1‑150 (MQ=255)

                                                                                                                                                  ccATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGTCTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAgg     >  2:3076552/1‑150 (MQ=255)

                                                                                                                                                   cATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGa    >  1:1903687/1‑150 (MQ=255)

                                                                                                                                                   cATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAACAGCGCCTGGACCTTGGGCTTGAGGa    >  2:2929955/1‑150 (MQ=255)

                                                                                                                                                    aTGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGAt   >  1:4199932/1‑150 (MQ=255)

                                                                                                                                                    aTGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGAt   >  1:8335647/1‑150 (MQ=255)

                                                                                                                                                     tGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGAtt  >  2:1014089/1‑150 (MQ=255)

GCCCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGATT  >  NC_002947/5655526‑5655824

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