Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,420,554	C→T	7.1%	R12C (CGT→TGT)	PP_1245 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,420,554	0	C	T	7.1%	163.6 / 2.3	57	R12C (CGT→TGT)	PP_1245	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base T (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

CCCTGCGCCCCTCGCTGGACGGCTGAATTGTCGTCTGCGGCCAGTGCTGCACACTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCGGCTCGACCATGGGCAGTG  >  NC_002947/1420406‑1420695                                                                                                                                                      |                                                                                                                                              cccTGCGCCCCTCGCTGGACGGCTGAATTGTCGTCTGCGGCCAGTGCTGCACACTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGc                                                                                                                                               <  5:231141/149‑1 (MQ=255)      cgcCCCTCGCTGGACGGCTGAATTGTCGTCTGCGGCCAGTGCTGCACACTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGt                                                                                                                                             >  6:165777/1‑146 (MQ=255)                       cTGAATTGTCGTCTTCGGCCAGTGCTGCACACTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAATTTCGTCGCGTCCCt                                                                                                                         <  8:257144/149‑1 (MQ=255)                              gtcgtcTGCGGCCAGTGCTGAACAATTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCTAAAATGTAGTGCTTCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGc                                                                                                                   <  1:48837/148‑1 (MQ=255)                               tcgtcTGCGGCCAGTGCTGCACACTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTAGCGTAAGAAATTCGTCGCGTCCCTTCTGGCgg                                                                                                                 >  8:115248/1‑149 (MQ=255)                                  tcTGCGGCCAGTGCTGCACACTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCg                                                                                                              >  4:284520/1‑149 (MQ=255)                                     gCGGCCAGTGCTGCACACTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCa                                                                                                           >  6:168956/1‑149 (MQ=255)                                       ggCCAGTGCTGCACACTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATc                                                                                                         <  6:45818/149‑1 (MQ=255)                                                  cacaCTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTAcc                                                                                                <  6:33364/147‑1 (MQ=255)                                                  cacaCTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCAc                                                                                              >  3:204181/1‑149 (MQ=255)                                                   acacTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACg                                                                                             <  2:10338/149‑1 (MQ=255)                                                   acacTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGGAAGAACTtc                                                                                                                                    >  8:215165/1‑110 (MQ=255)                                                     acgttCTGAAAACCCCGCACAAACTCGAATTCTGCACCCAACCATTGGCTGCGGTAATCGGCTGTGC‑GCCACACCGAATAAGGAGTGCTACACTAGAGAAAGAAATTCCTCGCGTCCCTTCTGTCGGTCGCCATAGCCTCTACCACgga                                                                                           >  1:130700/4‑148 (MQ=255)                                                     acTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGc                                                                                           >  2:140650/1‑149 (MQ=255)                                                     acTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGc                                                                                           <  7:115248/149‑1 (MQ=255)                                                       tttCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGc                                                                                                       <  1:122131/135‑1 (MQ=255)                                                       tttCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGc                                                                                                       >  2:122131/1‑135 (MQ=255)                                                       tttCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCt                                                                                         >  1:323930/1‑149 (MQ=255)                                                             aaaaGCCCGCACAAACTCGAATGCTGCCGCCAACCATTTTCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCtttg                                                                                      <  6:2755/146‑1 (MQ=255)                                                             aaaaCCCCGCACAAACTCGAATGCTGCCGCAAACCATTGGCTGCGGTAAGCGGCTGTGC‑GCCACACAGAAAATGGAGTGATCAACTTGCGTAAGAAATTCGTCGCCTCCCTACTGGCGGTCGCCATCGCATCTACCACGGCCtgtg                                                                                      >  5:2755/1‑146 (MQ=255)                                                                  cccGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTc                                                                              >  8:277242/1‑149 (MQ=255)                                                                       acaTACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGAt                                                                         <  4:211480/149‑1 (MQ=255)                                                                         aaaCTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATtc                                                                                                                                    <  6:51143/88‑1 (MQ=255)                                                                         aaaCTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATtc                                                                                                                                    >  5:51143/1‑88 (MQ=255)                                                                          aaCTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATc                                                                        >  7:283171/1‑147 (MQ=255)                                                                          aaCTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATc                                                                        <  8:283171/147‑1 (MQ=255)                                                                           aCTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGc                                                                     <  5:168956/149‑1 (MQ=255)                                                                               gAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCTCAGCTCGGGATCAGCaagg                                                                 >  5:176846/1‑147 (MQ=255)                                                                               gAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCaaa                                                                  <  8:57617/148‑1 (MQ=255)                                                                                       ccgccAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGTGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCg                                                         >  8:71481/1‑149 (MQ=255)                                                                                       ccgccAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGTGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCg                                                         >  8:67156/1‑149 (MQ=255)                                                                                       ccgccAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGc                                                          <  3:284520/148‑1 (MQ=255)                                                                                        cgccAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCgg                                                        >  1:119318/1‑149 (MQ=255)                                                                                           cAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCa                                                     <  4:204181/149‑1 (MQ=255)                                                                                            aaCCATTGGCTGCGGTACGCGGCTGTGCGGCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGc                                                                     >  1:24082/1‑133 (MQ=255)                                                                                            aaCCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGc                                                                     <  2:24082/132‑1 (MQ=255)                                                                                               cATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGc                                                                                           >  4:202646/1‑107 (MQ=255)                                                                                                 ttGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCa                                                                       <  8:153007/125‑1 (MQ=255)                                                                                                    gCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCAcc                                                   >  4:201677/1‑142 (MQ=255)                                                                                                    gCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCAcc                                                   <  3:201677/142‑1 (MQ=255)                                                                                                     cTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCgg                                           <  7:277242/149‑1 (MQ=255)                                                                                                      tGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGt                                          >  1:51380/1‑149 (MQ=255)                                                                                                          gTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGtggcc                                      <  2:51380/149‑1 (MQ=255)                                                                                                                ggCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCg                                >  4:61532/1‑149 (MQ=255)                                                                                                                     tgC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGtgtg                           <  2:323930/149‑1 (MQ=255)                                                                                                                             cACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGtggc                                       >  7:5622/1‑130 (MQ=255)                                                                                                                             cACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGtggc                                       <  8:5622/130‑1 (MQ=255)                                                                                                                             cACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCg                    >  7:247163/1‑149 (MQ=255)                                                                                                                                  aaaaTGGAGTGCTCCACTTGCGTAAGAAATTCGTCGAGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCGGCTc                >  3:70523/1‑148 (MQ=255)                                                                                                                                        gAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCGGCTCGACCATg         <  3:61532/149‑1 (MQ=255)                                                                                                                                              tCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCGCGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCGGCTCGACCATg         >  2:230076/1‑143 (MQ=255)                                                                                                                                              tCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCGCGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCGGCTCGACCATg         <  1:230076/143‑1 (MQ=255)                                                                                                                                              tCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCGGCt                 >  2:19091/1‑135 (MQ=255)                                                                                                                                              tCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCGGCt                 <  1:19091/135‑1 (MQ=255)                                                                                                                                               ccACTTGTGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCGGCTCGACCATGGGCAGTg  <  7:67156/149‑1 (MQ=255)                                                                                                                                                cACTTGTGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCGGCTCGACCATGGGCAGTg  <  7:71481/148‑1 (MQ=255)                                                                                                                                                 aCTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGtgg                                        >  3:241214/1‑109 (MQ=255)                                                                                                                                                 aCTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGtgg                                        <  4:241214/109‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                              CCCTGCGCCCCTCGCTGGACGGCTGAATTGTCGTCTGCGGCCAGTGCTGCACACTTTCTGAAAACCCCGCACAAACTCGAATGCTGCCGCCAACCATTGGCTGCGGTACGCGGCTGTGC‑GCCCCACCGAAAATGGAGTGCTCCACTTGCGTAAGAAATTCGTCGCGTCCCTTCTGGCGGTCGCCATCGCCTCTACCACGGCCTGTGCACAGCTCGGGATCAGCAAAGAGCAAGCGGGCACCGTGATCGGTGGCCTGGCCGGTGTGGCCATCGGCTCGACCATGGGCAGTG  >  NC_002947/1420406‑1420695

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