Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	76,727	G→A	H191H (CAC→CAT)	sgrR ←	DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	76,727	0	G	A	100.0%	116.9 / NA	35	H191H (CAC→CAT)	sgrR	DNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/17);  total (18/17)

CCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAGTAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACATCGTCCATTTCCAGTTCACGACCATGGTGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGAAATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCGATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCG  >  NC_000913/76483‑76963                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             ccGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGaaaa                                                                                                                                                                                                                                        <  2:80356/251‑1 (MQ=255)  cGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAAt                                                                                                                                                                                                                                       >  1:174224/1‑251 (MQ=255)       aCGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACt                                                                                                                                                                                                                                  <  2:127266/251‑1 (MQ=255)       aCGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACt                                                                                                                                                                                                                                  >  1:87305/1‑251 (MQ=255)                            ggctggcAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACTACCATGATGAAAATGGACTCCTGGACTCAAAAAGAAACg                                                                                                                                                                                                              <  2:192818/250‑1 (MQ=255)                               gggcAAAGTTACTGAGGGTTTTCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCGGGTTGCGTGAGATGGATATCCAGCGTCGAGGGCGTCGGCGACACAATGTCAGCAATATTCTAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACGGTGATGAAAATGGAGTTCTGGACGCAAAAAGAAACGCCAg                                                                                                                                                                                                          <  2:113146/250‑1 (MQ=255)                                   aaaattACGGAGGGTTTCCCATTCGCGCGGCACGAGCATCGCCGGAACTTGCCCCTGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGAAGGATATCCAGCGTCCAGTGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACTTCTTCCATATCCAGTTCACGATCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGaa                                                                                                                                                                                                      <  2:231027/247‑1 (MQ=255)                                              gggTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATAtt                                                                                                                                                                                           >  1:48144/1‑251 (MQ=255)                                                   tCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTtgct                                                                                                                                                                                       <  1:11031/250‑1 (MQ=255)                                                   tCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTtgct                                                                                                                                                                                       <  1:12568/250‑1 (MQ=255)                                                                                      ccccAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGATGGGCGATGTCTGCTTCCAGTTCCCCAt                                                                                                                                                   >  2:107178/1‑251 (MQ=255)                                                                                           gcaaaaaCGGGAACCAGCGGGGCGGTTGCGTGAGAGGGATAGCCAGCGTCCGTGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGGAGCGTATTGATTCGGTGTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATTATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAGGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTcc                                                                                                                                              <  2:26204/246‑1 (MQ=255)                                                                                            caatAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCt                                                                                                                                             >  1:259743/4‑251 (MQ=255)                                                                                              atAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGCTATCCAGCGTCCAGGGCGTCGGCGACACACTGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTc                                                                                                                                            >  2:17900/2‑250 (MQ=255)                                                                                                   ggTAACCAGCGGTCCGGGTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTGCACGACCATGATGAAAATGGACTCCTGTACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGGCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTAt                                                                                                                                      <  1:14438/251‑1 (MQ=255)                                                                                                                ccGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGa                                                                                                                          <  1:167704/250‑1 (MQ=255)                                                                                                                 cGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAAc                                                                                                                        >  2:279053/1‑251 (MQ=255)                                                                                                                     tGCGTGAGATGGATATCCAGCGTCCAGGGCGTAGGCGACACAATTTCCGCACTATGCGAATAGAGCGGCCGCGTATTGATTCGTTTTACAGAGGCGATCACGTCGTCCATTTCCAGTTCCCGACAATGATGAAAATGGACTCCTGAACGCAAAAAGACACGCCAGTGAACCGGTGATCTTTGCTTCCAGTGGTGGGCGATGTCTGCTTCCCGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGACCTGa                                                                                                                     >  2:321081/1‑250 (MQ=255)                                                                                                                                            ccAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTTCCGGGCAATATGGGtt                                                                                              <  1:255823/250‑1 (MQ=255)                                                                                                                                              aGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCg                                                                                           <  2:249395/251‑1 (MQ=255)                                                                                                                                                   gTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACg                                                                                      >  2:95796/1‑251 (MQ=255)                                                                                                                                                   gTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACg                                                                                      >  2:187994/1‑251 (MQ=255)                                                                                                                                                      ggCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCAAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGgcg                                                                                   >  2:266650/1‑251 (MQ=255)                                                                                                                                                          acacAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGCGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAAt                                                                               >  1:241294/1‑251 (MQ=255)                                                                                                                                                               aTGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCt                                                                          >  2:67943/1‑251 (MQ=255)                                                                                                                                                                                aTAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGAGCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGTCAAGATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTAc                                                         <  1:200921/251‑1 (MQ=255)                                                                                                                                                                                 tAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACg                                                        >  1:113461/1‑251 (MQ=255)                                                                                                                                                                                         cAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGAc                                                >  2:107351/1‑251 (MQ=255)                                                                                                                                                                                               aTTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGAtagt                                          <  1:74038/251‑1 (MQ=255)                                                                                                                                                                                                     tCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGc                                     <  1:1479/250‑1 (MQ=255)                                                                                                                                                                                                          ttAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCAcgcg                                <  2:63929/250‑1 (MQ=255)                                                                                                                                                                                                                    gATCACGTCGTCCATTTCCCGTTCACGCCCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGAAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTgcc                     >  2:264042/1‑251 (MQ=255)                                                                                                                                                                                                                          gtCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCATGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCt               <  1:200542/250‑1 (MQ=255)                                                                                                                                                                                                                             gTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGGGCCGCCCCTGGc            >  2:113674/1‑251 (MQ=255)                                                                                                                                                                                                                                       aGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCg  >  1:305309/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             CCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAGTAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACATCGTCCATTTCCAGTTCACGACCATGGTGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGAAATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCGATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCG  >  NC_000913/76483‑76963

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