Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,759,481	C→G	7.0%	I38M (ATC→ATG)	folEB →	GTP cyclohydrolase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,759,481	0	C	G	7.0%	179.6 / 3.9	57	I38M (ATC→ATG)	folEB	GTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base G (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CCCGGCTGGAACTGATTCATACGCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGCTGTATCGCCTGCTGGACGCCTGGGGG  >  NC_002947/3759333‑3759619                                                                                                                                                      |                                                                                                                                           cccGGCTGGAACTGATTCATACGCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATc                                                                                                                                            <  1:342644/149‑1 (MQ=255)  ccGGCTGGAACTGATTCATACGCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATc                                                                                                                                            >  6:70455/1‑148 (MQ=255)  ccGGCTGGAACTGATTCATACGCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCg                                                                                                                                           >  6:9057/1‑149 (MQ=255)   cGGCTGGAACTGATTCATACGCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGa                                                                                                                                          <  1:232273/149‑1 (MQ=255)             tGATTCATACGCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGt                                                                                                                                 >  3:174142/1‑148 (MQ=255)                 tCATACGCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAgg                                                                                                                             >  4:56945/1‑148 (MQ=255)                    tACGCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATcc                                                                                                                         <  7:171370/149‑1 (MQ=255)                      cGCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCaa                                                                                                                       <  5:325712/149‑1 (MQ=255)                       gCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGCCCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATcc                                                                                                                         >  6:370698/1‑146 (MQ=255)                       gCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAg                                                                                                                      <  5:291078/149‑1 (MQ=255)                          ttGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAAc                                                                                                                   >  1:244599/1‑149 (MQ=255)                                    gTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCa                                                                                                         <  4:155213/149‑1 (MQ=255)                                        aaGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGt                                                                                                     <  3:56945/149‑1 (MQ=255)                                                 accacTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCAc                                                                                            >  5:26268/1‑149 (MQ=255)                                                  ccacTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGc                                                                                                       <  4:203263/137‑1 (MQ=255)                                                  ccacTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGc                                                                                                       >  3:203263/1‑137 (MQ=255)                                                        gATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCgg                                                                                      <  5:9057/148‑1 (MQ=255)                                                          tataCATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCg                                                                                   <  5:70455/149‑1 (MQ=255)                                                             aCATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAg                                                                                                                      <  6:252574/111‑1 (MQ=255)                                                             aCATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAg                                                                                                                      >  5:252574/1‑111 (MQ=255)                                                             aCATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGAt                                                                                >  8:155286/1‑149 (MQ=255)                                                              cATGATGATCTCAGCATTACCCGACGTTTCCGTCACCGACCGGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGAt                                                                                <  7:99153/148‑1 (MQ=255)                                                              cATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGAt                                                                                <  2:244599/148‑1 (MQ=255)                                                                      gCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTc                                                                       >  8:123595/1‑149 (MQ=255)                                                                        tCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTCCCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGa                                                                     >  3:42204/1‑149 (MQ=255)                                                                           gCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCt                                                                  >  7:81398/1‑149 (MQ=255)                                                                            cATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTg                                                                 >  5:318565/1‑149 (MQ=255)                                                                             aTTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGTCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGc                                                                <  4:252037/149‑1 (MQ=255)                                                                             aTTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAgg                                                                         <  8:63336/140‑1 (MQ=255)                                                                             aTTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAgg                                                                         >  7:63336/1‑140 (MQ=255)                                                                                   cGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTaa                                                                                                                               >  6:200800/1‑80 (MQ=255)                                                                                   cGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTaa                                                                                                                               <  5:200800/80‑1 (MQ=255)                                                                                   cGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCAc                                                          <  7:252197/149‑1 (MQ=255)                                                                                        tttCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCAt                                                     <  1:156010/149‑1 (MQ=255)                                                                                         ttCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATg                                                    >  4:55516/1‑149 (MQ=255)                                                                                            cGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATGGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCa                                                           <  7:336020/139‑1 (MQ=255)                                                                                            cGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATGGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCa                                                           <  7:331281/139‑1 (MQ=255)                                                                                            cGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATGGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCa                                                           >  8:336020/1‑139 (MQ=255)                                                                                            cGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATGGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCa                                                           >  8:331281/1‑139 (MQ=255)                                                                                            cGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTc                                                  <  7:123595/148‑1 (MQ=255)                                                                                              tCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCCGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCaaa                                               <  5:90485/149‑1 (MQ=255)                                                                                               cACCGCCCCGGCCCCCACCCCCCACCCCGC‑GGAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAg                                              <  4:128942/149‑1 (MQ=255)                                                                                               cACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGc                                                                <  6:345984/131‑1 (MQ=255)                                                                                               cACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGc                                                                >  5:345984/1‑131 (MQ=255)                                                                                                  cGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAg                                                                          >  6:308853/1‑118 (MQ=255)                                                                                                      cTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTca                                       <  6:26268/149‑1 (MQ=255)                                                                                                         gCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACAt                                    >  6:166781/1‑149 (MQ=255)                                                                                                          cccccA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATg                                   >  6:368217/1‑149 (MQ=255)                                                                                                           ccccA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGt                                  >  7:319305/1‑149 (MQ=255)                                                                                                                 gCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAgg                             <  7:155286/149‑1 (MQ=255)                                                                                                                        cccGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGCTGTATc                      <  5:370698/149‑1 (MQ=255)                                                                                                                        cccGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGCTGTATc                      <  4:42204/149‑1 (MQ=255)                                                                                                                               cAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGCTGTATCGCctgct               <  3:211497/149‑1 (MQ=255)                                                                                                                               cAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGCTGTATCGCctgc                <  3:55516/148‑1 (MQ=255)                                                                                                                                  tGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGCTGTATCGCCTGCTGGa            <  6:318565/149‑1 (MQ=255)                                                                                                                                     gTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGCTGTATCGCCTGCTGGACGc         >  2:151163/1‑149 (MQ=255)                                                                                                                                        ggCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGCTGTATCGCCTGCTGGACGCCTg      <  5:368217/149‑1 (MQ=255)                                                                                                                                         gCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGCTGTATCGCCTGCTGGACGCCTgg     >  1:175524/1‑149 (MQ=255)                                                                                                                                            tGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGTTGTATCGCCTGCTGGACGCCTggggg  >  5:347231/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                           CCCGGCTGGAACTGATTCATACGCCTTGCTGCCACGTGCAAGGCTCCCACCACTGGATATACATGATGAGCTCAGCATTACCCGACGTTTCCGTCACCGACCTGGCCCCCA‑GCCTCAGCCCGCTGCAATGGGTCGGCATGCAGGGTATCGACCTGCCGGTAAGGATCCAAGAACCGACCTACCAGCGTGAGCTGCACGCCCGGGCGGATGTCCAGGTCGACCTGCCGGCACCGCATGTCAAAGGCATTCACATGTCCAGGCTGTATCGCCTGCTGGACGCCTGGGGG  >  NC_002947/3759333‑3759619

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