Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,579,424	T→C	intergenic (‑87/+75)	mcrB ← / ← symE	5‑methylcytosine‑specific restriction enzyme subunit McrB/toxic protein SymE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,579,424	0	T	C	100.0%	224.5 / NA	63	intergenic (‑87/+75)	mcrB/symE	5‑methylcytosine‑specific restriction enzyme subunit McrB/toxic protein SymE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/29);  total (34/29)

GCTTGCTTAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAC  >  NC_000913/4579294‑4579546                                                                                                                                   |                                                                                                                           gcttgcttAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTa                                                                                                                   <  2:666493/140‑1 (MQ=255)  cttgcttAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTAc                                                                                                                  >  1:974181/1‑140 (MQ=255)  cttgcttAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTAc                                                                                                                  >  1:60619/1‑140 (MQ=255)  cttgcttAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTAc                                                                                                                  >  1:496887/1‑140 (MQ=255)        taataaATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAggg                                                                                                            <  1:1033823/140‑1 (MQ=255)          ataaATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTg                                                                                                                >  1:1082334/1‑134 (MQ=255)          ataaATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTg                                                                                                                <  2:1082334/134‑1 (MQ=255)          ataaATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAggggg                                                                                                          <  2:433262/140‑1 (MQ=255)          ataaATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAggggg                                                                                                          <  1:933447/140‑1 (MQ=255)          ataaATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAggggg                                                                                                          <  1:460823/140‑1 (MQ=255)          ataaATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAggggg                                                                                                          >  1:597133/1‑140 (MQ=255)            aaaTTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGaa                                                                                                        <  2:430433/140‑1 (MQ=255)               tttttCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGt                                                                                                     <  2:571178/140‑1 (MQ=255)               tttttCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGt                                                                                                     >  1:968938/1‑140 (MQ=255)                     aaTCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTAt                                                                                               <  2:974181/140‑1 (MQ=255)                      aTCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTAtt                                                                                              <  2:597133/140‑1 (MQ=255)                         cAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAg                                                                                           >  2:906417/1‑140 (MQ=255)                                        tCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTcc                                                                            <  2:796487/140‑1 (MQ=255)                                        tCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTAATGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTcc                                                                            <  2:1028948/140‑1 (MQ=255)                                         ccATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTccc                                                                           <  1:17498/140‑1 (MQ=255)                                                  atGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTaa                                                                  >  2:344662/1‑140 (MQ=255)                                                  atGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTaa                                                                  >  1:872700/1‑140 (MQ=255)                                                    gCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAgg                                                                                <  1:75430/124‑1 (MQ=255)                                                    gCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAgg                                                                                >  2:75430/1‑124 (MQ=255)                                                    gCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAcc                                                                >  1:421313/1‑140 (MQ=255)                                                         ccTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAg                                                           >  2:157931/1‑140 (MQ=255)                                                         ccTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAg                                                           >  2:1128039/1‑140 (MQ=255)                                                         ccTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAg                                                           >  1:810848/1‑140 (MQ=255)                                                          cTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGt                                                          <  2:810848/140‑1 (MQ=255)                                                               tGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACg                                                     >  1:629546/1‑140 (MQ=255)                                                                    aGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGtt                                                >  1:194970/1‑140 (MQ=255)                                                                       ctctATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTAcg                                             <  2:896411/140‑1 (MQ=255)                                                                         ctATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTAcgcg                                           >  1:822501/1‑140 (MQ=255)                                                                         ctATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTAcgcg                                           >  2:445324/1‑140 (MQ=255)                                                                           attatCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGAc                                         <  1:157931/140‑1 (MQ=255)                                                                                cGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCt                                    >  1:274982/1‑140 (MQ=255)                                                                                  cTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGt                                  <  2:60619/140‑1 (MQ=255)                                                                                    ttCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGttt                                >  2:839554/1‑140 (MQ=255)                                                                                    ttCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGttt                                >  1:135076/1‑140 (MQ=255)                                                                                    ttCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGttt                                >  2:1005982/1‑140 (MQ=255)                                                                                      cGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTAc                              <  1:344662/140‑1 (MQ=255)                                                                                      cGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTAc                              >  1:69505/1‑140 (MQ=255)                                                                                        cAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCg                            >  2:962630/1‑140 (MQ=255)                                                                                             tACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                       <  2:421313/140‑1 (MQ=255)                                                                                             tACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                       <  1:962630/140‑1 (MQ=255)                                                                                             tACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                       <  2:69505/140‑1 (MQ=255)                                                                                             tACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                       >  1:766374/1‑140 (MQ=255)                                                                                             tACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                       <  2:629546/140‑1 (MQ=255)                                                                                                  aaaCAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTc                  <  2:822501/140‑1 (MQ=255)                                                                                                  aaaCAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTc                  <  1:1128039/140‑1 (MQ=255)                                                                                                   aaCAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTcc                 >  2:79201/1‑140 (MQ=255)                                                                                                     cAATAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCaa               >  1:577704/1‑140 (MQ=255)                                                                                                       aTAGATTTTTACTGCAAAATCAGACTGGCAGATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACACTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATa             >  2:581817/1‑140 (MQ=255)                                                                                                       aTAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGGGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                       >  2:548986/1‑130 (MQ=255)                                                                                                       aTAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGGGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                       <  1:548986/130‑1 (MQ=255)                                                                                                       aTAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                       <  2:226120/130‑1 (MQ=255)                                                                                                       aTAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                       >  1:226120/1‑130 (MQ=255)                                                                                                       aTAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATa             <  2:661854/140‑1 (MQ=255)                                                                                                       aTAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATa             >  2:581747/1‑140 (MQ=255)                                                                                                       aTAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATa             >  2:446556/1‑140 (MQ=255)                                                                                                       aTAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATa             >  1:386107/1‑140 (MQ=255)                                                                                                        tAGATTTTTACTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTGTACCGGCAATCACTCCAATaa            >  2:1066143/1‑140 (MQ=255)                                                                                                                  cTGCAAAATCAGACTGGCAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAc  <  2:872700/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                           GCTTGCTTAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAC  >  NC_000913/4579294‑4579546

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