Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,377,184	T→C	T137A (ACA→GCA)	menH ←	2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,377,184	0	T	C	100.0%	204.5 / NA	58	T137A (ACA→GCA)	menH	2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/22);  total (36/22)

CTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAGCCCCTGGCAAGCCGCCA  >  NC_000913/2377050‑2377312                                                                                                                                        |                                                                                                                                 cTCCAGCATGGCGGCAAGCGTTGCGCCTTTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATCCCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGaa                                                                                                                             <  2:1321618/140‑1 (MQ=255)       cagggcggcAAGCGTTGCGCCTTTTTTTTTGCTGCGCAGCCCCACCAGCCCCCGGCGTTGATCGTCATTGAGTGAGGAAAAAAC‑AGGCTGTTGATCCCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCt                                                                                                                       <  2:1275565/137‑1 (MQ=255)       cATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCt                                                                                                                       <  1:625270/140‑1 (MQ=255)           gagacaggCGTTGCGCCATTATTGTTGCTGCGCAGCGCCCCCAGCCCCCGGCGTTGACCGTCATTGAGTGAGGAAAAAAC‑GGGCTGTTGATCCCAGTCGGAAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCAc                                                                                                                   <  2:32502/133‑1 (MQ=255)            cggcAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCAcc                                                                                                                  >  1:452788/1‑140 (MQ=255)                 aGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTg                                                                                                             >  1:280164/1‑140 (MQ=255)                       gccccATTATTTTTGCTGCGCAGCGCCACCAGCCCCCGGCGTTGATCGTCATTGAGTGAGGAAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATc                                                                                                       <  2:452788/137‑1 (MQ=255)                             ttattGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATctgtctc                                                                                                 >  1:32502/1‑134 (MQ=255)                             ttattGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACg                                                                                                 >  2:449724/1‑140 (MQ=255)                             ttattGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACg                                                                                                 >  2:283576/1‑140 (MQ=255)                                     gctgcGCAGCGCCACCAGCTCCGGGCGTTGATCGTCATTGAGGGAGGAAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGtt                                                                                         <  2:670937/140‑1 (MQ=255)                                     gctgcGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGtt                                                                                         >  1:29264/1‑140 (MQ=255)                                         cgcAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGc                                                                                     >  2:1370349/1‑140 (MQ=255)                                           cAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCAc                                                                                   >  2:1174914/1‑140 (MQ=255)                                           cAGCCCCCCCAGCTCCCGGCTTTGATCGTCTTGGAGTGA‑GCAAAAACAAGGTTGTTGATCCAAGTCGGAAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCAc                                                                                   <  2:676781/140‑1 (MQ=255)                                              cgcCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGGAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACgtt                                                                                >  2:1460534/1‑140 (MQ=255)                                                ccaccaGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACgttgt                                                                              >  1:1564693/1‑140 (MQ=255)                                                      gCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCa                                                                        >  2:468240/1‑140 (MQ=255)                                                       cTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCAt                                                                       >  1:1341159/1‑140 (MQ=255)                                                        tCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCTGCACGTTGTTCAGCAtt                                                                      >  1:177161/1‑140 (MQ=255)                                                                         aTTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGcc                                                     >  2:1469422/1‑140 (MQ=255)                                                                           tgagtgagGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCg                                                   <  2:177161/140‑1 (MQ=255)                                                                           tgagtgagGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCg                                                   <  1:1469422/140‑1 (MQ=255)                                                                           ggagtgagGAAAAAAC‑AGGTTGTTGACCCCAGTCGAAAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCg                                                   <  2:1211034/139‑1 (MQ=255)                                                                             agtgagGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGcc                                                 >  1:1211034/1‑140 (MQ=255)                                                                             agtgagGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTACACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGcc                                                 >  2:792281/1‑140 (MQ=255)                                                                              gtgagGAAAAAAC‑AGGCTGTTGATCCCAGTCGGAAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCt                                                <  2:280164/140‑1 (MQ=255)                                                                               tgagGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCtt                                               >  1:152526/1‑140 (MQ=255)                                                                               ggagGAAAAAAC‑GGGCTGTTGATACCAGTCGGCAAATCCCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCtt                                               <  2:29264/139‑1 (MQ=255)                                                                                  ggCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCg                                                   <  1:150487/133‑1 (MQ=255)                                                                                  ggCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCg                                                   >  2:150487/1‑133 (MQ=255)                                                                                  ggCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGa                                            >  2:990664/1‑140 (MQ=255)                                                                                  ggCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGa                                            >  1:954205/1‑140 (MQ=255)                                                                                  ggCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGa                                            >  1:1213195/1‑140 (MQ=255)                                                                                  ggCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGa                                            >  1:996525/1‑140 (MQ=255)                                                                                  ggCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGa                                            >  1:658738/1‑140 (MQ=255)                                                                                   gCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGAc                                           >  1:179583/1‑140 (MQ=255)                                                                                   aaaaaaaC‑AGGCTGTTGATCCCAGTCGGAAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGAc                                           <  2:954205/138‑1 (MQ=255)                                                                                     aaaaaC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACaa                                         >  2:1515965/1‑140 (MQ=255)                                                                                     aaaaaC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACaa                                         <  1:1601412/140‑1 (MQ=255)                                                                                     aaaaaC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTACACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACaa                                         <  1:792281/140‑1 (MQ=255)                                                                                      aaaaC‑AGGCTTTTGATCCCAGTCGGCAAATCCCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAAt                                        <  2:1564693/140‑1 (MQ=255)                                                                                        aaC‑AGGTTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATaa                                      <  2:109507/140‑1 (MQ=255)                                                                                        aaC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATaa                                      <  1:283576/140‑1 (MQ=255)                                                                                                 gTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCCATCGGAACGCTGACGTTCCGCACGGTGTTCAGCATTTTGCAGCCCGGGATGCCCGCCTTAGACAATAACCCCACaa                              >  2:423891/1‑140 (MQ=255)                                                                                                    gATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGc                           >  2:939736/1‑140 (MQ=255)                                                                                                           tccGGAAAATCCCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCa                    <  1:719188/138‑1 (MQ=255)                                                                                                            tCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGcc                     >  2:1638310/1‑138 (MQ=255)                                                                                                            tCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGcc                     <  1:1638310/138‑1 (MQ=255)                                                                                                            tCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAg                   >  2:794710/1‑140 (MQ=255)                                                                                                            tCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAg                   >  1:967427/1‑140 (MQ=255)                                                                                                            tCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAg                   >  1:812788/1‑140 (MQ=255)                                                                                                            tCGGCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAg                   >  2:1357347/1‑140 (MQ=255)                                                                                                              ggCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAGcc                 >  2:56669/1‑140 (MQ=255)                                                                                                              ggCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACACATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAGcc                 >  1:382698/1‑140 (MQ=255)                                                                                                               gCAAATACCGCCGTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAGccc                >  1:574426/1‑140 (MQ=255)                                                                                                                     cccgccgTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAGCCCCTGGCa          <  1:1357347/139‑1 (MQ=255)                                                                                                                        gccgTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAGCCCCTGGCAAgc       >  1:1128527/1‑140 (MQ=255)                                                                                                                         ccgTTAACGGTTCTGCGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAGCCCCTGGCAAgcc      >  2:557159/1‑140 (MQ=255)                                                                                                                             taccGGTTCTGCGAGAAACCGCTGCACCCATTGGCGATGGGAACGCTGCCGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATCCCCCCCTTCGACAATAACCCCACAAACCCCCGCCAGCCCCTGGCAAGCCGCca  <  2:812788/137‑1 (MQ=255)                                                                                                                                        |                                                                                                                                 CTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAAC‑AGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGTTGTTCAGCATTTTGCAGCCCCGGATGCCCGCCTTCGACAATAACCCCACAAAGCCCCGCCAGCCCCTGGCAAGCCGCCA  >  NC_000913/2377050‑2377312

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑