BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,994,004	A→C	5.1%	H510P (CAC→CCC)	oplA →	5‑oxoprolinase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,994,004	0	A	C	5.1%	1205.1 / 33.6	357	H510P (CAC→CCC)	oplA	5‑oxoprolinase A
Reads supporting (aligned to +/- strand): ref base A (151/187); new base C (12/7); total (163/194)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 360 total aligned reads displayed.

CATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCAGTG  >  NC_002947/3993855‑3994153

cATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCa                                                                                                                                                       <  2:2864604/150‑1 (MQ=255)

cATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCAGGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCa                                                                                                                                                       <  2:8032835/150‑1 (MQ=255)

    cACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTc                                                                                                                                                   <  1:4196973/150‑1 (MQ=255)

    cACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTc                                                                                                                                                   <  1:3391332/150‑1 (MQ=255)

       gtggtgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGAc                                                                                                                                                >  1:8295559/1‑150 (MQ=255)

       gtggtgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGAc                                                                                                                                                <  2:1115104/150‑1 (MQ=255)

       gtggtgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGAc                                                                                                                                                <  2:4091986/150‑1 (MQ=255)

           tgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTgg                                                                                                                                            <  2:8348039/150‑1 (MQ=255)

           tgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAgg                                                                                                                                            <  2:6647410/150‑1 (MQ=255)

             cgcCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGgcg                                                                                                                                          <  1:4596600/150‑1 (MQ=255)

                cTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCaa                                                                                                                                       <  1:6691960/150‑1 (MQ=255)

                cTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCaa                                                                                                                                       >  2:6586439/1‑150 (MQ=255)

                 tGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAAc                                                                                                                                      >  1:2616569/1‑150 (MQ=255)

                 tGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAAc                                                                                                                                      >  2:8290708/1‑150 (MQ=255)

                 tGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAAc                                                                                                                                      >  1:7615450/1‑150 (MQ=255)

                 tGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAAc                                                                                                                                      >  1:1589761/1‑150 (MQ=255)

                      gaagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAAc                                                                                                                                 >  1:7568665/1‑150 (MQ=255)

                      gaagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAAc                                                                                                                                 >  1:4801430/1‑150 (MQ=255)

                       aagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAAcc                                                                                                                                >  2:6823102/1‑150 (MQ=255)

                        agaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCa                                                                                                                               >  1:6839110/1‑150 (MQ=255)

                        agaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCa                                                                                                                               >  1:7164641/1‑150 (MQ=255)

                          aaaTTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAgg                                                                                                                             <  1:4022906/150‑1 (MQ=255)

                               ggCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAgggtagt                                                                                                                        <  2:2043772/150‑6 (MQ=255)

                                gCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCCCGTCGACCTGGCGCAACTCAACCAGGCCCTCg                                                                                                                       >  1:1722434/1‑150 (MQ=255)

                                  aTGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGcc                                                                                                                     >  1:8495520/1‑150 (MQ=255)

                                     cacaCCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGaa                                                                                                                  <  2:7787399/150‑1 (MQ=255)

                                     cacaCCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGaa                                                                                                                  <  2:691892/150‑1 (MQ=255)

                                      acacCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAAc                                                                                                                 <  2:738980/150‑1 (MQ=255)

                                      acacCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCCACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAAc                                                                                                                 <  1:4304562/150‑1 (MQ=255)

                                       cacCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCCGCCTGTCGATGCGCCTGGACGCCGGGCCCGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACt                                                                                                                >  2:5530822/1‑150 (MQ=255)

                                        acCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTa                                                                                                               >  1:1385908/1‑150 (MQ=255)

                                         ccGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTAc                                                                                                              <  2:5838791/150‑1 (MQ=255)

                                           gTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGt                                                                                                            <  1:5141773/150‑1 (MQ=255)

                                               tGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGc                                                                                                        <  2:5203701/150‑1 (MQ=255)

                                                gATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCa                                                                                                       >  1:5836902/1‑150 (MQ=255)

                                                   ccccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCATGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAgg                                                                                                    >  1:7056897/1‑150 (MQ=255)

                                                    cccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGc                                                                                                   >  1:4831707/1‑150 (MQ=255)

                                                     ccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGcc                                                                                                  >  2:7197278/1‑150 (MQ=255)

                                                     ccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGcc                                                                                                  >  1:7953903/1‑150 (MQ=255)

                                                     ccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAgg                                                                                                                                            <  2:2764267/108‑1 (MQ=37)

                                                     ccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAgg                                                                                                                                            >  1:2764267/1‑108 (MQ=37)

                                                     ccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGAAGCGCCAGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAATGcc                                                                                                  >  1:8133426/1‑150 (MQ=255)

                                                      cAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCt                                                                                                 <  1:5814846/150‑1 (MQ=255)

                                                       aaGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCtt                                                                                                <  1:3247081/150‑1 (MQ=255)

                                                        aGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTg                                                                                               <  1:2153647/150‑1 (MQ=255)

                                                             gCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAAc                                                                                          <  1:5607883/150‑1 (MQ=255)

                                                               cTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCCCGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCGTGCCAACCt                                                                                        >  1:29036/1‑150 (MQ=255)

                                                               cTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCt                                                                                        >  2:6459143/1‑150 (MQ=255)

                                                                     cTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGa                                                                                  >  1:4760622/1‑150 (MQ=255)

                                                                      ttGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGaa                                                                                 <  1:4487699/150‑1 (MQ=255)

                                                                      ttGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGaa                                                                                 <  2:1203053/150‑1 (MQ=255)

                                                                      ttGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGaa                                                                                 <  2:4661755/150‑1 (MQ=255)

                                                                      ttGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGaa                                                                                 <  2:227218/150‑1 (MQ=255)

                                                                         tgtgCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGAc                                                                              <  1:6351049/150‑1 (MQ=255)

                                                                              gTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGACGACGGCtt                                                                         >  2:6264858/1‑150 (MQ=255)

                                                                               ttCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCttt                                                                        <  1:908035/150‑1 (MQ=255)

                                                                                tCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTg                                                                       >  2:1363237/1‑150 (MQ=255)

                                                                                tCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTg                                                                       <  1:8275873/150‑1 (MQ=255)

                                                                                  ggCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGc                                                                     <  1:1079559/150‑1 (MQ=255)

                                                                                  ggCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGc                                                                     <  2:2403402/150‑1 (MQ=255)

                                                                                     cAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGAc                                                                  >  1:8582461/1‑150 (MQ=255)

                                                                                      aGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGAcc                                                                 >  1:5753695/1‑150 (MQ=255)

                                                                                       gATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCa                                                                <  2:6148110/150‑1 (MQ=255)

                                                                                       gATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCa                                                                <  1:749222/150‑1 (MQ=255)

                                                                                        aTCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAg                                                               <  2:4323255/150‑1 (MQ=255)

                                                                                            tGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCca                                                           <  2:3198229/150‑1 (MQ=255)

                                                                                                 gACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACACCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGt                                                      >  1:2415665/1‑150 (MQ=255)

                                                                                                  aCGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGtt                                                     >  1:5420583/1‑150 (MQ=255)

                                                                                                  aCGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGtt                                                     >  2:5075354/1‑150 (MQ=255)

                                                                                                    gTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGc                                                   <  1:6597496/150‑1 (MQ=255)

                                                                                                    gTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGc                                                   <  2:7028534/150‑1 (MQ=255)

                                                                                                          tACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTac                                             >  2:4912805/1‑150 (MQ=255)

                                                                                                          tACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTac                                             >  2:7891204/1‑150 (MQ=255)

                                                                                                          tACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTac                                             >  1:6082633/1‑150 (MQ=255)

                                                                                                             gACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACAcc                                          >  1:6668506/1‑150 (MQ=255)

                                                                                                              aCCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACGACAccc                                         >  2:2597853/1‑150 (MQ=255)

                                                                                                                cAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTg                                       >  2:1721500/1‑150 (MQ=255)

                                                                                                                 aGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTgg                                      <  2:8534087/150‑1 (MQ=255)

                                                                                                                  gCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCCGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGga                                     <  1:3643589/150‑1 (MQ=255)

                                                                                                                  gCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGga                                     <  2:5143974/150‑1 (MQ=255)

                                                                                                                  gCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGga                                     <  2:6291892/150‑1 (MQ=255)

                                                                                                                  gCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGga                                     <  1:6200075/150‑1 (MQ=255)

                                                                                                                      aCCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATg                                 <  1:4914093/150‑1 (MQ=255)

                                                                                                                          gCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCATGCCAGCAGTTGCCACTACACCCTGGAGATGCGCt                             <  1:2378943/150‑1 (MQ=255)

                                                                                                                           ccTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTa                            <  2:879267/150‑1 (MQ=255)

                                                                                                                           ccTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTa                            >  1:6437757/1‑150 (MQ=255)

                                                                                                                              gTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCt                         >  2:559005/1‑150 (MQ=255)

                                                                                                                                 gATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTggg                      <  2:4831707/150‑1 (MQ=255)

                                                                                                                                  aTGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTgggg                     <  1:7351428/150‑1 (MQ=255)

                                                                                                                                    gcgcCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCa                   >  2:3749530/1‑150 (MQ=255)

                                                                                                                                    gcgcCAGGTCGACGTGCCCGGCGTCCAGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCa                   >  1:4082343/1‑150 (MQ=255)

                                                                                                                                      gcCTGGACGCCGGGCACGTCGACCTGGCTCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAACGCCTTCCCAAACTGCGCGAAGACGGCTTAGGCGACCAGGCCAGCAGTTGACACAACACCCTGGAGATGCGCTACCTGGGGCAGa                 >  2:3308922/1‑150 (MQ=255)

                                                                                                                                        cTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGAtt               >  2:4562803/1‑150 (MQ=255)

                                                                                                                                          ggACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCa             <  2:4714637/150‑1 (MQ=255)

                                                                                                                                           gACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCAt            <  2:5314237/150‑1 (MQ=255)

                                                                                                                                             cGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGa          >  2:3965947/1‑150 (MQ=255)

                                                                                                                                             cGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGa          <  2:6729756/150‑1 (MQ=255)

                                                                                                                                             cGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGa          <  1:3850530/150‑1 (MQ=255)

                                                                                                                                                cGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTg       >  2:3887416/1‑150 (MQ=255)

                                                                                                                                                     aCGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCAgtg  >  1:5273607/1‑150 (MQ=255)

CATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCAGTG  >  NC_002947/3993855‑3994153

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