Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,653,773 T→G 10.7% I206S (ATT→AGT)  U712_08495 → GTP‑sensing transcriptional pleiotropic repressor CodY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,653,7730TG10.7% 186.4 / 13.9 64I206S (ATT→AGT) U712_08495GTP‑sensing transcriptional pleiotropic repressor CodY
Reads supporting (aligned to +/- strand):  ref base T (20/37);  new base G (2/5);  total (22/42)
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

GAAGAGGAAGCAAGAAGCAAAGCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTAGGAATGAAAGGTACTTATATCAAGGTACTAAACAACAAATTCCTAAT  >  CP006881/1653643‑1653908
                                                                                                                                  |                                                                                                                                       
gAAGAGGAAGCAAGAAGCAAAGCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGAc                                                                                                                                  <  1:95999/138‑1 (MQ=255)
gAAGAGGAAGCAAGAAGCAAAGCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGAc                                                                                                                                  <  1:8183/138‑1 (MQ=255)
   gagGAAGCAAGAAGCAAAGCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCgtg                                                                                                                              <  5:95297/139‑1 (MQ=255)
    agGAAGCAAGAAGCAAAGCTGTCGTACAAATGGCTATCAGCTCTCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCgtgt                                                                                                                             <  3:224132/139‑1 (MQ=255)
    agGAAGCAAGAAGCAAAGCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCgtgt                                                                                                                             <  7:63122/139‑1 (MQ=255)
      gaagcaagaagcaaAGCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCg                                                                                                                           >  2:222646/1‑139 (MQ=255)
          caagaagcaaAGCTGTCGTACATATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCAt                                                                                                                       <  4:291544/139‑1 (MQ=255)
             gaagcaaAGCTTTCGTACAAATTGCTATCAGCTCGCTTTGTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGTGGAGCTTGACGGAATTGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTAc                                                                                                                    <  8:5653/139‑1 (MQ=255)
                gcaaAGCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCg                                                                                                                 <  3:149126/139‑1 (MQ=255)
                gcaaAGCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAAGTGCTGACCGTGTCGGCATTACCCg                                                                                                                 <  3:94908/139‑1 (MQ=255)
                gcaaAGCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAAGTGCTGACCGTGTCGGCATTACCCg                                                                                                                 <  3:152344/139‑1 (MQ=255)
                  aaaGCTGTCGTACAAATGGCTATCAGCGCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGtt                                                                                                               <  1:141785/139‑1 (MQ=255)
                    aGCTGTCGTACAAATGGTTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCt                                                                                                             >  1:99324/1‑139 (MQ=255)
                     gCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTg                                                                                                            <  4:296982/139‑1 (MQ=255)
                     gCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTg                                                                                                            <  6:304721/139‑1 (MQ=255)
                         tCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCGTTACCCGTTCTGTTAt                                                                                                        >  7:287447/1‑139 (MQ=255)
                         tCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTCTTAt                                                                                                        >  4:144786/1‑139 (MQ=255)
                           gTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTTTTGCAAGTAAAATTGCTGACCGTGTCGGCGTTACCCGTTCTGTTATtg                                                                                                      >  6:141526/1‑139 (MQ=255)
                              cAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGa                                                                                                                                   <  1:126885/107‑1 (MQ=255)
                              cAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGa                                                                                                                                   >  2:126885/1‑107 (MQ=255)
                               aaaTGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGaa                                                                                                  <  8:213505/139‑1 (MQ=255)
                                aaTGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAAc                                                                                                 >  5:5718/1‑139 (MQ=255)
                                  tGGCTATCAGCTCGCTTTCTTACAGAGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAAGTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGc                                                                                               <  8:99714/139‑1 (MQ=255)
                                        tCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCa                                                                                          <  6:5718/138‑1 (MQ=255)
                                        tCAGCTAGCTTTCATACAGGGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCa                                                                                          <  6:98334/138‑1 (MQ=255)
                                                        cAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAgcgc                                                                          <  4:99874/138‑1 (MQ=255)
                                                           tGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCgct                                                                      >  7:56409/1‑137 (MQ=255)
                                                           tGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAAGTGCTGACCGTGTCGGCATTACCCGTTCTGt                                                                                                           >  1:209016/1‑102 (MQ=255)
                                                           tGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAAGTGCTGACCGTGTCGGCATTACCCGTTCTGt                                                                                                           <  2:209016/102‑1 (MQ=255)
                                                             aGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGATGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGgtg                                                                     >  1:110839/1‑138 (MQ=255)
                                                             aGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGATGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGgtg                                                                     <  2:110839/138‑1 (MQ=255)
                                                             aGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCgtgt                                                                                                                             <  6:179444/82‑1 (MQ=255)
                                                             aGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCgtgt                                                                                                                             >  5:179444/1‑82 (MQ=255)
                                                             aGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGtt                                                                                                               >  5:211645/1‑96 (MQ=255)
                                                             aGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGtt                                                                                                               <  6:211645/96‑1 (MQ=255)
                                                             aGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGgtgt                                                                    >  5:213082/1‑139 (MQ=255)
                                                             aGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAAGTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAgcgc                                                                          <  4:86079/133‑1 (MQ=255)
                                                             aGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAAGTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAgcgc                                                                          >  3:86079/1‑133 (MQ=255)
                                                             aGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTAGACGGAAATGAAGGTCTTCTTGTTGCAATTAAAATTGCTGACAGTGTCGGCATTACCCGTTCTGTTATTGTCAACGCACTCAGAAAGCTGGAGAGCGCCGatgt                                                                    >  2:7988/1‑139 (MQ=255)
                                                               cTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGtt                                                                                                               >  1:77536/1‑94 (MQ=255)
                                                               cTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGtt                                                                                                               <  2:77536/94‑1 (MQ=255)
                                                                  gAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTAt                                                                 <  1:222646/137‑1 (MQ=255)
                                                                  gAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATc                                                                <  5:123740/138‑1 (MQ=255)
                                                                         ttGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCa                                                                                          >  8:21984/1‑105 (MQ=255)
                                                                         ttGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCa                                                                                          <  7:21984/105‑1 (MQ=255)
                                                                         ttGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAgcgc                                                                          >  2:278961/1‑121 (MQ=255)
                                                                         ttGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAgcgc                                                                          <  1:278961/121‑1 (MQ=255)
                                                                              cacaTTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCAt                                                   >  8:105052/1‑139 (MQ=255)
                                                                               acaTTTTTGAGCAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGTCATTACCCGTGCTGTTATTGTGAACGCACTCAGAAAGCTGGCGAGCGCCGGTGTTATCGAGTCTAGATCAtt                                                  <  4:118897/139‑1 (MQ=255)
                                                                                caTTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTa                                                 <  6:87762/139‑1 (MQ=255)
                                                                                      tGAGGAGCTTGTCGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGCCTAGATCATTAGGAAt                                            <  3:293452/138‑1 (MQ=255)
                                                                                      tGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGtt                                                                                                               <  4:13894/71‑1 (MQ=255)
                                                                                      tGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGtt                                                                                                               >  3:13894/1‑71 (MQ=255)
                                                                                       gaggagCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTAGGAATg                                           <  4:39639/138‑1 (MQ=255)
                                                                                        atgagGTTGACGGAAATGAAGTTCTTCTTGTTGCAAGTAAAATTTCTGACCTTGTCGGCATTACCCGTTGTGTTATTGTGTACGCACTCAGAAAGCTGGAGATCGCCTGTGTTATCGAGGCTAGATCATTatgaatga                                          <  1:7988/136‑1 (MQ=255)
                                                                                        aggggCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTAGGAATGaa                                         <  6:128653/139‑1 (MQ=255)
                                                                                            gCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCa                                                                                          >  7:44993/1‑86 (MQ=255)
                                                                                            gCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCa                                                                                          <  8:44993/86‑1 (MQ=255)
                                                                                             cTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTAGGAATGAAAGGTa                                    >  4:72942/1‑139 (MQ=255)
                                                                                                gACGGAAATGAAGGTCTTCTTGTTGCAAGTAATATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTAGGAATGAAAGGTACt                                  <  1:262411/138‑1 (MQ=255)
                                                                                                gACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCGTTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTAGGAAAGAAAGGTACt                                  <  8:287447/138‑1 (MQ=255)
                                                                                                      aaTGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTAGGAATGAAAGGTACTTATATCa                           <  3:72942/139‑1 (MQ=255)
                                                                                                              tcttctTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTAGGAATGAAAGGTACTTATATCAAGGTACTa                   <  7:105052/139‑1 (MQ=255)
                                                                                                                               aaaTTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTAGGAATGAAAGGTACTTATATCAAGGTACTAAACAACAAATTCCTAAt  <  2:99324/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
GAAGAGGAAGCAAGAAGCAAAGCTGTCGTACAAATGGCTATCAGCTCGCTTTCTTACAGTGAGCTTGAAGCAATTGAGCACATTTTTGAGGAGCTTGACGGAAATGAAGGTCTTCTTGTTGCAAGTAAAATTGCTGACCGTGTCGGCATTACCCGTTCTGTTATTGTGAACGCACTCAGAAAGCTGGAGAGCGCCGGTGTTATCGAGTCTAGATCATTAGGAATGAAAGGTACTTATATCAAGGTACTAAACAACAAATTCCTAAT  >  CP006881/1653643‑1653908

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