Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	646,861	G→A	L217F (CTC→TTC)	citG ←	triphosphoribosyl‑dephospho‑CoA synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,861	0	G	A	100.0%	120.0 / NA	35	L217F (CTC→TTC)	citG	triphosphoribosyl‑dephospho‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/16);  total (19/16)

CTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGG  >  NC_000913/646622‑647107                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATACCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGCCAGCGCAGGACCCCCCCCCCCCCCCGCGACG‑GAACGGTGGGATCGCCGGTGGTAGCCCTCAGTAGGAACAAgggat                                                                                                                                                                                                                                                >  1:309618/1‑245 (MQ=255)             aaTCTGTGCTAAAACACAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGtt                                                                                                                                                                                                                                 <  2:8003/251‑1 (MQ=255)                            ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCt                                                                                                                                                                                                                  <  1:125789/251‑1 (MQ=255)                                        cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGTAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCaga                                                                                                                                                                                                       <  2:144200/250‑1 (MQ=255)                                        cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCagag                                                                                                                                                                                                      <  1:67863/251‑1 (MQ=255)                                         aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGt                                                                                                                                                                                                     <  1:240623/251‑1 (MQ=255)                                                          gccTTGACTGAGATTGCGTGCGATACACTCCCTGTCGCACGGCCGGATATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCGCCTCGCCACCGCGCGATG‑CAACGTTTGTATCGTCGATGAGCGCCATCAGTATGAACAAGGTATCGAGCAAGGCCAGGGCAGGATTTAGCCCCTGACCCAGCAGATTGATGTAATGccaaaaga                                                                                                                                                                                    <  2:50332/251‑8 (MQ=255)                                                              ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAAt                                                                                                                                                                                             <  1:243930/238‑1 (MQ=255)                                                                aCTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTTAATAATTTTTTCTCCCCGCGCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGATGGCAA‑GTTGTGATCCGCGGTGATCGCCCTCATGAGGAAAAAAGTATCGAGCAATGCCAGTTTAGGATCTAAACCCTGATCCAGCagagggggggaagtggggcagggggggt                                                                                                                                                                              >  1:100505/1‑227 (MQ=255)                                                                                                                      ggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCa                                                                                                                         >  1:188671/1‑250 (MQ=255)                                                                                                                              ttCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGCAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTgg                                                                                                                >  1:148049/1‑251 (MQ=255)                                                                                                                               tCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCg                                                                                                                                                         <  2:291223/209‑1 (MQ=255)                                                                                                                               tCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCg                                                                                                                                                         >  1:291223/1‑209 (MQ=255)                                                                                                                                      cccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCg                                                                                                        <  1:20226/251‑1 (MQ=255)                                                                                                                                      cccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCg                                                                                                        >  1:353419/1‑251 (MQ=255)                                                                                                                                                gCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGcc                                                                                              <  2:100505/251‑1 (MQ=255)                                                                                                                                                 cAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCg                                                                                             <  2:309618/251‑1 (MQ=255)                                                                                                                                                   ataatGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTc                                                                                           >  1:102623/1‑251 (MQ=255)                                                                                                                                                      atGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGACACCGCGCGATG‑CAACGTTGGTATCGCCGTGGATCGCCATCAGTAGGAAAAAGGGATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGCGGTAACGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGt                                                                                        >  2:188878/1‑251 (MQ=255)                                                                                                                                                          tttcgcCTCGGGCTGTAGCCAGCGGAGGCCACCCTCGCCCCCGAGCGATG‑ACACGTTGGTATCGCCGCTTATCGCCATCAGCAGGAACAAGGCATCGAGCAATGCACGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGAAAGACGTGATTGCTCACCAGCGCATAACCCGCTTCGGCTTCACCGCGGGAGCCGGTAAGGCAAAGCTGTTGGTAAAACCGTTGACCTGCCGTCAGTtgtg                                                                                    >  2:60895/4‑251 (MQ=255)                                                                                                                                                           tGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTtgtg                                                                                    <  2:83020/250‑1 (MQ=255)                                                                                                                                                            gcgcCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGaa                                                                                  <  1:402640/251‑1 (MQ=255)                                                                                                                                                             cgcCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAat                                                                                 >  1:128527/1‑251 (MQ=255)                                                                                                                                                                cTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAattat                                                                              >  1:254601/1‑251 (MQ=255)                                                                                                                                                                    cgcTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGCATCGCCGTTGATCGCCATCAGTAGGAACAAGGCATCGAGCAATGCCCGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCAGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGt                                                                          >  2:67305/1‑251 (MQ=255)                                                                                                                                                                    cgcTGTAGCCAGCGCAGGCCCACCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAATGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGt                                                                          <  2:240597/251‑1 (MQ=255)                                                                                                                                                                                 gcAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCg                                                                                                                                  <  2:49851/182‑1 (MQ=255)                                                                                                                                                                                 gcAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCg                                                                                                                                  >  1:49851/1‑182 (MQ=255)                                                                                                                                                                                              gCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCAc                                                >  1:306908/1‑251 (MQ=255)                                                                                                                                                                                              ctccccctctctcttcaacttttgTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGcc                                                  >  1:221352/24‑250 (MQ=255)                                                                                                                                                                                                    gcgcgATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGc                                             >  1:324637/1‑248 (MQ=255)                                                                                                                                                                                                              aaCGTTGGTATCGCCGTTGATCGCCATCAGTAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTgccg                                 >  1:26981/1‑251 (MQ=255)                                                                                                                                                                                                                                            tAGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCg   >  1:59706/1‑251 (MQ=255)                                                                                                                                                                                                                                             aGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCgg  <  2:181025/251‑1 (MQ=255)                                                                                                                                                                                                                                             aGGAACAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCgg  >  1:256993/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       CTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGG  >  NC_000913/646622‑647107

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