Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,481,795	G→A	Q128* (CAG→TAG)	emrY ←	tripartite efflux pump membrane subunit EmrY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,481,795	0	G	A	100.0%	135.5 / NA	37	Q128* (CAG→TAG)	emrY	tripartite efflux pump membrane subunit EmrY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/4);  total (33/4)

CCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGAAAGCC  >  NC_000913/2481551‑2482041                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       ccTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGGTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATAGCGGAGCTATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAg                                                                                                                                                                                                                                                  <  2:285644/251‑1 (MQ=255)    ggTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTg                                                                                                                                                                                                                                                <  1:745163/250‑1 (MQ=255)             aTTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACgg                                                                                                                                                                                                                                      >  1:520049/1‑251 (MQ=255)              ttttGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGc                                                                                                                                                                                                                                     >  2:130625/1‑251 (MQ=255)                         gAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTa                                                                                                                                                                                                                           >  1:139356/1‑250 (MQ=255)                            gTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAAccc                                                                                                                                                                                                                       >  2:217298/1‑251 (MQ=255)                                      ctctTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGCTAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGACCCGCCATTAACCCCTGAACGAcc                                                                                                                                                                                                             >  2:668045/1‑250 (MQ=255)                                                         tAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACACCTGAACGACTCTAAAAAATATCAGCACAt                                                                                                                                                                                          >  1:145875/1‑251 (MQ=255)                                                          aaGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTAATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGAATCTAAAAAATATCAGCACATc                                                                                                                                                                                         >  1:143812/1‑251 (MQ=255)                                                                aaTGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAAc                                                                                                                                                                                                                 >  1:111634/1‑221 (MQ=255)                                                                            aTAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGACCATCAAGATTGGTCGATAGGCt                                                                                                                                                                       >  2:636427/1‑251 (MQ=255)                                                                                ccccaaaaggaacatttaaaaaaaaTACCCACCCCCAACAAAAATTATCAAAAAAATACCCCCCCAAAACCCGCCCACATATCGGGGCGATAATCACGGTCATTGACAATAATGCCAGAAAAAATTTTTTTTTTTCTGGCGAATA‑TTCCTTTAATAAAAAACTCTAGGACAGTGAATTTGACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTacac                                                                                                                                                                   <  1:247447/231‑1 (MQ=255)                                                                                    cATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                                                                                                                >  1:343744/1‑250 (MQ=255)                                                                                      tAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAAt                                                                                                                                                             >  1:164079/1‑251 (MQ=255)                                                                                             aTTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCAAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGAACATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGCTAGGCTACAAATTAATGAAGAc                                                                                                                                                       >  1:700172/1‑250 (MQ=255)                                                                                                  ttAAAAATATCCAACCCCCGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACACCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGa                                                                                                                                                  >  2:632195/1‑250 (MQ=255)                                                                                                   tAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaa                                                                                                                                                 >  2:689154/1‑250 (MQ=255)                                                                                                     aaaaTATCCAACCCCAGCTAACGTTATCACCAATATAACCGCACAATATCGGCCCACATATCGGAGCGATAATCACGATCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAACTATCAGCACATCACGATTGGTCGTTAGGCTACACATCAATGAAGACCGCGaaaa                                                                                                                                               >  2:556286/1‑250 (MQ=255)                                                                                                         attCCAACCCCAGATAAATATATCACAAAAAAACCCGCCCAATATCGGCCCACATATGGGGGCGATAATCACGGTCATTAACAATACTGCAAGAAAAAATGTTCTTTTTTCTGGCGGATAATTCCTA‑AATAATAAACTCTATAAAAGTGAAATTACCGCCCCCCCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATAAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAg                                                                                                                                          <  2:247928/249‑1 (MQ=255)                                                                                                               aCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTg                                                                                                                                    >  1:467634/1‑251 (MQ=255)                                                                                                                             tATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAAt                                                                                                                                                             >  1:589491/1‑212 (MQ=255)                                                                                                                               tCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGCTTGGTCGATAGGCTACACATTAATGAAGACCGCGAAAAAAAAGTGACTGAAAGTAAAaataatat                                                                                                                    >  2:65426/1‑249 (MQ=255)                                                                                                                                         aaCCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCt                                                                                                          >  2:484710/1‑251 (MQ=255)                                                                                                                                                  atatCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtgt                                                                                                 >  1:39381/1‑251 (MQ=255)                                                                                                                                                                 tCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAcc                                                                                          >  1:2381/1‑243 (MQ=255)                                                                                                                                                                   ggAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAAc                                                                                >  1:283240/1‑251 (MQ=255)                                                                                                                                                                            aTCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGc                                                                       >  1:413984/1‑251 (MQ=255)                                                                                                                                                                                            ctAATGGCAGCGCAACTGCTCTTTTTTCTGGCGGATAATTACTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAAACCCTGAACGACTCGAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGCCTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTac                                                        >  2:235403/2‑250 (MQ=255)                                                                                                                                                                                               aTGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACAcc                                                     >  2:715606/1‑250 (MQ=255)                                                                                                                                                                                                 gCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCaa                                                   >  1:560510/1‑250 (MQ=255)                                                                                                                                                                                                          aaTGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCAttt                                                             >  1:370939/1‑231 (MQ=255)                                                                                                                                                                                                             gTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAAc                                      >  1:413195/1‑251 (MQ=255)                                                                                                                                                                                                                                    ccTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAACAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATg               >  1:444725/1‑251 (MQ=255)                                                                                                                                                                                                                                    ccTT‑AATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGtt                   >  2:298587/1‑247 (MQ=255)                                                                                                                                                                                                                                               aaaCTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGc                                                                 >  1:285253/1‑191 (MQ=255)                                                                                                                                                                                                                                               aaaCTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGaa      >  1:111813/1‑250 (MQ=255)                                                                                                                                                                                                                                                  ctctATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGAAAGcc  >  1:495485/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       CCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGAAAGCC  >  NC_000913/2481551‑2482041

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