BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,285,033	T→C	5.0%	G86G (GGT→GGC)	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,033	0	T	C	5.0%	659.9 / 14.3	199	G86G (GGT→GGC)	PP_1122	OmpA family protein
Reads supporting (aligned to +/- strand): ref base T (98/90); new base C (7/3); total (106/93)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 223 total aligned reads displayed.

TTGGTGGCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACA  >  NC_002947/1284884‑1285182

 tggtggCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCggtg                                                                                                                                                      <  2:194408/150‑1 (MQ=255)

   gtggCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCggtggt                                                                                                                                                    >  1:2827497/1‑150 (MQ=255)

     ggCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGtggtgg                                                                                                                                                  <  1:2033091/150‑1 (MQ=255)

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

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

               cTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGgc                                                                                                                                        >  1:2794506/1‑150 (MQ=255)

               cTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGgc                                                                                                                                        <  1:1447072/150‑1 (MQ=255)

                 gCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGgcgc                                                                                                                                      >  1:2622835/1‑150 (MQ=255)

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

                   tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCa                                                                                                                                    >  1:3780007/1‑150 (MQ=255)

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

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

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

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

                                tACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCt                                                                                                                       <  2:2666544/150‑1 (MQ=255)

                                tACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCt                                                                                                                       >  2:3987292/1‑150 (MQ=255)

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

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

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

                                    cacaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAtt                                                                                                                   <  2:2377340/150‑1 (MQ=255)

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

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

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

                                      caGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTg                                                                                                                 >  2:3781667/1‑150 (MQ=255)

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

                                           aTGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACg                                                                                                            >  2:3445970/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                        ccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCGGATTTGTTACGCCCAGGACGGTGa                                                                                               >  2:1030908/1‑150 (MQ=255)

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

                                                          tcaatcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCATGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACg                                                                                             >  1:1541287/1‑150 (MQ=255)

                                                            aatcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGt                                                                                                 <  1:2801584/144‑1 (MQ=255)

                                                              tcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCGGATTTGTTACGCCCAGGACGGTGACGAAGa                                                                                         >  1:396867/1‑150 (MQ=255)

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

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

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

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

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

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

                                                                        tGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGc                                                                               >  1:4030669/1‑150 (MQ=255)

                                                                            ccATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAttt                                                                           >  1:2118571/1‑150 (MQ=255)

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

                                                                             cATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCGTGGGCGCCATTGCCGGCGGGCTGATTTGTTAAGCCCAGGACGGTGACGAAGATGGTGACGGCAtttt                                                                          >  1:3790383/1‑150 (MQ=255)

                                                                               tCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCg                                                                        >  1:3672686/1‑150 (MQ=255)

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

                                                                                      ccTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGgcg                                                                 >  1:1425495/1‑150 (MQ=255)

                                                                                         gctgGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGa                                                              <  1:272752/150‑1 (MQ=255)

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

                                                                                               cggcTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTg                                                        >  1:2557896/1‑150 (MQ=255)

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

                                                                                                   ttGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTTGGCGCCAATGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAATTTAGACCGGCGCGATAACTGcccc                                                    >  1:800973/1‑150 (MQ=255)

                                                                                                       tCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTTCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCgaga                                                >  1:3594313/1‑150 (MQ=255)

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

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

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

                                                                                                             gtggCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCgg                                          <  1:3634063/150‑1 (MQ=255)

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

                                                                                                                  cTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAac                                     >  2:1888629/1‑150 (MQ=255)

                                                                                                                   ttGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAaca                                    <  1:3467484/150‑1 (MQ=255)

                                                                                                                     ggTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCCGCCAACAcg                                  >  1:2435249/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                                                                   ttcttcCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACAt                    >  2:2474265/1‑150 (MQ=255)

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

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

                                                                                                                                       tcCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGc                <  2:2661693/150‑1 (MQ=255)

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

                                                                                                                                         cTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTg              <  2:598338/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                                                gccgGTGGTGGCGCCTTGGGCGCCATTTCCGGCGGGCTGATTTGTTACGCCCAGGAAGGTGACGAAGATGGTGAAGGCATTTTAGACCGGCGCGATCACTGCCCCCAGACCCCGGCCAACACGGCGGTTGAACACATTGGCTGCCCACTg       >  1:469636/1‑150 (MQ=255)

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

                                                                                                                                                 ccgGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGa      <  2:2827497/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)

                                                                                                                                                    gtggtgGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCTATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAAc   <  1:3695561/150‑1 (MQ=255)

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

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

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

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

TTGGTGGCTGCCAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACA  >  NC_002947/1284884‑1285182

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