Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,083,929	C→T	P392S (CCG→TCG)	sbcB →	exodeoxyribonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,083,929	0	C	T	91.4%	84.8 / ‑3.2	36	P392S (CCG→TCG)	sbcB	exodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (12/20);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

AATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAG  >  NC_000913/2083692‑2084172                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     aaTACGCTACGCCCTGAAGAGTCCGACAGACTGGGAACTAAGCGACAGCATATCCTCGATAGCCTGAAAATTCTGCGTGAAAATCCACAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTTCGCCTTCAGATAACGTGGATGCAAAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGatat                                                                                                                                                                                                                                         <  2:188382/251‑1 (MQ=255)          cgcCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTttgt                                                                                                                                                                                                                                >  2:114716/1‑251 (MQ=255)                                   gAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCtgttg                                                                                                                                                                                                       <  1:116274/251‑1 (MQ=255)                                        aaTCGTCAGCATGGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCaa                                                                                                                                                                                                  <  1:466818/251‑1 (MQ=255)                                                         cGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAAc                                                                                                                                                                                  <  1:46429/250‑1 (MQ=255)                                                                           gCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATACCGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCATTTATCGGGCACGCAAATTCCCGGGGACGCTGGAtt                                                                                                                                                               >  1:318852/1‑251 (MQ=255)                                                                                         aaGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAgca‑gca                                                                                                                                                 >  1:147328/1‑250 (MQ=255)                                                                                                      aGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGc                                                                                                                                   >  1:178035/1‑251 (MQ=255)                                                                                                            ggtgatgatcccgGAAGCCGTACCGTTTACGCCATCAGATAACGTGGATGCACAGCTTTATAACGGCGTTTGCAGTGACGCCGACCGTGCAGCCACTAAAATTGTGCTGGCCACCGAGCCGCGTAAGTTATCGGCGCTGGCTAACACGTTTGTTGACAAACGGATTGCAGAGCTGTTGCTCAATTATCGCGCTCGTAACTTCCCGGGGACCCTGGATTATGCCGTGCAAGCA‑CACTGGCTGGTGCACCGTc                                                                                                                             <  2:397632/240‑1 (MQ=255)                                                                                                                    tCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGt                                                                                                                      <  1:544775/250‑1 (MQ=255)                                                                                                                       gagAAGCCTAACCGTTTTCGCAGTCAGATAACGTGGATGCACAGCTGTATAACGGCTTGTTCCATGACGCAGATCGTACAGCAATGAAAATTGTGCTGGAAACAGAGCAGCGTAATTTACCGGCACTGGATATCACGTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGTATTATGCCTTGCA‑GCAACGGTGGTTAGATCACCGTCGCCAGGTCTTc                                                                                                                  <  2:233037/249‑1 (MQ=255)                                                                                                                                   cGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGttt                                                                                                      >  1:193775/1‑251 (MQ=255)                                                                                                                                    gTTTACGCCTTCAGCTAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGCCGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCt                                                                                                                    >  2:261307/1‑236 (MQ=255)                                                                                                                                    gTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTAGCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCt                                                                                                                    <  1:261307/236‑1 (MQ=255)                                                                                                                                                 gATAACGTGGATGTACAGGTTTATAACGGCGTATTCAGGGACGCAGATCGTGAAGCAATGAAATATGTGATGGAAGACGAGCCGCGTAATTTAGCGGAAATGGATATCACTTTTGTTGCAAAGAGGATTGAAAAGCTGTTGGTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGCTTATGCAGGGCA‑GCAACGCTGGCTGGAGCCCCGGCGTCAGGTCTACAGGCAGGAGTTTTTGAAGGGTTaggg                                                                                        >  2:70187/1‑248 (MQ=255)                                                                                                                                                    aaCGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCtga                                                                                     <  1:114716/251‑1 (MQ=255)                                                                                                                                                     ccGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTTGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCtgaa                                                                                    <  2:91676/250‑2 (MQ=255)                                                                                                                                                      cGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGc                                                                                                  <  1:58237/236‑1 (MQ=255)                                                                                                                                                      cGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGATGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGc                                                                                                  >  2:58237/1‑236 (MQ=255)                                                                                                                                                          gATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAAtt                                                                               <  1:259487/251‑1 (MQ=255)                                                                                                                                                           aTGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTg                                                                              <  1:279163/251‑1 (MQ=255)                                                                                                                                                                  gCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGAAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGACGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCCGAGTTTTTGCAGGGTTATGATGCTGAATTGCCGATGc                                                                       >  1:560547/1‑251 (MQ=255)                                                                                                                                                                  gCTTTATAACGGCTTTTTAAGTGACGCAGATCGTGCCGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTCTGCCGAGCA‑GCAACGCTGGCTGGAGCACCGCCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGACTTGCAGATGc                                                                       >  2:212738/1‑251 (MQ=255)                                                                                                                                                                        taACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGAtt                                                                                                                                                                                                                   <  1:482931/106‑1 (MQ=255)                                                                                                                                                                        taACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCAATGGATATCACTTTTGTTGATAAACGGAtt                                                                                                                                                                                                                   >  2:482931/1‑106 (MQ=255)                                                                                                                                                                           cGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTacaac                                                              <  1:53253/251‑1 (MQ=255)                                                                                                                                                                                ttttCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAGCCGAGCCTAGTAATTTATCGGCTCTGGATATCACGTTTGTTGATAAACCGATTGTAAAGCTGTAGTCCAATTATCGTGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑TCACCGCTGGTTGGAGCACATTCTTCATGTCTGCACGCCAGAGTTTTTTCAGGGTTATTCTGCTGAATTGCAGATGCTGTTACAACAata                                                          <  2:485020/250‑1 (MQ=255)                                                                                                                                                                                  ttCAGTGACGCAGATCGTGCAGCAATGATAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACACTATGc                                                       >  1:74556/1‑251 (MQ=255)                                                                                                                                                                                                 cGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGAGATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTTGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAgaga                                         <  2:244517/250‑1 (MQ=255)                                                                                                                                                                                                  gTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGaa                                        >  1:182569/1‑250 (MQ=255)                                                                                                                                                                                                       gcaATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGgc                                  <  1:429131/251‑1 (MQ=255)                                                                                                                                                                                                          ttGAAAAATGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGTGCA‑GCAAGGCTGGCTGGAGCACCGTCGCCATGTCTTCACGCACGAGTTTTTGCAGGGTTATCTTGATGAATTGCAGTTGCTGGCACCACAATATGCCGGTGACaaggataacgtgtgcct                               >  2:135183/2‑236 (MQ=255)                                                                                                                                                                                                          aTGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCt                               <  2:193775/251‑1 (MQ=255)                                                                                                                                                                                                                         gAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTg                <  1:197821/251‑1 (MQ=255)                                                                                                                                                                                                                              ccaaccGCGTAATTTATCGGCACTGGATATCTCTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAAGTATCGGGCCCGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGTAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAATAATATGCCGATGACAAAGAGAATGTGGCGCTGTTAAAAgtctttataaagc           <  2:560547/247‑13 (MQ=255)                                                                                                                                                                                                                                         aTTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTTAAAAGCTGTTGTTCAATTATCGGGAACGCAACTTCCCGGGGACGCTGGATTATACCGAGCA‑GCAACGCTGGATGGAGCAACGTCTCCAGGTCTTCACGCCAGAGTTTTTGCAAGGTTATGCTGATGAATTGCAGATGCTGGTACAACCATATGCCGCTGACAAAGCGAAATTGTAGCTGTTAAAAGCACTTTggtcgacagcggcag  >  2:83320/1‑236 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     AATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCA‑GCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAG  >  NC_000913/2083692‑2084172

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