Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,478,324	0	T	G	54.7%	52.1 / 88.2	53	G33G (GGT→GGG)	ynbB	putative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base T (3/21);  new base G (29/0);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 6.37e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-10

GGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTT  >  NC_000913/1478174‑1478474                                                                                                                                                        |                                                                                                                                                       ggACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGt                                                                                                                                                        <  2:658151/151‑1 (MQ=255)   aCTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGa                                                                                                                                                      <  1:1130987/151‑1 (MQ=255)        gATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTa                                                                                                                                                 <  1:533289/151‑1 (MQ=255)                      ggAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTT‑CCGGAGGGGAAAGGGGGGGAGGGTACATTACGTAATCg                                                                                                                                    >  2:222728/1‑150 (MQ=255)                      ggAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGt                                                                                                                                   <  1:146108/151‑1 (MQ=255)                       gAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTa                                                                                                                                  <  2:910465/151‑1 (MQ=255)                          tgtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑TCGGGGAGAAAAGGGGGGGGGGGTTTATTATGTATTTGttact                                                                                                                               >  1:21869/1‑147 (MQ=255)                                      gAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATc                                                                                                                    <  1:736141/150‑1 (MQ=255)                                      gAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATc                                                                                                                    <  2:343408/150‑1 (MQ=255)                                       aGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTT‑TCGGGGAGAAAAGGGGGGGAGGGTACCTTACGTATTCGTATCTGGTggggggt                                                                                                                     >  1:901913/1‑145 (MQ=255)                                        gcgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTc                                                                                                                 <  1:457832/151‑1 (MQ=255)                                         cgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTC‑CTGGGGGGGAAAGAGGGGGAGGGTATATTTCGTATTCGTTCCTGGGGGGGGAttg                                                                                                                   >  1:953115/1‑146 (MQ=255)                                               aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGtt                                                                                                          <  2:681279/151‑1 (MQ=255)                                                  ataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCTGAGCGAAAAAGGGGGGAGGTTTAATTTCGTATTTGTTCCTTGTGGGGGGTCGTCACtttttttt                                                                                                       >  1:653159/1‑151 (MQ=255)                                                  ataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGtttt                                                                                                        <  2:817077/150‑1 (MQ=255)                                                           gAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCCGGGGGGAAAAGGGGGGGGGGTTAC‑TTACGTATTTCTAACTGGTGGGGGGTCGTCATTTGttttt                                                                                                       >  1:191822/1‑142 (MQ=255)                                                               aaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGGGAGGGTAAATTACGTATTTGTACCTGGTGGGGGATCGTCATTTTTTTTTTTATGGGGGttt                                                                                           >  2:803879/1‑150 (MQ=255)                                                                 tctcTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGGGGGGAAAGGGGGGGAGGGTTCATTTTGTATTTGCTCCTTGTGGGGGCACGTCATTTTTTTTTTAATGGGGtttttc                                                                                         >  1:392376/1‑148 (MQ=255)                                                                    ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCa                                                                                         <  2:511414/147‑1 (MQ=255)                                                                         aaCACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTT‑CCGGGGAGAAAAGGGGGGGGGGGTACATTTCTTATTTGTACCCGGGGGGGCATCGTCATTTGTTTTTTTATGGGGATTTTAGGGCCACGc                                                                                >  1:415390/1‑151 (MQ=255)                                                                            acTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑TCCGAGAGAAAAAGGGGGGAGGGTTCATTACGTATTTCTACCTGGTGGGGCGTCGTCCTTTGTTTTTTTATGGTGGTTTTAGGTCCACGCTgg                                                                             >  1:433537/1‑151 (MQ=255)                                                                               tttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAt                                                                           <  1:244634/150‑1 (MQ=255)                                                                                ttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑ACGGGGGGAAAAGGGGGGGGGGGTACATTACGCATTTGTAACTGGGGGGGGATCGGCATTTGTTTTTTAAATGTGGTTTTAggggccccctgg                                                                             >  2:595524/1‑137 (MQ=255)                                                                                   gcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGGGGGAAAAGGGGGGGAGGGTACATTACGGATTCGTACCTGGTGGGGGATCGTCATTTTTTTTTCAATTGTGGTTTTAGGGCCACGCCGGGtgacg                                                                       >  1:844503/1‑150 (MQ=255)                                                                                     gttgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttg                                                                    <  2:957278/151‑1 (MQ=255)                                                                                       tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGGGAGAAAAAGGGGGGGGGGTACATTTCGTATTTGTACCTGGGGGGGGATCGGCCTTTTTTTTTTAATGGTGATTTTAGGGCCCCGCTGGGtgacg                                                                       >  2:434733/1‑146 (MQ=255)                                                                                           aaTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGttt                                                              <  1:240716/151‑1 (MQ=255)                                                                                              tttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGGGGGAAAAAGGGGGGAGGGTACATTACGTATTTGTACCTGGGGGGGGATCGTCATTTGTTTTTTAATGGTGGTTTCAGGTCCACGCTGGATGACGTTGACGtttttt                                                           >  2:871935/1‑151 (MQ=255)                                                                                                 aGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑TCGGGGAGAAAAGGGGGGGAGGGTACATTAAGTATTTGTACCTGGGGGGGCATCGGGCTTTGTTTTTTTATGGGGATTTTAGGTGCACGCCGGGTGAGGTTTTCGgttttt                                                           >  1:1001064/1‑148 (MQ=255)                                                                                                 aGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTgc                                                         <  1:226678/150‑1 (MQ=255)                                                                                                    aaCGCTGATAAACCGTTTTCTGCTATGGCGGTTT‑CCGGAGGGAAAAAGGGGGGGGGGTACATTACGTATTCGTACCTGGGGGGGCATCGTCCTTTGTTTTTTAATGGGGATTTTAGGGCCCCCCTGGGTGACGTTGACGGTTTTTgtg                                                        >  2:945620/1‑146 (MQ=255)                                                                                                     aCGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑ACGGGGAGAAAAGGGGGGGAGGGTTAATTTCGTATTCGTACCTGGGGGGGGATCGTCATTTGTTTTTCAATGGTGATTTTCGGTCCACGCTGGGTGGCGTTGACGtttttt                                                           >  1:885126/1‑144 (MQ=255)                                                                                                     aCGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑ACGGAGGGAAAAGGGGGGGAGGGTACATATCGTATTTGGACCTGGGGGGGGATTGTTATTTGTTTTTTTATGGTGATTTTAGGGCCACGCCGGATGAGGTTTACGTTTTTTgcgc                                                       >  1:281371/1‑148 (MQ=255)                                                                                                      cGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGGGAGAAAAGGGGGGGAGGGTATATTACGTATTTGTTGCTTGGGGGGGGTTGTTATTTTTTTTTTAAATGTGATTTCAGGGCCACGCGGGGTGACGTTGTCGtttttt                                                           >  2:316714/1‑143 (MQ=255)                                                                                                      cGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGGGAGGGTACATTACGTATTCGTACCTGGGGGGGGGTCGTCATTTGTTTTTTAATGGGGGTTTTAGGGCCACGCTGGGTGAAGTTGACGgttttt                                                           >  2:925282/1‑143 (MQ=255)                                                                                                      cGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCCGGGGGTAAAGGGGGGGAGGGTACATTACTTATTTGTACCTGGTGGGGGCTCGTTATTTGTTTTTTTATTGTGATTTTAGGTCCCCGCTGGGTGACGTTGGCGTTTTTTTCGCTg                                                     >  1:791081/1‑149 (MQ=255)                                                                                                      cGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑ACGGGGGGAAAAGGGGGGGGGGGTAAATTACGGATTTGTACCCGGGGGGGGATAGTCATTTGTTTTTTAATGGTGATTTCTAGGCCACGCCGGATGGCGTTTAAGTTTTTTTTGCTggt                                                   >  1:530765/1‑149 (MQ=255)                                                                                                          gATAAACCGTTTTCTGCTATGGCGGTTA‑CCCGGGAGAAAAAGGGGGGAGGGTACATTAAGTATTTGTACCCGGGGGGGGATCGGCATTTGTTTTTTTATGGTGATTTCAGGGCCACGCCGGATGACGGTGACGGTTTTTTCGCTGATAAGt                                               >  1:231796/1‑151 (MQ=255)                                                                                                             aaaCCGTTTTCTGCTATGGCGGTTT‑CCGGGGAGAAAAAGGGGGGAGGTTACATTATGTATTCGCACCTGGGGGGGGATCGTCCTTTGTTTTTTAATGGGGGTTTTAGGGCCACGCTGGAAGACGTTTACGTTTTTTTAGCTGTTAAGgttt                                            >  2:330957/1‑151 (MQ=255)                                                                                                              aaCCGTTTTCTGCTATGGCGGTTA‑GCGGAGGGAGAAGGGGGGGGTGGTACCTTTTTTAGTCGTAGCTGGTGGGGGACCTGCAATTTTTTTTCAAAGGGGAGTTTATGGCCACGGTTGGCGAGCGTGAGTGTTTTTTCGCTGATAAGttttt                                           >  1:356765/1‑150 (MQ=255)                                                                                                               aCCGTTTTCTGCTATGGCGGTTA‑CCCGGGAGAAAAAGGGGGGAGGTTACATTACGTATTCGTACCTGGGGGGGGATCGTCATTTGTTTTTCAATGGGGATTTCAGGTCCCCGCTGGATGACGTTGACGTTTTTTTCGCTGATAAGTTTTCt                                          >  2:258524/1‑151 (MQ=255)                                                                                                                ccGTTTTCTGCTATGGCGGTTT‑CCGGAGGGGAAACGGGGGGGGGTTTCATTACGTTTTTGTTCCTGGGGGGGGCGCGTCATTTTTTTTTTTATGGGGATTTTTGGGCCACGCTGGGTGACGGTGTCGTTTTTTTTGCTGATTAGtttttt                                          >  1:63650/1‑148 (MQ=255)                                                                                                                   ttttCTGCTATGGCGGTTA‑ACGGGGGGGAAAAGGGGGGGGGTTTCATTACGTATTCGTACCTGGGGGGGGATCGGCCTTTGTTTTTTAATGGGGATTTTAGGGTCCCGCCGGGTGGCGTTGGCGGTTTTTTTGCTGGTAAGGTTTTTgg                                        >  1:515448/1‑149 (MQ=255)                                                                                                                    tttCTGCTATGGCGGTTA‑CCGGGGAGAAAAGGGGGGGGGGGTACATTAAGTATTTGTACCTGGGGGGGCATCGTCATTTGTTTTTTTATGGGGATTTTAGGTCCACGCTGGGTGAAGGTGGCGTTTTTTTTGCTGGTTAGTTTTTTTGCa                                      >  1:250022/1‑150 (MQ=255)                                                                                                                      tCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCAtt                                    <  1:675259/150‑1 (MQ=255)                                                                                                                         gCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGaaa                                <  2:791374/151‑1 (MQ=255)                                                                                                                            aTGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAg                               <  1:335293/149‑1 (MQ=255)                                                                                                                                cGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATAtt                         <  1:403773/151‑1 (MQ=255)                                                                                                                                 ggTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTg                        <  1:862821/151‑1 (MQ=255)                                                                                                                                   ttA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGt                       >  2:31165/1‑150 (MQ=255)                                                                                                                                     a‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACg                    >  1:314651/1‑151 (MQ=255)                                                                                                                                          gagagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCtt                 >  1:767372/1‑150 (MQ=255)                                                                                                                                                   gggggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTAc       <  1:125885/151‑1 (MQ=255)                                                                                                                                                      ggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAtt    <  1:45188/151‑1 (MQ=255)                                                                                                                                                        tGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAtttt  >  1:82176/1‑151 (MQ=255)                                                                                                                                                        tGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAtttt  >  1:9985/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       GGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTT  >  NC_000913/1478174‑1478474

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