Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,057,368:1 +G 7.3% coding (294/2406 nt) PP_3575 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,057,3681.G 7.3% 209.5 / 11.6 55coding (294/2406 nt)PP_3575ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base . (25/26);  new base G (2/2);  total (27/28)
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.70e-01

GGCGTGCCACTTGGCCACCGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCG  >  NC_002947/4057226‑4057515
                                                                                                                                                |                                                                                                                                                    
ggCGTGCCACTTGGCCACCGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTg                                                                                                                                                 <  7:236522/148‑1 (MQ=255)
 gCGTGCCACTTGGCCACCGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTg                                                                                                                                                 >  8:254546/1‑147 (MQ=255)
  cGTGCCACTTGGCCACCGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAg                                                                                                                                              <  7:254546/149‑1 (MQ=255)
   gTGCCACTTGGCCACCGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAgg                                                                                                                                             >  2:138763/1‑149 (MQ=255)
    tGCCACTTGGCCACCGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGc                                                                                                                                            >  1:39825/1‑149 (MQ=255)
        aCTTGGCCACCGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGc                                                                                                                                         >  6:285778/1‑148 (MQ=255)
            ggCCACCGGCTTACGGGTGTGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACt                                                                                                                                     <  8:250409/148‑1 (MQ=255)
            ggCCACCGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACt                                                                                                                                     <  1:138763/148‑1 (MQ=255)
                 ccGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCCACTGGTAGCGCCTGCAGGCCGCTACTGCgca                                                                                                                                <  1:46406/149‑1 (MQ=255)
                 ccGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCgcag                                                                                                                               <  2:39825/149‑1 (MQ=255)
                    gCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCaa                                                                                                                            >  4:338708/1‑149 (MQ=255)
                      aTACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACt                                                                                                                          >  7:194187/1‑149 (MQ=255)
                        aCGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCg                                                                                                                                          >  8:211560/1‑131 (MQ=255)
                        aCGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCg                                                                                                                                          <  7:211560/131‑1 (MQ=255)
                               cgcTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGaa                                                                                                                   <  3:99797/147‑1 (MQ=255)
                                  ttCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCtt                                                                                                              >  2:240109/1‑149 (MQ=255)
                                          cTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGaca                                                                                                     >  2:209305/1‑149 (MQ=255)
                                                ccAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTTAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTaa                                                                                               <  5:181918/149‑1 (MQ=255)
                                                       gccgACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCg                                                                                         >  3:234537/1‑148 (MQ=255)
                                                         cgacgaTGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGa                                                                                       >  2:315504/1‑148 (MQ=255)
                                                            cgaTGCGGGTGGGGGCGATTTCAAGCGCCCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAg                                                                                   <  3:338708/149‑1 (MQ=255)
                                                                  gggTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTCACGACCGGACAAGCCGTTa                                                                             >  4:202707/1‑149 (MQ=255)
                                                                      gggggCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCAGTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACttt                                                                         <  8:194187/149‑1 (MQ=255)
                                                                       ggggCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACttt                                                                         >  4:323356/1‑148 (MQ=255)
                                                                       ggggCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACttt                                                                         <  3:323356/148‑1 (MQ=255)
                                                                              tttCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCg                                                                     >  8:228535/1‑145 (MQ=255)
                                                                              tttCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCg                                                                     <  7:228535/145‑1 (MQ=255)
                                                                               ttCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCGATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGcc                                                                  <  1:240109/148‑1 (MQ=255)
                                                                                   aGCGCGCCACGGGCCTGCGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTAc                                                            <  2:92940/149‑1 (MQ=255)
                                                                                           aCTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCAACCAATGc                                                      >  3:311504/1‑147 (MQ=255)
                                                                                           aCTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCAACCAATGc                                                      <  4:311504/147‑1 (MQ=255)
                                                                                             tGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCggg                                                  >  1:28866/1‑149 (MQ=255)
                                                                                               gCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTc                                                <  5:285778/149‑1 (MQ=255)
                                                                                                 cTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTc                                                >  4:27533/1‑147 (MQ=255)
                                                                                                  tGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAg                                             >  1:154476/1‑149 (MQ=255)
                                                                                                         cTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAgg                                      >  2:39455/1‑149 (MQ=255)
                                                                                                          tGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGc                                     >  5:112835/1‑149 (MQ=255)
                                                                                                            gCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCt                                    <  5:271589/148‑1 (MQ=255)
                                                                                                              aTCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGa                                  >  2:33332/1‑148 (MQ=255)
                                                                                                               tCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGaca                                >  6:140675/1‑149 (MQ=255)
                                                                                                                 aGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGacac                               <  1:33332/148‑1 (MQ=255)
                                                                                                                     tGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGcc                          <  3:27533/149‑1 (MQ=255)
                                                                                                                       aaGCCCCCCTACCCCACTGGTAGGCCCTGCAGGCCGCTACTTCGCAGCAACTGGTTGAAGCCTTCATCGAAACTGTACCGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGACAGGCTGACACCGGCCTg                        <  3:237890/149‑1 (MQ=255)
                                                                                                                           cccccTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTCACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCg                    <  3:202707/149‑1 (MQ=255)
                                                                                                                            ccccTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCgg                   >  2:176495/1‑149 (MQ=255)
                                                                                                                                  aCCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAAt             <  4:234537/149‑1 (MQ=255)
                                                                                                                                     tgcgcagTAGCGGCCTGCAGGCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCg                 >  5:155970/7‑143 (MQ=255)
                                                                                                                                     tgcgcagTAGCGGCCTGCAGGCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCg                 <  6:155970/137‑1 (MQ=255)
                                                                                                                                     tgcgcagTAGCGGCCTGCAGGCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATgg           >  3:290806/7‑149 (MQ=255)
                                                                                                                                     ccACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATgg           <  1:315504/149‑1 (MQ=255)
                                                                                                                                     ccACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATgg           >  2:297493/1‑149 (MQ=255)
                                                                                                                                     ccACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATgg           <  5:140675/149‑1 (MQ=255)
                                                                                                                                     ccACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATgg           <  1:209305/149‑1 (MQ=255)
                                                                                                                                       cgcagTAGCGGCCTGCAGGCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATgg           <  4:290806/143‑1 (MQ=255)
                                                                                                                                         tGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTc        >  3:159534/1‑148 (MQ=255)
                                                                                                                                              gcgcCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAcgcg  >  6:350526/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGCGTGCCACTTGGCCACCGGCATACGGCTGCGCTTCGTTGCCTGCGGCCAGGCCGCCGACGATGCGGGTGGGGGCGATTTCAAGCGCGCCACTGGCCTGAGGGGCTGGCATCAGGGTGAAGGCCCCCTCACCCACTGGTAGCGCCTGCAGGCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCG  >  NC_002947/4057226‑4057515

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: