BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	418,298	C→G	8.9%	L26L (CTC→CTG)	PP_0345 →	heptose kinase WapP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	418,298	0	C	G	8.9%	360.8 / 21.6	120	L26L (CTC→CTG)	PP_0345	heptose kinase WapP
Reads supporting (aligned to +/- strand): ref base C (53/56); new base G (6/5); total (59/61)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

Only 100 of 120 total aligned reads displayed.

GCTCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCAAA  >  NC_002947/418152‑418446

gctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAg                                                                                                                                                    <  2:83503/150‑1 (MQ=255)

 ctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGCGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                   >  1:2147061/1‑150 (MQ=255)

  tCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAggg                                                                                                                                                  <  2:1283753/150‑1 (MQ=255)

   cGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGc                                                                                                                                                 <  2:1855549/150‑1 (MQ=255)

   cGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGc                                                                                                                                                 <  1:1249368/150‑1 (MQ=255)

      ctggctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGaa                                                                                                                                              >  2:2979533/1‑150 (MQ=255)

        ggctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACg                                                                                                                                            <  2:1947696/150‑1 (MQ=255)

         gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTGGAGGGCGAAcac                                                                                                                                           <  1:791282/150‑3 (MQ=255)

         gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGCGTTCGCCCTGGAGGGCGAAcac                                                                                                                                           <  1:125754/150‑3 (MQ=255)

          ctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGc                                                                                                                                          <  1:2001146/150‑1 (MQ=255)

          ctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGAGGCGCAACGCTTGTGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGc                                                                                                                                          <  2:81115/150‑1 (MQ=255)

             gAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTga                                                                                                                                       <  1:2324480/150‑1 (MQ=255)

              ggcgcAAGGCCCAGACGCGCTACGATCGCAAGAAACACAATCAGCACCAGGTCGCATGGCGAGTTGGACCCTGGCCCCTAGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgac                                                                                                                                      <  1:2983902/149‑1 (MQ=255)

              aGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCAAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCCCGCTGCGCCCTCCAGGGCGAACGGCTgac                                                                                                                                      <  1:2385711/150‑1 (MQ=255)

                      gCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacc                                                                                                                              >  2:2403361/1‑150 (MQ=255)

                      gCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacc                                                                                                                              >  1:162547/1‑150 (MQ=255)

                        ccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacccg                                                                                                                            >  2:1500429/1‑150 (MQ=255)

                        ccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacccg                                                                                                                            >  1:826087/1‑150 (MQ=255)

                        ccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGAACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacgcg                                                                                                                            >  2:291436/1‑150 (MQ=255)

                        ccAGAAACTCTACGATCGCAAGAAACGAACCGGGGATGCACAGTAATGGCGAGTTGGACACGGGCGACGGGGAGCGAGCACCTGGGAGCGGACTTCGGCAGCCGCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacccg                                                                                                                            >  2:350772/1‑150 (MQ=255)

                          aGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTTACCCGCGACCCGCt                                                                                                                          <  1:2332634/150‑1 (MQ=255)

                           gAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTc                                                                                                                         >  1:575399/1‑150 (MQ=255)

                            aaaCTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCa                                                                                                                        <  1:2979533/150‑1 (MQ=255)

                                tctACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGa                                                                                                                    >  2:577253/1‑150 (MQ=255)

                                 ctACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGaa                                                                                                                   <  1:1500429/150‑1 (MQ=255)

                                    cGATGGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTg                                                                                                                <  1:1499596/150‑1 (MQ=255)

                                      aTCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGt                                                                                                              >  1:1112367/1‑150 (MQ=255)

                                        cGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTAc                                                                                                            <  1:2963833/150‑1 (MQ=255)

                                         gcaagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCt                                                                                                                          <  2:827494/135‑1 (MQ=255)

                                            agcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCAt                                                                                                        >  2:2551935/1‑150 (MQ=255)

                                             gcaaCGCTATGGGGATGCACTCTGATGGCGAGATGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATc                                                                                                       >  2:349955/1‑150 (MQ=255)

                                              caaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTTGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCg                                                                                                      >  2:2437917/1‑150 (MQ=255)

                                               aaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGa                                                                                                     >  2:1505212/1‑150 (MQ=255)

                                                   cTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTg                                                                                                                >  2:1435335/1‑135 (MQ=255)

                                                   cTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGt                                                                                                >  2:1234806/1‑150 (MQ=255)

                                                       ggggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACg                                                                                            <  1:2403361/150‑1 (MQ=255)

                                                         ggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACggg                                                                                          >  1:612339/1‑150 (MQ=255)

                                                            tGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATc                                                                                       <  2:2894444/150‑1 (MQ=255)

                                                               aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGt                                                                                                <  2:2467838/138‑1 (MQ=255)

                                                               aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGt                                                                                                >  1:2467838/1‑138 (MQ=255)

                                                               aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAac                                                                                    >  1:1094250/1‑150 (MQ=255)

                                                                 tctGATGGCGAGTTTGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATc                                                                                       <  1:154144/145‑1 (MQ=255)

                                                                 tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGTCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACTCGCTCAGTGTAGTGGTACGCATCGA‑GCGTGACGGGATCAacta                                                                                  <  1:1424911/150‑1 (MQ=255)

                                                                 tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATc                                                                                       >  2:154144/1‑145 (MQ=255)

                                                                          gAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCaa                                                                         >  2:2391180/1‑150 (MQ=255)

                                                                            gTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAgc                                                                       <  1:2170139/150‑1 (MQ=255)

                                                                                 acacTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTac                                                                  >  1:2800070/1‑150 (MQ=255)

                                                                                  cacTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTAACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTaca                                                                 >  1:1138180/1‑150 (MQ=255)

                                                                                   acTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACggg                                                                                          <  1:2586566/124‑1 (MQ=255)

                                                                                   acTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACggg                                                                                          >  2:2586566/1‑124 (MQ=255)

                                                                                   acTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTCTCGATCTACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTacac                                                                >  1:191545/1‑150 (MQ=255)

                                                                                   acTGGCGCCGGGGTACGAACACCGGGCGGCGGCCTTCGGCAGCCCCGACGCGGTGTTCAGCCTGCAGGACGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGAACAACTACTACGTCAAGCGCTacac                                                                >  2:1190731/1‑150 (MQ=255)

                                                                                      ggCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCgg                                                             <  2:309099/150‑1 (MQ=255)

                                                                                      ggCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCgg                                                             <  2:575399/150‑1 (MQ=255)

                                                                                       gcgcCGGGGTACGAACACCTGGCGGCGGACTTCGGTAGCCTCGACGCGGTGTTCGCCCGCCAGGTCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCggg                                                            <  1:777801/150‑1 (MQ=255)

                                                                                       gcgcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTTCTTTAAGCGCTACACCggg                                                            <  2:1809757/150‑1 (MQ=255)

                                                                                       gcgcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCggg                                                            <  2:162547/150‑1 (MQ=255)

                                                                                          ccGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGcc                                                         <  2:1854205/150‑1 (MQ=255)

                                                                                           cGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCg                                                        <  2:826087/150‑1 (MQ=255)

                                                                                              ggTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAgc                                                                       >  2:345518/1‑132 (MQ=255)

                                                                                              ggTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAgc                                                                       <  1:345518/132‑1 (MQ=255)

                                                                                              ggTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCa                                                     >  1:1580699/1‑150 (MQ=255)

                                                                                               gTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCaa                                                    <  2:2800070/150‑1 (MQ=255)

                                                                                                tACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAg                                                   >  2:144505/1‑150 (MQ=255)

                                                                                                tACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAg                                                   <  1:577253/150‑1 (MQ=255)

                                                                                                   gAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGcac                                                <  2:2670165/150‑1 (MQ=255)

                                                                                                   gAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGcac                                                >  2:411583/1‑150 (MQ=255)

                                                                                                    aaCACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGcaca                                               >  2:2137324/1‑150 (MQ=255)

                                                                                                       acCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGaa                                                                                                                   >  2:1742796/1‑80 (MQ=255)

                                                                                                       acCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGaa                                                                                                                   <  1:1742796/80‑1 (MQ=255)

                                                                                                       acCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGC‑GCCTTACTGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATgc                                            <  1:2521077/150‑1 (MQ=255)

                                                                                                          tGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGcc                                         <  1:2391180/150‑1 (MQ=255)

                                                                                                               gcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAac                                                                                    >  1:1843317/1‑102 (MQ=255)

                                                                                                                     ttCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGGGAGCGGCTGACCCGCGAGCCGCTCAGTGACGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTgggg                              >  1:1397072/1‑150 (MQ=255)

                                                                                                                     ttCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTgggg                              >  2:414634/1‑150 (MQ=255)

                                                                                                                      tCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggc                             <  1:350772/150‑1 (MQ=255)

                                                                                                                       cGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGAAAACACACGT‑ACGAGATCAACTACTACGTCAAGCGCTCCACCAGGGCCGGCAAGCACAGGCGCGGGTACCTGGggcg                            >  1:1301851/1‑150 (MQ=255)

                                                                                                                         gCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCCTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcggc                          <  1:1690846/150‑1 (MQ=255)

                                                                                                                          cAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGGGCCGGTACCTGGGGCGGcc                         >  2:2019113/1‑150 (MQ=255)

                                                                                                                           aGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCa                        <  1:1505212/150‑1 (MQ=255)

                                                                                                                               tCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCa                    >  2:2917006/1‑150 (MQ=255)

                                                                                                                                 gACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATc                  >  2:1785969/1‑150 (MQ=255)

                                                                                                                                  aCGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCa                 >  1:1339682/1‑150 (MQ=255)

                                                                                                                                   cgcgGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGa                                                                                         <  1:3066939/77‑1 (MQ=255)

                                                                                                                                   cgcgGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGa                                                                                         >  2:3066939/1‑77 (MQ=255)

                                                                                                                                    gcgGTGTTCGCCCTCCGGGGCGAACGGCTGACCCGCGACCCGCGCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAg               >  2:1722412/1‑150 (MQ=255)

                                                                                                                                      ggTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGc             <  2:1339682/150‑1 (MQ=255)

                                                                                                                                      gccgtTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggc                             <  2:1334316/131‑1 (MQ=255)

                                                                                                                                      gccgtTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggc                             >  1:1334316/4‑134 (MQ=255)

                                                                                                                                       ccgtTCGCCCTGGAGGGCGAACAGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCgggg                                                           >  1:2533000/3‑103 (MQ=255)

                                                                                                                                       ccgtTCGCCCTGGAGGGCGAACAGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCgggg                                                           <  2:2533000/101‑1 (MQ=255)

                                                                                                                                        cgtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTgg                                <  1:2448749/128‑1 (MQ=255)

                                                                                                                                        cgtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTgg                                >  2:2448749/2‑129 (MQ=255)

                                                                                                                                         gtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAac                                                                                    >  2:781897/1‑76 (MQ=255)

                                                                                                                                         gtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAac                                                                                    <  1:781897/76‑1 (MQ=255)

                                                                                                                                          ttCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACCCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGaa         >  2:2678005/1‑150 (MQ=255)

                                                                                                                                           tCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAAt        <  1:2019113/150‑1 (MQ=255)

                                                                                                                                           tCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAAt        <  1:1785969/150‑1 (MQ=255)

                                                                                                                                             gCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcggc                          <  2:555514/130‑1 (MQ=255)

                                                                                                                                             gCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcggc                          >  1:555514/1‑130 (MQ=255)

GCTCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGA‑GCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCAAA  >  NC_002947/418152‑418446

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`