Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,070,233	T→C	12.0%	intergenic (+33/+109)	aroP‑I → / ← PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,070,233	0	T	C	12.0%	123.9 / 8.0	50	intergenic (+33/+109)	aroP‑I/PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base T (19/25);  new base C (2/4);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

TGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTA  >  NC_002947/1070085‑1070377                                                                                                                                                      |                                                                                                                                                     tGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGt                                                                                                                                                      >  8:6300/1‑149 (MQ=255)    tAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATCCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGtttt                                                                                                                                                   <  6:211987/149‑1 (MQ=255)     aGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGtt                                                                                                                                                     <  4:48231/146‑1 (MQ=255)     aGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGtt                                                                                                                                                     >  3:48231/1‑146 (MQ=255)         aTGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGgcg                                                                                                                                              >  7:359441/1‑149 (MQ=255)          tGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCTGCATGGATCGCTGTTCTGTGGGTGGCTT‑TCAAGCTGACGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGGATGCCTTTAGGCTTTcatcgc                                                                                                                                             <  7:278397/149‑6 (MQ=255)          tGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGgcgc                                                                                                                                             <  6:150355/149‑1 (MQ=255)                tGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt                                                                                                                                       >  6:118282/1‑147 (MQ=255)                  aCGGCGCCGATTCGCATCTCTGTTATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Ga                                                                                                                                   <  7:299402/147‑1 (MQ=255)                  aCGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt                                                                                                                                       <  6:361793/145‑1 (MQ=255)                  aCGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt                                                                                                                                       >  5:361793/1‑145 (MQ=255)                  aCGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Ga                                                                                                                                   <  4:154174/147‑1 (MQ=255)                    gccgccGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGCTTTCATCGCCt‑‑gt                                                                                                                                       >  1:240014/1‑143 (MQ=255)                      cgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Gagatcga                                                                                                                             <  5:235533/149‑7 (MQ=255)                      cgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Gagatcga                                                                                                                             <  5:156864/149‑7 (MQ=255)                       gccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Gagatcgag                                                                                                                            <  2:201679/149‑8 (MQ=255)                            ttCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGACAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGcgccg                                                                                                                         >  3:342999/1‑149 (MQ=255)                               gCATCTCGGT‑ATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCcgcgcg                                                                                                                     >  2:120837/1‑149 (MQ=255)                                cATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCcgcgcg                                                                                                                     <  7:6300/149‑1 (MQ=255)                                   ctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTGTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcggg                                                                                                                   >  1:85633/1‑149 (MQ=255)                                   ctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTA‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcgggc                                                                                                                  >  5:269779/1‑149 (MQ=255)                                      ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcggg                                                                                                                   <  7:52790/145‑1 (MQ=255)                                         aaTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGACGCGCGGGCGGCGCt                                                                                                            >  4:341747/1‑149 (MQ=255)                                          aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                            <  3:289936/148‑1 (MQ=37)                                           tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                          <  6:35883/149‑1 (MQ=38)                                             cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGTGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                        <  6:225152/149‑1 (MQ=38)                                             cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                         <  2:240014/148‑1 (MQ=25)                                             cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                        >  5:225152/1‑149 (MQ=38)                                                aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                      <  5:118282/148‑1 (MQ=255)                                                     ggCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                                <  3:341747/149‑1 (MQ=255)                                                         tGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGGGA‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                            <  8:359441/149‑1 (MQ=255)                                                          ggATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                            <  4:342999/148‑1 (MQ=255)                                                               ggtgtgctgtgGGTGGCTT‑TCATGCTTAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                      <  2:382338/147‑1 (MQ=38)                                                                   tgctgtgGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTGGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                                    <  1:120837/147‑1 (MQ=255)                                                                       gtgggtggCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcggg                                                                                                                   >  8:145075/1‑112 (MQ=255)                                                                          ggtggCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                                           <  1:241300/149‑1 (MQ=255)                                                                              gCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTACGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGACGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                       <  4:246354/149‑1 (MQ=38)                                                                                   tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGCTCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggc                                                                   >  4:212940/1‑149 (MQ=255)                                                                                   tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggc                                                                   >  4:140349/1‑149 (MQ=255)                                                                                         tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                                             >  5:349853/1‑149 (MQ=255)                                                                                         tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                            >  1:262359/1‑102 (MQ=18)                                                                                         tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                            <  2:262359/102‑1 (MQ=18)                                                                                            agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                          <  3:140349/149‑1 (MQ=255)                                                                                                 gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                     >  5:171662/1‑149 (MQ=255)                                                                                                   tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                           <  1:173075/93‑1 (MQ=25)                                                                                                   tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                           >  2:173075/1‑93 (MQ=25)                                                                                                       cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                               >  6:85975/1‑149 (MQ=255)                                                                                                           ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                      <  4:343187/106‑1 (MQ=255)                                                                                                           ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                      >  3:343187/1‑106 (MQ=255)                                                                                                            cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                          <  2:85633/149‑1 (MQ=255)                                                                                                                                    aaGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg                  >  1:56303/1‑149 (MQ=255)                                                                                                                                     aGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg                  <  6:269779/148‑1 (MQ=255)                                                                                                                                                gTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg      >  6:325707/1‑149 (MQ=255)                                                                                                                                                    ggTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa  <  2:56303/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTT‑TCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTA  >  NC_002947/1070085‑1070377

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: ‑