BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,285,044	T→A	5.9%	L90* (TTG→TAG) ‡	PP_1122 →	OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,285,044	0	T	A	5.9%	756.2 / 24.7	220	L90* (TTG→TAG) ‡	PP_1122	OmpA family protein
Reads supporting (aligned to +/- strand): ref base T (103/103); new base A (9/4); total (113/107)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 225 total aligned reads displayed.

CAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTACCCGGCCG  >  NC_002947/1284895‑1285193

cAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCtt                                                                                                                                                       <  1:2112633/150‑1 (MQ=255)

  gCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTgg                                                                                                                                                     <  1:1190845/150‑1 (MQ=255)

  gCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTgg                                                                                                                                                     <  1:2655214/150‑1 (MQ=255)

      gCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGgcgc                                                                                                                                                 <  1:3864720/150‑1 (MQ=255)

        tCACGGGCTGCGCTACGCACAGCGATGGCAGAGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCTCCTTGGGCGCCa                                                                                                                                               >  2:2453128/1‑150 (MQ=255)

          aCGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCAtt                                                                                                                                             <  2:1037886/150‑1 (MQ=255)

          aCGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCAtt                                                                                                                                             <  1:4252099/150‑1 (MQ=255)

            gggCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGc                                                                                                                                           >  1:4420162/1‑150 (MQ=255)

              gCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCcg                                                                                                                                         <  1:3677457/150‑1 (MQ=255)

               cTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCcgg                                                                                                                                        >  1:2805459/1‑150 (MQ=255)

               cTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGCAATGGCGCCCAAGGCGCCaccaccgg                                                                                                                                        <  1:4129694/150‑8 (MQ=255)

                  cgcTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCcggcgg                                                                                                                                     >  1:505745/1‑150 (MQ=255)

                  cgcTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCcggcgg                                                                                                                                     <  1:4075294/150‑1 (MQ=255)

                    cTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGc                                                                                                                                   >  2:710749/1‑150 (MQ=255)

                     tACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCGTGGGCGCCCTTGCCGGCGGGcg                                                                                                                                  >  1:427950/1‑149 (MQ=255)

                      aCGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGCCCGCCGGCAATGGCGCCCAAGGCGCCACCACcggcggccca                                                                                                                                 <  2:795950/150‑5 (MQ=37)

                        gCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAt                                                                                                                               >  2:598246/1‑150 (MQ=255)

                        gCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAt                                                                                                                               >  1:946475/1‑150 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                             >  1:2259344/1‑150 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                             >  2:3658946/1‑150 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                             >  1:374600/1‑150 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                             >  1:4244236/1‑150 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                             <  1:1859138/150‑1 (MQ=255)

                               gATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCGGATTTGTTAc                                                                                                                        >  1:449579/1‑150 (MQ=255)

                                  ggCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCATGGGCGCCATTGCCGGCGGGCTGATTTGTTACGcc                                                                                                                     >  1:2860462/1‑150 (MQ=255)

                                   gCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGccc                                                                                                                    <  1:2401842/150‑1 (MQ=255)

                                    cAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCa                                                                                                                   <  2:2259344/150‑1 (MQ=255)

                                    cAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCa                                                                                                                   <  1:987321/150‑1 (MQ=255)

                                     aGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAg                                                                                                                  <  2:2823923/150‑1 (MQ=255)

                                     aGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAg                                                                                                                  <  2:3248098/150‑1 (MQ=255)

                                       tgcgcCCCACAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCGTGGGCGCCATTGCcggcg                                                                                                                                      >  1:3668786/2‑128 (MQ=255)

                                            cccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTg                                                                                                           <  1:2121825/150‑1 (MQ=255)

                                             ccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGa                                                                                                          <  1:1974230/150‑1 (MQ=255)

                                             ccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGa                                                                                                          <  1:2169912/150‑1 (MQ=255)

                                             ccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGa                                                                                                          <  2:12117/150‑1 (MQ=255)

                                               tcaatcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACg                                                                                                        >  2:3523894/1‑150 (MQ=255)

                                                  atcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAg                                                                                                     >  1:2598899/1‑150 (MQ=255)

                                                  atcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAg                                                                                                     >  2:3990925/1‑150 (MQ=255)

                                                   tcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCCTTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGa                                                                                                    >  1:396726/1‑150 (MQ=255)

                                                   tcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGa                                                                                                    <  2:3780007/150‑1 (MQ=255)

                                                     aaaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATg                                                                                                  <  1:3219827/150‑1 (MQ=255)

                                                     aaaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATg                                                                                                  >  2:2758997/1‑150 (MQ=255)

                                                      aaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATgg                                                                                                 >  2:1475219/1‑150 (MQ=255)

                                                      aaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCAGTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATgg                                                                                                 >  1:3106590/1‑150 (MQ=255)

                                                           ccTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACg                                                                                            <  2:2134694/150‑1 (MQ=255)

                                                           ccTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACg                                                                                            <  1:829644/150‑1 (MQ=255)

                                                            cTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACgg                                                                                           >  1:4403275/1‑150 (MQ=255)

                                                               gCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAt                                                                                        >  1:2796966/1‑150 (MQ=255)

                                                                 ccATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCAATGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAttt                                                                                      >  1:575579/1‑150 (MQ=255)

                                                                  cATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAtttt                                                                                     <  1:161070/150‑1 (MQ=255)

                                                                  cATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAtttt                                                                                     <  1:3781667/150‑1 (MQ=255)

                                                                  cATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAtttt                                                                                     >  1:1367283/1‑150 (MQ=255)

                                                                      tGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGAc                                                                                 >  1:1790661/1‑150 (MQ=255)

                                                                         aGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCgg                                                                              <  2:3438806/150‑1 (MQ=255)

                                                                           ccTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGgcg                                                                            <  2:2846659/150‑1 (MQ=255)

                                                                                 gggCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCa                                                                      <  2:4215877/150‑1 (MQ=255)

                                                                                     ggcTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGc                                                                  <  2:3767416/150‑1 (MQ=255)

                                                                                            tCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCgaga                                                           >  1:3594388/1‑150 (MQ=255)

                                                                                                 ggtggCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCg                                                      <  1:1329872/150‑1 (MQ=255)

                                                                                                 ggtggCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCg                                                      <  1:230892/150‑1 (MQ=255)

                                                                                                      ccTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCaa                                                 <  1:1474105/150‑1 (MQ=255)

                                                                                                           gTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACAcgg                                            >  1:671207/1‑150 (MQ=255)

                                                                                                             gCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACAcggcg                                          <  1:2947907/150‑1 (MQ=255)

                                                                                                              ccATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACAcggcgg                                         <  2:26639/150‑1 (MQ=255)

                                                                                                              ccATCGAGAGTTCGTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAATAcggcgg                                         <  2:4020290/150‑1 (MQ=255)

                                                                                                               cATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTTGGCGCCAATGCCGGCGGGCTGATTTGTTAAGCCCAGGACGGTGACGAAGATGGTTACGGCATTTTAGACCGGCGCGATAACAGCCCCCAGACCCCCGCCAACACCGCGGt                                        >  1:2287481/1‑150 (MQ=255)

                                                                                                                aTCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGtt                                       >  2:2051871/1‑150 (MQ=255)

                                                                                                                aTCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGtt                                       <  2:2658114/150‑1 (MQ=255)

                                                                                                                 tCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTTCCGGCGGGGTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTAGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTg                                      >  1:254178/1‑150 (MQ=255)

                                                                                                                  cGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGa                                     <  1:856316/150‑1 (MQ=255)

                                                                                                                   gagaGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGAc                                    >  2:2727847/1‑150 (MQ=255)

                                                                                                                   gagaGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGAc                                    >  2:1098174/1‑150 (MQ=255)

                                                                                                                    agagTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGAcc                                   <  2:4144350/150‑1 (MQ=255)

                                                                                                                        ttcttcCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGCCGGCATTTTCGACCGGCGCGATAACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACAt                               >  2:3964058/1‑150 (MQ=255)

                                                                                                                         tcttcCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATg                              <  2:440432/150‑1 (MQ=255)

                                                                                                                          cttcCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATgg                             <  1:3478225/150‑1 (MQ=255)

                                                                                                                             ccTGGGCCGCCGGTGGTGGCGCCGTGGGCGCCCATGCCGGCGGGCTGATTTGTTACGCCCCGGACGGGGACGAAGATGGTGACGGCATTTTCGACCCGCGCGATCACAGCCCCCAGACCCCCGCCAACACCGCGGTTGACCCCCAGGGCt                          >  1:598338/1‑150 (MQ=255)

                                                                                                                                gggCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGcc                       <  2:2805459/150‑1 (MQ=255)

                                                                                                                                agcccgccgGCAATGGCGCCCAAGGCGCCACCACCGGCGGCCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGcc                       >  2:470275/4‑150 (MQ=37)

                                                                                                                                 ggCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGccc                      >  2:4305677/1‑150 (MQ=255)

                                                                                                                                  gccgccGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCa                     <  1:1169900/150‑1 (MQ=255)

                                                                                                                                  gccgccGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCa                     <  2:1882918/150‑1 (MQ=255)

                                                                                                                                  gccgccGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATGTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCGGCCCa                     >  1:1102868/1‑150 (MQ=255)

                                                                                                                                   ccgccgGCAATGGCGCCCAAGGCGCCACCACCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCAc                    >  2:2595451/1‑150 (MQ=38)

                                                                                                                                   ccgccgGCAATGGCGCCCAAGGCGCCACCACCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCAc                    >  2:1901726/1‑150 (MQ=38)

                                                                                                                                      ccgGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGa                 <  2:3109688/150‑1 (MQ=255)

                                                                                                                                      ccgGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGa                 <  2:4035623/150‑1 (MQ=255)

                                                                                                                                        ggtggtGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaaa               <  1:1067051/150‑1 (MQ=255)

                                                                                                                                        ggtggtGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaaa               <  2:1617223/150‑1 (MQ=255)

                                                                                                                                         gtggtgGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAAc              <  1:1427592/150‑1 (MQ=255)

                                                                                                                                         gcaatggCGCCCAAGGCGCCACCACCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAAc              >  1:3460032/5‑150 (MQ=255)

                                                                                                                                          tggtggCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACa             >  1:2214354/1‑150 (MQ=255)

                                                                                                                                          tggtggCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACa             <  1:3983529/150‑1 (MQ=255)

                                                                                                                                          tggtggCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACa             <  1:4093752/150‑1 (MQ=255)

                                                                                                                                          tggtggCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACa             <  1:4196219/150‑1 (MQ=255)

                                                                                                                                             tggCGCCCAAGGCGCCACCACCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTa          <  1:358312/150‑1 (MQ=255)

                                                                                                                                                cgcCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTAccc       <  2:224325/150‑1 (MQ=255)

                                                                                                                                                 gcCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTACCCg      >  2:3039859/1‑150 (MQ=255)

                                                                                                                                                     tGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTACCCGGCcg  <  1:2794121/150‑1 (MQ=255)

                                                                                                                                                     tGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTACCCGGCcg  <  1:2431134/150‑1 (MQ=255)

CAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTACCCGGCCG  >  NC_002947/1284895‑1285193

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