Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,249,662 A→C 7.1% K428Q (AAG→CAG)  ureC → urease subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,249,6620AC 7.1% 173.1 / 4.4 56K428Q (AAG→CAG) ureCurease subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (24/28);  new base C (2/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCC  >  NC_002947/3249518‑3249804
                                                                                                                                                  |                                                                                                                                                
cAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCt                                                                                                                                              <  3:315810/149‑1 (MQ=255)
  gATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTgg                                                                                                                                            >  2:276439/1‑149 (MQ=255)
    tGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGc                                                                                                                                           <  3:79350/148‑1 (MQ=255)
       aGCTGCGTCGCGGGCCGCTGGCACCGGAGACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGAc                                                                                                                                       <  7:257263/149‑1 (MQ=255)
        gCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGgacgaca                                                                                                                                      >  6:235708/1‑148 (MQ=255)
         cTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCCAGCTTGCCGACCt                                                                                                                                     <  1:242231/149‑1 (MQ=255)
                 cgGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCACACGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtg                                                                                                                              <  1:276439/149‑1 (MQ=255)
                          ggCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAg                                                                                                                    >  5:255374/1‑149 (MQ=255)
                                      cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGc                                                                                                                       >  1:251041/1‑134 (MQ=255)
                                      cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGc                                                                                                                       <  2:251041/134‑1 (MQ=255)
                                      cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCgg                                                                                                        >  5:93205/1‑149 (MQ=255)
                                      cccGTAAAGCTACAACTTCCGGGTCAACCGCTACATCGCCAAGTACACCATCAACCCGCCAATCACCCATCGCATAGGTCACGAAGTCGGCTCCGAGCAGGTCGGCAACCTCGCCGACCTGCTCATCTGGTCGCCACCGATATTTGCg                                                                                                         >  8:193256/1‑148 (MQ=255)
                                       ccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTc                                                                                                                         >  8:298600/1‑131 (MQ=255)
                                       ccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTc                                                                                                                         <  7:298600/131‑1 (MQ=255)
                                         gTACAGCGACAACTTCCGGTTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCa                                                                                                     <  7:269961/149‑1 (MQ=255)
                                            cAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGc                                                                                                  <  1:240961/149‑1 (MQ=255)
                                                gACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGc                                                                                              <  5:62056/149‑1 (MQ=255)
                                                 aCAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACATCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGc                                                                                              >  1:88857/1‑149 (MQ=255)
                                                    aCTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTg                                                                                          <  6:255374/149‑1 (MQ=255)
                                                      ttCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGt                                                                                        >  2:107697/1‑149 (MQ=255)
                                                      ttCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGt                                                                                        >  4:25873/1‑149 (MQ=255)
                                                                aGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGc                                                                            >  3:177167/1‑149 (MQ=255)
                                                                 gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGc                                                                                      <  3:78935/140‑1 (MQ=255)
                                                                 gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGc                                                                                      >  4:78935/1‑140 (MQ=255)
                                                                 gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCg                                                                           <  5:308936/149‑1 (MQ=255)
                                                                 gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCg                                                                           >  8:60125/1‑149 (MQ=255)
                                                                       cATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGgatga                                                                     >  5:155452/1‑149 (MQ=255)
                                                                       cATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGgatga                                                                     >  3:121181/1‑149 (MQ=255)
                                                                        tcccccAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGgatgat                                                                    <  5:130671/145‑1 (MQ=255)
                                                                         tCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGTGATGATc                                                                   <  2:54325/149‑1 (MQ=255)
                                                                         tCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATc                                                                   >  5:255378/1‑149 (MQ=255)
                                                                             cAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCa                                                               <  6:255378/149‑1 (MQ=255)
                                                                              aaGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCAc                                                              <  6:93205/149‑1 (MQ=255)
                                                                                    acCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGcc                                                        >  4:40806/1‑149 (MQ=255)
                                                                                          aaCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATgg                                                   <  3:25873/148‑1 (MQ=255)
                                                                                          aaCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATgg                                                   <  4:121181/148‑1 (MQ=255)
                                                                                           aCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGc                                                 <  2:88857/149‑1 (MQ=255)
                                                                                            cccGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGTTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCg                                                >  7:259882/1‑149 (MQ=255)
                                                                                             ccGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGa                                               >  5:98993/1‑149 (MQ=255)
                                                                                              cGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGAc                                              <  7:60125/149‑1 (MQ=255)
                                                                                               ggCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACa                                             >  5:349024/1‑149 (MQ=255)
                                                                                                 cACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATc                                           <  6:349024/149‑1 (MQ=255)
                                                                                                 cACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGGGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATc                                           <  4:310036/149‑1 (MQ=255)
                                                                                                       aCCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAAc                                        <  6:155452/147‑1 (MQ=255)
                                                                                                       aCCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACgg                                      >  1:236743/1‑149 (MQ=255)
                                                                                                                                  cTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCgcgc                                                         >  1:127446/1‑104 (MQ=255)
                                                                                                                                  cTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCgcgc                                                         <  2:127446/104‑1 (MQ=255)
                                                                                                                                  cTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGc            <  6:98993/149‑1 (MQ=255)
                                                                                                                                  cTCGGTGGAGGTAGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGGTCAACCTGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGt              <  1:107697/149‑1 (MQ=255)
                                                                                                                                      gtgtAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTa        >  3:24254/1‑149 (MQ=255)
                                                                                                                                      gTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTa        >  7:32927/1‑149 (MQ=255)
                                                                                                                                           ggTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGc            <  1:6415/140‑1 (MQ=255)
                                                                                                                                           ggTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGc            >  2:6415/1‑140 (MQ=255)
                                                                                                                                           ggTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGc   <  7:193256/149‑1 (MQ=255)
                                                                                                                                            gTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGcc  >  7:189838/1‑149 (MQ=255)
                                                                                                                                            gTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGcc  <  3:40806/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCC  >  NC_002947/3249518‑3249804

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