Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	531,370	A→C	D284E (GAT→GAG)	allS ←	allD operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	531,370	0	A	C	100.0%	222.5 / NA	57	D284E (GAT→GAG)	allS	allD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/30);  total (27/30)

CATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAAC  >  NC_000913/531224‑531516                                                                                                                                                   |                                                                                                                                                   cATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCtctt                                                                                                                                                >  1:100043/1‑151 (MQ=255)       tAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCg                                                                                                                                          >  1:715047/1‑151 (MQ=255)       tAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCg                                                                                                                                          >  1:767132/1‑151 (MQ=255)       tAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCg                                                                                                                                          >  1:919400/1‑151 (MQ=255)        aaaaGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCg                                                                                                                                          <  1:874934/150‑1 (MQ=255)          aaGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCtt                                                                                                                                       >  2:166122/1‑151 (MQ=255)              gAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGc                                                                                                                                    <  1:612129/150‑1 (MQ=255)              gAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGc                                                                                                                                    <  1:579522/150‑1 (MQ=255)               aaCCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCg                                                                                                                                  <  2:595068/151‑1 (MQ=255)                   tCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGc                                                                                                                              >  1:22419/1‑151 (MQ=255)                    cAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGc                                                                                                                              <  2:478716/150‑1 (MQ=255)                         cAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAAtt                                                                                                                        >  2:861329/1‑151 (MQ=255)                          aGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAAttt                                                                                                                       >  1:711363/1‑151 (MQ=255)                                 gcgcTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCAt                                                                                                                <  1:709125/151‑1 (MQ=255)                                    cTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGcc                                                                                                             >  2:387349/1‑151 (MQ=255)                                      aTTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCa                                                                                                            <  1:591125/150‑1 (MQ=255)                                        tttACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGa                                                                                                          <  1:452579/150‑1 (MQ=255)                                         ttACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACt                                                                                                        >  2:523100/1‑151 (MQ=255)                                              cgcCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATg                                                                                                   <  1:823111/151‑1 (MQ=255)                                               gcCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATgg                                                                                                  <  1:652244/151‑1 (MQ=255)                                                 cAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCg                                                                                                <  1:13785/151‑1 (MQ=255)                                                    cgcgCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGa                                                                                               >  1:103800/1‑149 (MQ=255)                                                      cgcAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAagg                                                                                            >  1:16711/1‑150 (MQ=255)                                                       gcAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAagga                                                                                           <  1:456004/150‑1 (MQ=255)                                                        cAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAaggagg                                                                                         >  2:727589/1‑151 (MQ=255)                                                             aaaCTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCa                                                                                    <  2:553674/151‑1 (MQ=255)                                                                  tGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGtt                                                                               >  2:534631/1‑151 (MQ=255)                                                                     tttttAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTggg                                                                            >  1:174314/1‑151 (MQ=255)                                                                         tAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGAtt                                                                         <  2:626151/150‑1 (MQ=255)                                                                              gCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCgg                                                                   >  2:720714/1‑151 (MQ=255)                                                                              gCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCgg                                                                   >  2:670709/1‑151 (MQ=255)                                                                               cGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGc                                                                  >  1:493004/1‑151 (MQ=255)                                                                                  gggTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGa                                                               <  2:900009/151‑1 (MQ=255)                                                                                       gCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAg                                                           <  1:380049/150‑1 (MQ=255)                                                                                         cGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCAGCGTTGGGATTGCCCGGCTGACCAgttg                                                        <  2:24512/151‑1 (MQ=255)                                                                                           aaaaTTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGAAAACAAGGCCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGAGCAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAgttgtt                                                      >  2:669017/1‑151 (MQ=255)                                                                                             aaTTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGa                                                    <  2:658060/151‑1 (MQ=255)                                                                                                ttCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGAtt                                                  <  2:526675/150‑1 (MQ=255)                                                                                                    aaaaaTCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGa                                             <  1:156901/151‑1 (MQ=255)                                                                                                      aaaTCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATc                                           <  1:330071/151‑1 (MQ=255)                                                                                                      aaaTCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATc                                           <  2:839439/151‑1 (MQ=255)                                                                                                       aaTCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCa                                          <  1:580175/151‑1 (MQ=255)                                                                                                          ccGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCAtt                                        <  1:292219/150‑1 (MQ=255)                                                                                                             ttGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACt                                    >  2:506513/2‑151 (MQ=255)                                                                                                               gATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTgg                                  >  1:858768/1‑151 (MQ=255)                                                                                                                aTTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGc                                 >  1:133117/1‑151 (MQ=255)                                                                                                                  ttCCGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCa                                >  1:752010/1‑150 (MQ=255)                                                                                                                     cGGCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGc                            <  2:445478/151‑1 (MQ=255)                                                                                                                       gCCTGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCg                           >  1:540280/1‑150 (MQ=255)                                                                                                                          tGCGCTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGAttt                       >  2:220138/1‑151 (MQ=255)                                                                                                                             gcTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTgg                    <  1:745355/151‑1 (MQ=255)                                                                                                                              cTGGGTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGc                   >  1:362230/1‑151 (MQ=255)                                                                                                                                 ggTAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCaaa                <  2:725121/151‑1 (MQ=255)                                                                                                                                   tAAACAAGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCaaaa               >  2:867180/1‑150 (MQ=255)                                                                                                                                        aaGGTCACAATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCaa         <  1:166122/151‑1 (MQ=255)                                                                                                                                             cacaATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGcc     <  1:727607/150‑1 (MQ=255)                                                                                                                                               caATCTCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAAc  <  2:612158/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CATAGTTAAAAGTGAACCTCAAATCAGCTTGCGCGCTATTTTACGCGCCAACGCGCAAGGAAACTTGATTTTTAACTGCGTGGGTTGCCGAAAATTTCTAAAAATCCGCTGATTTCCGGCCTGCGCTGGGTAAACAAGGTCACAATATCTTCTACCGCTTTGCCGCTGCCAAATTTGCGCCATGCCAGACTCAATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAAC  >  NC_000913/531224‑531516

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