Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	387,624	G→A	E218E (GAG→GAA)	tauC →	taurine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	387,624	0	G	A	97.1%	108.4 / ‑4.5	34	E218E (GAG→GAA)	tauC	taurine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (17/16);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GGTATTGGTGAAACCTCGAAGATCTTACTGATCTATTTAGCGATTTTTGCACCGGTGGCGATGTCGGCGCTGGCGGGGGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTCTGCGCGCGTTACAGCGCCGCCTGACGCCCTGGCATGGAGAAGTACAATGAGTG‑AACGTCTGAGCATTACCCCGCTGGGGCCGTATATCGGCGCACAAATTTCGGGTGCCGACCTG  >  NC_000913/387385‑387863                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 ggTATTGGTGAAACCTCGAAGAGCTTACTGATCTATTTAGCGATTGTGGCGCCGGTGGCGATGGCGGCGCGGGCGGGGGTGAAAAGCGTGCAGCAGGTGCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGACGGAAATCCTCACCGGATTGGGGATTGGTCTGGGGGGGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCcgc                                                                                                                                                                                                                                       <  2:211604/251‑1 (MQ=255)        gTGAAACCTCGAAGATCTTACTGATCTATTTAGCGATTTTTGCACCGGTGGCGATGTCGGCGCTGGCGGGGGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGAcgcgc                                                                                                                                                                                                                                <  1:155860/251‑1 (MQ=255)           aaaCCTCGAAGATCTTACTGATCTATTTAGCGATTTTTGCACCGGTGGCGATGTCGGCGCTGGCGGGGGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGGGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGt                                                                                                                                                                                                                             >  1:3932/1‑251 (MQ=255)                              gATCTATTTAGCGATTTTTGCACCGGTGGCGATGTCGGCGCTGGCGGGGGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTtca                                                                                                                                                                                                            <  1:128242/249‑1 (MQ=255)                                   attgaGCGAGTTTTGCACCGGTGGCGATGTCGGCGCTGGCGGGGGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGACCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTGGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCGGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTtcagtcag                                                                                                                                                                                                       <  1:66162/245‑1 (MQ=255)                                                tGCACCGGTGGCGATGTCGGCGCTGGCGGGGGTGAAAAGCGTGCAGCAGGGTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATCTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGGGCGGGTTAGGATTTATTGTTCAGTCAGCGGGTGAATTtctac                                                                                                                                                                                        >  2:369954/1‑249 (MQ=255)                                                            gATGTCGGCGCTGGCGGGGGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGttt                                                                                                                                                                                                                           <  1:331008/204‑1 (MQ=255)                                                            gATGTCGGCGCTGGCGGGGGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGGGAATTTCTCGCAACTGACgtgg                                                                                                                                                                            >  1:7015/1‑251 (MQ=255)                                                            gATGTCGGCGCTGGCGGGGGGGAAAAGCGTGCAGCAGGTTCGCATTCGCGCCGCCCCGTCGCTGGGGGCCAGCCGTGCGCTGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCCGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGGTGGTCTACGCTGGTGGCGGCGGAAATGGTTGCCGGGACGCGCGGttt                                                                                                                                                                                                                           >  2:331008/1‑204 (MQ=255)                                                             aTGTGGGCGCTGGCGGGTGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGTGGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGGGGAACTGATTGCCGCGACGCGCGTTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACgtggt                                                                                                                                                                           <  2:357015/251‑1 (MQ=255)                                                                  ggCGCTGGCGGGGGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGGCAGCGGGGGAATTTCTCGCAACTGACGTGGTGCTgg                                                                                                                                                                      >  2:389652/1‑251 (MQ=255)                                                                      cTGGCGGGGGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCggg                                                                                                                                                                  >  1:90326/1‑251 (MQ=255)                                                                            ggggTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATcgc                                                                                                                                                            <  2:315067/251‑1 (MQ=255)                                                                              ggTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGGTCGCgg                                                                                                                                                          >  1:289035/1‑251 (MQ=255)                                                                                       cGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCgg                                                                                                                                                 >  1:171839/1‑250 (MQ=255)                                                                                                     gCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAAt                                                                                                                                                                                              <  1:302694/192‑1 (MQ=255)                                                                                                     gCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAAt                                                                                                                                                                                              >  2:302694/1‑192 (MQ=255)                                                                                                             gccgccCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAAc                                                                                                                            <  1:203488/250‑1 (MQ=255)                                                                                                              ccgccCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTg                                                                                                                          <  1:216876/251‑1 (MQ=255)                                                                                                                  ccAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTc                                                                                                                      >  1:260536/1‑251 (MQ=255)                                                                                                                                      gTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTCTGCGCGCGTTACAGcgccgc                                                                                                  >  1:205704/1‑251 (MQ=255)                                                                                                                                         cgcAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGAGCGCCGAGATTGCGATTATCGCCTTTCTTTTAAAAATGGGTCTTCGCGCGTTTCAACGCCACCTg                                                                                               >  1:372556/1‑251 (MQ=255)                                                                                                                                            aGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTCTGCGCGCGTTACAGCGCCGCCTGACg                                                                                            <  1:63642/251‑1 (MQ=255)                                                                                                                                                    tgGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGGGATTGCGATTATCGCCTTTCTTTTTGAACTTGGTTTGCGCGCGGTACaacccccccatgcgccccgcg                                                                                     >  1:187564/1‑230 (MQ=255)                                                                                                                                                     ggTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTCTGCGCGCGTTACAGCGCCGCCTGACGCCCTGGCa                                                                                    <  2:205704/250‑1 (MQ=255)                                                                                                                                                      gTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGGTCGCGGTGATTGCGGGTATCGCCTTTCTTTTTGAACTGGGTTCGCGCGCGTTACAGCGCCGCCTGACGCCCTGGCATg                                                                                  >  2:363717/1‑251 (MQ=255)                                                                                                                                                                ttGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTCTGCGCGCGTTACAGCGCCGCCTGACGCCCTGGCATGGAGAAGTAc                                                                         <  1:236895/250‑1 (MQ=255)                                                                                                                                                                           gcTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTCTGCGCGCGTTACAGCGCCGCCTGACGCCCTGGCATGGAGAAGTACAATGAGTG‑AACg                                                            <  2:111538/251‑1 (MQ=255)                                                                                                                                                                                        cTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCCTTCTTTTTGAACTGGGTCTTCGCGCGTTTCAGCGCCGCCTGCCGCCCTGGCATGGAGAAGTACAATGGGTG‑AACGCCTGGGCAtttcc                                               >  2:140635/1‑248 (MQ=255)                                                                                                                                                                                                  tACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTCTGCGCGCGTTACAGCGCCGCCTGACGCCCTGGCATGGAGAAGTACAATGAGTG‑AACGTCTGAGCATTACCCCGCTgggg                                      <  1:140635/250‑1 (MQ=255)                                                                                                                                                                                                      tATTGGTCTGGTGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTCGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCCACTGAAGTGGTGCTGGCGGGGATCGCGGGGGTTTCGGATAACTCCTTTTTTTTAGCACTGGGGTTGGGAGCGTTTCACCGCCGCCTGGCGTCCTGGCATGGAGAAGTACAAGGAGTGCAACGttcggagttaccccccctgggccgtt                                  >  2:367861/1‑226 (MQ=255)                                                                                                                                                                                                       aTTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTCTGCGCGCGTTACAGCGCCGCCTGACGCCCTGGCATGGAGAAGTACAATGAGTG‑AACGTCTGAGCATTACCCCGCTGGGGCCGta                                 <  2:289035/250‑1 (MQ=255)                                                                                                                                                                                                                                 gCTGGTGGCGGCGGAACTGCTTGCCGCGACGCGAGGTTTAGGATTTATGGTTCAGTCAGCGGGGGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACAGGGGCTTCGCGCGTTACAGCGCCGCCTGACGCCCAGGCATGGAGAAGTACAAAGAGAG‑AACGTCTGGACCTTAACCCGCGGGGGCCGGATATTGGAGCACCAATTTTGGGTGCCgc      >  1:17824/1‑250 (MQ=255)                                                                                                                                                                                                                                     gtggCGGCGGAACTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTCTGCCCGCGTTACAGCGCCGCCTGACGCCCTGGCATGGGGAAGTACAATGAGTG‑AAAGTCTGAGCATTACCCCGCCGGGGCCGTATATCGGCGCACCAATTTCGGGGGCCGACCTg  >  1:310481/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 GGTATTGGTGAAACCTCGAAGATCTTACTGATCTATTTAGCGATTTTTGCACCGGTGGCGATGTCGGCGCTGGCGGGGGTGAAAAGCGTGCAGCAGGTTCGCATTCGTGCCGCCCAGTCGCTGGGTGCCAGCCGTGCGCAGGTGCTGTGGTTTGTCATTTTGCCCGGTGCGCTGCCGGAAATCCTCACCGGATTACGTATTGGTCTGGGGGTGGGCTGGTCTACGCTGGTGGCGGCGGAGCTGATTGCCGCGACGCGCGGTTTAGGATTTATGGTTCAGTCAGCGGGTGAATTTCTCGCAACTGACGTGGTGCTGGCGGGGATCGCGGTGATTGCGATTATCGCCTTTCTTTTAGAACTGGGTCTGCGCGCGTTACAGCGCCGCCTGACGCCCTGGCATGGAGAAGTACAATGAGTG‑AACGTCTGAGCATTACCCCGCTGGGGCCGTATATCGGCGCACAAATTTCGGGTGCCGACCTG  >  NC_000913/387385‑387863

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