BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,285,045	G→A	5.9%	L90L (TTG→TTA) ‡	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,045	0	G	A	5.9%	675.0 / 21.4	219	L90L (TTG→TTA) ‡	PP_1122	OmpA family protein
Reads supporting (aligned to +/- strand): ref base G (100/103); new base A (9/4); total (112/107)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

AGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTACCCGGCCGC  >  NC_002947/1284896‑1285194

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

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

 gCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCGGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTgg                                                                                                                                                      <  1:4261597/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

                     aCGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTg                                                                                                                                  <  1:4140254/150‑1 (MQ=255)

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

                      cGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGa                                                                                                                                 >  1:2390446/1‑150 (MQ=255)

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

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

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

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

                         acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCAACAGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCCTTGACGGCGGGCTGAttt                                                                                                                              >  1:4231791/1‑150 (MQ=255)

                           aGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCCTTGCCGGCGGGCTGATTTGt                                                                                                                            >  1:58349/1‑150 (MQ=255)

                             cGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTa                                                                                                                          >  2:2000364/1‑150 (MQ=255)

                                  gCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGccc                                                                                                                     >  1:440432/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                    aaaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGGTGTTTTGTTAAGCCCAGGACGGTGACGAAGATg                                                                                                   >  1:692994/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                                 cATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCGTGGGCGCCATTGCCGGCGGGCTGATTTGTTAAGCCCAGGACGGTGACGAAGATGGTGACGGCAtttt                                                                                      >  1:3790383/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)

                                                                           ctgctgGGCGGCTTGGTCGGTGGTGACCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGgcgc                                                                            <  1:311549/150‑1 (MQ=255)

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

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

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

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

                                                                                  gcggcTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACt                                                                     <  1:599918/150‑1 (MQ=255)

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

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

                                                                                       ttGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGcccc                                                                <  1:2203130/150‑1 (MQ=255)

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

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

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

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

                                                                                                ggtggCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCg                                                       <  1:1329872/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)

                                                                                                         ggTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACAcg                                              <  1:4481053/150‑1 (MQ=255)

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

                                                                                                             ccATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTTGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTAGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACAcggcgg                                          >  1:2836117/1‑150 (MQ=255)

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

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

                                                                                                              cATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCGGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTAGACCGGCGCGATCACTGCCCCCAGACCCCGGCCAACACGGCGGt                                         >  1:2287484/1‑150 (MQ=255)

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

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

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

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

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

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

                                                                                                                    gagTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACca                                   <  1:4088040/150‑1 (MQ=255)

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

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

                                                                                                                         cttcCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATgg                              <  1:3478307/150‑1 (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:877381/150‑1 (MQ=255)

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

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

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

                                                                                                                                   cgccgGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACt                    <  1:3280924/150‑1 (MQ=255)

                                                                                                                                   cgccgGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACt                    <  1:996236/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                                              gcgcCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTAcc         >  1:2218879/1‑150 (MQ=255)

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

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

AGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTACCCGGCCGC  >  NC_002947/1284896‑1285194

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