Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,086,153	C→T	G123R (GGG→AGG)	pgaD ←	poly‑N‑acetyl‑D‑glucosamine synthase subunit PgaD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,086,153	0	C	T	97.1%	96.9 / ‑4.0	34	G123R (GGG→AGG)	pgaD	poly‑N‑acetyl‑D‑glucosamine synthase subunit PgaD
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (13/20);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGCCCAGACAATTAACACGACGGCATTCGCCACTGCCAGCAAAAAATAGAACTGAAGTCGGCTTCTGGCCTCGCTTTGCCAGTAAT  >  NC_000913/1085907‑1086389                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             cAATGGGATGTGCCTCGCGGCGTGTGCTCCGTATTAAGTGACGCGCTGGAACAGTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGTTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGTCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTctgg                                                                                                                                                                                                                                          <  2:375504/251‑1 (MQ=255)          ttGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGaa                                                                                                                                                                                                                                  <  1:416508/250‑1 (MQ=255)            gccgcgcgGCGTGTGCTCCTGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAgt                                                                                                                                                                                                                                <  1:223678/246‑1 (MQ=255)                     ggtgcgCTCCGGATTAAGTGGTGCTCTGGATAGTTTCGAAGAAGATGAAATGGCTGGGATCGTCCGGTTCTGTAAAGAGGAATGCGACCGTTCGGCACTTTGTCAATGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGATCCGGGCTTTGTCTTGGGTGTTTATGCCCGGACTAGCGCTTGTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCa                                                                                                                                                                                                                       <  2:467601/245‑1 (MQ=255)                               ggATTAAGTGACGCGCTGGAACGTTTCGCAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCtt                                                                                                                                                                                                            >  1:409124/1‑251 (MQ=255)                                        gACGCGCTGGAACGTTGCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGt                                                                                                                                                                                                   <  1:68408/251‑1 (MQ=255)                                                  aaCGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGAt                                                                                                                                                                                            >  2:262866/1‑248 (MQ=255)                                                                 aTGATATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTg                                                                                                                                                                          <  2:405790/251‑1 (MQ=255)                                                                            gggATCGTCCGGTTCGGTAAAGAGGACCGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGctctct                                                                                                                                                               <  1:49502/251‑1 (MQ=255)                                                                                     cGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCt                                                                                                                                                      >  2:88153/1‑251 (MQ=255)                                                                                         tCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTTCAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGTCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTgc                                                                                                                                                   <  2:65546/250‑1 (MQ=255)                                                                                                            aCGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTg                                                                                                                                <  2:39282/250‑1 (MQ=255)                                                                                                                    cACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGGGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTTGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGAtgct                                                                                                                       <  2:191032/251‑1 (MQ=255)                                                                                                                    cACTTTGCCAACGTACGCTGCATGACTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTAATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGAtgct                                                                                                                       >  1:360032/1‑251 (MQ=255)                                                                                                                          gCCAACGTACGCTGCATGCTTACAGTTATGGTTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCGTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTTCTAAGCTCTCTGCATATTCGTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGtt                                                                                                                 <  2:480136/251‑1 (MQ=255)                                                                                                                                  aCGCTGCATGGCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGaaa                                                                                                          <  2:190920/250‑1 (MQ=255)                                                                                                                                       gCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGGTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAg                                                                                                    <  1:337961/251‑1 (MQ=255)                                                                                                                                            ccTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTAt                                                                                               >  1:118856/1‑251 (MQ=255)                                                                                                                                               aCAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAg                                                                                                                                    <  1:214/211‑1 (MQ=255)                                                                                                                                               aCAGTTAAGTGAGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAg                                                                                                                                    >  2:214/1‑211 (MQ=255)                                                                                                                                                 ggTTAAGTGTGTTATCGGTGCAGAGGCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGGGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGGATTGCTAAGCGCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTAc                                                                                          <  2:251707/250‑1 (MQ=255)                                                                                                                                                       gtgtgtTATCGGTGCAGAGCCCGGGCGAACCGGGCTGTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTGATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGcc                                                                                    <  1:37671/251‑1 (MQ=255)                                                                                                                                                       gtgtgtTAGCGGTGCAGAGCCCGGGCGAACCGGGCTTTTTCTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACATCCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGGATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGcc                                                                                    <  1:470523/251‑1 (MQ=255)                                                                                                                                                           gtTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGccc                                                                                   >  2:227557/1‑248 (MQ=255)                                                                                                                                                              aTCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCGCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAAAGCAGCTTATTGTACAGCGCCCAGACa                                                                              >  2:334796/1‑250 (MQ=255)                                                                                                                                                                  gTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGCCCAGACAATTaa                                                                         <  1:227557/251‑1 (MQ=255)                                                                                                                                                                                           ttGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGATGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGCCCAGACAATTAACACGACGGCATTCGCCACTGCCa                                                  <  1:320820/249‑1 (MQ=255)                                                                                                                                                                                                   ggTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGCCCAGACAATTAACACGACGGCATTCGCCACTGCCAGCAAAAAAt                                         >  1:354397/1‑250 (MQ=255)                                                                                                                                                                                                       tttATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGCCCAGACAATTAACACGACGGCATTCGCCACTGCCAGCAAAAAATAGAAc                                    >  2:106513/1‑251 (MQ=255)                                                                                                                                                                                                                 gACTAGCGCTTTGTCTGAAACAACCATTTTTATTGTCCTCTGGCTGGGGAAGTGTACGCTCATCCTGTGGCTTTGTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGCCCAGACAATTAACACGACGGCATTCGCCACTGCCAGCAAAAAATAGAACTGAAGTCGGc                          <  1:260729/251‑1 (MQ=255)                                                                                                                                                                                                                       cgcTTTTTCTGAAACAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCCGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGCCCAGACAATTAACACGACGGCATTCGCCACTGCCAGCAAAAAATAGAACTGAAGTCGGCTTCTgg                    >  2:155332/1‑251 (MQ=255)                                                                                                                                                                                                                                  aaaCAACCATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGCCCAGACAATTAACACGACGGCATTCGCCACTGCCAGCAAAAAATAGAACTGAAGTCGGCTTCTGGCCTCGCTTTGc         >  2:110054/1‑251 (MQ=255)                                                                                                                                                                                                                                     gaacCATTTTTATTTGCCTCTGGCTGGAGCAGTGTACGCTCATCCTGAGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTACTGTTTTTGACAACGCAGCTTATTGTACAGCGCCCAGACAATTAACACGACGGCATTCGCCACTGCCAGCCAAAAAGAGAACGGAAGTCGGCTTCTGGCCTCGCTTAGcc        >  2:256025/2‑249 (MQ=255)                                                                                                                                                                                                                                         cATTTTTATTTGCCTCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGCCCAGACAATTAACACGACGGCATTCGCCACTGCCAGCAAAAAATAGAACTGAAGTCGGCTTCTGGCCTCGCTTTGCCAGtaat  <  1:165893/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             CAATGCGATTTGCCTCGCGGCGTGTGCTCCGGATTAAGTGACGCGCTGGAACGTTTCGAAGAAGATGAAATGGTTGGGATCGTCCGGTTCGGTAAAGAGGACTGCGTACGTTCGGCACTTTGCCAACGTACGCTGCATGCCTACAGTTAAGTGTGTTATCGGTGCAGAGCCCGGGCGAACCGGGCTTTGTTTTGGGTGTTTATGCCCGGACTAGCGCTTTTTCTGAAACAACCATTTTTATTTGCCCCTGGCTGGTGAAGTGTACGCTCATCCTGTGGCTTTTTTGTAGTTGCTGATAGAGCTCATCAGGTATTGCTAAGCTCTCTGCATATTCTTGCGGCGTATATTGGTAGGCTGCATGATGCTGCTGTTTTTGAAAACGCAGCTTATTGTACAGCGCCCAGACAATTAACACGACGGCATTCGCCACTGCCAGCAAAAAATAGAACTGAAGTCGGCTTCTGGCCTCGCTTTGCCAGTAAT  >  NC_000913/1085907‑1086389

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