Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,318,175	A→G	Y568H (TAC→CAC)	rcsC ←	histidine kinase RcsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,318,175	0	A	G	100.0%	120.0 / NA	36	Y568H (TAC→CAC)	rcsC	histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/19);  total (17/19)

TACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGCCGATCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATA  >  NC_000913/2317929‑2318421                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       atcgctgttcgccgttcccgcccgggAAAAAGGGATCAAACAAGCCCGCCACTTACTTGGCCGGGAGCCCCACGCCGGGATCGCTAACCCTAATGGGGAGATAACCGCCGTCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCGTCGGGTCGCCATTTAAGGCCACACGCGGATCCTGTTCAATAAAGCAGTGCAAg                                                                                                                                                                                                                                                    <  2:350574/227‑1 (MQ=255)                  ccGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTac                                                                                                                                                                                                                                   <  1:103853/251‑1 (MQ=255)                   cGACCTGGAAGAAGGGATCAAACAAGCGCCCCCCTGCTTCCGCCGGTACCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTacc                                                                                                                                                                                                                                  <  2:277871/251‑1 (MQ=255)                         gaaagacgggtgcAAACAAGCGCACCACTTCTTTCGCCGGTATCGCTACGCCGGTATCGCGAACACGGATAGTGAGATAATCGCCATCCGCTGGAACATGCAAATCTATACAGCCGGTATCGGTTAATTTTATGGCGTTACTCAACAGGTTGGAGTTGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTAccacca                                                                                                                                                                                                                              <  1:190622/237‑1 (MQ=255)                              aagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTa                                                                                                                                                                                                                        >  2:105902/1‑250 (MQ=255)                                gggATCAAACAAGCGCACCACTTCTTTCGCCGGTACCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAAt                                                                                                                                                                                                                     <  2:88823/251‑1 (MQ=255)                                     caaccAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGtt                                                                                                                                                                                                                 <  1:276351/246‑1 (MQ=255)                                      aaaCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTgg                                                                                                                                                                                                               >  2:273892/1‑251 (MQ=255)                                                     ttcttTCGCCGGTATCCCCACGCAGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAACCTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTtc                                                                                                                                                                                                 >  2:184045/1‑250 (MQ=255)                                                                                     gaacacggattaagaacataTCGCCATCCGCGCTAACATGGAAATCTATACAGCCGGTATCGGGTTATTTTAGGGGGTTTCTCAGCAGGTTGGAGATGTCCTGCTGTAAACGCATCGGTTCGCCATTTAAGGCCACTGGCTCATTCGGTTCAATAAAGCAGTGCGAGCCTAACTGCTTGCGTACCACCAGCGTTAAATATTTGTCGTTTATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCg                                                                                                                                                                <  1:348591/232‑1 (MQ=255)                                                                                     gAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAAGTATACAGCCGGTATCGGTGAATTTTATGGCGTAACTCAACAGGGTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAGGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCg                                                                                                                                                                >  1:160888/1‑251 (MQ=255)                                                                                          cGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCtt                                                                                                                                                           >  2:104381/1‑251 (MQ=255)                                                                                               agagagATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACt                                                                                                                                                      <  1:10958/251‑1 (MQ=255)                                                                                                    gaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTcc                                                                                                                                                 >  1:358688/1‑251 (MQ=255)                                                                                                        aTCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGAtt                                                                                                                                             <  1:94904/251‑1 (MQ=255)                                                                                                        aTCGCCACCCGAGCGAACATGCAAAACTATACAGCCGGTATCGGTGACTTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTAAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCCACTGTTCCGAtt                                                                                                                                             >  2:229256/1‑251 (MQ=255)                                                                                                                    gcgAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCgaga                                                                                                                                 <  1:317104/251‑1 (MQ=255)                                                                                                                    gcgAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCg                                                                                                                                    >  1:339063/1‑248 (MQ=255)                                                                                                                          aTGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCg                                                                                                                            >  1:124601/1‑250 (MQ=255)                                                                                                                              aaaaCTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAat                                                                                                                       >  1:341137/1‑251 (MQ=255)                                                                                                                                                           tATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTgg                                                                                          >  1:67103/1‑251 (MQ=255)                                                                                                                                                                   tactcaaccggtcggcaatgacgaactgTAAACGCTCCGAGTCGCCATATAAGGCCACTGGCACATCCGGTTCAATGAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTACGATTCCATCTTCGAGAAATCGAGAATATCGCTGATAATTTCAAAAAACAGGCTGGAAGAgttg                                                                                  >  2:179813/26‑251 (MQ=255)                                                                                                                                                                    aCTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAgttg                                                                                  <  1:9113/250‑1 (MQ=255)                                                                                                                                                                       cAACAGGTTGGAGATGACCTGCTGTAAACGCATAGGGTCGCCATTTAAGGCCACTGGCACATCCGGATCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTACGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTAAACAACAGGCTGGAAGAGTTGTTc                                                                               >  1:345041/1‑250 (MQ=255)                                                                                                                                                                             gTTGGAGATGACCTGCTGTAAACGCATCGCGTCGCCATTTAAGGCCACTGGCACCTCCGGTTCAATAAACCAGTGCAGGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATGCAATCTTCGAGAAATCGAGAATATCGCGGATAATTTGCAACAACGGGCTGGAAGATGTGGTCATTGCCg                                                                        <  2:217000/251‑1 (MQ=255)                                                                                                                                                                                      gACCTGCTGTAAACTCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGcc                                                               <  2:322174/251‑1 (MQ=255)                                                                                                                                                                                           gctgTAAACGCATCGGTTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTccccctcccatc                                                           <  1:290333/250‑12 (MQ=255)                                                                                                                                                                                           gctgTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGCCGATc                                                           <  2:142348/250‑1 (MQ=255)                                                                                                                                                                                            ctgTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGCCGATCg                                                          >  2:45974/1‑250 (MQ=255)                                                                                                                                                                                                         gggTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTc                                                                                                                                     >  2:105599/1‑162 (MQ=255)                                                                                                                                                                                                         gggTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTc                                                                                                                                     <  1:105599/162‑1 (MQ=255)                                                                                                                                                                                                                     aaGGCCACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGCCGATCGACGCCTTTCGGTAACTCTTTGGTTTg                                <  2:135312/251‑1 (MQ=255)                                                                                                                                                                                                                          cACTGGCACATCCGGTTCAATAAAGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGCCGATCGACGCCTTTCGGTAACTCTTTGGTTTGCAAc                            >  2:272116/1‑250 (MQ=255)                                                                                                                                                                                                                                                aaGCAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGCCGATCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGcc     >  1:284550/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gCAGTGCATGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGGCACCAGCCGATCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCAt   <  2:133112/251‑1 (MQ=255)                                                                                                                                                                                                                                                   cAGTGCAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGCCGATCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATa  <  1:130856/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       TACGCTGTACGCCCGTTCCGACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCACGCGGTGAAAACTCACGCGGTTCGATCTTCAACTGTTCCGATTCAATCTTCGAGAAATCGAGAATATCGCTGATAATTTTCAACAACAGGCTGGAAGAGTTGTTCATTGCCGTCACCAGCCGATCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATA  >  NC_000913/2317929‑2318421

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