Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,513,844	T→C	A155A (GCA→GCG)	fecB ←	ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,513,844	0	T	C	100.0%	198.5 / NA	56	A155A (GCA→GCG)	fecB	ferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/29);  total (27/29)

CAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGC  >  NC_000913/4513700‑4513989                                                                                                                                                 |                                                                                                                                                  cAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGc                                                                                                                                             <  1:131509/151‑1 (MQ=255) cAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGc                                                                                                                                             >  2:668835/1‑151 (MQ=255)   gagagGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCAt                                                                                                                                           >  1:324412/1‑151 (MQ=255)       ggCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCAtctc                                                                                                                                        <  1:493150/150‑1 (MQ=255)         ccAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCg                                                                                                                                     <  2:477789/151‑1 (MQ=255)          cAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGc                                                                                                                                    <  1:102063/151‑1 (MQ=255)          cAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGc                                                                                                                                    <  2:408521/151‑1 (MQ=255)            gCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCtt                                                                                                                                  <  2:130984/151‑1 (MQ=255)                gCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCttttt                                                                                                                               >  2:351373/1‑150 (MQ=255)                   gCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTAccc                                                                                                                           <  1:560601/151‑1 (MQ=255)                    ccGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTAccc                                                                                                                           >  1:485711/1‑150 (MQ=255)                      ggTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACccacc                                                                                                                        <  1:189846/151‑1 (MQ=255)                          cAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCAttt                                                                                                                    >  2:33663/1‑151 (MQ=255)                          cAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCAttt                                                                                                                    <  2:177482/151‑1 (MQ=255)                            ggTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCg                                                                                                                  <  2:634925/151‑1 (MQ=255)                              tctcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGcc                                                                                                                >  2:211866/1‑151 (MQ=255)                                 ccTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgat                                                                                                             >  1:423426/1‑151 (MQ=255)                                      gTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCACCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAg                                                                                                        <  1:668777/151‑1 (MQ=255)                                          gCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGc                                                                                                    >  1:666013/1‑151 (MQ=255)                                                tGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGAtt                                                                                               >  2:720720/1‑150 (MQ=255)                                                    ttGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCaa                                                                                           >  1:3401/1‑150 (MQ=255)                                                     tGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCaaa                                                                                          <  1:774746/150‑1 (MQ=255)                                                     tGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAAt                                                                                         >  1:559930/1‑151 (MQ=255)                                                      gCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAAtt                                                                                        >  2:440058/1‑151 (MQ=255)                                                       cTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAAtt                                                                                        <  1:209560/150‑1 (MQ=255)                                                           tCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCa                                                                                    <  1:171730/150‑1 (MQ=255)                                                            cgcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCa                                                                                    <  1:343274/149‑1 (MQ=255)                                                             gcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCg                                                                                 <  1:637525/151‑1 (MQ=255)                                                               gTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTa                                                                               <  2:423477/151‑1 (MQ=255)                                                                      gCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTc                                                                         <  1:651490/150‑1 (MQ=255)                                                                      gCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCg                                                                        <  1:332355/151‑1 (MQ=255)                                                                      gCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCg                                                                        >  2:194151/1‑151 (MQ=255)                                                                         aaaGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTg                                                                     <  1:513310/151‑1 (MQ=255)                                                                         aaaGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTg                                                                     >  2:199595/1‑151 (MQ=255)                                                                           aGGCCACTCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCg                                                                   >  1:160860/1‑151 (MQ=255)                                                                           aGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCATTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCg                                                                   >  2:715593/1‑151 (MQ=255)                                                                              ccACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCggg                                                                 >  2:128345/1‑150 (MQ=255)                                                                               cACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGAc                                                               <  1:211834/151‑1 (MQ=255)                                                                                aCGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACt                                                              >  1:277247/1‑151 (MQ=255)                                                                                aCGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACt                                                              >  2:459324/1‑151 (MQ=255)                                                                                    gtgtCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTaa                                                           <  1:33632/150‑1 (MQ=255)                                                                                            ttGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGcagcag                                                    <  1:677343/149‑1 (MQ=255)                                                                                            ttGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGt                                                   >  2:153464/1‑150 (MQ=255)                                                                                               ggAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTa                                                  <  1:246085/148‑1 (MQ=255)                                                                                                gAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCgg                                              >  2:753612/1‑151 (MQ=255)                                                                                                  aGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGgcg                                            <  2:225331/151‑1 (MQ=255)                                                                                                    ctggctggCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGcgcg                                          >  1:453274/1‑151 (MQ=255)                                                                                                                  gcgcCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAgg                             <  2:494256/150‑1 (MQ=255)                                                                                                                     ccATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGAt                         <  2:704655/151‑1 (MQ=255)                                                                                                                         ccTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTa                      <  2:547609/150‑1 (MQ=255)                                                                                                                          ctctctTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTaa                     >  2:72202/1‑150 (MQ=255)                                                                                                                                   tgttgtTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCAt           >  2:24647/1‑151 (MQ=255)                                                                                                                                   tgttgtTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCAt           >  1:487637/1‑151 (MQ=255)                                                                                                                                      tgttCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATgg         <  2:43694/150‑1 (MQ=255)                                                                                                                                        ttCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCg       >  1:646559/1‑150 (MQ=255)                                                                                                                                            aGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGActgc  >  1:93493/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGC  >  NC_000913/4513700‑4513989

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