BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,637,084	C→T	5.4%	P283S (CCT→TCT) ‡	lepA →	elongation factor 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,637,084	0	C	T	5.4%	1008.9 / 30.3	314	P283S (CCT→TCT) ‡	lepA	elongation factor 4
Reads supporting (aligned to +/- strand): ref base C (148/148); new base T (3/14); total (152/162)
Fisher's exact test for biased strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 323 total aligned reads displayed.

AGGTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAA  >  NC_002947/1636935‑1637233

  gTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCt                                                                                                                                                     <  2:3025174/150‑1 (MQ=255)

  gTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCt                                                                                                                                                     >  2:1860686/1‑150 (MQ=255)

   tGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTg                                                                                                                                                    >  2:2744576/1‑150 (MQ=255)

    gCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGa                                                                                                                                                   <  1:1251377/150‑1 (MQ=255)

    gCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGa                                                                                                                                                   <  1:5017726/150‑1 (MQ=255)

        ttGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTc                                                                                                                                               >  1:4367998/1‑150 (MQ=255)

          ggTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGa                                                                                                                                             <  2:1549310/150‑1 (MQ=255)

            tCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAgg                                                                                                                                           <  2:6752140/150‑1 (MQ=255)

            tCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAgg                                                                                                                                           >  1:1831392/1‑150 (MQ=255)

             cGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGgt                                                                                                                                          <  2:1902588/150‑1 (MQ=255)

              gACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGAggac                                                                                                                                         <  2:1492933/150‑3 (MQ=255)

              gACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGgtg                                                                                                                                         <  2:4619201/150‑1 (MQ=255)

               aCAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGACCTCGACCTCAGGGGTCGAggaca                                                                                                                                        <  1:1178584/150‑4 (MQ=255)

                cAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGtt                                                                                                                                       <  2:4644146/150‑1 (MQ=255)

                cAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGtt                                                                                                                                       <  2:351660/150‑1 (MQ=255)

                   cGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGcc                                                                                                                                    >  2:6302057/1‑150 (MQ=255)

                    gTCGGTGTCTTCACCCCGATGCACACCCAGACCGCTGCTCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACCTCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCa                                                                                                                                   <  2:2033439/150‑1 (MQ=255)

                     tCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAg                                                                                                                                  >  1:6843129/1‑150 (MQ=255)

                     tCGGTGTCTTAAACCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGTCTTCATCATCGCAAGCATCAAGGACATCCAAGCTGCCCCGGTTGGTGAAACCATGACACTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAg                                                                                                                                  >  2:1363060/1‑150 (MQ=255)

                        gtgtCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCt                                                                                                                               >  1:1188658/1‑150 (MQ=255)

                           tCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAg                                                                                                                                            <  1:3496459/134‑1 (MQ=255)

                            cTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGTCTTCaa                                                                                                                           <  1:1363060/150‑1 (MQ=255)

                            cTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                           <  2:3814224/150‑1 (MQ=255)

                            cTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                           <  1:7234281/150‑1 (MQ=255)

                            cTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                           <  2:50725/150‑1 (MQ=255)

                            cTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                           <  1:2975649/150‑1 (MQ=255)

                            cTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                           <  1:3119588/150‑1 (MQ=255)

                             ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                          >  2:6965111/1‑150 (MQ=255)

                               cACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCGGAGGTCGAGGTGTTGCCAGGCGTCaagac                                                                                                                        >  2:5537140/1‑149 (MQ=255)

                                 ccccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAagaaga                                                                                                                      <  1:4155923/150‑1 (MQ=255)

                                 ccccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAagaaga                                                                                                                      <  2:4853172/150‑1 (MQ=255)

                                 ccccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAagaaga                                                                                                                      <  2:2665519/150‑1 (MQ=255)

                                      aaGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAg                                                                                                                 <  1:769225/150‑1 (MQ=255)

                                          acacCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGc                                                                                                             <  2:3555574/150‑1 (MQ=255)

                                           cacCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCa                                                                                                            <  2:4079074/150‑1 (MQ=255)

                                             cccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAgg                                                                                                          <  2:525770/150‑1 (MQ=255)

                                             cccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAgg                                                                                                          >  1:3681761/1‑150 (MQ=255)

                                              ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCg                                                                                                                                              <  1:1733986/113‑1 (MQ=255)

                                              ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCg                                                                                                                                              >  2:1733986/1‑113 (MQ=255)

                                              ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                         <  1:615379/150‑1 (MQ=255)

                                              ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                         >  1:4785526/1‑150 (MQ=255)

                                              ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                         <  1:4335559/150‑1 (MQ=255)

                                               cAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTc                                                                                                        >  2:5647814/1‑150 (MQ=255)

                                                  aCCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTAc                                                                                                     <  2:943848/150‑1 (MQ=255)

                                                    cGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGc                                                                                                   >  1:5109175/1‑150 (MQ=255)

                                                     gCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGcc                                                                                                  >  2:5148745/1‑150 (MQ=255)

                                                      cTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCg                                                                                                 <  1:6722081/150‑1 (MQ=255)

                                                       tGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGtt                                                                                                                                       <  2:3489059/111‑1 (MQ=255)

                                                       tGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCgg                                                                                                <  1:2911888/150‑1 (MQ=255)

                                                             gAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAg                                                                                                                                            <  1:6622076/100‑1 (MQ=255)

                                                             gAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGtt                                                                                          <  2:5471002/150‑1 (MQ=255)

                                                                aGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTccc                                                                                       >  2:2524780/1‑150 (MQ=255)

                                                                 gCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCg                                                                                      <  2:6612649/150‑1 (MQ=255)

                                                                 gCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCg                                                                                      <  2:1093446/150‑1 (MQ=255)

                                                                 gCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCg                                                                                      <  1:311351/150‑1 (MQ=255)

                                                                 gCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCg                                                                                      <  2:6396021/150‑1 (MQ=255)

                                                                     gCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCa                                                                                  <  1:2191792/150‑1 (MQ=255)

                                                                          gTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTcc                                                                             >  1:1929422/1‑150 (MQ=255)

                                                                             ggCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCcgac                                                                          <  2:840650/150‑1 (MQ=255)

                                                                                   atcatcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTc                                                                    <  2:7196608/150‑1 (MQ=255)

                                                                                    tcatcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCg                                                                   >  1:2677950/1‑150 (MQ=255)

                                                                                     catcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                  <  1:5948988/150‑1 (MQ=255)

                                                                                      atcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAg                                                                 <  2:312718/150‑1 (MQ=255)

                                                                                      atcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAg                                                                 <  2:5776007/150‑1 (MQ=255)

                                                                                       tcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACCTCGAgg                                                                >  1:3284619/1‑150 (MQ=255)

                                                                                           cAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGAttt                                                            >  1:1981258/1‑150 (MQ=255)

                                                                                               aTCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCcgc                                                        <  2:3900117/150‑1 (MQ=255)

                                                                                                tCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCcgcg                                                       <  1:4779008/150‑1 (MQ=255)

                                                                                                  aaGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGAc                                                     >  2:374822/1‑150 (MQ=255)

                                                                                                    ggACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGGAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGc                                                   <  2:5875507/150‑1 (MQ=255)

                                                                                                    ggACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGc                                                   >  1:5310121/1‑150 (MQ=255)

                                                                                                        aTCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTg                                               <  1:623466/150‑1 (MQ=255)

                                                                                                             cGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGaa                                          <  2:3702264/150‑1 (MQ=255)

                                                                                                              ggTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAg                                         <  1:735016/150‑1 (MQ=255)

                                                                                                               gTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGc                                        >  1:1000755/1‑150 (MQ=255)

                                                                                                                  ccccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAgg                                                                                                                                           <  2:4561145/48‑1 (MQ=255)

                                                                                                                  ccccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGGTTTCCGCGACGCTCTGCAGAAGCTGa                                     <  2:1188658/150‑1 (MQ=255)

                                                                                                                    ccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCCGGTCTACGCCGGCCTGTTCCCGGTCCGCTCCGACGACTTCGAGGATTTCCGCGACGCGCTGCAGAAGCTGAcc                                   >  2:4994046/1‑150 (MQ=255)

                                                                                                                    ccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAcc                                   <  1:1860686/150‑1 (MQ=255)

                                                                                                                       gTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTg                                <  1:5490440/150‑1 (MQ=255)

                                                                                                                          ggTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAAc                             <  2:176159/150‑1 (MQ=255)

                                                                                                                                accctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTc                                                                                         >  1:4569357/1‑84 (MQ=255)

                                                                                                                                  cctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcctc                     >  2:6457387/1‑150 (MQ=255)

                                                                                                                                  cctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcctc                     >  1:2505571/1‑150 (MQ=255)

                                                                                                                                  cctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcctc                     >  1:4204955/1‑150 (MQ=255)

                                                                                                                                      accctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTg                 <  2:4318759/150‑1 (MQ=255)

                                                                                                                                        cctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa               <  2:4117190/150‑1 (MQ=255)

                                                                                                                                        cctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa               >  2:2286293/1‑150 (MQ=255)

                                                                                                                                        cctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa               >  1:6317886/1‑150 (MQ=255)

                                                                                                                                         ctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAg              >  2:7382578/1‑150 (MQ=255)

                                                                                                                                             ccTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACa          >  2:3324095/1‑150 (MQ=255)

                                                                                                                                              cTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACAt         >  1:1552504/1‑150 (MQ=255)

                                                                                                                                               tCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATg        >  2:3278025/1‑150 (MQ=255)

                                                                                                                                               tCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATg        >  2:2117890/1‑150 (MQ=255)

                                                                                                                                               tCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATg        >  2:6667717/1‑150 (MQ=255)

                                                                                                                                                 gACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGcc      >  1:4118787/1‑150 (MQ=255)

                                                                                                                                                    cccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGa   >  1:5995444/1‑150 (MQ=255)

                                                                                                                                                    cccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGa   >  1:2300899/1‑150 (MQ=255)

                                                                                                                                                     ccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGaa  >  2:6077709/1‑150 (MQ=255)

                                                                                                                                                     ccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGaa  >  1:2092775/1‑150 (MQ=255)

AGGTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAA  >  NC_002947/1636935‑1637233

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