Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 257,908 Δ776 bp [crl] [crl]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 257908–258675 258683 9–776 63 [1] [1] 64 [crl] [crl]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2579071 (0.020)60 (1.210) 53/268 0.0 99.2% intergenic (+8/‑769) crl/crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator,regulator, Surface structures, transcriptional regulator of cryptic csgA gene for curli surface fibers/pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator,regulator, Surface structures, transcriptional regulator of cryptic csgA gene for curli surface fibers
?NC_000913 258684 = 0 (0.000)pseudogene (9/331 nt) crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator,regulator, Surface structures, transcriptional regulator of cryptic csgA gene for curli surface fibers

TTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCC  >  NC_000913/257768‑258047
                                                                                                                                           |                                                                                                                                            
ttAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                                                                                                           <  2:863815‑M1/143‑4 (MQ=255)
 tAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGCTTGATCAAAACATTTACCGCACTAGGCCCGTATCTCCGTGAAGCTAagtg                                                                                                                                          >  1:15997‑M1/1‑139 (MQ=255)
 tAAAAATTTGCTAAATTTAGCCAATTTGGTAGAAAAGTTGAATAACAACAGGAGATAGCTGTGAAGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAATGACTTTAAAGATTAAGGCCCGTATATTCGTGAaggtaagtg                                                                                                                                          <  1:17692/143‑2 (MQ=255)
   aaaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGAAGATTGATCAAAAAATTTACGGCCCAAGGCCCGAGTATAGGAGAAGGAAagtgca                                                                                                                                        >  1:186190‑M1/1‑137 (MQ=255)
   aaaaTTTGCTAAATTTTGCCAATTTGGAAAAACAGTTGCATCACAACAGGAGATAGCAATGCCGTTACCGCGTGGACACCCGAAGAGCAGATTGCTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGCAGGTAagtgca                                                                                                                                        >  2:436029‑M1/1‑137 (MQ=255)
       tttGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcacaga                                                                                                                                    >  1:4507‑M1/1‑133 (MQ=255)
           cTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATCCGTGAAGGTAagtgcaaagataat                                                                                                                                >  2:390559‑M1/1‑129 (MQ=255)
             aaaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagatactcc                                                                                                                              >  1:367884‑M1/1‑127 (MQ=255)
              aaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcga                                                                                                                             <  1:257653‑M1/143‑18 (MQ=255)
               aTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagggcaaagataatcgac                                                                                                                            >  2:190687‑M1/1‑125 (MQ=255)
                   tGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctt                                                                                                                        >  2:913040‑M1/1‑121 (MQ=255)
                   tGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATGTACCGCACTAGGCCCGTATATTGGTGAAGGTAagtgcaaagataatcgattctt                                                                                                                        >  2:500637‑M1/1‑121 (MQ=255)
                    gCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGCATATTCGTGAAGGTAagtgcaaagataatcgattcttt                                                                                                                       >  1:777056‑M1/1‑120 (MQ=255)
                     ccAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATGAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGGAagtgcaaagataatcgattctttt                                                                                                                      >  1:738647‑M1/1‑119 (MQ=255)
                      cAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttttt                                                                                                                     <  2:169816‑M1/143‑26 (MQ=255)
                      cAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttttt                                                                                                                     <  1:505699‑M1/143‑26 (MQ=255)
                       aaTTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                                                                                                                    <  2:482683‑M1/143‑27 (MQ=255)
                         tctggTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcga                                                                                                                  <  2:562959‑M1/141‑29 (MQ=255)
                             gTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctctgacgaatgt                                                                                                              >  1:354136‑M1/1‑111 (MQ=255)
                              tAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACCCCCGAAGCGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGATGGTAagtgcaaagataatcgattctttttctattgtc                                                                                                             >  2:586824‑M1/1‑110 (MQ=255)
                               aaaaCAGTTGCATCACAACAGGAGATAGCAATGATGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtct                                                                                                            <  2:562540‑M1/143‑35 (MQ=255)
                                aaaCAGTTGCATCACAACAGGAGATAGCACTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGGCTATTCGTGAGGGGAagagcaaagataatcgattccttttcgattgtctg                                                                                                           >  2:412768‑M1/1‑108 (MQ=255)
                                aaaCAGTGGCAGCACAACAGGAGATAGCAAGGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctg                                                                                                           <  1:256157‑M1/143‑36 (MQ=255)
                                 aaCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                                                                                                          <  2:518354‑M1/143‑37 (MQ=255)
                                      ttGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataagcgattctttttcgattgtctggctgta                                                                                                     >  1:177067‑M1/1‑102 (MQ=255)
                                       tGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtat                                                                                                    <  2:836276‑M1/143‑43 (MQ=255)
                                           tCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgt                                                                                                >  2:337060‑M1/1‑97 (MQ=255)
                                            cacaACAGGATATAGCAATGACGTTACCGAGTGTAAACCAGTTGGGGGGATTAAACAAAAAATTTACCGCACTAGGCGCGAATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtc                                                                                               <  2:419543‑M1/143‑48 (MQ=255)
                                            cacaACAGGAGATAGCAATGACGTTACCGAGTGGACACACGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGATAagtgcaaagatcatcgaatcttcttcgattgtctggctgtatgcgtc                                                                                               >  1:149306‑M1/1‑96 (MQ=255)
                                                 ccggagATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataaccgattctttttcgattgtctggctgtatgcgtcaacgt                                                                                          <  1:746679‑M1/141‑53 (MQ=255)
                                                       aTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataagcgattcttttgcgattggctggctgtatgcgtcaacgtgaaacc                                                                                    <  1:869619‑M1/143‑59 (MQ=255)
                                                           cAATGACGTTACCGAGTGGACACCTGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggca                                                                                <  2:152646‑M1/143‑63 (MQ=255)
                                                                 cGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaa                                                                          <  2:844648‑M1/143‑69 (MQ=255)
                                                                  gTTACCGATTGGACACCCGCAGAGGAGATGGATTAAATAATTTACCGCACTAGGCCCGTATATTCGTGAAGGCAagtgcgctgatactcgcttctttttcgagtgtctggctgtatgcgtcaacgtcaaaccggcaccggaag                                                                         <  1:188039‑M1/143‑70 (MQ=255)
                                                                      ccGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcg                                                                     >  1:877806‑M1/1‑70 (MQ=255)
                                                                         aGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTGTATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtga                                                                  >  1:563151‑M1/1‑67 (MQ=255)
                                                                          gcggacACCCGAAGGGCAGATTGGTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaa                                                                 <  1:665185‑M1/141‑78 (MQ=255)
                                                                             gACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattc                                                              <  2:527541‑M1/143‑81 (MQ=255)
                                                                                acgcgaaGAGCAGATTGAGCAAAAAAGTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgg                                                           <  1:859180‑M1/140‑84 (MQ=255)
                                                                                 cccGAAGAGCAGATTGATCAAAAAAGTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggg                                                          <  2:367884‑M1/143‑85 (MQ=255)
                                                                                   cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaggatattagattctttttcgattgtctggctgtatgcgtcaacgtgtaaccggcaccggaaggtcgtgaattcgggggc                                                        >  2:252810‑M1/1‑57 (MQ=255)
                                                                                      agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctgg                                                     <  2:503335‑M1/143‑90 (MQ=255)
                                                                                             aTTGATCAAAAAATTTACCGCACTAGGCCCGCATATTCGTGAAGGTAagtgcacagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggcagtgcgtgaactctggggccggtagcggg                                              >  1:82281‑M1/1‑47 (MQ=255)
                                                                                                 aTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgcattctggggctggtggatggagct                                          >  2:635906‑M1/1‑43 (MQ=255)
                                                                                                 aTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagct                                          >  2:850945‑M1/1‑43 (MQ=255)
                                                                                                 aTCAAAAAATATACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagct                                          <  2:8344‑M1/143‑101 (MQ=255)
                                                                                                  tCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgggaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctg                                         >  2:82197‑M1/1‑42 (MQ=255)
                                                                                                    aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgggaattctggggctggtggatggggctggg                                       >  1:543587‑M1/1‑40 (MQ=255)
                                                                                                     aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagcgggag                                      >  1:805191‑M1/1‑39 (MQ=255)
                                                                                                               cGCACTAGGCCCTTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaat                            <  1:70186‑M1/143‑115 (MQ=255)
                                                                                                               cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaggcgcgggaat                            >  2:797935‑M1/1‑29 (MQ=255)
                                                                                                               cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggcggctggagctggaagcgccggaat                            >  2:598871‑M1/1‑29 (MQ=255)
                                                                                                                gCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatc                           >  2:886130‑M1/1‑28 (MQ=255)
                                                                                                                  aCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatctattctttttctattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggattgagctggaagcgcaggaatccc                         <  2:594076‑M1/143‑118 (MQ=255)
                                                                                                                    tcggccCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaattgcgtgaattctggggctggtggatggagctggaagcgcaggaatccagt                       <  1:873376‑M1/141‑120 (MQ=255)
                                                                                                                    tAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtggcgccgccaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgt                       <  2:250391‑M1/143‑120 (MQ=255)
                                                                                                                     gggcCCGTATATTCGTGAAGGTAagtgcaaagatcatcgattctccttcgattgtctggctgtatgcggcaacgtgagacctgcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtt                      <  1:412768‑M1/142‑121 (MQ=255)
                                                                                                                         ccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttac                  <  2:777056‑M1/143‑125 (MQ=255)
                                                                                                                               atTCGTGCAGGGAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacag            <  1:287656‑M1/143‑131 (MQ=255)
                                                                                                                                   gTGAAGGTAagtgcaaatttaatcgattttttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttac        <  2:683480‑M1/143‑135 (MQ=255)
                                                                                                                                   gTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttac        >  2:152965‑M1/1‑9 (MQ=255)
                                                                                                                                     gAAGGGAagtggaaagatgagcgatgctgtgtcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttacca      <  2:39745‑M1/143‑137 (MQ=255)
                                                                                                                                         gTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtgttacctacagttaccagtgt  >  2:371877‑M1/1‑3 (MQ=255)
                                                                                                                                           |                                                                                                                                            
TTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCC  >  NC_000913/257768‑258047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: