Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 627,738 C→T 6.9% H15H (CAC→CAT entB → isochorismatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913627,7380CT 6.9% 183.6 / 3.0 58H15H (CAC→CATentBisochorismatase
Reads supporting (aligned to +/- strand):  ref base C (32/22);  new base T (2/2);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATTGCTGCGCTGCGC  >  NC_000913/627593‑627885
                                                                                                                                                 |                                                                                                                                                   
gtgtgtGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGCTGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGCGTCTCATGa                                                                                                                                                   >  6:289933/1‑148 (MQ=255)
gtgtgtGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCATGa                                                                                                                                                   >  6:289929/1‑148 (MQ=255)
 tgtgtgGATTCACTTCAGCTGTCGGCGGTCGGGAAAGTCTATAAAAAACAATTACGTCATTGGCTGGCGTCTCGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGat                                                                                                                                                  <  7:29051/148‑1 (MQ=255)
    gtgGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATAtt                                                                                                                                               >  1:368573/1‑148 (MQ=255)
    gtgGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATAtt                                                                                                                                               >  3:33017/1‑148 (MQ=255)
     tgGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTc                                                                                                                                              >  7:43991/1‑148 (MQ=255)
           cACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAg                                                                                                                                         >  5:96332/1‑147 (MQ=255)
             cTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAAt                                                                                                                                      <  8:43991/148‑1 (MQ=255)
                      aCGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTTGCGTCACGCGCATCAGCCTGAAGGAGAGAACACTATGGCTATTCCAAAATTACAGGCTTACGCACTTCCGGAGTCTCATGATATTCCGCAGAATAAAGTTg                                                                                                                               <  5:289933/146‑1 (MQ=255)
                      aCGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCATGATATTCCGCAGAATAAAGTTg                                                                                                                               <  5:289929/146‑1 (MQ=255)
                      aCGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTg                                                                                                                               <  2:368573/146‑1 (MQ=255)
                      aCGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTg                                                                                                                               <  2:35369/146‑1 (MQ=255)
                      aCGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTg                                                                                                                               <  6:249463/146‑1 (MQ=255)
                      aCGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTg                                                                                                                               <  4:33017/146‑1 (MQ=255)
                                             aaaaCAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGc                                                                                                      <  6:237100/148‑1 (MQ=255)
                                                               gctggcGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGa                                                                                    >  8:211175/1‑148 (MQ=255)
                                                                                cAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGCCTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTg                                                                   >  8:375172/1‑148 (MQ=255)
                                                                                  gCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTc                                                                 <  6:96332/148‑1 (MQ=255)
                                                                                   ccTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGt                                                                  >  8:43098/1‑146 (MQ=255)
                                                                                      gAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCt                                                             >  7:225339/1‑148 (MQ=255)
                                                                                              gaACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGat                                                                                   <  5:266793/118‑1 (MQ=255)
                                                                                              gaACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGat                                                                                   >  6:266793/1‑118 (MQ=255)
                                                                                              gaACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGc                                                     >  3:115509/1‑148 (MQ=255)
                                                                                               aaCACGATTGCTATTCCAAAATTAGAGGCTTACGCAGTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGat                                                                                   <  6:357164/117‑1 (MQ=255)
                                                                                               aaCACGATGGCTATTCCTAAATTACAGGCTTCCGCACTGGCGTAGTCTCACGATATTCCGCAGAATATAGTTGACTGGGCCTTTGAACCGCAACGTGCGGCGTTGTTAATCCATGat                                                                                   <  1:69164/117‑1 (MQ=255)
                                                                                                acacGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGTAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGCACCGCAACGTGCCGCGTTGTTAATCCATGata                                                                                  >  2:69164/1‑117 (MQ=255)
                                                                                                  acGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCgaga                                                 <  7:211175/148‑1 (MQ=255)
                                                                                                   cGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGaa                                                <  7:43098/148‑1 (MQ=255)
                                                                                                        gCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCt                                                             >  2:197959/1‑130 (MQ=255)
                                                                                                        gCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCt                                                             <  1:197959/130‑1 (MQ=255)
                                                                                                            ttCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCgaga                                                 <  5:124217/138‑1 (MQ=255)
                                                                                                            ttCAAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCgaga                                                 >  6:124217/1‑138 (MQ=255)
                                                                                                             tCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCgat                                       >  1:258044/1‑147 (MQ=255)
                                                                                                              ccAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCgatga                                     >  2:394250/1‑148 (MQ=255)
                                                                                                               cAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCgatgat                                    <  7:375172/148‑1 (MQ=255)
                                                                                                               cAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCgatg                                      <  5:330297/146‑1 (MQ=255)
                                                                                                               cAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCgatg                                      >  6:330297/1‑146 (MQ=255)
                                                                                                                     tACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTACTTTGTCAGCTTCTGGGGCg                                                    <  3:331740/126‑1 (MQ=255)
                                                                                                                     tACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTAATTTGTCAGCTTCTGGGGCg                                                    >  4:331740/1‑126 (MQ=255)
                                                                                                                               cGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCtt                                                            <  8:3423/108‑1 (MQ=255)
                                                                                                                               cGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCtt                                                            >  7:3423/1‑108 (MQ=255)
                                                                                                                               cGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATcgc                     >  4:215074/1‑147 (MQ=255)
                                                                                                                                     gCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGatgatg                                   >  2:177832/1‑127 (MQ=255)
                                                                                                                                     gCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGatgatg                                   <  1:177832/127‑1 (MQ=255)
                                                                                                                                     gCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTg                          >  8:325981/1‑136 (MQ=255)
                                                                                                                                     gCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTg                          <  7:325981/136‑1 (MQ=255)
                                                                                                                                     gCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATAtt              >  5:195439/1‑148 (MQ=255)
                                                                                                                                          aGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATtgctgc         >  4:321527/1‑148 (MQ=255)
                                                                                                                                           gTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATTgctgcg        >  4:194748/1‑148 (MQ=255)
                                                                                                                                           gTCACACGATATTCCGAAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTCCCAGATGATGGAGCAGGTGATCGCGATTATTgctccg        >  4:199523/1‑148 (MQ=255)
                                                                                                                                             ctcACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGATCAGGAGATCGCGAATATTgctgcgct      >  2:136566/1‑148 (MQ=255)
                                                                                                                                             ctcACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAgg                            >  6:204892/1‑126 (MQ=255)
                                                                                                                                             ctcACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAgg                            <  5:204892/126‑1 (MQ=255)
                                                                                                                                             ctcACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATTgctgcgct      >  6:218360/1‑148 (MQ=255)
                                                                                                                                               cACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGACCCGCAACGTGACGCGTTGTTAATCCATGATATGCAGGACGATTTTGTCATATTCTGGGGCGAGAACTGCCCGATGATGGa                                 >  7:278974/1‑119 (MQ=255)
                                                                                                                                               cACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATTgctgcgctgc    >  3:348798/1‑148 (MQ=255)
                                                                                                                                               cACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATTgctgcgctgc    >  1:207402/1‑148 (MQ=255)
                                                                                                                                               cACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATTgctgcgctgc    >  2:251574/1‑148 (MQ=255)
                                                                                                                                                 cGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATTgctgcgctgc    >  4:145334/1‑146 (MQ=255)
                                                                                                                                                 cGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATTgctgcgctgc    <  3:145334/146‑1 (MQ=255)
                                                                                                                                                 cGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATTGCtgcgctgcgc  <  6:195439/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTGTGTGGATTCACTTCCGCTGACGGCGGTCGGGAAAGTCGATAAAAAACAATTACGTCAGTGGCTGGCGTCACGCGCATCAGCCTGAAGGAGAGAACACGATGGCTATTCCAAAATTACAGGCTTACGCACTGCCGGAGTCTCACGATATTCCGCAGAATAAAGTTGACTGGGCCTTTGAACCGCAACGTGCCGCGTTGTTAATCCATGATATGCAGGACTATTTTGTCAGCTTCTGGGGCGAGAACTGCCCGATGATGGAGCAGGTGATCGCGAATATTGCTGCGCTGCGC  >  NC_000913/627593‑627885

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