BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,863,399	T→C	6.1%	F348S (TTC→TCC) ‡	PP_1666 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,863,399	0	T	C	6.1%	515.6 / 14.3	163	F348S (TTC→TCC) ‡	PP_1666	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (76/77); new base C (4/6); total (80/83)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 167 total aligned reads displayed.

AATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTGG  >  NC_002947/1863251‑1863545

  tCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGa                                                                                                                                                 >  1:4821251/1‑150 (MQ=255)

     tGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACgg                                                                                                                                              >  1:4573861/1‑150 (MQ=255)

        aCAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCTCAGCCAGCACGCAAGCCGTCGAACGGCtt                                                                                                                                           <  1:3830303/150‑1 (MQ=255)

         cAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTCGAACGGCTTg                                                                                                                                          <  2:52333/150‑1 (MQ=255)

          aGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCACAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTCGAACGGCTTgg                                                                                                                                         <  1:2386479/150‑2 (MQ=255)

            cTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGt                                                                                                                                       >  1:4349045/1‑150 (MQ=255)

                gTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTc                                                                                                                                   <  1:1035542/150‑1 (MQ=255)

                gTGCCCAGCTTGGCCTGGCCAAGCTGAAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTc                                                                                                                                   <  1:2243726/150‑1 (MQ=255)

                  gCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGtgctggc                                                                                                                                 <  1:4616228/150‑7 (MQ=255)

                  ctgcAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGc                                                                                                                                 >  2:2699338/4‑150 (MQ=255)

                   tgcAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCt                                                                                                                                >  2:3314922/3‑150 (MQ=255)

                   cccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGcg                                                                                                                                >  2:1268401/1‑149 (MQ=255)

                   cccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCt                                                                                                                                <  1:154107/150‑1 (MQ=255)

                     cAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGACGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTgg                                                                                                                              >  2:4201076/1‑150 (MQ=255)

                      aGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGt                                                                                                                             >  2:2237194/1‑150 (MQ=255)

                        cTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTaa                                                                                                                           <  1:2788087/150‑1 (MQ=255)

                         ttGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAg                                                                                                                          >  1:5235852/1‑150 (MQ=255)

                         ttGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAg                                                                                                                          <  1:3425181/150‑1 (MQ=255)

                           ggcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGa                                                                                                                        <  1:4684367/150‑1 (MQ=255)

                             cctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAgc                                                                                                                      >  1:139395/1‑150 (MQ=255)

                              ctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAgcg                                                                                                                     >  2:1414337/1‑150 (MQ=255)

                                 gccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCtt                                                                                                                  <  2:1029910/150‑1 (MQ=255)

                                  ccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTg                                                                                                                 >  2:1342193/1‑150 (MQ=255)

                                   cAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGc                                                                                                                >  2:4240726/1‑150 (MQ=255)

                                   cAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGc                                                                                                                >  2:5053152/1‑150 (MQ=255)

                                    aaGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGc                                                                                                                                             >  2:4557942/1‑120 (MQ=255)

                                      gCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGt                                                                                                             <  2:3198830/150‑1 (MQ=255)

                                       cTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtg                                                                                                            >  2:1063925/1‑150 (MQ=255)

                                        tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGTGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                                                                                                           >  2:865034/1‑150 (MQ=255)

                                        tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                                                                                                           <  1:3353971/150‑1 (MQ=255)

                                        tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                                                                                                           >  2:2599281/1‑150 (MQ=255)

                                        tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                                                                                                           <  1:1951001/150‑1 (MQ=255)

                                           aGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgactg                                                                                                        <  1:408539/150‑1 (MQ=255)

                                           aGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgactg                                                                                                        <  1:408588/150‑1 (MQ=255)

                                                  gCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATg                                                                                                   >  2:3479249/1‑148 (MQ=255)

                                                  gCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATg                                                                                                   <  1:3479249/148‑1 (MQ=255)

                                                  gCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCg                                                                                                 >  1:971117/1‑150 (MQ=255)

                                                   ccGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGt                                                                                                >  1:1469294/1‑150 (MQ=255)

                                                      gTCGAGCAGGCACCGGCGGTGCCGGCCATGCAAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGc                                                                                             >  2:4108109/1‑150 (MQ=255)

                                                         gAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTgg                                                                                          <  2:3631905/150‑1 (MQ=255)

                                                               gCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCtt                                                                                                                                           <  1:4310739/95‑1 (MQ=255)

                                                                 aCCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCGTGCCGTGACTGACATGCGTCGCTGGATctggct                                                                                  >  2:3825822/1‑150 (MQ=255)

                                                                         tGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTg                                                                          <  2:1837823/150‑1 (MQ=255)

                                                                            gtccactGCCGGTACTGGCA‑TGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCAct                                                                      <  1:3959844/144‑1 (MQ=255)

                                                                            gtccactGCCGGTACTGGCA‑TGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCAct                                                                      <  1:1661025/144‑1 (MQ=255)

                                                                            cGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCAc                                                                       >  1:1413248/1‑150 (MQ=255)

                                                                               ccATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgc                                                                    <  1:1662200/150‑1 (MQ=255)

                                                                               cactGCCGGTACTGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgc                                                                    >  1:1441190/4‑150 (MQ=255)

                                                                                cATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgct                                                                   >  1:554590/1‑150 (MQ=255)

                                                                                cATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgct                                                                   <  2:1179336/150‑1 (MQ=255)

                                                                                 aTGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgctg                                                                  >  2:1152818/1‑150 (MQ=255)

                                                                                  tGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGa                                                                 >  2:2099214/1‑150 (MQ=255)

                                                                                  tGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGATTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGa                                                                 <  1:677785/150‑1 (MQ=255)

                                                                                     cAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCg                                                              <  1:1152818/150‑1 (MQ=255)

                                                                                       gTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActg                                                                     >  2:436540/1‑141 (MQ=255)

                                                                                        tACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCtg                                                           <  1:2598275/150‑1 (MQ=255)

                                                                                          ccGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCtgtg                                                         >  2:1486258/1‑150 (MQ=255)

                                                                                           cGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGtgc                                                        >  2:3307981/1‑150 (MQ=255)

                                                                                             gCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGtgctg                                                      <  1:2730147/150‑1 (MQ=255)

                                                                                               aGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTgctgct                                                    <  1:696363/150‑1 (MQ=255)

                                                                                               aGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTgctgct                                                    <  1:5053152/150‑1 (MQ=255)

                                                                                                  ggACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTta                                                 >  2:1543733/1‑150 (MQ=255)

                                                                                                   gACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTtat                                                >  1:3264882/1‑150 (MQ=255)

                                                                                                   gACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTtat                                                >  1:3264954/1‑150 (MQ=255)

                                                                                                   gACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTtat                                                >  1:4519012/1‑150 (MQ=255)

                                                                                                    aCGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTtata                                               <  1:4530487/150‑1 (MQ=255)

                                                                                                      gccgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTaa                                                                                                                           >  1:919160/1‑72 (MQ=255)

                                                                                                       ccgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGcc                                            <  2:214659/150‑1 (MQ=255)

                                                                                                        cgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCt                                           >  2:1357463/1‑150 (MQ=255)

                                                                                                          cAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGc                                         >  2:3607802/1‑150 (MQ=255)

                                                                                                          cAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGc                                         >  2:2690422/1‑150 (MQ=255)

                                                                                                               gcgcAGCCCGCGCCTGCGCAGCCAGCACCCCAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAat                                    <  2:2756409/150‑1 (MQ=255)

                                                                                                                cgcAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAata                                   <  2:5225576/150‑1 (MQ=255)

                                                                                                                  cAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATAtt                                 <  2:5415852/150‑1 (MQ=255)

                                                                                                                  cAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATAtt                                 <  2:3040225/150‑1 (MQ=255)

                                                                                                                     cccGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTATGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCtt                              >  2:2888307/1‑150 (MQ=255)

                                                                                                                     cccGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCtt                              <  1:2570514/150‑1 (MQ=255)

                                                                                                                      ccGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTa                             <  2:4055386/150‑1 (MQ=255)

                                                                                                                      ccGCGCCGGCGCAGCCAGCACCCACGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTa                             <  1:5335970/150‑1 (MQ=255)

                                                                                                                       cgcgCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTAc                            <  2:4757718/150‑1 (MQ=255)

                                                                                                                        gcgcCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTg                                          <  1:4100977/135‑1 (MQ=255)

                                                                                                                        gcgcCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACg                           >  1:1560388/1‑150 (MQ=255)

                                                                                                                             tGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGtt                      <  2:4012338/150‑1 (MQ=255)

                                                                                                                              gcgcAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTc                     <  1:3020868/150‑1 (MQ=255)

                                                                                                                                     cagcACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCg              <  2:1560388/150‑1 (MQ=255)

                                                                                                                                     cagcACCCAAGCCGTTCGACGGCTTGCGTTGTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCt                   >  2:1813527/1‑145 (MQ=255)

                                                                                                                                       gcACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGc            <  1:3877174/150‑1 (MQ=255)

                                                                                                                                         aCCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTcc                    <  1:2655875/140‑1 (MQ=255)

                                                                                                                                         aCCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTcc                    >  2:2655875/1‑140 (MQ=255)

                                                                                                                                           gcAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGtgct                                                       <  2:4778779/102‑1 (MQ=255)

                                                                                                                                           gcAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTg                                          <  2:3088214/115‑1 (MQ=255)

                                                                                                                                           gcAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGcatca        >  2:1740769/2‑150 (MQ=255)

                                                                                                                                           ccAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGcatca        >  1:4684993/1‑150 (MQ=255)

                                                                                                                                            cAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGcatcat       >  1:3371798/1‑150 (MQ=255)

                                                                                                                                             aaGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATg      >  1:2326623/1‑150 (MQ=255)

                                                                                                                                               gCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCt    <  1:4429088/150‑1 (MQ=255)

                                                                                                                                                ccGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTg   <  1:1200048/150‑1 (MQ=255)

                                                                                                                                                 cGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTgg  <  2:4821251/150‑1 (MQ=255)

                                                                                                                                                   ttCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGc                                                             <  2:4816834/89‑1 (MQ=255)

                                                                                                                                                   ttCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGc                                                             >  1:4816834/1‑89 (MQ=255)

AATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTGG  >  NC_002947/1863251‑1863545

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