Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	4,348,705	A→G	61.0%	intergenic (‑70/+404)	PP_3820 ← / ← galU	Group II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,348,705	0	A	G	61.0%	1.2 / 44.8	41	intergenic (‑70/+404)	PP_3820/galU	Group II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (6/10);  new base G (9/16);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGT  >  AE015451/4348508‑4348927                                                                                                                                                                                                      |                                                                                                                                                                                                                               gccgccGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAg                                                                                                                                                                                          <  1:321111/236‑1 (MQ=17)    gccgGCACCTTTGTACTGTCAGGCATCACCCTCGGCTTCTGGCCGTGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                                                                                                <  1:237110/227‑1 (MQ=11)                   cTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                                                                                                <  2:451606/212‑1 (MQ=11)                         ggCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                                                                                                           >  2:511644/1‑227 (MQ=17)                              cACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                                                                                              <  1:417663/235‑1 (MQ=33)                                 ccTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAcc                                                                                                                                                          <  1:440075/236‑1 (MQ=18)                                   tCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCa                                                                                                                                                                                            >  2:503069/1‑200 (MQ=31)                                   tCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCCAGCGTACAGATCTCCCAGGGt                                                                                                                                                                       >  1:44536/1‑221 (MQ=17)                                   tCGGCTTCTGGCCGGGATTCGGGATAACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGc                                                                                                                                                                                             >  2:217464/1‑199 (MQ=14)                                              ccGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGc                                                                                                                                             >  1:292135/1‑236 (MQ=17)                                              ccGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGc                                                                                                                                             >  1:307992/1‑236 (MQ=17)                                              ccGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACCGATCTccc                                                                                                                                                                            >  2:352337/1‑205 (MQ=17)                                                         gATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCt                                                                                                                                  <  1:290054/236‑1 (MQ=33)                                                                         gtgGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCt                                                                                                                                  <  1:161169/220‑1 (MQ=18)                                                                                       cATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcg                                                                                                                                                                 >  1:431794/1‑175 (MQ=18)                                                                                          tACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTCCGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAg                                                                                                                                                                                                 <  2:213862/140‑1 (MQ=11)                                                                                           aCGGTTACCAACGAGACGGCGACGCCTCCTGGCGGCATACCCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                                                                                                <  1:79936/140‑1 (MQ=14)                                                                                           aCGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAcc                                                                                                                                                          <  1:380717/178‑1 (MQ=33)                                                                                                          aCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCAcc                                                                                 <  1:508925/236‑1 (MQ=33)                                                                                                            gtcGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgc                                                                                                                                                                <  2:309287/153‑1 (MQ=18)                                                                                                                 cGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCAccc                                                                           <  1:159571/235‑1 (MQ=18)                                                                                                                 cGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCAc                                                                             >  1:157254/1‑233 (MQ=17)                                                                                                                 cGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTCTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACCGCGTTCCg                                                                                                                  >  1:512922/1‑196 (MQ=18)                                                                                                                 cGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCCGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCg                                                                                             <  1:512350/217‑1 (MQ=17)                                                                                                                    cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCt                                                                       <  1:380650/236‑1 (MQ=18)                                                                                                                         ggcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACCGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGcc                                                                                                                                            <  2:408420/162‑1 (MQ=21)                                                                                                                         ggcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTgc                                                                   <  1:445165/235‑1 (MQ=17)                                                                                                                                      gTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTa                                                                                                                                 <  1:96959/160‑1 (MQ=33)                                                                                                                                       ttCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCt                                                     <  1:29009/235‑1 (MQ=17)                                                                                                                                           gCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTg                                                >  1:185874/1‑236 (MQ=18)                                                                                                                                               ttGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCg                                            <  1:477196/236‑1 (MQ=33)                                                                                                                                                      cgcgGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                                                                                                     <  2:30974/92‑1 (MQ=37)                                                                                                                                                           ttACCCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAAGTCGCGCCACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCa                                                                              <  2:485031/190‑1 (MQ=16)                                                                                                                                                            acaCCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGc                                                                                                                                                                                         <  2:58475/80‑1 (MQ=37)                                                                                                                                                                      ttACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGtt                                              <  1:318437/211‑1 (MQ=18)                                                                                                                                                                        actacGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGCATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTgg                                                                                               >  1:295754/1‑160 (MQ=18)                                                                                                                                                                             ggCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGcttc                 <  2:185870/233‑1 (MQ=17)                                                                                                                                                                                  cGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCAcccc                                                                          >  1:258077/1‑171 (MQ=18)                                                                                                                                                                                   ggCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGcttc                 >  1:265764/1‑227 (MQ=17)                                                                                                                                                                                     cTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGt                                           >  1:4404/1‑199 (MQ=33)                                                                                                                                                                                         cTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGt  <  1:503073/236‑1 (MQ=32)                                                                                                                                                                                                      |                                                                                                                                                                                                                               GCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGT  >  AE015451/4348508‑4348927

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