Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 715,435 C→T V8V (GTC→GTT ybfP → lipoprotein YbfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913715,4350CT96.7% 192.0 / ‑4.0 60V8V (GTC→GTTybfPlipoprotein YbfP
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (34/24);  total (35/25)
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

GGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTGCGGTAAAT  >  NC_000913/715299‑715570
                                                                                                                                        |                                                                                                                                       
ggTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGgttgtc                                                                                                                                      >  1:281523/1‑140 (MQ=255)
 gTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCa                                                                                                                                     >  1:1271217/1‑140 (MQ=255)
  tGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCAt                                                                                                                                    >  1:21262/1‑140 (MQ=255)
  tGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCAt                                                                                                                                    <  2:21262/140‑1 (MQ=255)
   gCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATa                                                                                                                                   <  2:467606/140‑1 (MQ=255)
      atatatTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAgtg                                                                                                                                <  2:543634/140‑1 (MQ=255)
           ttGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACt                                                                                                                           >  1:145622/1‑140 (MQ=255)
            tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTa                                                                                                                          >  2:1135348/1‑140 (MQ=255)
            tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTa                                                                                                                          >  1:698092/1‑140 (MQ=255)
                 ttAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTAc                                                                                                                     >  2:661711/1‑140 (MQ=255)
                   aaGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCg                                                                                                                   >  2:1394181/1‑140 (MQ=255)
                     gCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGcc                                                                                                                 >  2:286756/1‑140 (MQ=255)
                          tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACAtt                                                                                                            >  1:1314238/1‑140 (MQ=255)
                                   cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGAc                                                                                                   >  1:1051906/1‑140 (MQ=255)
                                   cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGAc                                                                                                   >  2:210542/1‑140 (MQ=255)
                                   cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGAc                                                                                                   >  1:1273732/1‑140 (MQ=255)
                                     gATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTg                                                                                                 >  2:819248/1‑140 (MQ=255)
                                     gATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTg                                                                                                 >  1:482858/1‑140 (MQ=255)
                                      aTAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGc                                                                                                >  1:773971/1‑140 (MQ=255)
                                       tAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCa                                                                                               <  1:661711/140‑1 (MQ=255)
                                            ggtccACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAACCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATgcgc                                                                                          <  2:935291/136‑1 (MQ=255)
                                            cGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATgcgc                                                                                          >  2:862645/1‑140 (MQ=255)
                                               tCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTActg                                                                                                        >  1:47225/1‑123 (MQ=255)
                                               tCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCa                                                                                       >  1:661479/1‑140 (MQ=255)
                                               tCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCa                                                                                       >  2:1191513/1‑140 (MQ=255)
                                               tCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCa                                                                                       >  1:696099/1‑140 (MQ=255)
                                               tCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAACCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTActg                                                                                                        <  2:47225/123‑1 (MQ=255)
                                                cACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGccc                                                                                        <  1:1184682/138‑1 (MQ=255)
                                                cACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGccc                                                                                        >  2:1184682/1‑138 (MQ=255)
                                                cACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCaa                                                                                      >  2:50487/1‑140 (MQ=255)
                                                cACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCaa                                                                                      >  1:824075/1‑140 (MQ=255)
                                                         tCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATgcgc                                                                                          <  2:671701/127‑1 (MQ=255)
                                                         tCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATgcgc                                                                                          >  1:671701/1‑127 (MQ=255)
                                                         tCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACaa                                                                             <  2:773971/140‑1 (MQ=255)
                                                          ccGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAat                                                                            <  1:1394181/140‑1 (MQ=255)
                                                          ccGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAACCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAat                                                                            <  1:199501/140‑1 (MQ=255)
                                                                  aacaacGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTa                                                                    >  1:38633/1‑140 (MQ=255)
                                                                      acGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTg                                                                >  2:1279512/1‑140 (MQ=255)
                                                                       cgaccccTTTTTCAACAACCTGATGAGCTAAAGGAAACCCAAAGAAAGACAAAGGGGGGTCTGTTTTTCATGGTTTCACTATTTCCCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTgg                                                               <  2:696099/134‑1 (MQ=255)
                                                                       cGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTgg                                                               >  1:688988/1‑140 (MQ=255)
                                                                          acacTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGa                                                            <  2:1051906/140‑1 (MQ=255)
                                                                          acacTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGa                                                            <  2:145622/140‑1 (MQ=255)
                                                                                   caacaaACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACt                                                   <  1:1135348/140‑1 (MQ=255)
                                                                                   caacaaACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACt                                                   <  1:1191513/140‑1 (MQ=255)
                                                                                      caaACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGAc                                                >  2:1256856/1‑140 (MQ=255)
                                                                                          cTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCt                                            <  1:1256856/140‑1 (MQ=255)
                                                                                           tGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTa                                           <  1:819248/140‑1 (MQ=255)
                                                                                            gAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTaa                                          >  1:209844/1‑140 (MQ=255)
                                                                                                   tCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATaaa                                   <  1:862645/140‑1 (MQ=255)
                                                                                                       ggAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTa                                           <  1:1189868/128‑1 (MQ=255)
                                                                                                       ggAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTa                                           >  2:1189868/1‑128 (MQ=255)
                                                                                                       ggAAATCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCa                               <  2:482858/140‑1 (MQ=255)
                                                                                                       ggAAATCCCAATGAAGACAAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACAGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCTTAAAGATAAAACCa                               >  2:858310/1‑140 (MQ=255)
                                                                                                            tCCCAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGa                          >  1:1102008/1‑140 (MQ=255)
                                                                                                               cAATGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGAtca                       <  1:736950/140‑3 (MQ=255)
                                                                                                               cAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTa                       <  1:286756/140‑1 (MQ=255)
                                                                                                                 aTGAAGACCAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTGAc                     >  2:1161593/1‑140 (MQ=255)
                                                                                                                         cAACAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCa             <  2:661479/140‑1 (MQ=255)
                                                                                                                           aCAGGAGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCAtt           <  2:633543/140‑1 (MQ=255)
                                                                                                                                aGTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTGCGGt      <  2:209402/140‑1 (MQ=255)
                                                                                                                                 gTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTGCg        <  2:1168111/137‑1 (MQ=255)
                                                                                                                                 gTCTGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTGCg        >  1:1168111/1‑137 (MQ=255)
                                                                                                                                    tGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTGCGGTAAat  >  1:1309348/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTGCGGTAAAT  >  NC_000913/715299‑715570

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