Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 710,595 G→T P18T (CCT→ACT)  fur ← DNA‑binding transcriptional dual regulator Fur

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,5950GT98.3% 173.7 / NA 59P18T (CCT→ACT) furDNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (27/31);  total (27/32)
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

GTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGA  >  NC_000913/710455‑710731
                                                                                                                                            |                                                                                                                                        
gtcgtcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTaa                                                                                                                                        <  1:971152/143‑1 (MQ=255)
ctcgtcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTaa                                                                                                                                        <  1:267048/142‑1 (MQ=255)
    tcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGt                                                                                                                                    <  1:383271/143‑1 (MQ=255)
          tGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACttt                                                                                                                              >  1:226001/1‑143 (MQ=255)
          tGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATCGTTGTCCGGCTCCTGACGAACTTACAGGCTTTTTAAACGAGTAAGCGTTACttt                                                                                                                              >  2:614444/1‑143 (MQ=255)
            gTTCAGTAGGCGATATGCCGTTGCCAGACCACTTTCTGCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCa                                                                                                                            <  2:830593/143‑1 (MQ=255)
                  gaCGCGATATGCCGGAGCCAGACCAAGTTCTTCACCCGTATCGATCAGACGTTTGTATAAATCTTCCGCACTGTCGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAg                                                                                                                      <  1:469917/142‑1 (MQ=255)
                    cgcgATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGct                                                                                                                    <  1:607557/143‑1 (MQ=255)
                     gcgATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGctt                                                                                                                   >  1:134020/1‑143 (MQ=255)
                          ataCCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGctttcttt                                                                                                              <  1:332131/143‑1 (MQ=255)
                           taCCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTa                                                                                                             <  2:134020/143‑1 (MQ=255)
                             ctgTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAgg                                                                                                           <  1:542530/141‑1 (MQ=255)
                              cGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTCCTTTCAGGCCAGctttctttcagg                                                                                                          >  2:209149/1‑139 (MQ=255)
                              cGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAggg                                                                                                          <  1:146659/143‑1 (MQ=255)
                                tAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACTAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCg                                                                                                        <  1:794766/143‑1 (MQ=255)
                                   ccAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGGTTTCGTTAGGGCGGTa                                                                                                     >  1:689891/1‑143 (MQ=255)
                                     aGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCAGTACTTTCAGGCCCGCTTTCTTTAGGGCCATAtt                                                                                                   >  2:33527/1‑143 (MQ=255)
                                      ggccAATTTCTTCCCCCAGAGCGGTCGGACGTTTGTATAAAGCTTNCGCCCTGACGTGATGGTTGGCCGGCTCGTNAAGATCTTCCAGGATTTTNAAGCGANTAAGCGNTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTg                                                                                                  <  2:606530/141‑1 (MQ=255)
                                            tttCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCa                                                                                            <  1:442192/143‑1 (MQ=255)
                                                  cACCCATATCTATCTTTCGTTTGTATATTTCTTCCGCACTGACGTGATGGGTGTCCGGCTCCTGAAGAACGTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATg                                                                                      <  1:433699/143‑1 (MQ=255)
                                                     ccATCTCGATCAGACGTTTGTATAAATCTTCCGCACCGACGTCATGGTTGGCCGGCTCCCGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCgg                                                                                   <  1:487029/143‑1 (MQ=255)
                                                        tatCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAAt                                                                                <  1:1013145/143‑1 (MQ=255)
                                                          tCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGCTTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCt                                                                              >  2:823925/1‑143 (MQ=255)
                                                             atcgcaCGTTTGCACACAGCTGCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTc                                                                           <  1:92261/138‑1 (MQ=255)
                                                              tCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTACACGAGTAAGCGTTACTTTAAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTTTcc                                                                          >  2:770007/1‑143 (MQ=255)
                                                                    gTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGAATTGTTATCAGTCATGCGGAATCTGTCCTGTTAc                                                                    >  1:1088970/1‑143 (MQ=255)
                                                                          ataAATCTTCCGCACTGACGTGATGGTTGTCCGGCTACTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGGTACTTTCAGGCCAGCTTTCTTCAGGGCGATATTGTTATCATTCATGCGGAATCTGTCCGGTTACTAAGct                                                              >  2:805021/1‑142 (MQ=255)
                                                                              aTCTTCCTCACTGACGTGATGGTTGCCCGGCTCCTGAAGACCTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTc                                                          <  2:202499/143‑1 (MQ=255)
                                                                                     gtattGACGTGATGGTTGCCCGGCTCCTGAAGTACTTCCCGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAt                                                   <  2:56058/139‑1 (MQ=255)
                                                                                      cACTGACGTGATGGTCGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAtt                                                  <  1:692401/143‑1 (MQ=255)
                                                                                          gACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTagaa                                              <  1:35591/143‑1 (MQ=255)
                                                                                             ttGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAgaagaa                                           >  1:736401/2‑143 (MQ=255)
                                                                                             gTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGATACTTTCAGGCCAGCTTTCTTTAGGGCCGAAAACACATCAGACAAACGCGGTCGGTCACGTTACAAAGCGGGTCAATTCATaagaagca                                           >  2:535300/1‑141 (MQ=255)
                                                                                               tttGGGTGTCCGGCTCTTGAAGAACTTCCAGGATTTGTAAACGGGTATGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTTTTGCTAAGCGGTTCACTTCATTAGAAGAAGt                                         <  2:34129/141‑1 (MQ=255)
                                                                                               gATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGt                                         <  1:1111224/143‑1 (MQ=255)
                                                                                                aTGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                        >  2:220201/1‑143 (MQ=255)
                                                                                                aTGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                        >  1:244784/1‑143 (MQ=255)
                                                                                                aTGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                        >  1:850737/1‑143 (MQ=255)
                                                                                                       tCCTGCTCCTGAAGAACTTCCAGGATTTTTAAACGATTAAGCGTTACTTTCAGACCAGCTCTCTTTAGGGCGGTATTGTGATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAt                                 <  1:252536/143‑1 (MQ=255)
                                                                                                        ccGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTCTTGGTATCAGTCACACCAAAATTTTCTTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAtt                                >  1:1056821/1‑143 (MQ=255)
                                                                                                        ccGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCACGCCAGCTCTCTTTAGGGCGGTACTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGCAta                                <  1:770007/143‑2 (MQ=255)
                                                                                                          ggCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTg                              <  1:135748/143‑1 (MQ=255)
                                                                                                             tCCTGAAGCACTTCCAGGATTTGTAAGCGAGTAAGCGTTACTTTCAGTCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctc                           <  2:853307/143‑1 (MQ=255)
                                                                                                             tCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctc                           <  1:488706/143‑1 (MQ=255)
                                                                                                             tCCTGAAGAACTTCCAGGATCTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctc                           <  1:823925/143‑1 (MQ=255)
                                                                                                               cTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCtt                         >  1:735111/1‑143 (MQ=255)
                                                                                                                tGAAGAACTTCCAGGATTTTTAACCGACTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGTGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTg                        <  2:980384/143‑1 (MQ=255)
                                                                                                                  aagaaCTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgag                      >  2:727001/1‑143 (MQ=255)
                                                                                                                        ttCCAGGATTTTTAAATGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATg                >  1:638126/1‑143 (MQ=255)
                                                                                                                          ccAGGATTTTTAAACGCGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGTGCGGAATTGTTATCAGTAATGCGATGTCTTACCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCg              >  2:335615/1‑143 (MQ=255)
                                                                                                                           cAGGATATTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCACGCGGAATCTGCCCTGCTACTACGCGGATCACTACAATAGAAGAAGTGACAGAATTTGCTCACGAGATAATGCGt             >  2:18927/1‑143 (MQ=255)
                                                                                                                            agtaTTTGTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTa            <  2:657190/140‑1 (MQ=255)
                                                                                                                            aGGATTTTTAAACTAGTAAGCGTTACTTTCATGCCAGCTTTCTTTAGGGCGTTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTa            <  2:658744/143‑1 (MQ=255)
                                                                                                                            aGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTCAGCGGTTCACTTCATTAGAAGAAGTGACAGCATTTGCTCTTGAGATAATGCGTa            >  1:631553/1‑143 (MQ=255)
                                                                                                                              gATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTACCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATc          >  2:552945/1‑143 (MQ=255)
                                                                                                                              gATTTTTAAACGAGTAAGCGTTACTTTCAGGCAAGCTTTCTGTAGGGCGGAATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTGGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATc          >  1:286708/1‑143 (MQ=255)
                                                                                                                                 ttttAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAATGACAGAATTTGCTCTTGAAATAATGCggagcagt       >  2:140418/1‑136 (MQ=255)
                                                                                                                                   ttAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTGTTGAGATAATGCGTATAATtat     >  2:642510/1‑143 (MQ=255)
                                                                                                                                      aaCGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTAATGGTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGa  >  2:905067/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                        
GTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGA  >  NC_000913/710455‑710731

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: