Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,086,996 Δ13,531 bp yiiGrhaT yiiG, frvR, frvX, frvB, frvA, rhaM, rhaD, rhaA, rhaB, rhaS, rhaR, rhaT

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4086996 4100526 13531 38 [0] [0] 40 yiiG–rhaT yiiG,frvR,frvX,frvB,frvA,rhaM,rhaD,rhaA,rhaB,rhaS,rhaR,rhaT

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40869950 (0.000)38 (0.560) 35/492 0.6 100% intergenic (+146/‑7) fdhD/yiiG sulfurtransferase for molybdenum cofactor sulfuration/DUF3829 domain‑containing lipoprotein YiiG
?NC_000913 4100527 = 0 (0.000)intergenic (‑2/‑283) rhaT/sodA rhamnose/lyxose:H(+) symporter/superoxide dismutase (Mn)

GCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCAACCTGCCGCACCAGCTCCTGAAGCGAAACCTACTGAAGCTCCCGTTGCAAAAGCAGAAGCTAAACCTGAAACACCTGCGCAACCGGTGGTCGATGAACAAGCGGTTTTCGACGAAAAAATGGATG  >  NC_000913/4086747‑4087239
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
gCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTTACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCtt                                                                                                                                                                                                                                                    <  2:1634‑M1/251‑3 (MQ=255)
gCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCtt                                                                                                                                                                                                                                                    <  2:303359‑M1/251‑3 (MQ=255)
   tGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttct                                                                                                                                                                                                                                                  <  2:84564‑M1/250‑5 (MQ=255)
    gCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATGAGGCttctta                                                                                                                                                                                                                                                <  1:311214‑M1/251‑7 (MQ=255)
       gTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcc                                                                                                                                                                                                                                             <  1:12314‑M1/251‑10 (MQ=255)
       gTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttat                                                                                                                                                                                                                                               <  1:391984‑M1/249‑8 (MQ=255)
                   gAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACTGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatgtttc                                                                                                                                                                                                                                 <  2:271411‑M1/251‑22 (MQ=255)
                                 gACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatct                                                                                                                                                                                                                   <  1:265075‑M1/251‑36 (MQ=255)
                                                                    cGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcattgttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcttgctgtcgatt                                                                                                                                                                                <  2:198329‑M1/251‑71 (MQ=255)
                                                                                aTCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAACCAAAATGCGTGTAGGGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTCTCTCTATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgactggatgagtgtaaaaatcgtgctgtcgattaacctttcgcct                                                                                                                                                                    >  1:77031‑M1/1‑169 (MQ=255)
                                                                                             gAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttg                                                                                                                                                       >  1:429905‑M1/1‑156 (MQ=255)
                                                                                                       ttaaAATAAGTCCTGCGAAATATATTAACTTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttcct                                                                                                                                             >  2:102014‑M1/1‑146 (MQ=255)
                                                                                                              aaGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatc                                                                                                                                      >  1:82387‑M1/1‑139 (MQ=255)
                                                                                                                     gCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcat                                                                                                                               >  1:362178‑M1/1‑132 (MQ=255)
                                                                                                                            atatTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtg                                                                                                                         <  1:313308‑M1/250‑126 (MQ=255)
                                                                                                                              atTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggt                                                                                                                      <  1:285342‑M1/251‑129 (MQ=255)
                                                                                                                                         ggACTTATTCCTTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGCAGTGGCGTTTTCTTTTTTCATTTATTAAGTTCACTCTATTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcc                                                                                                           >  2:430677‑M1/1‑112 (MQ=255)
                                                                                                                                            cTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacat                                                                                                        >  1:281749‑M1/1‑109 (MQ=255)
                                                                                                                                              tATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTTGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcagcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacat                                                                                                        <  2:509920‑M1/249‑143 (MQ=255)
                                                                                                                                                 tCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgg                                                                                                   >  2:113556‑M1/1‑104 (MQ=255)
                                                                                                                                                              tttattattttattTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggc                                                                                                                                 <  1:350383‑M1/208‑118 (MQ=255)
                                                                                                                                                                 attattttattTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcctttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaa                                                                                    >  1:390908‑M1/1‑88 (MQ=255)
                                                                                                                                                                          tttatAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgcgatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgcagccgttgtcgttttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccataccct                                                                          >  2:7721‑M1/1‑79 (MQ=255)
                                                                                                                                                                             ataAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttac                                                                       <  1:113556‑M1/251‑176 (MQ=255)
                                                                                                                                                                                 aCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaa                                                                    >  2:135962‑M1/1‑72 (MQ=255)
                                                                                                                                                                                    aaaTGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagt                                                                >  2:442245‑M1/1‑69 (MQ=255)
                                                                                                                                                                                        gCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacgg                                                            >  1:409623‑M1/1‑65 (MQ=255)
                                                                                                                                                                                          gtgtAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggc                                                           >  1:164526‑M1/1‑63 (MQ=255)
                                                                                                                                                                                          gagtagTGGCGTTTTCTTTTTTCCTTTATTATGTTCACTCTATTATCTATATATTTATAGTGCttcttatccttattatttttcgtcgcgtcacatctccaacgagatgagagtaataatcgtgctgtagattaacctttcgcccgttgcccccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaacaatacccttacgaaaagtacggct                                                          >  2:491177‑M1/3‑63 (MQ=255)
                                                                                                                                                                                                tGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgat                                                     >  1:184411‑M1/1‑57 (MQ=255)
                                                                                                                                                                                                  gCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataat                                                  >  1:174846‑M1/1‑55 (MQ=255)
                                                                                                                                                                                                                      ttatGTTCACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaat                              >  1:109249‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                                             cACTCTCTTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaattaactat                       >  1:342367‑M1/1‑28 (MQ=255)
                                                                                                                                                                                                                                tctctTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaattaactataat                    >  1:9643‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                                                tctcgtctctatataAATATAATGCttcatatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacattttgcctgttgccgccgttgtcgatttattggcaatcaccgtattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacctgtacgaaaactacggcattgataatcattatcaatatcatttaattaactataat                    >  2:394880‑M1/8‑25 (MQ=255)
                                                                                                                                                                                                                                 ctctTATCTATATATTTATAAGGCttcttatcctcatcatttttcgtcgcgtcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaattaactataatg                   >  1:293665‑M1/1‑24 (MQ=255)
                                                                                                                                                                                                                                            atatTTATAAGGTttcttatcctcatcatttttcgtcgcttcacatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaattaactataatgaaccaactgc         <  2:23336‑M1/250‑238 (MQ=255)
                                                                                                                                                                                                                                                  ataAGGCttcgtatcctcatcatttttcggcgcgtcaaatctccgacgagatgagtgtaaaaatcgtgctgtcgattaacctttcgcctgttgccgccgttgtcgatttactggcaatcacggcattaagtgggtgatttgcttcacatctcgggcattttcctgcaaaaccatacccttacgaaaagtacggcattgataatcattttcaatatcatttaattaactataatgaaccaactgcttacgcg  <  2:120569‑M1/251‑245 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
GCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCAACCTGCCGCACCAGCTCCTGAAGCGAAACCTACTGAAGCTCCCGTTGCAAAAGCAGAAGCTAAACCTGAAACACCTGCGCAACCGGTGGTCGATGAACAAGCGGTTTTCGACGAAAAAATGGATG  >  NC_000913/4086747‑4087239

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