Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	775,120	G→A	A123A (GCG→GCA)	ybgC →	acyl‑CoA thioester hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	775,120	0	G	A	100.0%	217.0 / NA	66	A123A (GCG→GCA)	ybgC	acyl‑CoA thioester hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/33);  total (33/33)

ACAGACTGAAATAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCGCTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCAT  >  NC_000913/774985‑775253                                                                                                                                        |                                                                                                                                      actgACTGAAAGAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTGTTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACtt                                                                                                                                    >  2:612197/4‑139 (MQ=255)        gAAATAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGt                                                                                                                             >  1:220698/1‑139 (MQ=255)        gAAATAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGt                                                                                                                             >  1:60718/1‑139 (MQ=255)         aaaTAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTccc                                                                                                                                 >  2:826531/1‑134 (MQ=255)         aaaTAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTccc                                                                                                                                 <  1:826531/134‑1 (MQ=255)          aaTAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCt                                                                                                                           >  1:527233/1‑139 (MQ=255)               cATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGt                                                                                                                      >  2:455874/1‑139 (MQ=255)               cATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGt                                                                                                                      >  2:712221/1‑139 (MQ=255)                 tCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcg                                                                                                                    <  1:322641/139‑1 (MQ=255)                  cAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   <  2:60718/139‑1 (MQ=255)                  cAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   <  2:220698/139‑1 (MQ=255)                  cAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   <  1:283156/139‑1 (MQ=255)                   aaTGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgcg                                                                                                                  >  2:529945/1‑139 (MQ=255)                    aTGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   >  2:355845/1‑137 (MQ=255)                    aTGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   <  1:355845/137‑1 (MQ=255)                          ggCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTa                                                                                                           >  1:908559/1‑139 (MQ=255)                          ggCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTa                                                                                                           >  2:971474/1‑139 (MQ=255)                              ccTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACTCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCa                                                                                                       <  1:375085/139‑1 (MQ=255)                               ctctTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   <  2:491792/126‑1 (MQ=255)                               ctctTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   >  1:491792/1‑126 (MQ=255)                                   ttGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtgac                                                                                                  <  2:80685/139‑1 (MQ=255)                                    tGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtgact                                                                                                 <  1:79019/139‑1 (MQ=255)                                         ttCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACAt                                                                                            >  1:46220/1‑139 (MQ=255)                                           cACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGa                                                                                          <  2:1001396/139‑1 (MQ=255)                                            aCGCAACGTATTGTCAACGCCGAGTATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGaa                                                                                         <  1:588830/139‑1 (MQ=255)                                                 aCGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGaa                                                                                         <  2:932462/134‑1 (MQ=255)                                                 aCGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGaa                                                                                         >  1:932462/1‑134 (MQ=255)                                                          cAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGt                                                                           >  1:941266/1‑139 (MQ=255)                                                            aCGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTc                                                                         >  1:451692/1‑139 (MQ=255)                                                            aCGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTc                                                                         >  1:214374/1‑139 (MQ=255)                                                            aCGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTc                                                                         >  2:24528/1‑139 (MQ=255)                                                              gCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCt                                                                       <  2:263795/139‑1 (MQ=255)                                                              gCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCt                                                                       <  1:529945/139‑1 (MQ=255)                                                                  agaATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATAAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAg                                                                   >  2:384771/1‑139 (MQ=255)                                                                   gaATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAgg                                                                  >  2:612564/1‑139 (MQ=255)                                                                        cTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   <  2:229727/85‑1 (MQ=255)                                                                        cTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   >  1:229727/1‑85 (MQ=255)                                                                          ttGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtgactgac                                                                                              >  1:1022950/1‑104 (MQ=255)                                                                          ttGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtgactgac                                                                                              <  2:1022950/104‑1 (MQ=255)                                                                          ttGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCt                                                                                   <  1:334852/115‑1 (MQ=255)                                                                          ttGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCt                                                                                   >  2:334852/1‑115 (MQ=255)                                                                          ttGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCt                                                           <  2:214374/139‑1 (MQ=255)                                                                                  tgaaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTa                                                   >  1:935055/1‑139 (MQ=255)                                                                                     aGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGt                                                     <  2:835400/134‑1 (MQ=255)                                                                                     aGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGt                                                     >  1:835400/1‑134 (MQ=255)                                                                                     aGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAAc                                                >  2:517239/1‑139 (MQ=255)                                                                                         gagGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTAt                                            <  1:24528/139‑1 (MQ=255)                                                                                          agGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGAt                                                                               >  2:532540/1‑103 (MQ=255)                                                                                          agGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGAt                                                                               <  1:532540/103‑1 (MQ=255)                                                                                                       ttGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgattttga                              <  1:612197/139‑1 (MQ=255)                                                                                                                    cTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGa                                                                                          >  1:844328/1‑66 (MQ=255)                                                                                                                    cTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGa                                                                                          <  2:844328/66‑1 (MQ=255)                                                                                                                    cTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCtt                                                          <  2:127569/98‑1 (MQ=255)                                                                                                                    cTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCtt                                                          >  1:127569/1‑98 (MQ=255)                                                                                                                           tGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCt          <  1:971474/139‑1 (MQ=255)                                                                                                                           tGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCt          <  2:46220/139‑1 (MQ=255)                                                                                                                                ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc     <  2:451692/139‑1 (MQ=255)                                                                                                                                ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc     <  1:455874/139‑1 (MQ=255)                                                                                                                                ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc     <  2:935055/139‑1 (MQ=255)                                                                                                                                  tCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCa   <  1:384771/139‑1 (MQ=255)                                                                                                                                   cGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtga                                                                                                   >  1:999283/1‑42 (MQ=255)                                                                                                                                   cGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtga                                                                                                   <  2:999283/42‑1 (MQ=255)                                                                                                                                   cGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTg                                                                            <  2:256213/65‑1 (MQ=255)                                                                                                                                   cGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTg                                                                            >  1:256213/1‑65 (MQ=255)                                                                                                                                   cGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCat  >  1:863200/1‑139 (MQ=255)                                                                                                                                   cGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCat  >  1:522405/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                      ACAGACTGAAATAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCGCTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCAT  >  NC_000913/774985‑775253

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