Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 124,241 C→T Q409* (CAG→TAG)  aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913124,2410CT94.8% 143.2 / NA 60Q409* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base T (19/36);  total (22/38)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTGGTCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCA  >  NC_000913/124101‑124380
                                                                                                                                            |                                                                                                                                           
gtggcgTTCACGATCCGAAGAAAATCTACTCTGCATTCAAGAAAGCGCAGGAATCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAg                                                                                                                                           <  2:207623/143‑1 (MQ=255)
  ggtggtAACGACCCGAATAAAATCTGCGCTGCACTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGt                                                                                                                                         <  2:400293/143‑1 (MQ=255)
    tgggCACGATCCGAAGAAAATCTACGCTGCGTGGCAGAAAGCGCAGGACACCAAACGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTa                                                                                                                                       <  2:279855/143‑1 (MQ=255)
       tCACGATCCGAAGAAAATCTACTCTGCATTCAAGAAAGCGCAAGAAACCAAAGGCAAAGCGACAGTAATGCCTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGa                                                                                                                                    <  2:857772/143‑1 (MQ=255)
       tCACGAGCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGa                                                                                                                                    <  2:796422/143‑1 (MQ=255)
        cACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCGAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGaa                                                                                                                                   >  1:590927/1‑143 (MQ=255)
         aggagccGAGGAAAGTCTACGCGGGAGTCAAGAAAGCGGAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCAGGAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGaaa                                                                                                                                  <  1:808013/138‑1 (MQ=255)
           gATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAGGGTAAAACCATCGCGCACCAGGTTCAGAAAAt                                                                                                                                >  1:351650/1‑143 (MQ=255)
           gATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAAt                                                                                                                                <  1:499674/143‑1 (MQ=255)
           gATCCGAAGAAAATCTACGCCGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAAt                                                                                                                                <  1:187638/143‑1 (MQ=255)
            aTCCGAAGAAAATCTACGCTGCATTCCAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATg                                                                                                                               >  1:1028690/1‑143 (MQ=255)
              gcaaataaaaTCTACGCTGCATACAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGaa                                                                                                                             <  1:1120796/137‑1 (MQ=255)
                  agaaAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATg                                                                                                                         <  1:697712/143‑1 (MQ=255)
                   gaaAATCTACGCTGCATGCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGATACGGCAGGGGCGACGCGGCTGAAGGTAAAAACATCGCGCAATAGGGAATGAAAGTGAACATgg                                                                                                                        >  2:177660/1‑143 (MQ=255)
                     aaaTCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGAc                                                                                                                      >  1:116724/1‑143 (MQ=255)
                       ctCTACTCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACgg                                                                                                                    <  2:1053045/142‑1 (MQ=255)
                         ctccGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTGATACCAATAAAGGTTACGGCATGGGCGCCGCCGCTGAAGGTCAAAACATCGCGAACCAGGTTACGACAATGAACATGGACgacg                                                                                                                  >  2:712392/4‑140 (MQ=255)
                          tACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGgtgt                                                                                                                 <  1:250749/143‑1 (MQ=255)
                                   ttCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCCTTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGCTAAGAAAATGAGCATGGAAGGTGTGCGTCatat                                                                                                        >  2:420118/1‑143 (MQ=255)
                                    tCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATc                                                                                                       >  1:841890/1‑143 (MQ=255)
                                      aaGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCcg                                                                                                     >  1:673535/1‑143 (MQ=255)
                                      aaGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAAAATAGTGCATTGTGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCcg                                                                                                     >  1:519709/1‑143 (MQ=255)
                                        gAAAGCGCAGGCAACCAAAGGCAAAGCGGCAGTAATCCTAGCTCATACCCTTAAAGGTTACGGCATGGGCGACTCGGCCGACGGTAACAACATCGAGGACTAGGTTGAGCAAATCAACATGGAGGGGGTGAGGCatatgagag                                                                                                   >  2:798267/1‑138 (MQ=255)
                                        gAAAGCGCAGCGAACCAAAGGCGAAGCGACAGTAATCCTTGCTCATACCCCTAAATGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCcgcg                                                                                                   <  2:1045426/143‑1 (MQ=255)
                                               cAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGttt                                                                                            <  2:895245/143‑1 (MQ=255)
                                                 ggAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCa                                                                                          <  2:18760/143‑1 (MQ=255)
                                                  gAACCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAATATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCaa                                                                                         <  1:603943/143‑1 (MQ=255)
                                                  gAAAGCAAAGGCAAAGCGACAGTCATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCaa                                                                                         <  1:800288/143‑1 (MQ=255)
                                                    aaCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATAAACATGGACGGTGTGCGTCATATCCGCGACCGTGTCCAtg                                                                                       >  2:656008/1‑143 (MQ=255)
                                                       cAAAGGCAAAGCGACGGGAATCCTGGCGCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGc                                                                                    <  1:875954/143‑1 (MQ=255)
                                                       cAAAGGCAAAGCGACAGTAATCCTTGATCATACCATTAAAGGTTACGGCATGGGCGACGCGGATGAAGGTAAAACGATTGCGCACTAGGTTAAGGAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTGAATGTGc                                                                                    >  1:830236/1‑143 (MQ=255)
                                                         aaGGCAAAGCGACAGTAATCCTCGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCg                                                                                  <  1:691681/143‑1 (MQ=255)
                                                              aaaGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTACGTCATAACCCCGACCGTTTCAAAGTACCAGTGTc                                                                             >  1:2867/1‑143 (MQ=255)
                                                                aGCGACAGTAATCCTTGCTCATACTATTTATGGTTACGGCATGTACAACGCGGCTGAAGGTAAAAACACCGAGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTg                                                                           <  1:1124009/143‑1 (MQ=255)
                                                                aGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTg                                                                           <  1:452419/143‑1 (MQ=255)
                                                                      aGTAATCCGTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAg                                                                     <  1:162147/143‑1 (MQ=255)
                                                                        tAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGat                                                                   <  2:590927/143‑1 (MQ=255)
                                                                          aTCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGatat                                                                 >  1:245494/1‑143 (MQ=255)
                                                                             cTTGCTCATACCATTAAAGGTTAGGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATTGa                                                              <  2:212998/143‑1 (MQ=255)
                                                                             cTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGa                                                              >  1:808858/1‑143 (MQ=255)
                                                                             cTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAACAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGa                                                              >  1:806557/1‑143 (MQ=255)
                                                                              ttgcgtaTACCATTTAAACTTACGGCATGTTCGACTCGGCTGATGGTAAAAACATCTCGCACTAGCTTAACAAAATTAACATGGACGTTGTGCGTCATATCCTCTACCGTTTCAATGTGCCGTTGTCTGATGCAGATATCGaa                                                             <  2:5686/137‑1 (MQ=255)
                                                                                       ccATGAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACGAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCg                                                    <  1:557551/143‑1 (MQ=255)
                                                                                        cATTAAAGGTTACGGCATGGGCGACGCGGCTGGAGGAAAAAACATCGCGCACTGGGGTAAAGGAATGAACAGGGACGGTGTGCGTCATACCCGCGAGCGTTTCAATGTGCCGGTGACTAATGCAGGGATGGAAAAAATGCAGt                                                   >  2:949424/1‑143 (MQ=255)
                                                                                                       cATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCgg                                    <  1:372895/143‑1 (MQ=255)
                                                                                                         tGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTGTCGATGTGACGGTGTCTGCTGCAGATATCGAAAAACTGCCGTACATGACCTTCCCGGaa                                  >  1:309140/1‑143 (MQ=255)
                                                                                                                 gcgGCTGAAGGAAAAAACATCGCACACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGa                          <  1:1060888/143‑1 (MQ=255)
                                                                                                                  cgGCTGAAGGTAAGAACATCGCGTACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGaa                         <  2:865554/143‑1 (MQ=255)
                                                                                                                  cgGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAACTGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGaa                         <  1:688466/143‑1 (MQ=255)
                                                                                                                     cTGAAGGTAAAACCATCGCGCACTAGGTTAGGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGAGGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAagag                      <  2:1091714/143‑1 (MQ=255)
                                                                                                                      taaaGCTAAAAACATCGCGCACGAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGc                     <  1:1027696/141‑1 (MQ=255)
                                                                                                                         aGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGTGGGTGTGTGTCAGATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATa                  <  1:1075880/143‑1 (MQ=255)
                                                                                                                              aaaaCATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTAt             <  1:959684/143‑1 (MQ=255)
                                                                                                                              aaaaCATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTAt             >  1:634196/1‑143 (MQ=255)
                                                                                                                               aaaCATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATc            >  1:118805/1‑143 (MQ=255)
                                                                                                                                aaCCTCGCGCAGTAGGCTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAACGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCt           <  1:316756/143‑1 (MQ=255)
                                                                                                                                  cATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTATGAAGAGCCTACCTATCTGc         >  2:318372/1‑143 (MQ=255)
                                                                                                                                   agcgcgCACTTGGTGAAGAGAATGAACAGGGGCGGGGCGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGGTATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCa        <  2:665106/141‑1 (MQ=255)
                                                                                                                                         cACTAGGTTGAGACAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAACAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCa  <  2:563908/143‑1 (MQ=255)
                                                                                                                                         cACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCa  <  2:430437/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GTGGTGGTCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCA  >  NC_000913/124101‑124380

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: