Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,484,500	A→C	K4T (AAA→ACA)	prkB →	putative phosphoribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,484,500	0	A	C	100.0%	223.8 / NA	59	K4T (AAA→ACA)	prkB	putative phosphoribulokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (31/28);  total (31/28)

ACAGTGCGGAGAAGCGGTGCTGGTATGGTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCAAACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTCACCGTTACACCCGCCCGGAAATGG  >  NC_000913/3484355‑3484646                                                                                                                                                  |                                                                                                                                                   aCAGTGCGGAGAAGCGGTGCTGGTATGGTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATc                                                                                                                                               <  1:114500/151‑1 (MQ=255)  cAGTGCGGAGAAGCGGTGCTGGTATGGTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATcc                                                                                                                                              <  1:17783/151‑1 (MQ=255)       cGGAGAAGCGGTGCTGGTATGGTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCa                                                                                                                                         <  2:227281/151‑1 (MQ=255)        ggAGAAGCGGTGCTGGTATGGTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCAt                                                                                                                                        <  2:841261/151‑1 (MQ=255)                ggTGCTGGTATGGTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCCAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGc                                                                                                                                     >  2:291138/1‑146 (MQ=255)                 gTGCTGGTATGGTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAAc                                                                                                                               >  2:928512/1‑151 (MQ=255)                     tGGAATGGTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCAGGTCATTGCGGTAACAgg                                                                                                                            >  1:371223/1‑150 (MQ=255)                        tATGGTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATcc                                                                                                                        >  2:366361/1‑151 (MQ=255)                            gTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCTGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAgcg                                                                                                                    >  1:310832/1‑151 (MQ=255)                              cGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAgcggc                                                                                                                  <  2:76511/151‑1 (MQ=255)                                 aGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGcgcg                                                                                                               <  2:909289/151‑1 (MQ=255)                                      cACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGac                                                                                                          <  2:515818/151‑1 (MQ=255)                                           aaCGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGaccacc                                                                                                      >  1:407840/1‑150 (MQ=255)                                             cGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGaccaccacc                                                                                                   >  1:765417/1‑151 (MQ=255)                                                 cAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGAccaccaccacca                                                                                               <  1:693938/151‑1 (MQ=255)                                                  aaTATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGaccaccaccagcct                                                                                              <  2:165489/151‑5 (MQ=255)                                                  aaTATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGaccaccaccacc                                                                                                >  1:640054/1‑149 (MQ=255)                                                   atatTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGc                                                                                             >  2:375927/1‑151 (MQ=255)                                                       tatGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTcg                                                                                         >  2:542694/1‑151 (MQ=255)                                                       tatGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTc                                                                                          <  1:388536/150‑1 (MQ=255)                                                       tatGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACAACCAGCCTcg                                                                                         >  1:442720/1‑151 (MQ=255)                                                         tGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTcgcg                                                                                       >  1:904425/1‑151 (MQ=255)                                                           ccGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGtt                                                                                     >  2:320994/1‑151 (MQ=255)                                                             gcgcAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTc                                                                                   <  1:388146/151‑1 (MQ=255)                                                              cgcAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCg                                                                                  <  2:838768/151‑1 (MQ=255)                                                                cagtcagtTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTa                                                                                <  2:275171/151‑1 (MQ=255)                                                                  gtcagtTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTaa                                                                               >  1:613027/1‑150 (MQ=255)                                                                    cagtTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTaaaa                                                                             <  1:601002/150‑1 (MQ=255)                                                                           gcgAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTcgcg                                                                     <  2:170616/151‑1 (MQ=255)                                                                             gAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCa                                                                   <  2:293533/151‑1 (MQ=255)                                                                             gAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGCAAAATATTCgcgc                                                                    >  1:159637/1‑150 (MQ=255)                                                                                tAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGtt                                                                >  2:914382/1‑151 (MQ=255)                                                                                  aCCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTa                                                               <  1:162957/150‑1 (MQ=255)                                                                                   ccTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTaaa                                                             <  1:362793/151‑1 (MQ=255)                                                                                    cTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAAt                                                            >  1:348598/1‑151 (MQ=255)                                                                                    cTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGATAAAt                                                            >  1:672374/1‑151 (MQ=255)                                                                                         aCCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCtgc                                                        >  2:243811/1‑150 (MQ=255)                                                                                          ccGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCtgcat                                                      <  2:74101/151‑1 (MQ=255)                                                                                            gTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCtgcatgc                                                    <  1:514748/151‑1 (MQ=255)                                                                                                gcggACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCt                                                <  1:418948/151‑1 (MQ=255)                                                                                                 cggACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTg                                               >  1:87665/1‑151 (MQ=255)                                                                                                 cggACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTg                                               >  1:564137/1‑151 (MQ=255)                                                                                                    aCGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAgg                                            <  1:479651/151‑1 (MQ=255)                                                                                                      gagaTAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTg                                          >  2:875468/1‑151 (MQ=255)                                                                                                       agaTAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTgg                                         >  2:420933/1‑151 (MQ=255)                                                                                                         ataataATGATAGCTAAGTCATGGAGTTACTGTGTCTGCCACACATCCGGTCATTACGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGCTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGCGGTg                                          >  1:379211/1‑148 (MQ=255)                                                                                                         ataataATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTg                                          >  2:277854/1‑148 (MQ=255)                                                                                                         ataataATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGaa                                       <  1:460809/151‑1 (MQ=255)                                                                                                         ataataATGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGa                                        >  1:833232/1‑150 (MQ=255)                                                                                                              aatGATAGCTAAGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGa                                  <  1:204764/151‑1 (MQ=255)                                                                                                                        aaGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTc                          <  2:250233/149‑1 (MQ=255)                                                                                                                        aaGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTCAc                        >  1:50744/1‑151 (MQ=255)                                                                                                                         aGTCATGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTCAcc                       <  1:471633/151‑1 (MQ=255)                                                                                                                             aTGGAGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTCACCGtt                    >  2:597483/1‑150 (MQ=255)                                                                                                                                 aGTTACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTCACCGTTACAcc               <  2:338641/151‑1 (MQ=255)                                                                                                                                    tACTATGTCTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTCACCGTTACAcccgc            >  2:287108/1‑151 (MQ=255)                                                                                                                                            cTGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTCACCGTTACACCCGCCCGGAAAt    >  2:273292/1‑151 (MQ=255)                                                                                                                                             tGCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTAGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTCACCGTTACACCCGCCCGGAAATg   >  1:622323/1‑151 (MQ=255)                                                                                                                                              gCCACACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTCACCGTTACACCCGCCCGGAAATgg  <  2:215251/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   ACAGTGCGGAGAAGCGGTGCTGGTATGGTCGGAGCTGCACGAAACGGTCAATATTATGCCGCGCAGTCAGTTTCGCGAATAACCTGCAACCGTCGGCGGACGAGATAATAATGATAGCTAAGTCATGGAGTTACTATGTCTGCCAAACATCCGGTCATTGCGGTAACAGGATCCAGCGGCGCGGGGACCACCACCACCAGCCTCGCGTTTCGTAAAATATTCGCGCAGTTAAATCTGCATGCAGCTGAGGTGGAAGGCGACAGTTTTCACCGTTACACCCGCCCGGAAATGG  >  NC_000913/3484355‑3484646

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