Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 568,195 C→T pseudogene (161/213 nt) renD → DLP12 prophage; protein RenD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913568,1950CT98.4% 208.7 / NA 63pseudogene (161/213 nt)renDDLP12 prophage; protein RenD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/40);  minor base G (1/0);  total (23/40)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAA  >  NC_000913/568066‑568319
                                                                                                                                 |                                                                                                                            
cgcgGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGAt                                                                                                                           <  2:400326/133‑1 (MQ=255)
           aTCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTAt                                                                                                                <  1:55230/133‑1 (MQ=255)
            tCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTAtt                                                                                                               <  2:33149/133‑1 (MQ=255)
             ccTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTg                                                                                                              <  2:335326/133‑1 (MQ=255)
              cTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGa                                                                                                             <  2:51807/133‑1 (MQ=255)
              cTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGa                                                                                                             <  1:340244/133‑1 (MQ=255)
               tGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGGGAGGAAACGGGTATTGAg                                                                                                            >  1:174014/1‑133 (MQ=255)
                ggTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAgg                                                                                                           <  1:94421/133‑1 (MQ=255)
                ggTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAgg                                                                                                           <  1:41521/133‑1 (MQ=255)
                ggTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAgg                                                                                                           <  1:415410/133‑1 (MQ=255)
                     ttGATGGTAAGGTCTGGCTAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATATGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGtata                                                                                                      <  1:189789/133‑1 (MQ=255)
                       gATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGtatata                                                                                                    >  1:281754/1‑133 (MQ=255)
                                   tGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCaa                                                                                        >  1:482347/1‑133 (MQ=255)
                                    ggCGAACGGTGTATTTCCGGTTTGCTACCATAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAAt                                                                                       <  1:173410/133‑1 (MQ=255)
                                    ggCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAAt                                                                                       >  2:316200/1‑133 (MQ=255)
                                      cGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATgg                                                                                     <  2:94251/133‑1 (MQ=255)
                                          cGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAAt                                                                                       <  2:284598/127‑1 (MQ=255)
                                          cGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAAt                                                                                       >  1:284598/1‑127 (MQ=255)
                                           ggTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAgca                                                                                <  2:473174/133‑1 (MQ=255)
                                                 ttACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAatga                                                                          <  1:316200/133‑1 (MQ=255)
                                                    ccGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAAtgatga                                                                       <  1:369517/133‑1 (MQ=255)
                                                      ggTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGaaa                                                                     <  1:401134/133‑1 (MQ=255)
                                                      ggTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGaaa                                                                     <  1:310385/133‑1 (MQ=255)
                                                       gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGGGGGTATATGAAAGAACGTCAATGGGGGCACAATGATGAAAc                                                                    >  2:135192/1‑133 (MQ=255)
                                                       gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGGGGGTATATCAAATAACTTCAATGGGTGCACAATGATGAAAc                                                                    >  2:175992/1‑133 (MQ=255)
                                                       gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATGAAAAAACAACAATGGGAGCACCTTGATGAAAc                                                                    >  2:251409/1‑133 (MQ=255)
                                                       gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                    >  2:315446/1‑133 (MQ=255)
                                                       gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                    >  1:29398/1‑133 (MQ=255)
                                                       gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                    >  2:72593/1‑133 (MQ=255)
                                                       gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                    >  1:494600/1‑133 (MQ=255)
                                                       gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                    >  1:256854/1‑133 (MQ=255)
                                                       gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAAc                                                                    >  2:403231/1‑133 (MQ=255)
                                                         tCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGTGTATATAAGAGAGCATGAATGGGAGGACAAGGATGAAACAg                                                                  >  2:132127/1‑133 (MQ=255)
                                                         tCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAg                                                                  >  2:349290/1‑133 (MQ=255)
                                                         tCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAg                                                                  >  1:53417/1‑133 (MQ=255)
                                                                         aTGGGAAGGAAAGGTGAGCATGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCaaa                                                  <  1:409492/133‑1 (MQ=255)
                                                                           gggAAGGAAACGTGATCACGAATCGGATTTTTGAGGAGTGTAGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACt                                                <  2:195634/133‑1 (MQ=255)
                                                                            ggaaggaaAGGTGAGTACGAATCTGATTTTCAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTg                                               <  1:169195/133‑1 (MQ=255)
                                                                            ggaaggaaAGGTGAGTACGAATCTCATTATTATGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTg                                               <  2:174014/133‑1 (MQ=255)
                                                                            ggaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTg                                               <  1:286165/133‑1 (MQ=255)
                                                                             gaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTgt                                              >  2:472814/1‑133 (MQ=255)
                                                                                 gaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtagta                                          <  2:281754/133‑1 (MQ=255)
                                                                                 gaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtagta                                          <  1:317995/133‑1 (MQ=255)
                                                                                          gCACGAATTTGATTTTTAAGGACTGTCGCCGGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTGCAATtctct                                 <  1:251409/133‑1 (MQ=255)
                                                                                          gCACGAATCTGATTTCTACGGACTGTCGCCTGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtctct                                 <  1:135192/133‑1 (MQ=255)
                                                                                           cACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtc                                    >  1:387115/1‑129 (MQ=255)
                                                                                           cACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtc                                    <  2:387115/129‑1 (MQ=255)
                                                                                           cACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTACAAACTGTAGTACAATtctctc                                >  1:204937/1‑133 (MQ=255)
                                                                                                    gATTTTTAAGCAGTGCCGCCAGAGTGCCTTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGaa                       <  2:138035/133‑1 (MQ=255)
                                                                                                     aTTTTTAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAAc                      <  2:315410/133‑1 (MQ=255)
                                                                                                        tttAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAgg                   <  2:29398/133‑1 (MQ=255)
                                                                                                        tttAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAgg                   <  1:72593/133‑1 (MQ=255)
                                                                                                         ttAAGGAGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGa                  <  1:315446/133‑1 (MQ=255)
                                                                                                               aGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAat            >  1:449733/1‑133 (MQ=255)
                                                                                                               aGTGTCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAat            >  1:505759/1‑133 (MQ=255)
                                                                                                                 tgtCGCCAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAatat          >  2:405634/1‑133 (MQ=255)
                                                                                                                       cAGAGTGCCGTGATGCAACGGGCATTGAGGCTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg    <  1:132127/133‑1 (MQ=255)
                                                                                                                       cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg    <  2:98320/133‑1 (MQ=255)
                                                                                                                       cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg    <  2:482347/133‑1 (MQ=255)
                                                                                                                       cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg    <  2:275363/133‑1 (MQ=255)
                                                                                                                       cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg    <  1:87884/133‑1 (MQ=255)
                                                                                                                       cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg    <  1:472814/133‑1 (MQ=255)
                                                                                                                       cAGAGTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTACAGTTTGAACAGGAAAGAATATGCTATg    <  2:445163/133‑1 (MQ=255)
                                                                                                                         gagTGCCGTGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAGCACAATGATGAAACAGGTGAGTTGAGTTCAAACTGGAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGat  <  2:53417/133‑2 (MQ=255)
                                                                                                                                 |                                                                                                                            
CGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAA  >  NC_000913/568066‑568319

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