Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,223,638	T→C	intergenic (‑10/‑190)	iclR ← / → metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,638	0	T	C	97.1%	113.5 / ‑5.2	34	intergenic (‑10/‑190)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (8/25);  minor base A (1/0);  total (9/25)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

CAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAG  >  NC_000913/4223400‑4223860                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   cAACCTGACGCACGAAACCCTTCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAGCCCGGCTTGTTGCGCCAGTTCCGTGATTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATc                                                                                                                                                                                                                        <  1:92361/251‑1 (MQ=255)          gCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGAGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCTACGCGTTAAAGACTGAACCAGTCCAGTCGCTGGGGCGGTGGCAACGGCGGGTGTTTTGTCGCGGTTCGGGGGGAAGGGGGGGACCATGGAAGTCTCCTTTTTTTGTATCgggggaat                                                                                                                                                                                                                >  1:122872/1‑243 (MQ=255)             cGAAACCCTGCTGTTGCAGCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtca‑tt                                                                                                                                                                                                           <  2:250807/251‑1 (MQ=255)                aaCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtca‑ttttc                                                                                                                                                                                                        <  2:354567/251‑1 (MQ=255)                      gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtca‑ttttcattt                                                                                                                                                                                                    <  2:94862/249‑1 (MQ=255)                      gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtca‑ttttcattt                                                                                                                                                                                                    <  1:437694/249‑1 (MQ=255)                           tGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTg                                                                                                                                                                                              <  1:387503/250‑1 (MQ=255)                           tGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTg                                                                                                                                                                                              <  1:156175/250‑1 (MQ=255)                           ggCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCTCGTGTCGCTGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGt                                                                                                                                                                                             <  1:215789/250‑1 (MQ=255)                                cGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCt                                                                                                                                                                                        <  1:11429/251‑1 (MQ=255)                                        gCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAt                                                                                                                                                                                <  1:6636/251‑1 (MQ=255)                                            gCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGtt                                                                                                                                                                           <  1:123987/251‑1 (MQ=255)                                               gTGGGTCGTTGAATTGGTTAAGCCTGCTTGTTGCGCCAGTTTCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATGTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGGCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑GTTTCAGTTTTATTGTTAGCTAATCCAAT‑AGTTact                                                                                                                                                                        <  2:11632/251‑1 (MQ=255)                                                     cGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGGGCGACCATGAAAGTCTCCTTTTTCTTTTTCGTGGGAATAC‑TTTTCATTTTTTTTTTTATCTAATGCCAT‑AGTTactggactg                                                                                                                                                                  >  1:167317/1‑251 (MQ=255)                                                        ggAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATc                                                                                                                                                               <  1:138337/251‑1 (MQ=255)                                                               ggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGGATGGGGGCGGCCATGGCAGTCTCCTTTTTTTGTATCGTGGGAATCA‑TTTTCATTTTTATTGTTAGGTTATTCAAT‑AGGTACTTAACTGAttcggtg                                                                                                                                                          >  1:46747/1‑243 (MQ=255)                                                                        gCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGAGGAG‑TTAATGTTg                                                                                                                                              <  1:148521/251‑1 (MQ=255)                                                                                         cGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCCCCGTTCGCGGGGATGGGGGGGGACAAGACAGCATCCTTTTTTTTTGTTGTGGGAAACA‑TTTTCATTTTTTTTTTTATGTAATGTCAT‑AGTTACTGCAATGATACGATTGG‑TTAATTTTG‑AACAAAACTCATGTTgg                                                                                                                            >  2:75000/1‑250 (MQ=255)                                                                                            gAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCgtg                                                                                                                         <  2:167317/251‑1 (MQ=255)                                                                                                     ccTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCt                                                                                                                <  2:200684/250‑1 (MQ=255)                                                                                                        gCCATTGGATTCGGCAATCCCCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCtttt                                                                                                             <  1:16868/251‑1 (MQ=255)                                                                                                           aTTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc                                                                                                          <  1:360191/251‑1 (MQ=255)                                                                                                           aTTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCCCGGTTCCCGGGAATGGGGGCGGGCATGAACGACCCCTTTTTCTGTTTCGGGGGAATCA‑TTTTTCTTTTTTTTTTTAGATTATGCAAT‑AGTTACTGAAATGATCCGATGGGATTAAT‑TTGTAACAA‑TCTCATGTTTCGttgtggggggtttttcc                                                                                                          >  1:288166/1‑240 (MQ=255)                                                                                                              ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGGAATGGGGGCGACCCTGAAAGTCTCCTTTTTCTGTATCGTGGGAAtc                                                                                                                                                                                                               >  1:297152/1‑151 (MQ=255)                                                                                                              ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGTGAATGGGTGCGACCATGACAGTCTCCTTTTTCGGTATCGTGGAAAtc                                                                                                                                                                                                               <  2:297152/151‑1 (MQ=255)                                                                                                                   cGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGc                                                                                                  <  2:289077/251‑1 (MQ=255)                                                                                                                      cAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTGTACCACAGATGCGtt                                                                                               <  1:129969/251‑1 (MQ=255)                                                                                                                                            gCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTCCGGGGAGGGGGGGGGCCCTGGCAGTTTCCTTTTTCTGTTTCGTGGAAATCA‑TTTTCATTTTTTTTTTTAGCTAATTGTAT‑AGTTATTGAACTTATTCGATGGG‑TTTATTTTG‑AACAAAATTCATGTTGCGGGGGGGGCGCTTTTTCTACAGAAGCGGTTTTGCCAGTATGGTttgttta                                                                          >  2:14369/1‑248 (MQ=255)                                                                                                                                                   gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCAGTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttat                                                                  <  1:460126/251‑1 (MQ=255)                                                                                                                                                                 cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCAGCGATTTGAGGGGATGGGTGCGACCCTGGCAATTTCATTTTTCTTTATTGTGGAAATAAATTTTC‑TTTTTTTTTTTAGCTAAAGAAATCAGTT‑CTTAACTGATCCGGAGAG‑TTAATTTTG‑AAAAAATCTAATGTTTCGTGGTGGTCGCTTTTACCACAGCTTCGTTTATTCCAGTATGGTTTGGTGTATTTTTATTAAATATggggtgg                                                    >  1:431077/1‑247 (MQ=255)                                                                                                                                                                               gTGCGGTGGCAACGGCTGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCa                                      <  1:128445/251‑1 (MQ=255)                                                                                                                                                                                             gcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAgcagca                         >  2:9742/1‑250 (MQ=255)                                                                                                                                                                                                ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATTTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTGTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                     <  1:469955/251‑1 (MQ=255)                                                                                                                                                                                                                   cgGGAATGGGTGCGACCATTACAGTCTCCGTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTTAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTCGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAATTGGAACAACTGCGTGCGCAg  <  2:453542/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   CAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAAT‑AGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAG  >  NC_000913/4223400‑4223860

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: ‑