Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,264,918 A→T 7.1% D44V (GAC→GTC)  vanB → vanillate O‑demethylase oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,264,9180AT 7.1% 157.5 / 3.3 56D44V (GAC→GTC) vanBvanillate O‑demethylase oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base T (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

TGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCAGTGCG  >  NC_002947/4264773‑4265066
                                                                                                                                                 |                                                                                                                                                      
tgAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACgg                                                                                                                                                     >  7:89865/1‑149 (MQ=255)
tgAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGCGCACCAGCCCGTCgg                                                                                                                                                     >  7:2653/1‑149 (MQ=255)
     aaCAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTg                                                                                                                                                 <  2:208832/148‑1 (MQ=255)
      aCAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCGGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGt                                                                                                                                               <  1:124132/149‑1 (MQ=255)
      aCAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGCGCACCAGCCCGTCGGGCAgg                                                                                                                                                <  8:2653/148‑1 (MQ=255)
       cAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTg                                                                                                                                              >  1:184784/1‑149 (MQ=255)
          ccGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGCGCACCAGCCCGTCGGGCAGGTgcac                                                                                                                                           <  3:230120/149‑3 (MQ=255)
            gACATGATCGATGCCGTAGTGGTAGCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCa                                                                                                                                         <  8:193944/149‑1 (MQ=255)
             aCATGATCGATGCCGTAGTGGTATCCTGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAg                                                                                                                                        >  3:193078/1‑149 (MQ=255)
               aTGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACg                                                                                                                                                      <  6:42551/133‑1 (MQ=255)
               aTGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACg                                                                                                                                                      >  5:42551/1‑133 (MQ=255)
                          cGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTc                                                                                                                                   >  7:243825/1‑141 (MQ=255)
                          cGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTc                                                                                                                                   <  8:243825/141‑1 (MQ=255)
                                        gTAACGATGAAGCGCAGGGTATCTGCAGCTNCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCgaagaa                                                                                                             >  1:26659/1‑149 (MQ=255)
                                         tAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAAc                                                                                                            <  6:146554/149‑1 (MQ=255)
                                          aaCGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCg                                                                                                                                  <  1:159994/126‑1 (MQ=255)
                                          aaCGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCCACGGGCTGGTGCGCCAGTATTCg                                                                                                                                  >  2:159994/1‑126 (MQ=255)
                                             gATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCAcccc                                                                                                                   <  8:114470/138‑1 (MQ=255)
                                             gATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCAcccc                                                                                                                   >  7:114470/1‑138 (MQ=255)
                                                         ggTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGAt                                                                                          >  2:206800/1‑149 (MQ=255)
                                                          gTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGAtt                                                                                         >  8:39611/1‑149 (MQ=255)
                                                               tGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCg                                                                                     >  2:235581/1‑148 (MQ=255)
                                                                gCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTAATCGGTGTGCAACCCCCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGt                                                                                    >  2:140880/1‑148 (MQ=255)
                                                                        gAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACg                                                                            <  1:206800/148‑1 (MQ=255)
                                                                          gCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGAcc                                                                         >  1:242718/1‑149 (MQ=255)
                                                                              gCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAAc                                                                             >  8:210415/1‑141 (MQ=255)
                                                                              gCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAAc                                                                             <  7:210415/141‑1 (MQ=255)
                                                                                cgcgGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCtt                                                                   <  1:235581/149‑1 (MQ=255)
                                                                                     cAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCgg                                                              <  2:26659/149‑1 (MQ=255)
                                                                                        aTGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGTTTCGCGGGGc                                                           <  8:108701/149‑1 (MQ=255)
                                                                                         tGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCg                                                          <  4:90578/149‑1 (MQ=255)
                                                                                            cAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCgg                                                                      <  1:291193/134‑1 (MQ=255)
                                                                                            cAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCgg                                                                      >  2:291193/1‑134 (MQ=255)
                                                                                                ctgctgCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCg                                                                                                   <  5:136633/101‑1 (MQ=255)
                                                                                                ctgctgCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCg                                                                                                   >  6:136633/1‑101 (MQ=255)
                                                                                                   ctgcCGGCGTTCAGCGCCGGCGCCCATATCGACGTGTACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAg                                                >  5:114287/1‑149 (MQ=255)
                                                                                                   ctgcCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAg                                                >  4:110745/1‑149 (MQ=255)
                                                                                                         gCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCa                                          >  5:283654/1‑149 (MQ=255)
                                                                                                             tCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCgg                                                         <  5:211811/130‑1 (MQ=255)
                                                                                                             tCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCgg                                                         >  6:211811/1‑130 (MQ=255)
                                                                                                              cAGCGCCGGCGCCCATATCGACGTGTACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGaa                                      <  6:114287/148‑1 (MQ=255)
                                                                                                               aGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAAc                                     <  4:193078/148‑1 (MQ=255)
                                                                                                                 cggcGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTAGCTGTGCAACCAACCCGAAGAACGCCACCACGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACa                                    <  1:257881/146‑1 (MQ=255)
                                                                                                                        gcCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGc                           <  7:39611/149‑1 (MQ=255)
                                                                                                                          ccATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCg                         <  2:151246/149‑1 (MQ=255)
                                                                                                                          ccATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCg                         <  1:299267/149‑1 (MQ=255)
                                                                                                                              atCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGtt                                                       <  2:42383/115‑1 (MQ=255)
                                                                                                                              atCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGtt                                                       >  1:42383/1‑115 (MQ=255)
                                                                                                                                 gACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCg                  >  4:185251/1‑149 (MQ=255)
                                                                                                                                  tggcgCACCAGCCCGTCGGGCAGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCgg                 >  8:107839/5‑149 (MQ=255)
                                                                                                                                    gTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCg                                                   <  2:28248/113‑1 (MQ=255)
                                                                                                                                    gTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCg                                                   >  1:28248/1‑113 (MQ=255)
                                                                                                                                    gTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCt               >  4:126064/1‑149 (MQ=255)
                                                                                                                                        aCCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCg            >  6:177398/1‑148 (MQ=255)
                                                                                                                                             cccGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑CCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCa       <  2:184784/148‑1 (MQ=255)
                                                                                                                                               cGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCAGTg    <  8:89865/149‑1 (MQ=255)
                                                                                                                                                 aCGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCAGTgcg  <  6:283654/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
TGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCA‑‑TCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCAGTGCG  >  NC_002947/4264773‑4265066

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: