Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,178,145 C→T 6.9% I127I (ATC→ATT PP_1931 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,178,1450CT 6.9% 167.5 / 2.1 58I127I (ATC→ATTPP_1931hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/24);  new base T (2/2);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

TCCAGCCGCCGGCTTGAGCGTGAAGCTGAACGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTAC  >  NC_002947/2177999‑2178292
                                                                                                                                                  |                                                                                                                                                   
tCCAGCCGCCGGCTTGAGCGTGAAGCTGAACGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCct                                                                                                                                                   >  7:47154/1‑149 (MQ=255)
   aGCCGCCGGCTTGAGCGTGAAGCTGAACGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTg                                                                                                                                                >  3:266580/1‑149 (MQ=255)
    gccgccGGCTTGAGCGTGAAGCTGAACGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGc                                                                                                                                               <  6:39397/149‑1 (MQ=255)
                     gAAGCTGAACGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCttt                                                                                                                               >  6:236562/1‑148 (MQ=255)
                      aaGCTGAACGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCtttt                                                                                                                              <  4:266580/148‑1 (MQ=255)
                         cTGAACGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTctct                                                                                                                          >  5:260767/1‑149 (MQ=255)
                           gAACGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGa                                                                                                                        >  5:98470/1‑149 (MQ=255)
                            aaCGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGa                                                                                                                        <  7:39995/148‑1 (MQ=255)
                             aCGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAg                                                                                                                       <  8:47154/148‑1 (MQ=255)
                             aCGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGt                                                                                                                      >  1:153129/1‑149 (MQ=255)
                                     tCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATTCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCAtt                                                                                                              >  8:204884/1‑149 (MQ=255)
                                     tCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATTCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCAtt                                                                                                              >  8:210334/1‑149 (MQ=255)
                                         gTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGcc                                                                                                           >  4:159118/1‑148 (MQ=255)
                                                tgGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTtc                                                                                                                             <  1:43590/123‑1 (MQ=255)
                                                tgGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTtc                                                                                                                             >  2:43590/1‑123 (MQ=255)
                                                    tCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGgcag                                                                                                >  2:286251/1‑148 (MQ=255)
                                                      acaGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACTGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATTCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCaa                                                                                             <  7:204884/149‑1 (MQ=255)
                                                      acaGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATTCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCaa                                                                                             <  7:210334/149‑1 (MQ=255)
                                                      acaGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTctct                                                                                                                          >  1:167966/1‑120 (MQ=255)
                                                      acaGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTctct                                                                                                                          <  2:167966/120‑1 (MQ=255)
                                                              tCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGGCGGCAGCAAGTTCAAAg                                                                                     <  6:98470/149‑1 (MQ=255)
                                                                  gCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCa                                                                                   <  5:236562/147‑1 (MQ=255)
                                                                    tCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTg                                                                                >  2:191587/1‑148 (MQ=255)
                                                                    tCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTg                                                                                >  3:318682/1‑148 (MQ=255)
                                                                     ccGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGa                                                                               >  2:23453/1‑148 (MQ=255)
                                                                       ggATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTaa                                                                                                                                          >  4:107697/1‑87 (MQ=255)
                                                                       ggATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTaa                                                                                                                                          <  3:107697/87‑1 (MQ=255)
                                                                         aTTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAAt                                                                          <  6:260767/149‑1 (MQ=255)
                                                                          tttCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGAc                                                                                                     <  2:286430/121‑1 (MQ=255)
                                                                           ttCAAAACCATCGCCGACTTTCGCAAGGACTACGGTAATGCCATTCGTAGCGTCTGCCGTCAGTTCGTTATCCTCTGCCGTATCCTCAACGTTTTCTCTGAGGCGATCATTGGCAGCGTCGGCAGCTAGTTCAAAGCAGTGAACAATc                                                                         <  4:317262/148‑1 (MQ=255)
                                                                              aaaaCCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGa                                                                     >  4:271527/1‑149 (MQ=255)
                                                                                    aTCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTaa                                                               >  1:233962/1‑149 (MQ=255)
                                                                                     tggCCGACTTTCGAAAGGACAACGGTAAAGCCATTCGTAGCGGCTGCCTTCAGTTCGTGATCCTCTGCCTTAACCTCAACCTTTTCTCTGAGTCTATCATTGCCATCGACGGCATCAAGTTCAAAGCAGTGAACAATCGCGACCGTAAc                                                              <  2:153129/147‑1 (MQ=255)
                                                                                        ccGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTc                                                           >  7:281185/1‑149 (MQ=255)
                                                                                                 gCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGcc                                                  <  1:23453/149‑1 (MQ=255)
                                                                                                 gCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCAGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGACGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGc                                                   >  1:241024/1‑148 (MQ=255)
                                                                                                   aaGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGaa                                                                              >  2:232705/1‑119 (MQ=255)
                                                                                                   aaGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGaa                                                                              <  1:232705/119‑1 (MQ=255)
                                                                                                   aaGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCAcc                                                        >  8:85190/1‑141 (MQ=255)
                                                                                                   aaGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCAcc                                                        <  7:85190/141‑1 (MQ=255)
                                                                                                     ggACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAgg                                              <  2:233962/149‑1 (MQ=255)
                                                                                                     ggACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAg                                               <  3:160296/148‑1 (MQ=255)
                                                                                                     ggACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAg                                               >  4:160296/1‑148 (MQ=255)
                                                                                                              tAAATCCATTCGTAGGGTCTGCCGTCAGTTCGTGATCATCTGCCGTATCCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAATGCTGTGTACTATCGCGTCCGTAACTTCACCGAGGCCAAGGTGAAGGccc                                     <  2:241024/149‑1 (MQ=255)
                                                                                                                   ccATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGa                                           <  7:314080/138‑1 (MQ=255)
                                                                                                                   ccATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAATTAAAAGCAGTTAACAATCGCAACCGTAACTTCACCCAGGCCAAGGTGaaatatcgc                                   >  8:314080/1‑139 (MQ=255)
                                                                                                                    cATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATg                                >  3:135251/1‑148 (MQ=255)
                                                                                                                    cATTCGTAGCGGCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATgc                               <  8:281185/149‑1 (MQ=255)
                                                                                                                       tCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATgcag                             >  7:129864/1‑148 (MQ=255)
                                                                                                                         gTAGCGTCTGCTGTAAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATgcagca                           >  4:167637/1‑148 (MQ=255)
                                                                                                                            gCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATc                       <  3:159118/149‑1 (MQ=255)
                                                                                                                               tCTGCTGTTAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgag                    <  3:167637/149‑1 (MQ=255)
                                                                                                                                      tCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGAGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagcagagc              >  8:73448/1‑148 (MQ=255)
                                                                                                                                          ttCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCAtt           >  8:57961/1‑147 (MQ=255)
                                                                                                                                          ttCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCAtt           >  7:102722/1‑147 (MQ=255)
                                                                                                                                            cGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAgc       >  8:300499/1‑149 (MQ=255)
                                                                                                                                              tGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAgcgc     >  3:40677/1‑149 (MQ=255)
                                                                                                                                                 tCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTAc  <  7:300499/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCCAGCCGCCGGCTTGAGCGTGAAGCTGAACGCAATGTCGAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTAC  >  NC_002947/2177999‑2178292

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: