Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,564,525	G→A	100%	M1M (ATG→ATA) †	yffS →	CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,564,525	0	G	A	100.0%	288.5 / NA	86	M1M (ATG→ATA) †	yffS	CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (52/34);  total (52/34)

GTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564428‑2564620                                                                                                  |                                                                                                gtgATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAg                                                                                               >  2:2329020/1‑100 (MQ=255)   gATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGtt                                                                                             >  2:363941/1‑100 (MQ=255)    aTCTGCGTCAATATTCCCTTCCGTTACGCTTACAATCCACCAATACGCGTGTTAACCTCCGAGTGGATTCTATGAGATGGGAGATAAAAATCATAAGTTa                                                                                            >  2:1969139/1‑100 (MQ=255)       tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                         >  1:158325/1‑100 (MQ=255)          gTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAAt                                                                                      <  1:1503151/100‑1 (MQ=255)           tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATc                                                                                     >  2:1361749/1‑100 (MQ=255)           tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATc                                                                                     >  2:461849/1‑100 (MQ=255)           tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATc                                                                                     >  1:1553611/1‑100 (MQ=255)            cAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCa                                                                                    >  1:1733948/1‑100 (MQ=255)             aaTATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCaa                                                                                   >  2:1972437/1‑100 (MQ=255)             aaTATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCaa                                                                                   >  1:1470349/1‑100 (MQ=255)                 ttCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAAtt                                                                               >  2:43442/1‑100 (MQ=255)                    ccTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGc                                                                            >  1:1241250/1‑100 (MQ=255)                     cTTCCGTTACGCTTGCCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           <  1:1821613/100‑1 (MQ=255)                     cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCTAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGc                                                                            >  2:1433781/1‑99 (MQ=255)                     cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           <  2:1499632/100‑1 (MQ=255)                      ttCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGa                                                                          <  2:45199/100‑1 (MQ=255)                       tCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGAc                                                                         <  1:1078343/100‑1 (MQ=255)                       tCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGAc                                                                         <  1:330277/100‑1 (MQ=255)                          gTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACAt                                                                       >  1:56209/1‑99 (MQ=255)                          gTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATc                                                                      >  1:1331121/1‑100 (MQ=255)                             aCGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTc                                                                   <  2:1775438/100‑1 (MQ=255)                              cgattACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  >  2:1046223/4‑100 (MQ=255)                              cGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  <  1:2275042/100‑1 (MQ=255)                                cTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTCTGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaaa                                                                >  2:1509211/1‑100 (MQ=255)                                cTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaaa                                                                <  2:1401468/100‑1 (MQ=255)                                cTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaaa                                                                >  2:1812357/1‑100 (MQ=255)                                   aCCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgc                                                             <  1:1651680/100‑1 (MQ=255)                                      cTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCt                                                          >  2:1705577/1‑100 (MQ=255)                                        ccaccaATACGCGTGTTAACATCCCAATGGATTCTAGGAGAAGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTAAAAgtcac                                                           >  2:1575179/1‑93 (MQ=255)                                         caccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCg                                                       >  2:2379672/1‑100 (MQ=255)                                          accaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGa                                                      >  1:1409331/1‑100 (MQ=255)                                          accaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGa                                                      <  2:1470333/100‑1 (MQ=255)                                           ccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAt                                                     >  1:223425/1‑100 (MQ=255)                                           ccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAt                                                     >  2:297836/1‑100 (MQ=255)                                            caATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATg                                                    <  1:1534489/100‑1 (MQ=255)                                              aTACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATgg                                                   >  1:1579868/1‑99 (MQ=255)                                               tACGCGTGTTAACCTCCGAATGGATTCGAGGAGAGGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGaa                                                 <  1:236694/100‑1 (MQ=255)                                                aCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAg                                                >  1:2364640/1‑100 (MQ=255)                                                 cgcgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAgg                                               >  1:2296117/1‑100 (MQ=255)                                                 cgcgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAgg                                               >  1:1010468/1‑100 (MQ=255)                                                     tgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACaa                                           >  1:1604340/1‑100 (MQ=255)                                                     tgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACaa                                           <  2:1362095/100‑1 (MQ=255)                                                      gtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAat                                          <  2:1679020/100‑1 (MQ=255)                                                      gtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTACGACATCCTCAAAGGCGCTGCGATGGAAGGACAat                                          <  2:2028143/100‑1 (MQ=255)                                                       ttAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAata                                         >  1:1654258/1‑100 (MQ=255)                                                       ttAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAat                                          <  1:162720/99‑1 (MQ=255)                                                         aaCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAata                                         >  2:1158289/1‑98 (MQ=255)                                                            cTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAgg                                    <  1:800289/100‑1 (MQ=255)                                                            cTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATAAAAATTTGCGACATCCTGAAAGGCACTGCGATGGAAGGACAATATAAg                                     >  2:1012076/1‑99 (MQ=255)                                                             tCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAggg                                   >  1:237393/1‑100 (MQ=255)                                                             tCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAggg                                   <  2:1698658/100‑1 (MQ=255)                                                              cccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgc                                  <  1:1178187/100‑1 (MQ=255)                                                              cccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgc                                  <  1:2379700/100‑1 (MQ=255)                                                                     gATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGAAATCCTCAAAGGCGATGAGATGGAAGGACAATATAAGGgcagt                               >  1:2221899/1‑93 (MQ=255)                                                                      atgcTATGAGATGGGAAATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCGAAGGCGCTGCTATGGAAGGACAATATAAGGGCGCtgcacg                          <  2:2363232/97‑6 (MQ=255)                                                                        tCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTg                        >  1:291324/1‑100 (MQ=255)                                                                          tATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTggg                      <  1:1345269/100‑1 (MQ=255)                                                                          tATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTggg                      <  1:1934588/100‑1 (MQ=255)                                                                          tATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTggg                      >  2:1925188/1‑100 (MQ=255)                                                                          tATGAGATAGGAGATAAAAATCATAAGTTACGATATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTggg                      <  1:54213/100‑1 (MQ=255)                                                                          naTGAGATGGGAGATAAAAATTATAAGATACGAAATCAAAATTTGAGACATCCGTAAAGGCGCTGCGATGGAAGGCAAATATAAGGGCGCTCAACg                          >  2:2201102/2‑96 (MQ=255)                                                                           aTGAGATGGGAGATACAACTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTgggg                     <  2:1428983/100‑1 (MQ=255)                                                                            tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    >  1:2113307/1‑100 (MQ=255)                                                                             gagaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCATCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGATCAACGTggggac                   >  2:1183996/1‑98 (MQ=255)                                                                              agaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCa                  <  1:1207402/100‑1 (MQ=255)                                                                              agaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCa                  <  2:291321/100‑1 (MQ=255)                                                                              agaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCa                  <  2:2299097/100‑1 (MQ=255)                                                                              agaAGGGAGACAAAAAGTATAAGGTACGAAACTAAAATTTACGCAGTCCTCAAAGGCGCTGCGCTGGATGTACAATATAACGGCGAGCACCGTGGAGc                    <  2:1495173/98‑1 (MQ=255)                                                                                aTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCaaa                <  2:1553593/100‑1 (MQ=255)                                                                                  gggAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATg              >  2:1504234/1‑100 (MQ=255)                                                                                  gggAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATg              >  1:1827470/1‑100 (MQ=255)                                                                                    gagaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCg            >  2:1265101/1‑100 (MQ=255)                                                                                    gagaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCg            >  1:686205/1‑100 (MQ=255)                                                                                      gaTAAAGACCATAAGTTAAGAAATAAAAATTTGCAACAACCTCAAAGAGACTGCGATGGCAGGAGAATATAAGGGCACTTAACGTGGGGCCAAATACgag          >  2:1093676/1‑100 (MQ=255)                                                                                      gaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          >  2:1582360/1‑100 (MQ=255)                                                                                      gaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          >  1:1778760/1‑100 (MQ=255)                                                                                      gaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          <  2:1602281/100‑1 (MQ=255)                                                                                       aTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagg         <  2:709968/100‑1 (MQ=255)                                                                                        tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  2:2357880/1‑100 (MQ=255)                                                                                        tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  2:220198/1‑100 (MQ=255)                                                                                        tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  1:2308763/1‑100 (MQ=255)                                                                                        tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  1:59173/1‑100 (MQ=255)                                                                                        tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAAAGAGGGGCCAAATAAgagga        >  1:1173644/1‑100 (MQ=255)                                                                                              tCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGc  <  1:911237/100‑1 (MQ=255)                                                                                              tCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGc  >  1:620189/1‑100 (MQ=255)                                                                                                  |                                                                                                GTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564428‑2564620

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