Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	157,491	A→C	7.3%	F79V (TTC→GTC)	citN ←	citrate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	157,491	0	A	C	7.3%	158.5 / 4.7	55	F79V (TTC→GTC)	citN	citrate transporter
Reads supporting (aligned to +/- strand):  ref base A (19/32);  new base C (2/2);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

TCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGATCGGTACCAG  >  NC_002947/157344‑157638                                                                                                                                                     |                                                                                                                                                    tCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGaa                                                                                                                                                     >  6:157751/1‑148 (MQ=255) tCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGaa                                                                                                                                                     >  4:357672/1‑148 (MQ=255)   ggCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACa                                                                                                                                                   <  5:373057/148‑1 (MQ=255)     cagcaTGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGc                                                                                                                                                <  2:193415/149‑1 (MQ=255)      agcaTGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCTGTTCGACCCGGcc                                                                                                                                               <  2:375974/149‑1 (MQ=255)       gcaTGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCg                                                                                                                                              <  6:30110/149‑1 (MQ=255)          tGGCGGCGCAGCAGATCATGTAGGTGGTCGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCTGTTCGACCCGGCCgttc                                                                                                                                           <  7:109193/149‑4 (MQ=255)           ggcggaGAAGCAGCTCATGTAGGTGGTGGCACCGTCACCTTCGAGCGAGACCACAAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGAGGATAGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAAt                                                                                                                                          >  4:329187/1‑149 (MQ=255)            gcggcgCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATc                                                                                                                                         <  5:157751/149‑1 (MQ=255)            gcggcgCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATc                                                                                                                                         <  3:357672/149‑1 (MQ=255)              ggcgCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCAGGTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATc                                                                                                                                         >  1:112090/1‑148 (MQ=255)                      agATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAgg                                                                                                                               >  5:231638/1‑149 (MQ=255)                              tAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGaaa                                                                                                                         >  8:260406/1‑147 (MQ=255)                                ggtggtGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                                     <  8:113276/149‑1 (MQ=255)                                ggtggtGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                                     <  6:76933/149‑1 (MQ=255)                                ggtggtGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                                     <  2:112090/149‑1 (MQ=255)                                  tggtggCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAgg                                                                                                                   <  1:78462/149‑1 (MQ=255)                                   ggtggCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGa                                                                                                                  >  5:108054/1‑149 (MQ=255)                                    gtggCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGag                                                                                                                 <  5:271827/149‑2 (MQ=255)                                      ggCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCCGAATCGATCATAAGGGCGAAATAGAGGATgc                                                                                                               >  7:333040/1‑148 (MQ=255)                                       gCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGc                                                                                                              <  2:195104/149‑1 (MQ=255)                                            gtcaccgtcGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACa                                                                                                         >  2:41930/1‑149 (MQ=255)                                             tcaccgtcGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACAt                                                                                                        <  2:382088/149‑1 (MQ=255)                                                   tcGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                                     >  2:14743/1‑130 (MQ=255)                                                   tcGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                                     <  1:14743/130‑1 (MQ=255)                                                     gagcgagACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCa                                                                                                >  8:322292/1‑149 (MQ=255)                                                          agaCCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAAc                                                                                                          >  6:237740/1‑134 (MQ=255)                                                          agaCCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAAc                                                                                                          <  5:237740/134‑1 (MQ=255)                                                           gaCCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCgg                                                                                          <  7:322292/149‑1 (MQ=255)                                                            accaccAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGAGAAATAGAGGATGGCGAACATCAGCATCACGCCgg                                                                                          >  5:265593/1‑148 (MQ=255)                                                                 caGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTg                                                                                    >  8:67759/1‑149 (MQ=255)                                                                    ggCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGa                                                                                 <  6:265593/149‑1 (MQ=255)                                                                           tACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGc                                                                           <  6:108054/149‑1 (MQ=255)                                                                            actggtgTGCCGACCGAGACTTTAAGCGGATCTCCCTTGACCAGCTTTAGGATTTTGCGTACGGCCGGGTCGAAAAGGCCGGAATCGATAATCAGTTCGAAATTGAGGATGGCGAACATCAGCATCACGCCGGTAGTTGCGTGCTTGc                                                                           <  8:333040/142‑1 (MQ=255)                                                                                  gTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGcc                                                                    >  4:70628/1‑149 (MQ=255)                                                                                   tGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGccc                                                                   <  1:41930/149‑1 (MQ=255)                                                                                     ccgaccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCctc                                                                 <  7:260406/149‑1 (MQ=255)                                                                                       gaccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTa                                                               <  3:372839/149‑1 (MQ=255)                                                                                       gaccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGACGGCTCGAACAGGCGGGAATCGATCATCAGGGCGAAAAAGAGGATGCCGAACATCAGCATCACGCCGGTTGGTGCGAGCTTCCTGATGCCctc                                                                 >  8:280521/1‑147 (MQ=255)                                                                                        accgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTac                                                              >  6:223415/1‑148 (MQ=255)                                                                                        accgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAg                                                              <  3:237590/149‑1 (MQ=255)                                                                                         ccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGAt                                                                                                                 >  1:322471/1‑97 (MQ=255)                                                                                         ccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGAt                                                                                                                 <  2:322471/97‑1 (MQ=255)                                                                                           gagaCTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGCTCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGcat                                                           <  2:256103/149‑1 (MQ=255)                                                                                            agaCTTTCAGCGGTTCGCACTTGACCAGCTTGAGGATCTTGCGTACGGCCGGGTCGAACAGGCCGGAATCGATCATCATGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGcatc                                                          <  6:195105/149‑1 (MQ=255)                                                                                                            gCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTg                                           <  7:67759/148‑1 (MQ=255)                                                                                                             cccTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGc                                         <  4:345714/149‑1 (MQ=255)                                                                                                                          tGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAaccacc                             >  8:306597/1‑148 (MQ=255)                                                                                                                           gAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCg                            >  2:158580/1‑148 (MQ=255)                                                                                                                                    gCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTTCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAg                       <  3:329187/144‑1 (MQ=255)                                                                                                                                    gCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAgg                      <  7:306597/145‑1 (MQ=255)                                                                                                                                      ttccGGCCTGTTCGACCCGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGaa                 >  3:242114/4‑148 (MQ=255)                                                                                                                                      ttccGGCCTGTTCGACCCGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGaa                 >  6:14011/4‑148 (MQ=255)                                                                                                                                       aaCGGCCGGGTCTAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACa               <  6:231638/149‑1 (MQ=255)                                                                                                                                        aCGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAg              <  7:16762/149‑1 (MQ=255)                                                                                                                                           gCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGAGCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGag           >  6:232061/1‑148 (MQ=255)                                                                                                                                                   gAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGATCGGTACCa   <  3:70628/149‑1 (MQ=255)                                                                                                                                                    aaCAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGATCAGtccagg  >  3:112251/1‑144 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGATCGGTACCAG  >  NC_002947/157344‑157638

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