Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,533,553	A→G	L78L (TTG→CTG)	yjhQ ←	GNAT family putative N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,533,553	0	A	G	100.0%	219.7 / NA	57	L78L (TTG→CTG)	yjhQ	GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/34);  total (23/34)

GGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAAGAGGCCACTAAATCGGCTTCTTTG  >  NC_000913/4533407‑4533695                                                                                                                                                   |                                                                                                                                               ggggTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGt                                                                                                                                            <  1:526551/151‑1 (MQ=255)  gggTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTg                                                                                                                                           <  2:879740/151‑1 (MQ=255)      aaTAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAAt                                                                                                                                        <  1:565803/150‑1 (MQ=255)      aaTAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATc                                                                                                                                       >  1:119565/1‑151 (MQ=255)      aaTAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATc                                                                                                                                       <  2:404009/151‑1 (MQ=255)       aTAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATcc                                                                                                                                      <  1:903489/151‑1 (MQ=255)        tAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATcc                                                                                                                                      >  1:288102/1‑150 (MQ=255)         aGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCAt                                                                                                                                    <  1:561419/151‑1 (MQ=255)          ggTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCAtt                                                                                                                                   <  1:826591/151‑1 (MQ=255)             ggcgtggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTtct                                                                                                                                >  2:42466/1‑151 (MQ=255)                 tggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCtt                                                                                                                            >  1:789345/1‑151 (MQ=255)                 tggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCt                                                                                                                             <  2:379233/150‑1 (MQ=255)                  ggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCttt                                                                                                                           >  2:582414/1‑151 (MQ=255)                  ggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCtt                                                                                                                            <  1:528363/150‑1 (MQ=255)                   gcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTg                                                                                                                          <  2:474613/151‑1 (MQ=255)                      gagaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAg                                                                                                                       <  2:547764/151‑1 (MQ=255)                        gaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGt                                                                                                                     >  2:921809/1‑151 (MQ=255)                           cAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGc                                                                                                                  <  1:44278/151‑1 (MQ=255)                            aaaaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGcc                                                                                                                 >  1:954480/1‑151 (MQ=255)                              aaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCg                                                                                                               <  1:725688/151‑1 (MQ=255)                              aaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCg                                                                                                               <  1:339848/151‑1 (MQ=255)                                      tGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaat                                                                                                        >  2:167842/1‑150 (MQ=255)                                       ggCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaata                                                                                                       >  1:451129/1‑150 (MQ=255)                                        gCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaataaa                                                                                                     <  1:370227/151‑1 (MQ=255)                                          aGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaat                                                                                                   >  1:783037/1‑151 (MQ=255)                                          aGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaat                                                                                                   >  2:109712/1‑151 (MQ=255)                                           gCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaata                                                                                                  <  1:79498/151‑1 (MQ=255)                                                 cAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGccc                                                                                             <  2:498037/150‑1 (MQ=255)                                                 cAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCg                                                                                            <  2:726635/151‑1 (MQ=255)                                                 cAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCg                                                                                            >  2:257357/1‑151 (MQ=255)                                                    ccTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACg                                                                                         <  2:338730/151‑1 (MQ=255)                                                       cAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACgg                                                                                        <  2:342348/149‑1 (MQ=255)                                                               cGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTc                                                                              <  1:925195/151‑1 (MQ=255)                                                                  tGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCAt                                                                            >  1:614209/1‑150 (MQ=255)                                                                   gCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATAc                                                                          >  2:724463/1‑151 (MQ=255)                                                                                aaCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCg                                                             >  1:579627/1‑151 (MQ=255)                                                                                  ccGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGAt                                                           <  2:939122/151‑1 (MQ=255)                                                                                        tccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAAcgcg                                                     >  1:70814/1‑151 (MQ=255)                                                                                         cccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAAcgc                                                      <  1:828679/149‑1 (MQ=255)                                                                                         cccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCgg                                                    <  2:582238/151‑1 (MQ=255)                                                                                             cGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGgcg                                                <  2:860920/151‑1 (MQ=255)                                                                                             cGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGgc                                                 <  1:585732/150‑1 (MQ=255)                                                                                                  aTCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGgcgcgcgc                                           >  2:820589/1‑151 (MQ=255)                                                                                                      cTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTc                                       >  1:945583/1‑151 (MQ=255)                                                                                                          gTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATcc                                   >  1:880783/1‑151 (MQ=255)                                                                                                                cGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGt                             <  1:463741/151‑1 (MQ=255)                                                                                                                 ggTATGACCGCCAGGGTAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTa                            <  1:820459/151‑1 (MQ=255)                                                                                                                    aTGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTaaa                          <  1:46524/150‑1 (MQ=255)                                                                                                                     tGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTaa                           >  1:270986/1‑148 (MQ=255)                                                                                                                     tGACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTaa                           >  1:393270/1‑148 (MQ=255)                                                                                                                      gACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAagag                       <  2:995412/151‑1 (MQ=255)                                                                                                                      gACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAaga                        <  2:451076/150‑1 (MQ=255)                                                                                                                      gACCGCCAGGGGAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAaga                        <  2:946052/150‑1 (MQ=255)                                                                                                                      gACCGCCAGGGGAGCGAGGATGTGCATCAGCGGAGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCCCAAGTAAagtg                       >  2:774534/1‑149 (MQ=255)                                                                                                                               gggAGCGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAAGAGGCCACTaaa              <  2:411658/151‑1 (MQ=255)                                                                                                                                    cGAGGATGTGCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAAGAGGCCACTAAATCGGc         >  2:660268/1‑151 (MQ=255)                                                                                                                                           gtgCATCAGCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAAGAGGCCACTAAATCGGCTTCTTTg  <  1:740630/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               GGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAAGAGGCCACTAAATCGGCTTCTTTG  >  NC_000913/4533407‑4533695

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