Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,424,381	G→C	7.0%	intergenic (‑206/‑9)	PP_2124 ← / → yegS	glycosyl transferase/lipid kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,424,381	0	G	C	7.0%	179.9 / 2.6	57	intergenic (‑206/‑9)	PP_2124/yegS	glycosyl transferase/lipid kinase
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base C (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

GGGGAGCCACTGCCGGGCGCGTCCATGGCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAGCGGCTGGGT  >  NC_002947/2424233‑2424529                                                                                                                                                     |                                                                                                                                                     ggggAGCCACTGCCGGGCGCGTCCATGGCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCgg                                                                                                                                                      >  4:375104/1‑149 (MQ=255) ggggAGCCACTGCCGGGCGCGTCCATGGCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCgg                                                                                                                                                      >  8:226689/1‑149 (MQ=255)   ggAGCCACTGCCGGGCGCGTCCATGGCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGc                                                                                                                                                    <  8:248877/149‑1 (MQ=255)             ccGGGCGCGTCCATGGCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTgg                                                                                                                                          >  6:389222/1‑149 (MQ=255)              cGGGCGCGTCCATGGCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGt                                                                                                                                         >  6:160495/1‑149 (MQ=255)                   gcgTCCATGGCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTc                                                                                                                                        >  5:154943/1‑145 (MQ=255)                      tCCATGGCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAg                                                                                                                                 <  5:389222/149‑1 (MQ=255)                            gCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGggtg                                                                                                                           <  8:17437/149‑1 (MQ=255)                            gCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGggtg                                                                                                                           <  6:368000/149‑1 (MQ=255)                            gCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCAAGGCCATGCAACGTTCCGGGCAGGCGAGGGGTCGAGATAGGGggtg                                                                                                                           >  7:378340/1‑149 (MQ=255)                                  ggggTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCAt                                                                                                                     >  5:437876/1‑149 (MQ=255)                                  ggggTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCAt                                                                                                                     <  6:346065/149‑1 (MQ=255)                                   gggTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGCCCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCAtt                                                                                                                    >  3:131884/1‑149 (MQ=255)                                    ggTATTTGCTTAGCATAGCCGGGATTGCCAAGTGTGTAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCAttt                                                                                                                   <  4:131884/149‑1 (MQ=255)                                     ggaTTTGCTGAGGATAGCCGGGCTTGGCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTGCAGGCCATGCAACGATCCGCGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCAttt                                                                                                                   >  6:1458/3‑148 (MQ=255)                                     gTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCa                                                                                                                                                   >  4:68950/1‑116 (MQ=255)                                     gTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCa                                                                                                                                                   <  3:68950/116‑1 (MQ=255)                                     gTATTTGCTGAGCATAGCCGGGATTGCCAAGGGTGGAAAAAATGACCGGCGATGCGATAAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGCGCAGGCCAGTGGTCGTGATAGGGGGTGGTCCAttt                                                                                                                   >  2:140168/1‑148 (MQ=255)                                        tttGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACg                                                                                                               <  6:437876/149‑1 (MQ=255)                                               aGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGAcgc                                                                                                        >  1:145417/1‑149 (MQ=255)                                                 cATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGAcgcg                                                                                                       >  8:278184/1‑148 (MQ=255)                                                     gCCGGGATTGCCAAGTGTGGAAAAATTGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACAt                                                                                                  >  2:279343/1‑149 (MQ=255)                                                       cGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATg                                                                                                 >  7:207804/1‑148 (MQ=255)                                                          gATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAg                                                                                             <  3:318106/149‑1 (MQ=255)                                                                 aaGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGgcgc                                                                                        >  2:392664/1‑147 (MQ=255)                                                                        gAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACg                                                                                                               >  1:365634/1‑117 (MQ=255)                                                                        gAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACg                                                                                                               <  2:365634/117‑1 (MQ=255)                                                                        gAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATg                                                                               >  7:206475/1‑149 (MQ=255)                                                                             aaTGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGt                                                                          <  8:43085/149‑1 (MQ=255)                                                                                     gcgatgcgatCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGtg                                                                         >  2:24008/1‑142 (MQ=255)                                                                                     gcgatgcgatCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGtg                                                                         <  1:24008/142‑1 (MQ=255)                                                                                        atgcgatCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCa                                                               <  3:375104/149‑1 (MQ=255)                                                                                             atCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCAGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAg                                                          >  4:151684/1‑149 (MQ=255)                                                                                                 gTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGa                                                      >  2:374169/1‑149 (MQ=255)                                                                                                  ttGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGAt                                                     <  2:281855/149‑1 (MQ=255)                                                                                                    gCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATg                                                    >  5:16544/1‑148 (MQ=255)                                                                                                      aGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAAc                                                 <  6:11103/149‑1 (MQ=255)                                                                                                      aGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAAc                                                 <  5:13293/149‑1 (MQ=255)                                                                                                       ggCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACg                                                <  5:160495/149‑1 (MQ=255)                                                                                                              cgcCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGt                                         <  8:207804/149‑1 (MQ=255)                                                                                                               gcCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTg                                        >  5:440129/1‑149 (MQ=255)                                                                                                                ccGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTgc                                       <  8:378340/149‑1 (MQ=255)                                                                                                                   gCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCa                                    <  1:279343/149‑1 (MQ=255)                                                                                                                     ccAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGGGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAgc                                  <  6:442323/149‑1 (MQ=255)                                                                                                                      cAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAgcg                                 >  2:291800/1‑149 (MQ=255)                                                                                                                      cAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAgcg                                 >  7:346251/1‑149 (MQ=255)                                                                                                                                 ccAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGAc                      >  4:134006/1‑149 (MQ=255)                                                                                                                                 ccAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGAc                      >  7:112147/1‑149 (MQ=255)                                                                                                                                 ccAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGAc                      >  8:352630/1‑149 (MQ=255)                                                                                                                                  cAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGAc                      <  6:16544/148‑1 (MQ=255)                                                                                                                                  cAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGAc                      <  2:145417/148‑1 (MQ=255)                                                                                                                                            aaCGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAg           >  5:274666/1‑149 (MQ=255)                                                                                                                                            aaCGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAg           >  5:83161/1‑149 (MQ=255)                                                                                                                                               gATCCGCGCAGGCCAGTGGTCGTGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAACGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAGCgg        <  1:140168/149‑1 (MQ=255)                                                                                                                                               gATCCGCGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAGCgg        <  5:1458/149‑1 (MQ=255)                                                                                                                                                 tCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAGCGGCt      >  2:267038/1‑149 (MQ=255)                                                                                                                                                   cGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAGCGGCTgg    <  7:352630/149‑1 (MQ=255)                                                                                                                                                   cGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAGCGGCTgg    <  1:267038/149‑1 (MQ=255)                                                                                                                                                   cGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAGCGGCTgg    <  6:440129/149‑1 (MQ=255)                                                                                                                                                     ggCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAGCGGCTGGGt  <  7:226689/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGGGAGCCACTGCCGGGCGCGTCCATGGCAATAGGGGTATTTGCTGAGCATAGCCGGGATTGCCAAGTGTGGAAAAAATGACCGGCGATGCGATCAGTTGCAGGCCTGGCGCCGGCCCAATGTTTTTTCCAGGCCATGCAACGATCCGGGCAGGCCAGTGGTCGAGATAGGGGGTGGTCCATTTTACGGAGACGCGCACATGCAAGGGCGCAAGGCAATGCTGGTGCTGCACGGCAAGCAGGCGATGAACGAAGACGTGCGCAGCGCCGTGGGCGACCTGCGCGATAGCGGCTGGGT  >  NC_002947/2424233‑2424529

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑