Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	679,874	G→A	G123R (GGG→AGG)	ybeU →	DUF1266 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	679,874	0	G	A	100.0%	215.7 / NA	67	G123R (GGG→AGG)	ybeU	DUF1266 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/34);  total (33/34)

CGCTCTTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAG  >  NC_000913/679739‑680009                                                                                                                                        |                                                                                                                                        cGCTCTTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGt                                                                                                                                      >  1:65304/1‑139 (MQ=255)   ctctTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGttt                                                                                                                                    <  2:65304/139‑1 (MQ=255)    tctTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTc                                                                                                                                   <  2:526624/139‑1 (MQ=255)             gAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACggg                                                                                                                          <  1:193985/139‑1 (MQ=255)              aaGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGt                                                                                                                         >  1:385143/1‑139 (MQ=255)                  gTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgt                                                                                                                     <  1:919601/139‑1 (MQ=255)                   tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                                    >  1:344329/1‑139 (MQ=255)                   tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                                    >  1:990317/1‑139 (MQ=255)                   tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                                    >  2:512723/1‑139 (MQ=255)                   tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                                    >  2:521124/1‑139 (MQ=255)                   tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                                                    >  1:270608/1‑139 (MQ=255)                     gAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACggg                                                                                                                          >  2:913757/1‑131 (MQ=255)                     gAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACggg                                                                                                                          <  1:913757/131‑1 (MQ=255)                                ttGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTg                                                                                                       >  2:362727/1‑139 (MQ=255)                                             ggAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGgaagaa                                                                                          <  2:990317/139‑1 (MQ=255)                                               agagagGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAagaaga                                                                                        <  2:385143/139‑1 (MQ=255)                                                     ggTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGttt                                                                                  >  1:893041/1‑139 (MQ=255)                                                                aCCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAAt                                                                       >  1:114446/1‑139 (MQ=255)                                                                  ccAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGACTGCAATcg                                                                     <  1:6297/139‑1 (MQ=255)                                                                   cAATACGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGCGTTTATGGCTGCAATcgc                                                                    <  1:707843/139‑1 (MQ=255)                                                                    aaTTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcg                                                                   >  1:988728/1‑139 (MQ=255)                                                                    aaTTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcg                                                                   >  1:458156/1‑139 (MQ=255)                                                                    aaTTCGTTGCACAAACTGCCATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGATGTTTCTGTCACGGGTCGGTGTGCTCAACACGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcg                                                                   <  2:344329/139‑1 (MQ=39)                                                                       tCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGt                                                                <  1:436442/139‑1 (MQ=255)                                                                       tCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGt                                                                <  2:988728/139‑1 (MQ=255)                                                                         gTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCt                                                              <  1:521124/139‑1 (MQ=255)                                                                           tGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTAt                                                                                <  2:121917/119‑1 (MQ=255)                                                                           tGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTAt                                                                                >  1:121917/1‑119 (MQ=255)                                                                           tGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTAt                                                            >  1:64580/1‑139 (MQ=255)                                                                           tGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTAt                                                            >  2:515573/1‑139 (MQ=255)                                                                                    cGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCt                                                   <  2:458156/139‑1 (MQ=255)                                                                                    cGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGGCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCt                                                   <  2:114446/139‑1 (MQ=255)                                                                                      cTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCa                                                 <  1:515573/139‑1 (MQ=255)                                                                                       tATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGt                                                                >  1:779838/1‑123 (MQ=255)                                                                                       tATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGt                                                                <  2:779838/123‑1 (MQ=255)                                                                                        aTGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTGCGAGCTCAcc                                               <  1:362727/139‑1 (MQ=255)                                                                                              tgtgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTc                                                                                                                               >  2:811704/1‑53 (MQ=255)                                                                                              tgtgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTc                                                                                                                               <  1:811704/53‑1 (MQ=255)                                                                                              tgtgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGttt                                                                                  >  2:702482/1‑98 (MQ=255)                                                                                              tgtgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGttt                                                                                  <  1:702482/98‑1 (MQ=255)                                                                                                     aaGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcg                                                                     <  2:696179/104‑1 (MQ=255)                                                                                                     aaGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcg                                                                     >  1:696179/1‑104 (MQ=255)                                                                                                     aaGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCa                                                 >  1:700336/1‑124 (MQ=255)                                                                                                     aaGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCa                                                 <  2:700336/124‑1 (MQ=255)                                                                                                       gggggAATCAAAGCCTGGGATTATGTTCTAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGt                                                                >  1:664733/1‑107 (MQ=255)                                                                                                       gggggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGt                                                                <  2:664733/107‑1 (MQ=255)                                                                                                       gggggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCa                                                 <  2:227250/122‑1 (MQ=255)                                                                                                       gggggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCa                                                 >  1:227250/1‑122 (MQ=255)                                                                                                       gggggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTAt                                       <  2:677345/132‑1 (MQ=255)                                                                                                       gggggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTAt                                       >  1:677345/1‑132 (MQ=255)                                                                                                          ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTc                                                               <  2:27393/105‑1 (MQ=255)                                                                                                          ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTc                                                               >  1:27393/1‑105 (MQ=255)                                                                                                          ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTa                                                         <  1:586073/111‑1 (MQ=255)                                                                                                          ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTa                                                         >  2:586073/1‑111 (MQ=255)                                                                                                            aaTCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcg                                                                     <  2:1038711/97‑1 (MQ=255)                                                                                                            aaTCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcg                                                                     >  1:1038711/1‑97 (MQ=255)                                                                                                                  aGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTcaacaa                                                                                                          <  2:694557/54‑1 (MQ=255)                                                                                                                  aGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTcaacaa                                                                                                          >  1:694557/1‑54 (MQ=255)                                                                                                                    ccTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCt                                                   >  2:436662/1‑107 (MQ=255)                                                                                                                    ccTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCt                                                   <  1:436662/107‑1 (MQ=255)                                                                                                                      tGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTAttttt                 >  1:612658/1‑139 (MQ=255)                                                                                                                        ggATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTcc               >  1:1025206/1‑139 (MQ=255)                                                                                                                              tGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcgc                                                                  <  2:776696/82‑1 (MQ=255)                                                                                                                              tGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcgc                                                                  >  1:776696/1‑82 (MQ=255)                                                                                                                                   gAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCt                                <  1:357489/111‑1 (MQ=255)                                                                                                                                   gAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCt                                >  2:357489/1‑111 (MQ=255)                                                                                                                                     aTGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAg  <  1:512723/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        CGCTCTTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAG  >  NC_000913/679739‑680009

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