Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,936,766	Δ1 bp	12.1%	coding (1185/1776 nt)	fucI →	L‑fucose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,936,766	0	T	.	12.1%	89.5 / 7.4	33	coding (1185/1776 nt)	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base T (13/15);  new base . (2/2);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

ATCAATATCCCAATGGTGACACCGCCGAAGCGATCCTCAACAGTTCATTTGACTGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAA  >  NC_000913/2936501‑2937027                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                      aTCATTATCCCAATGGTGACACCGCCGAAGCGATCCTCAACAGTTCATTTGACTGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCCCCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTg                                                                                                                                                                                                                                                              >  2:227340/1‑275 (MQ=255) aTCAATATCCCAATGGTGACACCGCCGAAGCGATCCTCAACAGTTCATTTGACTGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACCCGGCATCATCCATTTGATCAACTCCGGTTC‑GCTGCGCTgg                                                                                                                                                                                                                                                             >  2:211980/1‑275 (MQ=255) aTCAATATCCCAATGGTGACACCGCCGAAGCGATCCTCAACAGTTCATTTGACTGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTC‑GCTGCGCTgg                                                                                                                                                                                                                                                             >  2:106273/1‑275 (MQ=255) aTCAATATCCCAATGGTGACACCGCCGAAGCGATCCTCAACAGTTCATTTGACTGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCGCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACCGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGCATGGCTGAAGACGGCATCATCCATTTGATCAACCCCGGTTCTGCTGCGCTg                                                                                                                                                                                                                                                              >  2:211487/1‑275 (MQ=255)      tatCCCAATGGTGACACCGCCGAAGCGATCCTCACCAGTTCATTTGACTGGAATGGCGGGCGCGAACCCGTTGTCGTGGCGACCGAAAACGACAGTCTTTACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGGGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACgg                                                                                                                                                                                                                                                         <  2:228191/275‑1 (MQ=255)      tatCCCAATGGTGACACCGCCGAAGCGAGCCTCAACAGTTCATTGGACTGGAATGGCGGGCGCGAACCCGTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGGGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACgg                                                                                                                                                                                                                                                         <  1:178810/275‑1 (MQ=255)              tgggtacgccgccaAAGCGATCCCCAACAGTTCATTTGACTGGAAGGGCGTGCGCGAACCCTTTGTCGTGGCCACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTa                                                                                                                                                                                                                                                 <  1:159624/267‑1 (MQ=255)                gTGACACCGCCGAAGCGATCCTCAACAGTTCATTTGACTGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTAACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAACCTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTCGTGCGCTGGCCGGTTCCTGTaaa                                                                                                                                                                                                                                               >  1:202653/1‑275 (MQ=255)                   acacCGCCGAAGCGATCCTCAACAGTTCATTTGACTGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTgcgc                                                                                                                                                                                                                                                                >  2:126975/1‑255 (MQ=255)                     cccgccgAATCGATCCTCAACAGTTCCTTTGACGGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACGTGCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAAcg                                                                                                                                                                                                                                          <  1:203432/274‑1 (MQ=255)                       cgccgAAGCGATCCTCAACAGTTCATTTGACTGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGa                                                                                                                                                                                                                                                            <  1:126974/255‑1 (MQ=255)                                            ttCATTTGACTGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGTCACCGCTCAGGTATTTGCCGATGTGCGTACATACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCGTAAACTGGATGGACTGGCAGAACACGGCATCATGCACTGGATCAACTCCGGTTCTGCTGCGCTGGACGGTTACTGTAAACAACGCGACCGCGCAGGTAACCcgacga                                                                                                                                                                                                                   >  2:46577/1‑275 (MQ=255)                                                                        cccTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCCACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCATAACACGGCATCATCCATTTGATCAACGCCGGTTCTGCTGGGCTGGAAGGTTCCTGTAAACAACGCGACAGCGAAGTTAACCCGACGATGAGGCCACACGGGGAAATCTAtcatca                                                                                                                                                                                       >  2:49389/1‑275 (MQ=255)                                                                                             aacaacacgcttaaagaggTGGCAATGCTGATGGTTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGGGCGTACCTACTGGTAACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTcgc                                                                                                                                                                  <  1:233621/257‑1 (MQ=255)                                                                                                                 ggCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTgg                                                                                                                                                                                                                                                             >  2:72940/1‑164 (MQ=255)                                                                                                                     aTGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACgg                                                                                                                                                                                                                                                         <  1:72939/164‑1 (MQ=255)                                                                                                                      tGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTC‑GCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATcc                                                                                                                                        <  1:211979/275‑1 (MQ=255)                                                                                                                      tGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTC‑GCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATcc                                                                                                                                        <  1:106272/275‑1 (MQ=255)                                                                                                                                            cGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGCACTCTCTCAGCAAGAGCCTGACGCTTGCCTCCCCCCTACCGAATGGTGCCCTGCGATCCCCGCATCCTTCCGCtgcggc                                                                                                                    >  1:162740/1‑274 (MQ=255)                                                                                                                                                              tgctagtgcgcagacggactgtaaccaagaagaAATTTAGAGTGGAACGGGGCCTAAACGGGGTGGCCTGGCAGAACGCTGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGGAAACAACGCGACAGCGAAGGTAACCCGGCGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTcc                                                                                                 <  1:11587/248‑1 (MQ=255)                                                                                                                                                                          aaCCTACTGGTAACCAAAAGCAATTGGGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAggcg                                                                                     <  1:227885/274‑1 (MQ=255)                                                                                                                                                                                  ggTCAGCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAACGTAACCCGACGATGAAGCCACACTGGGAAATCTCACAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATAGTGCCAGGCGATCCACGACTACTTCAGTGGCGGCGGCTACTATTCCCTCTTCCTTACCGCAGGCGGCGTCCCg                                                                             >  2:172559/1‑275 (MQ=255)                                                                                                                                                                                                      cGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCGACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCCACCAGTTCACCATCACCCGTGTCa                                                          >  2:5873/1‑274 (MQ=255)                                                                                                                                                                                                               tgtctggaacgggttgggcttgacgAACAGGGCTACATCCGTGGGAACAATTCCGGTTCGCGTGCGCTGGACGTGTCCTGTAAACAACGCGACAGCCAAGGGAGCCCGACGATGAAGCCACACTGGGAAATCTCTCGGCACGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCaa                                                <  1:44439/251‑1 (MQ=255)                                                                                                                                                                                                                      aaccgtagggaATGGCAGAACCCGGCATATTCCATGTGAGCAACACCGTTTCTTGTGCGCTGTGCGGTTCCTGTAAACCACTCTACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCGCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGCATTGTGCCCGTCGATCCACGAATATTTCCGTGGCGGCGGTTACTCGTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCt                                          <  1:211486/266‑1 (MQ=255)                                                                                                                                                                                                                      aaCTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCt                                          <  2:42317/274‑1 (MQ=255)                                                                                                                                                                                                                          ggatggaCGGGCAGAACACGGCATCATCCATTTGATCAACTCGGGTTCTGCTGCGCGGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGAc                                     <  1:49389/275‑1 (MQ=255)                                                                                                                                                                                                                               gaCGGGCAGAACACGGCATCATCCATTTGAGCAACTCCGGTTCTGCGGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGt                                 <  2:202654/274‑1 (MQ=255)                                                                                                                                                                                                                                   ggCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCAACATGACCAGTGTCAACATCATCAAAGGCCTGGGCCCGGTACTg                             >  1:121612/1‑274 (MQ=255)                                                                                                                                                                                                                                                    tcCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACACGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTgg           >  1:51912/1‑275 (MQ=255)                                                                                                                                                                                                                                                            aTCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAg                                              >  1:161057/1‑232 (MQ=255)                                                                                                                                                                                                                                                             tCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACACCGCGACAGCGAAGGTAGCCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGaa  <  1:46577/275‑1 (MQ=255)                                                                                                                                                                                                                                                                aCTCCGGTTCTGCTGCGCTTGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCt                                          <  2:161058/232‑1 (MQ=255)                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                      ATCAATATCCCAATGGTGACACCGCCGAAGCGATCCTCAACAGTTCATTTGACTGGAATGGCGTGCGCGAACCCTTTGTCGTGGCGACCGAAAACGACAGTCTTAACGGCGTGGCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAA  >  NC_000913/2936501‑2937027

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