Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,869,891 C→G 18.2% D79H (GAC→CAC)  U712_19940 ← Deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,869,8910CG18.2% 117.9 / 6.0 61D79H (GAC→CAC) U712_19940Deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (19/27);  new base G (8/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 3.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTGACAAGCGCTTTTCCGGCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTG  >  CP006881/3869755‑3870022
                                                                                                                                         |                                                                                                                                    
aCTTTGACAAGCGCTTTTCCGGCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCtt                                                                                                                                    <  7:359420/139‑1 (MQ=255)
aCTTTGACAAGCGCGTTTCCGGCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCtt                                                                                                                                    <  1:27073/139‑1 (MQ=255)
    tGACAAGCGCTTTTCCGGCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCGGTTTTAAATGGGGTGCttttt                                                                                                                                 >  2:146420/1‑138 (MQ=255)
     gACAAGCGCTTTTCCGGCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCGTTTTAAATGGCGTGTTTTGtt                                                                                                                               >  6:321919/1‑139 (MQ=255)
          gcgcTTTTCCGGCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGaa                                                                                                                          <  1:279872/139‑1 (MQ=255)
                   cGGCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGt                                                                                                                 <  1:272214/139‑1 (MQ=255)
                     gCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGttt                                                                                                               <  1:146420/139‑1 (MQ=255)
                     gCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGttt                                                                                                               <  1:211771/139‑1 (MQ=255)
                     gCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGttt                                                                                                               <  1:353085/139‑1 (MQ=255)
                     gCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGttt                                                                                                               <  5:287808/139‑1 (MQ=255)
                         cAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGGAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGt                                                                                                           <  5:117500/139‑1 (MQ=255)
                         cAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGt                                                                                                           <  1:108492/139‑1 (MQ=255)
                         cAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGt                                                                                                           <  7:45333/139‑1 (MQ=255)
                         cAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGt                                                                                                           <  6:286544/139‑1 (MQ=255)
                         cAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGt                                                                                                           <  4:68047/139‑1 (MQ=255)
                              tCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCATGTCTTTTAAAGCGGCAATATTAATGACCATATCCCCTTCAGTGGCGCCTTTTTAAATGGCGTCTTTTTTTTTGCAGGCGTTTGTTTACGTTGtta                                                                                                      >  8:167292/1‑137 (MQ=255)
                               ccACTACACCGCGGATATCAGCTTCCACCCCATCGTCTTCCTTGTCTGTTAAAGCGACAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGATATGTCGTCGTTTGTTTCGCACGCGTTTGTTTCAGTTGTagt                                                                                                     >  6:143702/1‑137 (MQ=255)
                                        cgcgGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCGTTTGAAATGGCGTCTTTTGTTTCGAAAGCGTTTGTTTGAGTTTTTTTGGCACCGAg                                                                                            >  4:197146/1‑139 (MQ=255)
                                          cgGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATTGCGTGTGTTTTTTCGCACGCGTTTGTTTGATGTGTATTTGGACCGAGCg                                                                                          >  3:247794/1‑139 (MQ=255)
                                          cgGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTGTTTTGTTTCGAACGCGTTTGTTTCAGGTGTATTGGCACCGAGCg                                                                                          >  8:60933/1‑139 (MQ=255)
                                          cgGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCg                                                                                          <  7:36213/139‑1 (MQ=255)
                                          cgGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGAGGCAATAATAATGACCATATCCACTTCAGTGGCGCCTATTTAAATGGCGTGGTTTTTTTGGTCAGCCTTTGTTTAAGTTTTTATTGCACCAAGCg                                                                                          >  6:1409/1‑139 (MQ=255)
                                            gATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTATAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCggg                                                                                        <  4:149096/139‑1 (MQ=255)
                                            gATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCggg                                                                                        <  8:62647/139‑1 (MQ=255)
                                               atCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTAAGTTTTATTTGCACCGCGCGGGAAg                                                                                     >  8:105165/1‑139 (MQ=255)
                                                  aGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTTAAATGGCGGCTTTTTGTTGGGCCGCGTTTGTTTTAGTTTTATTGGCCCCCGGCGGGGCGCCg                                                                                  >  8:186001/1‑139 (MQ=255)
                                                       ccaccaCATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGa                                                                              <  7:105165/138‑1 (MQ=255)
                                                       ccaccaCATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTAGAACGCTTTTGTTTCAGTTGTATTTGCACCCAGCGGGAATCCGAttaa                                                                             >  8:45539/1‑136 (MQ=255)
                                                       ccaccaCATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCGTTTTAAATTGCGTGGTTTTTTTTGTACGCGTTTTTTTTGTTTTTTTTTGGACCCAGCCGGTAGCCCATGaa                                                                             >  4:203045/1‑138 (MQ=255)
                                                             caTCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGACGCCGATGACCGTACa                                                                       <  7:356693/139‑1 (MQ=255)
                                                              aTCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGGGTTTTTTTTTGCACGCTTTTTGTTAAGTTGTTGTTGCCCCCACCCCGCAGCCGCTGACCCTACa                                                                       >  7:75427/1‑138 (MQ=255)
                                                                  tCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTGAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGt                                                                 >  8:106376/1‑139 (MQ=255)
                                                                  tCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACAATATCCACTTCAGTGGCGACATTTGAAATTGCGTGGTTTGTTTAGAACCCTTTTGTTAGATATGTAATGGCAACCAGCGGGAAGCAGATGACCCGACAAACga                                                                  >  2:116750/1‑138 (MQ=255)
                                                                   cTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCTGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTc                                                                <  5:321919/139‑1 (MQ=255)
                                                                     tCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTTTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATTAACGTACAAACGTCGa                                                              >  8:356693/1‑139 (MQ=255)
                                                                      ccTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTc                                                                <  7:45539/136‑1 (MQ=255)
                                                                           tCTTTTAAAGCGGCAAATTTAATGACCAAATCCACTGCTCTGGCGCCTTTTGAAATGGCGTATTTTGTTTCTAACGCTTTGGTTTCTGGTGTAGTGGCCCCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAg                                                        <  5:120685/139‑1 (MQ=255)
                                                                           gcTTTTAAAGCGGAAATATTAATGACCATATCCACTTCAGGGGCGCCTTTTTAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCTGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAg                                                        <  8:223647/138‑1 (MQ=255)
                                                                            cTTTTAAAGCGGAAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGGATTTTTTTTCGAACGCTTTTGTTTCAGTTGTTTTGCCACCGAGCGGGAAGCCTATGACCGTACAAACGTCGACTCCAg                                                        <  5:1409/138‑1 (MQ=255)
                                                                                tAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTc                                                                <  7:22358/126‑1 (MQ=255)
                                                                                tAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTc                                                                >  8:22358/1‑126 (MQ=255)
                                                                                tAAAGCGGCAATATTAATGACCATATCAACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGTAAGATCTTGTTTAAGTTGTATTGTCACCCCGCCGGATTCCGAATTACGGACAAAAGTCGACACCAGTTccc                                                   >  2:44466/1‑139 (MQ=255)
                                                                                     cGGCAATATTAATGACCATATCCACTTCCGTGGCGCCTTTTGAAATGGCGGCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGCGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGTTCCCTTAAg                                              <  3:271476/139‑1 (MQ=255)
                                                                                       gCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGGTCCCGTAAGcc                                            >  7:223647/1‑138 (MQ=255)
                                                                                                attcccATATACAATTCAGTGGCGCCTTTGGATAGGGAGTCTTTTGTAGGTACCGCTTTTGTTTCTGTTGTCTTGCCCCCGCGCGGGACGCCCTTGACCGTACAACCGGCGACTCCAGTTCCCTTAAGCTCTTTTgcag                                   <  2:177560/135‑1 (MQ=255)
                                                                                                     cATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGc                              >  7:366546/1‑139 (MQ=255)
                                                                                                     cATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGc                              >  8:205465/1‑139 (MQ=255)
                                                                                                        atCCACTTCAGTGGCGCCTTTTGAACTGGCGTGTTTTGTTTCGAACGCTTTTGTTTCTGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTc                            <  1:116750/138‑1 (MQ=255)
                                                                                                        atCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGt                                                       <  8:121152/111‑1 (MQ=255)
                                                                                                        atCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGt                                                       >  7:121152/1‑111 (MQ=255)
                                                                                                              ttCAGTGGCGCCTTTTGAAATGGCGTCGTTTGTTTGGCACGCTTTTGTTTGAGTTGTTTTGGCACCGCGCGGGGCTCCGATGACCGTACAAACCGAGCCACCAGTTCCCGTAAGCGCTTTTTCAGCCAGCACCCCCCaa                     >  8:176265/1‑138 (MQ=255)
                                                                                                               tCAGTGGCGCCTTTTGAAATGGCGTGTTTTGTTTCGAAAGCGTTTGTTTGAGGTGTTTTGGCACCGAGCGGGGATCCGCAGACCGGTCAAACGTCGGCTGCAGTGGCCGTCAGCGCTGTTGCTGCCAGGGCCCCCCagg                    >  4:174142/1‑137 (MQ=255)
                                                                                                                 aGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTc                  >  1:218076/1‑139 (MQ=255)
                                                                                                                 aGTGGCGCCTTTTGAAATGGCGTCATTTGTTTCGAACGCATTTGTTTCAGTTGTATTGGCACCGAGCGGCAAGCCGAAGACCGTACAAACGTCGCCACCAGTTCCCATAAGCTCTTTTGCAGCAACCTCCACCCATGTc                  >  7:44426/1‑139 (MQ=255)
                                                                                                                   tGGCGCTTTTTTAATTGGGGTTTTTTTTTTGGAAGGCTTTTGTTTCTGTTGTTTTGGCACGTCGGGGTACGCCTATGACCGTAAAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGt                   <  6:49052/136‑1 (MQ=255)
                                                                                                                   tGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTGAGTTGTATTTGCACCGCGCGGGTAGCCGATGACCGTACAAACGTCGACTCCAGTTGCCGTTAGCTCTTTTGCAGCAAGCAGCCCCCATGTCg                 >  4:338913/1‑138 (MQ=255)
                                                                                                                    ggCGCCTTTTTAAATGGGGTTTTTTGTTTCCAACGCTTTTGTTTCTGTTGTTTTGGCACCGAGCGGTAAGCCGATGACCTTACAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGa               <  1:245802/139‑1 (MQ=255)
                                                                                                                    ggCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGa               <  7:205465/139‑1 (MQ=255)
                                                                                                                    ggCGCATTTTGAAATGGCGTGTTTTGTTTCGAACGCTTTTGTTTGAGTTGTTTTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGa               <  8:75427/139‑1 (MQ=255)
                                                                                                                     gcgcCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCCCCGAGCGGGAAGCCGATGACCGTACAAACGCTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGa               <  2:33133/139‑1 (MQ=255)
                                                                                                                                 aTGGCGTATTTTGTTTCCACCGCTTTTGTTGCAGTTGTATTGGCACCGTGCGGTAAGCCGATAACCGTACAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTg  <  1:44466/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
ACTTTGACAAGCGCTTTTCCGGCTACAGCTTCCACTACACCGCGGATATCAGCTTCCACCACATCGTCTTCCTTGTCTTTTAAAGCGGCAATATTAATGACCATATCCACTTCAGTGGCGCCTTTTGAAATGGCGTCTTTTGTTTCGAACGCTTTTGTTTCAGTTGTATTGGCACCGAGCGGGAAGCCGATGACCGTACAAACGTCGACTCCAGTTCCCTTAAGCTCTTTTGCAGCAAGCTCCACCCATGTCGGATTGACACATACTG  >  CP006881/3869755‑3870022

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