Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,353,673	G→A	L282L (CTC→CTT)	sapC ←	putrescine ABC exporter membrane protein SapC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,353,673	0	G	A	100.0%	128.9 / NA	37	L282L (CTC→CTT)	sapC	putrescine ABC exporter membrane protein SapC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/14);  total (23/14)

AGTAACACGCCAGTTATCTTTATTCACCCCACAAATCGCTTTCGCAATCAAACTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCATCACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGC  >  NC_000913/1353437‑1353900                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                    agTAACACGCCAGTTATCTTTATTCACCCCACAAATCGCTTTCGCAATCAAACTCTTGACGGAAACGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACcaacaa                                                                                                                                                                                                                       <  2:250/251‑1 (MQ=255)     acacTCCAGTTAGCTTGATTCACCCCACAAATCGCTTTCGCAATCAAACTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTAGTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACg                                                                                                                                                                                                                   <  1:275972/251‑1 (MQ=255)             gTTATCTTTATTCACCCCACAAATCGCTTTCGCAATCAAACTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCAt                                                                                                                                                                                                           >  1:414805/1‑251 (MQ=255)                          aCCCCACAAATCGCTTTCGCAATCAAACTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCt                                                                                                                                                                                              >  1:224133/1‑251 (MQ=255)                           ccccACAAATCGCTTTCGCAATCAAACTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTg                                                                                                                                                                                             <  1:301311/251‑1 (MQ=255)                                 aaaTCGCTTTCGCAATCAAACTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGgcagc                                                                                                                                                                                        >  2:190284/1‑250 (MQ=255)                                  aaTCGCTTTCGCAATCAAACTCTTGCAGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCCTTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCAt                                                                                                                                                                                      >  1:127440/1‑251 (MQ=255)                                     cGCTTTCGCAATCAAACTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAGTCATAATTGCCGCACCTGGCAGCATGAc                                                                                                                                                                                   >  1:244661/1‑251 (MQ=255)                                         ttCGCAATCAACCTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAg                                                                                                                                                                                 <  1:103788/249‑1 (MQ=255)                                                tCAAACTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGATACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACgggg                                                                                                                                                                        >  1:380358/1‑251 (MQ=255)                                                   aaCTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCg                                                                                                                                                                      >  1:310520/1‑250 (MQ=255)                                                                              caaGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCGGGCAGCATGACAGGCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGa                                                                                                                                          <  2:151237/251‑1 (MQ=255)                                                                                aGACCGCGGATTTCACATTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGc                                                                                                                                        >  2:249017/1‑251 (MQ=255)                                                                                 aaCCGCTGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCa                                                                                                                                       <  2:321112/250‑1 (MQ=255)                                                                                  aCCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCAt                                                                                                                                      <  1:93988/251‑1 (MQ=255)                                                                                   ccGCGGATTTCACCTTCTGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCCCCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACTACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTGCCAg                                                                                                                                                     <  2:416202/235‑1 (MQ=255)                                                                                   ccGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAACTTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGGGCCATAAATCCGTTCCAg                                                                                                                                                     >  1:416202/1‑235 (MQ=255)                                                                                       ggATTTCACCTTCGGTTAAGGACATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATACTCGCACGACGAACGCCATCGCCAAGAAGGTTCACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTc                                                                                                                                 >  2:74735/1‑251 (MQ=255)                                                                                                          ggTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGATTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCAACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAAAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGAATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCATCACCGAGCATCGCTCCCCATTCAGGCGAgggg                                                                                                                >  1:139817/1‑249 (MQ=255)                                                                                                                            gACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCg                                                                                              >  2:415638/1‑249 (MQ=255)                                                                                                                             aCGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATAGCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCgaga                                                                                           <  2:94539/251‑1 (MQ=255)                                                                                                                               ggCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGaa                                                                                          >  2:48687/1‑250 (MQ=255)                                                                                                                                     aaNCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACAGAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCa                                                                                    <  2:380358/250‑1 (MQ=255)                                                                                                                                      aCCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCNAgc                                                                                  >  2:379299/1‑251 (MQ=255)                                                                                                                                          aCTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCgcgg                                                                              >  1:67699/1‑251 (MQ=255)                                                                                                                                                  ccGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCg                                                                      <  1:415638/251‑1 (MQ=255)                                                                                                                                                    ggTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCg                                                                      >  2:234049/1‑249 (MQ=255)                                                                                                                                                       tttAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAAt                                                                 >  1:302846/1‑251 (MQ=255)                                                                                                                                                             ttCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCAt                                                           >  1:92725/1‑251 (MQ=255)                                                                                                                                                                          ttACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCgg                                              >  1:416897/1‑251 (MQ=255)                                                                                                                                                                                  atCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGAtctc                                      >  1:320580/1‑251 (MQ=255)                                                                                                                                                                                       gTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGCACGCCATAGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGTCGACATACATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGCGACAGCAAAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTg                                  >  1:172007/1‑250 (MQ=255)                                                                                                                                                                                                  attCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCcggc                       <  1:197241/250‑1 (MQ=255)                                                                                                                                                                                                   ttCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCcggcgg                     >  2:216879/1‑251 (MQ=255)                                                                                                                                                                                                                  aatCGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATAGCGAGAATTGCCATCGACAGTGCGCGGGGGATCTCGGTGACCAGCCCGGCGGTGAAGTTTGGCATca      <  2:127440/251‑1 (MQ=255)                                                                                                                                                                                                                     cGCACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAGTCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATcaca    <  2:244661/250‑1 (MQ=255)                                                                                                                                                                                                                       cACGACGAACGCCATCGCCAAGAAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGc  >  1:314638/1‑250 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                    AGTAACACGCCAGTTATCTTTATTCACCCCACAAATCGCTTTCGCAATCAAACTCTTGCCGGAACCGGATTCACCAACAAGACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCATCACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGC  >  NC_000913/1353437‑1353900

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