Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,486,576	0	C	G	42.9%	87.6 / 19.4	56	A120G (GCG→GGG)	lptF	lipopolysaccharide export ABC permease
Reads supporting (aligned to +/- strand):  ref base C (9/23);  new base G (24/0);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 9.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCA  >  NC_000913/4486433‑4486725                                                                                                                                                |                                                                                                                                                       gctgATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGGGGGGcc                                                                                                                                                  >  2:900367/1‑150 (MQ=255) gctgATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGAcc                                                                                                                                                  <  1:985801/150‑1 (MQ=255)     atgaCGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAACGTGATGTGGGGGGGGCCgggg                                                                                                                                              >  2:840281/1‑150 (MQ=255)        aCGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCAt                                                                                                                                          <  1:420484/151‑1 (MQ=255)                 aaaCTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAAAGGGATGTGGGGGGGGGCGCGGGCATcgcggcct                                                                                                                                  >  1:81533/1‑146 (MQ=255)                  aaCTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGTGGTGTGGGGGGGGCCGGGGGCATcgcggccac                                                                                                                                 >  2:752908/1‑145 (MQ=255)                        ataCCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGTTCGCGACCCTGGTCATCACTTCATCAGTATg                                                                                                                           <  1:228877/150‑1 (MQ=255)                        ataCCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATg                                                                                                                           <  1:982064/150‑1 (MQ=255)                         taCCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGGGGTTTGGGCGGGGGCGTGGGCATcgcggcctccgggt                                                                                                                            >  2:565629/1‑138 (MQ=255)                          aCCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGaa                                                                                                                         <  1:491413/150‑1 (MQ=255)                               aaGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGGGATGGGGGGGGGGCCGGGGGGCTCGCGTCGACAGGGTGAAggggt                                                                                                                    >  1:476660/1‑146 (MQ=255)                                  tGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGGTGTGGGGGGGGCCGGGGGCATCGCGGCCTCAGGGTGGAggggt                                                                                                                    >  2:354257/1‑143 (MQ=255)                                   gAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGATGTGGGGGGGGCCGGGGGCATCGCGGCCTCCGGATGAAGAGtt                                                                                                                    >  2:901632/1‑146 (MQ=255)                                     aaTTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGGTGTGGGGGGGGGCGGGGGCATCGCGTCATCAGGGTGAAatggttgga                                                                                                                >  2:304454/1‑142 (MQ=255)                                     aaTTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTTACCTGATGTGGGGGGGGCCGTGGGCCTCGCGGCATCCGGGTGAAGTGTTTGCAGa                                                                                                              >  2:287908/1‑150 (MQ=255)                                     aaTTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAACGTGATGTGGGGGGGACCGTGGGCGTCCCCTCATAAGGagggaagggt                                                                                                                    >  1:621869/1‑135 (MQ=255)                                        tACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAACGGGATGTGGGGGGGGCCGGGGTCATCGCGTCATCAGGGTGAAGGGTTAGCAGaa                                                                                                             >  2:644728/1‑148 (MQ=255)                                         aCGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGGGGGGCCGCGGGCATCGCGTCCTCAGGGTGAAGGGTTTGCAGaaggg                                                                                                          >  2:392767/1‑148 (MQ=255)                                         aCGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGaagcga                                                                                                         <  2:839196/151‑1 (MQ=255)                                           ggTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGaagcgaaa                                                                                                       <  2:642435/151‑1 (MQ=255)                                            gTAATGCATGCCTGCGGCCTGAGTAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGaagcgaaat                                                                                                      <  1:222673/151‑2 (MQ=255)                                            gTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGGGGTGTGGGGGGGGGCGGGGGCATCGCGGCATCAGGGTGAAGGGTTTGCAGaa                                                                                                             >  2:285050/1‑144 (MQ=255)                                              aaTGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGaagcgaaagc                                                                                                     <  1:213731/150‑1 (MQ=255)                                               atgcatgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGTGGTGTGGGGGGGGCCGTGGGCATCGCGGCCTCAGGGTGAAGGGGTAGCAGaagaggaa                                                                                                       >  2:258725/1‑142 (MQ=255)                                               atgcatgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGCTGTGGGGGGGGCCGCGGGCATCCCGTCCTCCGGATGAAGGGGTAGCAGaagcggaa                                                                                                       >  2:769781/1‑147 (MQ=255)                                                 gcatgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAagcgaaagcgaa                                                                                                  <  1:673342/150‑1 (MQ=255)                                                  catgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGGCCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAAc                                                                                                 <  1:930167/150‑1 (MQ=255)                                                  catgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGGTGTGGGGGGGGCCGGGGGGATCGCGGCATCAGGAGGGAGGGTTAGCAGAagcgaaagagaaac                                                                                                >  2:424803/1‑149 (MQ=255)                                                    tgcCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGATGTGGGGGGGGGCGGGGGCATCGCGCCCTCCGGGGGAAGAGGTAGCAGAAGCGAAAAAGAAccc                                                                                               >  2:863905/1‑150 (MQ=255)                                                             ccTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCCGAAGCGAAAGCGAACCCTGGCATggcg                                                                                     <  2:889978/151‑1 (MQ=255)                                                              cTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTACGGGATGTGGGGGGGGGCGGGGGCATCGCGGCCTCAGGGTGAAGGGGTTGGAGAAACGGAAAAGAACCCCGGCATggcgg                                                                                    >  2:613882/1‑151 (MQ=255)                                                                gAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGGCCGGGGGCATCGCGGCATCAGGGTGAAGGGTTTGCAGAAACGAAAGCGAAACCCGGCCTgg                                                                                       >  2:900504/1‑146 (MQ=255)                                                                 aGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGGGGGGCCGGGGGGCTCGCGGCCTCAGGATGAAGGGTTTGCAGAAGCGAAAACGAAACCCGGAATGgcggcg                                                                                  >  2:779662/1‑150 (MQ=255)                                                                     aaGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTTAAGTGGTTTGGGGGGGGCCGTGGGGATCGCGGCCTCAGGGTGGAGGGGTTGGAGAAAAGAAAACGAACACCTGCATGGGGGGGCggg                                                                              >  2:854027/1‑150 (MQ=255)                                                                        cGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑Ca                                                                         <  2:993759/151‑1 (MQ=255)                                                                          gTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAg                                                                       <  1:119269/151‑1 (MQ=255)                                                                               ggTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGc                                                                    <  1:165881/149‑1 (MQ=255)                                                                                    aaGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGGGATGTGGGGGGGACCGTGGGCATCGCGGCATCAGGGTGGAGGGGTTGCAGAAACGAAAGCGAACCCCCGCATGGGGGGGCTGGCG‑GCAGGGGAAtt                                                                >  2:308006/1‑148 (MQ=255)                                                                                     aGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGGGGGGGCGGGGGCATCGCGGCATCAGGGTGAAGGGTTAGCAGAAGCGGAAGCGAAACCTGGCATGGCGGCGCCGGCG‑CCAGGGGaaattca                                                             >  1:942237/1‑145 (MQ=255)                                                                                      gCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGc                                                           <  1:752741/151‑1 (MQ=255)                                                                                            aTGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGTGATGTGGGCGGGGCCGCGGGCATCGCGGCATCAGGAGGAAGTGTTAGCAGAAGAGAAAGCGAAACCCGGCATGGGGGGGCGGGGG‑CGAGGGGAATTCCAGCAAGCg                                                      >  2:38954/1‑150 (MQ=255)                                                                                             tGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGGTAAAGTGCTGTGGGCGGGGCCCTGGGGATCGCGGCATCAGGGTGGGGGGGTAGCAGAAGGGAAAGAGAACCCCGGCCGGGCGGGGCTGGCG‑CAAGGGGAATTTCa                                                             >  1:286603/1‑142 (MQ=255)                                                                                                 ccTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTaa                                                 <  1:982996/150‑1 (MQ=255)                                                                                                    tGCGGTATTCACGGCAATCGTAGCGGCGGTTAAAGCGATGTGGGGGGGGCCGTGGGCATCGCGTCGTCAGGGTGAAGGGGTTGCAGAAACGAAAACGAAACCCCGCATGGGGGCGCCGGCG‑CCAGGGGGAATTCAGCCAGCGACTAAAgg                                              >  2:58711/1‑150 (MQ=255)                                                                                                         tATTCACGGCAATCGTAGCGGCGGTTAACGTGGTTTGGGGGGGACCCCGGTGATCCCGTCATCCGGGTGAAGGGGTAGCAGAAGGGGAAACGAACACCGGCCTGGCGGGGCCGGGGCCAGGGGGAATTCC‑GCAAACGACTaa                                                 >  1:502581/1‑142 (MQ=255)                                                                                                                ggCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTg                                 <  2:441941/151‑1 (MQ=255)                                                                                                                         aGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGaa                        <  2:286550/151‑1 (MQ=255)                                                                                                                           cggcggTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAg                      <  2:354197/151‑1 (MQ=255)                                                                                                                            ggcggTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGc                     >  2:165112/1‑151 (MQ=255)                                                                                                                                gTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTg                 >  2:845390/1‑151 (MQ=255)                                                                                                                                  tAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGAc               >  2:137456/1‑151 (MQ=255)                                                                                                                                   aaCGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGAc               <  2:81541/150‑1 (MQ=255)                                                                                                                                   aaCGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACg              >  1:987750/1‑151 (MQ=255)                                                                                                                                      gTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCa           <  2:829367/151‑1 (MQ=255)                                                                                                                                        gATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGc         >  1:88629/1‑151 (MQ=255)                                                                                                                                               gCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCa  <  1:565513/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                       GCTGATGACGCTGGGCAAACTGTATACCGAAAGTGAAATTACGGTAATGCATGCCTGCGGCCTGAGCAAAGCGGTTCTGGTGAAAGCGGCAATGATCCTTGCGGTATTCACGGCAATCGTAGCGGCGGTTAACGTGATGTGGGCGGGACCGTGGTCATCGCGTCATCAGGATGAAGTGTTAGCAGAAGCGAAAGCGAACCCTGGCATGGCGGCGCTGGCG‑CAAGGGCAATTCCAGCAAGCGACTAATGGCAGCTCGGTGCTGTTCATCGAAAGCGTTGACGGCAGCGATTTCA  >  NC_000913/4486433‑4486725

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