BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	508,778	A→T	6.4%	V79D (GTC→GAC)	estP ←	esterase EstP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	508,778	0	A	T	6.4%	725.2 / 27.5	217	V79D (GTC→GAC)	estP	esterase EstP
Reads supporting (aligned to +/- strand): ref base A (123/80); new base T (9/5); total (132/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 221 total aligned reads displayed.

CGGCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTT  >  NC_002947/508629‑508927

     ttGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGc‑cc                                                                                                                                                  <  1:1161443/150‑1 (MQ=255)

         ggTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggc                                                                                                                                              <  2:3462487/150‑1 (MQ=255)

         ggTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggc                                                                                                                                              <  1:3939778/150‑1 (MQ=255)

          gTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggcg                                                                                                                                             >  1:3408234/1‑150 (MQ=255)

          gTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggcg                                                                                                                                             <  2:2915597/150‑1 (MQ=255)

           tttCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGG‑GCGTCGCGCCAggcga                                                                                                                                            <  1:4295564/150‑2 (MQ=255)

               cTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGG‑GCGTCGCGCCAggcg                                                                                                                                             <  1:1527723/145‑1 (MQ=255)

                 ggggggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCg                                                                                                                                      >  2:1490320/1‑150 (MQ=255)

                  gggggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGa                                                                                                                                     <  1:5430596/150‑1 (MQ=255)

                  gggggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGa                                                                                                                                     <  1:1637816/150‑1 (MQ=255)

                           cGGACTCGGACCGTGTCGTGCTCGAGTCCAGGCTCCGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATCACCATCGGGCTGGATACCTACGGATGTCGCCGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCg                                                                                                                            <  2:3263024/150‑1 (MQ=255)

                             gTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCgac                                                                                                                          <  1:90268/150‑3 (MQ=255)

                             gTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGc                                                                                                                          <  1:2326518/150‑1 (MQ=255)

                               ctcGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGa                                                                                                                        <  1:4757604/150‑1 (MQ=255)

                                  ggACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGacac                                                                                                                     >  1:2950440/1‑150 (MQ=255)

                                        tGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGc                                                                                                               >  2:4914068/1‑150 (MQ=255)

                                          tCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGa                                                                                                             >  1:4460122/1‑150 (MQ=255)

                                           cGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGaa                                                                                                            <  1:3538704/150‑1 (MQ=255)

                                            gTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAg                                                                                                           >  2:3173920/1‑150 (MQ=255)

                                            gTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAg                                                                                                           >  1:3242798/1‑150 (MQ=255)

                                              gATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGtt                                                                                                         <  2:3364280/150‑1 (MQ=255)

                                               aTCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTg                                                                                                        >  2:88628/1‑150 (MQ=255)

                                               aTCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTg                                                                                                        >  1:1984056/1‑150 (MQ=255)

                                                   aGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGt                                                                                                    >  1:822114/1‑150 (MQ=255)

                                                     tCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCg                                                                                                  >  1:1226671/1‑150 (MQ=255)

                                                     tCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCg                                                                                                  >  2:4236179/1‑150 (MQ=255)

                                                     tCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCg                                                                                                  >  1:378931/1‑150 (MQ=255)

                                                       cAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGc                                                                                                <  2:1615120/150‑1 (MQ=255)

                                                               gctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGt                                                                                        >  2:1088113/1‑150 (MQ=255)

                                                               gctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGt                                                                                        >  1:4853427/1‑150 (MQ=255)

                                                               gctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGt                                                                                        >  2:5356599/1‑150 (MQ=255)

                                                               gctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGt                                                                                        >  1:3671299/1‑150 (MQ=255)

                                                                 tgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTg                                                                                      >  1:5223161/1‑150 (MQ=255)

                                                                 tgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTg                                                                                      <  1:2458242/150‑1 (MQ=255)

                                                                   ggtggtGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGt                                                                                    >  2:4735462/1‑150 (MQ=255)

                                                                       gtgtACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGaaa                                                                                >  1:93656/1‑150 (MQ=255)

                                                                         gtACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcg                                                                              >  2:1537324/1‑150 (MQ=255)

                                                                          tACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgc                                                                             >  2:3502062/1‑150 (MQ=255)

                                                                            cccGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                           <  2:5223495/150‑1 (MQ=255)

                                                                            cccGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                           >  1:1897220/1‑150 (MQ=255)

                                                                              cgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACgcgcgc                                                                         >  1:3966207/1‑150 (MQ=255)

                                                                                  gACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTgg                                                                     >  1:2257674/1‑150 (MQ=255)

                                                                                  gACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTgg                                                                     >  2:4323793/1‑150 (MQ=255)

                                                                                    cTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTgggg                                                                   <  2:1044886/150‑1 (MQ=255)

                                                                                       cccAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGTGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAc                                                                >  1:4192121/1‑150 (MQ=255)

                                                                                        ccAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAcc                                                               <  2:2257674/150‑1 (MQ=255)

                                                                                        ccAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAcc                                                               >  2:4025151/1‑150 (MQ=255)

                                                                                        ccAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAcc                                                               >  1:2928254/1‑150 (MQ=255)

                                                                                             tattaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaa                                                          <  2:5393904/150‑1 (MQ=255)

                                                                                              attaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaac                                                         >  1:940469/1‑150 (MQ=255)

                                                                                                taCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaccaa                                                       >  2:2964981/1‑150 (MQ=255)

                                                                                                taCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaacaa                                                       <  1:2163120/150‑1 (MQ=255)

                                                                                                 aCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTGGGTGAAACGCGCGCCTGGGGTACCGCCAACAAg                                                      >  2:4599780/1‑150 (MQ=255)

                                                                                                   cATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGaaa                                                                                <  2:1101443/122‑1 (MQ=255)

                                                                                                    aTCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                                   <  1:1003706/150‑1 (MQ=255)

                                                                                                    aTCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                                   <  2:521038/150‑1 (MQ=255)

                                                                                                     tCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCg                                                  <  1:600272/150‑1 (MQ=255)

                                                                                                     tCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCg                                                  >  1:300459/1‑150 (MQ=255)

                                                                                                      cGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCgg                                                 <  1:3025430/150‑1 (MQ=255)

                                                                                                      cGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCgg                                                 <  2:4985032/150‑1 (MQ=255)

                                                                                                        ggCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGa                                               <  1:113864/150‑1 (MQ=255)

                                                                                                          cTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                             <  1:3121028/150‑1 (MQ=255)

                                                                                                           tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAAc                                            >  1:4614122/1‑150 (MQ=255)

                                                                                                           tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAAc                                            >  2:2231156/1‑150 (MQ=255)

                                                                                                           tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAAc                                            >  2:2450168/1‑150 (MQ=255)

                                                                                                              aTACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGc                                         <  1:2544279/150‑1 (MQ=255)

                                                                                                                aCCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCAGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGccc                                       <  2:3584453/150‑1 (MQ=255)

                                                                                                                  cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCgg                                     >  2:674624/1‑150 (MQ=255)

                                                                                                                   tACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGc                                    >  1:2722462/1‑150 (MQ=255)

                                                                                                                      ggATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATc                                 >  2:125372/1‑150 (MQ=255)

                                                                                                                      ggATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATc                                 >  1:283560/1‑150 (MQ=255)

                                                                                                                       gATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCg                                <  2:3697959/150‑1 (MQ=255)

                                                                                                                           tCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGa                            <  1:4961700/150‑1 (MQ=255)

                                                                                                                           tCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGa                            >  2:5118165/1‑150 (MQ=255)

                                                                                                                           tCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGa                            <  1:2231156/150‑1 (MQ=255)

                                                                                                                            cGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAg                           <  1:2794169/150‑1 (MQ=255)

                                                                                                                            cGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAg                           <  2:360778/150‑1 (MQ=255)

                                                                                                                             gACTGGTTAAGGTCGCCTGGCTCGACGA‑CCAGGCGGCCAACGAGGATCATCGGCGACACCGGGGCGAATTTGCCGTCGGCGTCACTTTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCAGGCATCGTCGAgg                          <  2:2478146/150‑1 (MQ=255)

                                                                                                                             gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAgg                          >  1:2213932/1‑150 (MQ=255)

                                                                                                                             gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAgg                          >  1:4376485/1‑150 (MQ=255)

                                                                                                                             gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAgg                          >  2:1839982/1‑150 (MQ=255)

                                                                                                                                gggTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTg                       >  1:2638275/1‑150 (MQ=255)

                                                                                                                                  gTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTc                     >  1:3843694/1‑150 (MQ=255)

                                                                                                                                   ttAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCg                    >  2:2868538/1‑150 (MQ=255)

                                                                                                                                      aGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCa                 >  1:5379983/1‑150 (MQ=255)

                                                                                                                                      aGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCa                 >  1:2614714/1‑150 (MQ=255)

                                                                                                                                        gccGCCTGG‑GCGTCGCGCCAGGCGACCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCaaa               >  1:4157628/3‑150 (MQ=255)

                                                                                                                                        gTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCaaa               >  2:1340902/1‑150 (MQ=255)

                                                                                                                                        gTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCaaa               >  1:5083906/1‑150 (MQ=255)

                                                                                                                                         tCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAac              <  2:1997962/150‑1 (MQ=255)

                                                                                                                                         tCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAac              <  2:1267131/150‑1 (MQ=255)

                                                                                                                                          cGCCTGG‑GCGTCGCGCCAGGCGGTCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAaca             >  2:2842258/1‑150 (MQ=255)

                                                                                                                                            ccTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACg           >  1:4537149/1‑150 (MQ=255)

                                                                                                                                             cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGa          >  2:1783280/1‑150 (MQ=255)

                                                                                                                                                  gcgACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGt                                                                                    >  2:2643647/1‑71 (MQ=255)

                                                                                                                                                  gcgACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAgg     <  1:3284422/150‑1 (MQ=255)

                                                                                                                                                  gcgACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAgg     <  2:4411228/150‑1 (MQ=255)

                                                                                                                                                   cgACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAggg    <  1:1088113/150‑1 (MQ=255)

                                                                                                                                                    gACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGt   <  2:278714/150‑1 (MQ=255)

                                                                                                                                                     aCGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGtt  >  1:1740363/1‑150 (MQ=255)

CGGCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTT  >  NC_002947/508629‑508927

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