Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,574,051	0	G	T	38.9%	78.7 / 20.2	56	V384V (GTC→GTA)	pqqL	putative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base G (6/27);  new base T (21/0);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.14e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCG  >  NC_000913/1573906‑1574200                                                                                                                                                   |                                                                                                                                                        gTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGAcatc                                                                                                                                                    <  1:652349/151‑1 (MQ=255) gTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTAca                                                                                                                                                      >  1:733017/1‑149 (MQ=255)    tCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTAcatca                                                                                                                                                   >  2:832644/1‑149 (MQ=255)     cGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACatcatc                                                                                                                                                 <  2:1058683/150‑1 (MQ=255)      gAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGa                                                                                                                                               <  2:848611/151‑1 (MQ=255)      gAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTAcatca                                                                                                                                                   >  2:385932/1‑147 (MQ=255)       aaaGTTTATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAg                                                                                                                                              <  2:777680/151‑1 (MQ=255)          gTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GAAGATTTGACATCATCGAGttt                                                                                                                                           >  1:1116928/1‑150 (MQ=255)              atatGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTAAGGCGG‑GTTGAGTTAGcatccacatgttctttt                                                                                                                                       >  2:190137/1‑138 (MQ=255)                  gTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTcag‑ca                                                                                                                                   <  1:968621/150‑1 (MQ=255)                  gTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTcag‑ca                                                                                                                                   <  1:425054/150‑1 (MQ=255)                     tctTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTATATCATCGAGTTTTTCAG‑CAGaa                                                                                                                                >  1:1052671/1‑150 (MQ=255)                      ctTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGaa                                                                                                                                <  1:771986/149‑1 (MQ=255)                       ttCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGTGATTTTAAAATAACGAGGTCTTcagaccaaa                                                                                                                                >  1:858875/1‑143 (MQ=255)                        tCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGaaaa                                                                                                                              <  2:658440/149‑1 (MQ=255)                         ccGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGAGATCATCCAGGTTTTTTG‑CCGAAAAAc                                                                                                                            >  1:816779/1‑150 (MQ=255)                              gACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCt                                                                                                                      <  1:296108/151‑1 (MQ=255)                                 aaGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTACATCAACGAGTTTTTTCG‑CAGAAAAAACATGc                                                                                                                       >  1:1222257/1‑147 (MQ=255)                                  aGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATc                                                                                                                  <  2:1246458/151‑1 (MQ=255)                                   gAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGAATTTTCATCAACGAGTTTTTTAG‑CCGAAAAAACAAGCTg                                                                                                                     >  1:1041501/1‑147 (MQ=255)                                           tattatTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTACATCATCGAGGTCTTCAG‑CAGAAAAAACATGCTGGTCAATGGtt                                                                                                           >  1:111866/1‑149 (MQ=255)                                                   aaTGAGCTGGATGCCAGGCGACTGGTCAGCATTCGTAAATCTCGTTCAGCTTGCTCATCAATCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CTGAAAAACCTTGCTGATCAATGGTTGCCAActct                                                                                                  <  2:226192/150‑1 (MQ=255)                                                      gAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTACATCATCGAGGTTTTCTT‑CAGAAAAACAAAACTGGTCAAAGGTTGTCAActcttc                                                                                                >  1:13375/1‑147 (MQ=255)                                                       aGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTTGGTTTTACATCCTCGAATTCTTCTG‑CAGAAAAAACATGCTGATCAATGGGTTCCAACTCTG‑CCa                                                                                             >  1:705388/1‑150 (MQ=255)                                                             aTGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATg                                                                                       <  1:1023938/150‑1 (MQ=255)                                                               gCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCAt                                                                                    <  1:1177036/151‑1 (MQ=255)                                                                 cAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGAATTTATATCATCGAGGTCTTTTG‑CCGAAAAAACAAGCTGATCAAAGGGTGCCAACTCTT‑CCCTTTATGCAtt                                                                                   >  1:578438/1‑150 (MQ=255)                                                                       aCTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCgcagca                                                                            <  2:862518/151‑1 (MQ=255)                                                                            tCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTTACATCATCGAGGTTTTCAG‑CAGAAAAAACATGCTGATCAATGGGTGCCACCTCTT‑CCATTAATGCATTCGCAGCATcc                                                                         >  1:718022/1‑149 (MQ=255)                                                                            tCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTTCATCATCGAGTTCTTTCA‑CAGAAAAAACAAGCTGATCAAAGGTTGGCAACTCTT‑CCCTTAATTCATTCgcagca                                                                            >  2:787217/1‑146 (MQ=255)                                                                               gCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTATATCATCGAGTTTTTTTA‑CCGAAAAAACCCGCTGATCAAAGGTTTTCAACTCTG‑CCCTTTATGTATTTGCCACATCCTGc                                                                      >  2:70696/1‑149 (MQ=255)                                                                                cATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAACCCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCat                                                                    <  1:124592/150‑1 (MQ=255)                                                                                  tACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGCGTAGATTT‑TCAGCCTCGAGGTTTTTAGCCACAAAAACAAACCGGTTAAAAGGTTCCCAACCCT‑‑CCCTTTATGTATTCACAGCCACCTGCaa                                                                    >  1:922800/1‑147 (MQ=255)                                                                                   aCGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCAtatt                                                                 <  2:647802/150‑1 (MQ=255)                                                                                       aaaTCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCg                                                             <  2:500728/150‑1 (MQ=255)                                                                                           cACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTTACAACATCGAGGTCTTTAG‑CAGAAAAAACATACTGATCAAAGGTTGCCCACTCTT‑CCATTAATTCATTTGCAGCATCCTGCATATTTTTGtctc                                                         >  1:352701/1‑150 (MQ=255)                                                                                             cGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTg                                                      <  2:1073637/151‑1 (MQ=255)                                                                                             cGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTTTATTCTCTAGTTTTTCAG‑CCGAAAAAACCAGCTGGTTAAAGGTTGTCCAATCTGCCCTTTATTGAATTC‑CAACATCCTGCctat                                                                  >  2:148593/1‑139 (MQ=255)                                                                                              gTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTACATAATCGAGTTTTTTTG‑CCGAAAAAACAAGCTGATCAAATGTTTTCAACCCTG‑CCCTTAATGCATTTGCAGCAACCTGCCTATTATTGt                                                            >  2:941953/1‑144 (MQ=255)                                                                                                  aGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCAttt                                                 <  1:752535/151‑1 (MQ=255)                                                                                                        cTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTGTAATCATCGAGTTCTTTAG‑CAGAAAAAACATGATGGTCAATAGTTGCCAACTCTG‑CCCTTAATGCATTAGCCACATCCTGCATATTATCTTCTCGTGCATTTTCACgg                                           >  1:831071/1‑151 (MQ=255)                                                                                                              aaCCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTTATATCACTGCGTTTTTTCG‑CAAAAAAAACATGCTGGACAAAGGGTGCCAACTCTT‑CCCTTAATGCATTCTCACGATCCTGCATATTTTTGTCCCCGGCCTTTTCACAGgaaaa                                      >  2:512982/1‑150 (MQ=255)                                                                                                               aCCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTCCTGCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCGTCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACa                                    <  2:498658/151‑1 (MQ=255)                                                                                                                 cGCATTTTTCAGCCAGGTGAGGCGG‑GGAGATTTTTCATCATCGAGGTTTTTAG‑CCAAAAAAACAAGCTGGTTAATGGGTGGCAACTCCG‑CCATTAATGCATTTGCAGCCTCCCGCATATTATCGTTCCGTGCATTTAAACGGaaaaaaa                                    >  1:13063/1‑147 (MQ=255)                                                                                                                      ttttCAGCCAGGTGAGGCGG‑GGAGATTTTACATCCTCCAGGTTTTCAG‑CCGGAAAAACATGGTGGTCAAAGGGTGGCAACTCCT‑CCCTTAATTCATTTGCAGCATCCCTCATATTATTGTCTCGGGCATTTTCAAGGAAAAAAAACGAccg                             >  2:220060/1‑149 (MQ=255)                                                                                                                        ttCAGCCAGGTGAGGCGG‑GTAGATTTTTCAATATCGAGTTCTTCTG‑CAGAAAAAACCCGCTGGTTAAAGGTGGCCCACACTT‑CCCTTTATTCATTTGCCACCCCCCTCATATTTTCATCTCGTGCCTTTTCAACGaaaaaaaaagaacg                             >  2:798027/1‑140 (MQ=255)                                                                                                                         tCAGCCAGGTGAGGCGG‑GTAGATTTGGAATCAACGGATTCTTCAT‑CAGAAAAAACATGCTTGTTAATGGGTGCCAACTCTT‑CCCTTTAAAAAATCGCAGGATCCTGCATATTATCGTCTCGTTCCATTATACCGAAAAACa                                    >  1:985418/1‑141 (MQ=255)                                                                                                                            gCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATc                       <  1:1206242/151‑1 (MQ=255)                                                                                                                             ccAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACATGGAAAAACAGCGACTGATAATCg                      <  1:46268/151‑1 (MQ=255)                                                                                                                                 gTGAGGCGG‑GTAGATTTGACATCATCGAGTGCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTg                  <  1:722130/151‑1 (MQ=255)                                                                                                                                        gg‑gTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAtttt           <  1:349248/151‑1 (MQ=255)                                                                                                                                         g‑gTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTa          <  2:495498/151‑1 (MQ=255)                                                                                                                                           gTAGATTTGTCATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATGATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTaa         <  2:525411/151‑1 (MQ=255)                                                                                                                                           gTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTaa         <  2:816751/151‑1 (MQ=255)                                                                                                                                              gATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAACATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGc      >  2:741080/1‑151 (MQ=255)                                                                                                                                               aTTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCt     >  1:151963/1‑151 (MQ=255)                                                                                                                                                   gACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcg  <  1:385945/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                        GTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGG‑GTAGATTTGACATCATCGAGTTCTTCAG‑CAGAAAAACCATGCTGATCAATGGTTGCCAACTCTG‑CCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCG  >  NC_000913/1573906‑1574200

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