Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,422,240	T→G	L178R (CTG→CGG)	yfcH →	putative NAD‑dependent nucleotide‑sugar epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,422,240	0	T	G	100.0%	202.2 / NA	56	L178R (CTG→CGG)	yfcH	putative NAD‑dependent nucleotide‑sugar epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/31);  total (25/31)

AGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTA  >  NC_000913/2422093‑2422384                                                                                                                                                    |                                                                                                                                                 agGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGgcc                                                                                                                                               <  1:591175/151‑1 (MQ=255)    aaCCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGgccgcc                                                                                                                                            <  1:209159/151‑1 (MQ=255)          cgcATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCg                                                                                                                                      >  2:337339/1‑151 (MQ=255)            cATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGAACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTc                                                                                                                                    >  1:48469/1‑151 (MQ=255)                      ttACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTgg                                                                                                                           <  2:501553/150‑1 (MQ=255)                       tACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTgggc                                                                                                                         >  2:44768/1‑151 (MQ=255)                            aTAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGcc                                                                                                                    <  2:842123/151‑1 (MQ=255)                               aaCTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGAt                                                                                                                 >  1:162507/1‑151 (MQ=255)                                   ctGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTCGCCTGGGCGGGCCGATTGGc                                                                                                             >  1:693030/1‑151 (MQ=255)                                   ctGCGCTCGCTGGGAAGAAATTGCCTGACGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTg                                                                                                               <  2:480524/149‑1 (MQ=255)                                     gcgcTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTgg                                                                                                              <  2:144806/148‑1 (MQ=255)                                          cgctGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGt                                                                                                       <  1:872624/150‑1 (MQ=255)                                           gctGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCg                                                                                                     >  2:669227/1‑151 (MQ=255)                                              gggAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGc                                                                                                   >  2:781026/1‑150 (MQ=255)                                              gggAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCa                                                                                                  <  2:876536/151‑1 (MQ=255)                                              gggAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCa                                                                                                  <  1:559626/151‑1 (MQ=255)                                              gggAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCAGGTCGGCa                                                                                                  <  1:49549/151‑1 (MQ=255)                                               ggAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGc                                                                                                   >  1:297458/1‑149 (MQ=255)                                                      aTTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTgg                                                                                          <  2:784501/151‑1 (MQ=255)                                                      aTTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTgg                                                                                          <  2:291700/151‑1 (MQ=255)                                                        tgcctgccGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGTCGTTTCGGCTTGGCCTGGGCGGGCCGATTGGCTCTGTTCGGCAGTATCTGGcc                                                                                        >  1:567770/1‑151 (MQ=255)                                                          cctgccGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGcc                                                                                        <  1:869232/149‑1 (MQ=255)                                                                gggCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTc                                                                                 >  2:751076/1‑150 (MQ=255)                                                                gggCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTc                                                                                 >  2:489318/1‑150 (MQ=255)                                                                  gcgcAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCata                                                                              <  1:122884/151‑1 (MQ=255)                                                                        aGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATtgatg                                                                        <  2:198947/151‑1 (MQ=255)                                                                             cAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATg                                                                   >  1:144698/1‑151 (MQ=255)                                                                              aaaaCGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATgg                                                                  >  1:601471/1‑151 (MQ=255)                                                                                aaCGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTc                                                                >  2:577416/1‑151 (MQ=255)                                                                                  cgcgAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCa                                                               <  1:452337/150‑1 (MQ=255)                                                                                  cgcgAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCCGGCCTGGATTCATATTGATGATATGGTCa                                                               >  2:366953/1‑150 (MQ=255)                                                                                       gtgtGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTa                                                         <  1:1712/151‑1 (MQ=255)                                                                                        tgtgCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTa                                                         <  1:374760/150‑1 (MQ=255)                                                                                         gtgCCTGCTGCGTACCGGTGTAGTGCTGGCCCCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTAt                                                        <  1:686491/150‑1 (MQ=255)                                                                                           gCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATtc                                                      <  2:897150/150‑1 (MQ=255)                                                                                                      aCCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTgg                                          >  1:60642/1‑151 (MQ=255)                                                                                                        cGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGAt                                        >  2:406110/1‑151 (MQ=255)                                                                                                          gtgtAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATaa                                      <  2:17591/151‑1 (MQ=255)                                                                                                          gtgtAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATa                                       >  2:51713/1‑150 (MQ=255)                                                                                                            gtAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACg                                    >  2:623341/1‑151 (MQ=255)                                                                                                               gTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGc                                 >  2:62342/1‑151 (MQ=255)                                                                                                                    ggCACCGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTgc                              <  2:342202/149‑1 (MQ=255)                                                                                                                        ccGGATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGc                        <  1:891533/151‑1 (MQ=255)                                                                                                                          ggATGGCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGcc                       >  1:300331/1‑150 (MQ=255)                                                                                                                               gCGGTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTaa                 <  1:104778/151‑1 (MQ=255)                                                                                                                                 ggTATTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAat                >  1:689569/1‑150 (MQ=255)                                                                                                                                   tATGCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAatat              <  1:862194/150‑1 (MQ=255)                                                                                                                                    aTTCTCGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATgg            >  1:510838/1‑151 (MQ=255)                                                                                                                                        tcGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTc        >  2:563052/1‑151 (MQ=255)                                                                                                                                         cGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGtt          <  2:848108/148‑1 (MQ=255)                                                                                                                                         cGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTc        <  2:612983/150‑1 (MQ=255)                                                                                                                                         cGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTc        <  1:365807/150‑1 (MQ=255)                                                                                                                                         cGGTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTc        <  2:267029/150‑1 (MQ=255)                                                                                                                                          ggTAAAATGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTc        <  1:340669/149‑1 (MQ=255)                                                                                                                                              aaaTGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCt   >  1:569863/1‑150 (MQ=255)                                                                                                                                              aaaTGCGGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTa  <  1:673067/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 AGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATTGGCTCTGGTCGGCAGTATCTGGCCTGGATTCATATTGATGATATGGTCAACGGTATTCTCTGGCTGCTGGATAACGAGCTGCGCGGGCCATTTAATATGGTTTCGCCCTA  >  NC_000913/2422093‑2422384

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