Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,508,973	C→A	E281* (GAA→TAA)	yhfW ←	putative mutase YhfW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,973	0	C	A	100.0%	221.7 / NA	63	E281* (GAA→TAA)	yhfW	putative mutase YhfW
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (34/29);  total (34/29)

CTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCACGCCTGCTTCATACAGTTTTTGTG  >  NC_000913/3508839‑3509104                                                                                                                                       |                                                                                                                                    cTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGa                                                                                                                                >  1:427818/1‑140 (MQ=255) cTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGa                                                                                                                                >  1:360744/1‑140 (MQ=255)   caacaaGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAgg                                                                                                                              <  1:648509/140‑1 (MQ=255)   caacaaGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGCTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAgg                                                                                                                              <  1:596150/140‑1 (MQ=255)     acaaGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGt                                                                                                                            >  2:135690/1‑140 (MQ=255)       aaGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGa                                                                                                                          <  2:427818/140‑1 (MQ=255)       aaGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGa                                                                                                                          <  1:515824/140‑1 (MQ=255)             ggCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCa                                                                                                                    >  2:295036/1‑140 (MQ=255)             ggCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCa                                                                                                                    >  2:766212/1‑140 (MQ=255)              gCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCAt                                                                                                                   <  2:360744/140‑1 (MQ=255)              gCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCAt                                                                                                                   >  1:519317/1‑140 (MQ=255)                            aaCGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTAt                                                                                                     >  2:789887/1‑140 (MQ=255)                            aaCGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTAt                                                                                                     >  1:134461/1‑140 (MQ=255)                            aaCGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTAt                                                                                                     >  2:650235/1‑140 (MQ=255)                              cGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATcc                                                                                                   >  1:496190/1‑140 (MQ=255)                                aCCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATccac                                                                                                 <  1:210106/140‑1 (MQ=255)                                    gCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGa                                                                                             >  1:349761/1‑140 (MQ=255)                                    gCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGa                                                                                             >  2:367981/1‑140 (MQ=255)                                           ttCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGcc                                                                                      >  1:29732/1‑140 (MQ=255)                                             cGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAg                                                                                    >  2:662665/1‑140 (MQ=255)                                               gCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCt                                                                                  <  1:789887/140‑1 (MQ=255)                                                cATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTc                                                                                 <  1:650235/140‑1 (MQ=255)                                                 aTAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCa                                                                                <  2:779209/140‑1 (MQ=255)                                                 aTAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCa                                                                                <  2:29732/140‑1 (MQ=255)                                                 aTAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCa                                                                                <  2:349761/140‑1 (MQ=255)                                                 aTAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCa                                                                                <  1:295036/140‑1 (MQ=255)                                                 aTAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCa                                                                                <  1:766212/140‑1 (MQ=255)                                                 aTAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCa                                                                                <  2:793608/140‑1 (MQ=255)                                                 aTAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCa                                                                                <  2:412935/140‑1 (MQ=255)                                                    aCGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGc                                                                             <  2:490236/140‑1 (MQ=255)                                                            cGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGa                                                                     <  2:496190/140‑1 (MQ=255)                                                                  ctGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGAttgttg                                                               >  1:3535/1‑140 (MQ=255)                                                                    gCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCAATAATCCGCTGGCTACCCACCAGATTTGGCCAGCTCACGCCATAAGGATTGTTGAc                                                             <  2:776209/140‑1 (MQ=255)                                                                        gACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGAtt                                                                   >  1:530350/1‑130 (MQ=255)                                                                        gACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGAtt                                                                   <  2:530350/130‑1 (MQ=255)                                                                        gACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGa                                                              <  2:888846/135‑1 (MQ=255)                                                                        gACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGa                                                              >  1:888846/1‑135 (MQ=255)                                                                         aCCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGatat                                                        >  2:280706/1‑140 (MQ=255)                                                                          ccAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATc                                                       >  1:266506/1‑140 (MQ=255)                                                                           cAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCt                                                      >  1:261570/1‑140 (MQ=255)                                                                            aGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTg                                                     <  2:519317/140‑1 (MQ=255)                                                                              cGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGcc                                                   >  1:487313/1‑140 (MQ=255)                                                                               gAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCAGGCCATAAGGATTGTTGACGATATCTGCCa                                                  >  1:787400/1‑140 (MQ=255)                                                                               gAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCa                                                  <  2:487313/140‑1 (MQ=255)                                                                                     ggTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTAc                                            >  2:342669/1‑140 (MQ=255)                                                                                      gTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTAcc                                           <  2:301810/140‑1 (MQ=255)                                                                                            tGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAg                                     >  1:681575/1‑140 (MQ=255)                                                                                                 gTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACgg                                >  2:875621/1‑140 (MQ=255)                                                                                                   tGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGtt                              >  2:122696/1‑140 (MQ=255)                                                                                                   tGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGtt                              >  1:19916/1‑140 (MQ=255)                                                                                                    ggTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTg                             >  1:743728/1‑140 (MQ=255)                                                                                                    ggTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTg                             >  2:38588/1‑140 (MQ=255)                                                                                                     gTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTgg                            >  2:632384/1‑140 (MQ=255)                                                                                                     gTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTgg                            <  1:280706/140‑1 (MQ=255)                                                                                                     gTGCAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTgg                            <  1:875621/140‑1 (MQ=255)                                                                                                        cAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTgg                            <  1:851685/137‑1 (MQ=255)                                                                                                        cAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTgg                            >  2:851685/1‑137 (MQ=255)                                                                                                        cAAATAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCAc                         >  1:134160/1‑140 (MQ=255)                                                                                                          aaTAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCACGc                       <  1:36323/140‑1 (MQ=255)                                                                                                            tAAACGCCGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCACGCCt                     <  1:38588/140‑1 (MQ=255)                                                                                                                  ccGTCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCACGCCTGCTTCa               <  1:833824/140‑1 (MQ=255)                                                                                                                     tCGGATGGGTGTTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCACGCCTGCTTCATAc            >  1:491592/1‑140 (MQ=255)                                                                                                                               gtTAAATTAGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCACGCCTGCTTCATACAGTTTTtgtg  >  2:382658/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                    CTCAACAAGCCGGGCAAGGTTACGGTCAACGACCTGCAAACGTTCGGCATAACGTGCGACGTCTTCTGCATGACCAGCGAGGTCGGTTTCCTGAATGTTGGTGCAAATAAACGCCGTCGGATGGGTGTTAAATTCGTTGAGGGTGATATCCATAATCCGCTGGCTATCCACCAGATTTTGCCAGCTCACGCCATAAGGATTGTTGACGATATCTGCCACCTTACCCACCAGCACGGTTGGCACGCCTGCTTCATACAGTTTTTGTG  >  NC_000913/3508839‑3509104

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