Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,175,263 C→T E267E (GAG→GAA U712_10880 ← putative glycosyltransferase ypjH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,175,2630CT96.8% 184.2 / ‑3.1 64E267E (GAG→GAAU712_10880putative glycosyltransferase ypjH
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (38/23);  total (39/25)
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

CCTCCGGAATACCGCCAATGTTTGTTCCGATACAAGGCACTCCACAAGCCATCGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGG  >  CP006881/2175136‑2175390
                                                                                                                               |                                                                                                                               
ccTCCGGAATACCGCCAATGTTTGTTCCGATACAAGGCACTCCACAAGCCATCGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAaca                                                                                                                           <  6:63682/134‑1 (MQ=255)
ccTCCGGAATACCGCCAATGTTTGTTCCGATACAAGGCACTCCACAAGCCATCGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAaca                                                                                                                           >  5:63682/1‑134 (MQ=255)
ccTCAGGGATACCGCCAATGTTTGTTCCAATACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTc                                                                                                                      >  1:262877/1‑139 (MQ=255)
    agggaTACCGCCAATGTTTGTTCCAATACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCtt                                                                                                                    <  1:365596/133‑1 (MQ=255)
    agggaTACCGCCAATGTTTGTTCCAATACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGa                                                                                                                  >  2:407726/5‑139 (MQ=255)
           ccGCCAATGTTTGTTCCAATACAAGGCCCTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCaa                                                                                                           >  3:62832/1‑139 (MQ=255)
           ccGCCAATGTTTGTTCCAATACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAaca                                                                                                                           >  3:145276/1‑123 (MQ=255)
           ccGCCAATGTTTGTTCCAATACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTcc                                                                                                             >  8:387860/1‑137 (MQ=255)
           ccGCCAATGTTTGTTCCAATACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCaa                                                                                                           >  8:158078/1‑139 (MQ=255)
           ccGCCAATGTTTGTTCCAATACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCaa                                                                                                           >  8:261649/1‑139 (MQ=255)
                      tgttCCAATACAAGGCACTCCGCACGCCCTCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACt                                                                                                >  4:211595/1‑139 (MQ=255)
                      tgttCCAATACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAATCTCTTCATCACGGTCTTGATTTCCAAGCATTAAGACt                                                                                                >  7:239500/1‑139 (MQ=255)
                      tgttCCAATACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACt                                                                                                >  1:39864/1‑139 (MQ=255)
                         tCCAATACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGAACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTgg                                                                                             >  1:208108/1‑139 (MQ=255)
                            aaTACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGtct                                                                                          >  8:179793/2‑139 (MQ=255)
                             aTACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTcc                                                                                                             >  1:107232/1‑119 (MQ=255)
                             aTACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTcc                                                                                                             <  2:107232/119‑1 (MQ=255)
                              tACAAGGCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGtcttc                                                                                        >  5:55305/1‑139 (MQ=255)
                               aCAAGGCACTCCGCGAGCCATCGCGTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCATCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGtcttcc                                                                                       <  6:138452/139‑2 (MQ=255)
                                cAAGGCACTCCGCAAGCCATCGCTTCAAGCAGAACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGtcttcc                                                                                       <  2:208108/138‑2 (MQ=255)
                                   ggCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTg                                                                                              >  3:400499/1‑128 (MQ=255)
                                   ggCACTCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTg                                                                                              <  4:400499/128‑1 (MQ=255)
                                        tCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCa                                                                                                            <  8:74951/109‑1 (MQ=255)
                                        tCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCa                                                                                                            >  7:74951/1‑109 (MQ=255)
                                        tCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCatat                                                                              >  6:284017/1‑139 (MQ=255)
                                        tCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCatat                                                                              >  4:4755/1‑139 (MQ=255)
                                        tCCGCAAGCCATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCatat                                                                              >  8:151358/1‑139 (MQ=255)
                                                ccATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGat                                                                      <  2:262877/139‑2 (MQ=255)
                                                 cATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGata                                                                     <  6:143029/139‑3 (MQ=255)
                                                 cATCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTATCAGATAACAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGata                                                                     <  4:162694/139‑3 (MQ=255)
                                                  aTCGCTTCAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGat                                                                      >  4:182342/1‑136 (MQ=255)
                                                           agcagGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAgg                                                           <  7:179793/136‑1 (MQ=255)
                                                                gACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGat                                                                      <  8:137068/123‑2 (MQ=255)
                                                                gACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGat                                                                      >  7:137068/1‑122 (MQ=255)
                                                                gACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTc                                                      >  7:76905/1‑139 (MQ=255)
                                                                 aCAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTc                                                      >  1:99942/1‑138 (MQ=255)
                                                                  cAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGa                                                    <  2:39864/139‑1 (MQ=255)
                                                                     ggCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGAttt                                                 >  8:410469/1‑139 (MQ=255)
                                                                     ggCCAAAGCTATCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGAttt                                                 >  6:155535/1‑139 (MQ=255)
                                                                          aaGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTtctccg                                            >  3:423758/1‑137 (MQ=255)
                                                                          aaGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTtctccg                                            >  2:160103/1‑137 (MQ=255)
                                                                            gCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATCAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGa                                          <  1:160103/139‑1 (MQ=255)
                                                                                    tttCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGTTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGtctc                                   <  6:55305/138‑1 (MQ=255)
                                                                                    tttCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGtctc                                   <  3:4755/138‑1 (MQ=255)
                                                                                    tttCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGtctc                                   <  7:387860/138‑1 (MQ=255)
                                                                                      tCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAATACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCAccgc                                >  5:45931/1‑137 (MQ=255)
                                                                                       cAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGAc                               >  1:160818/1‑139 (MQ=255)
                                                                                              cagcagcTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCa                                                                                                            <  4:97930/55‑1 (MQ=255)
                                                                                              cagcagcTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCa                                                                                                            >  3:97930/1‑55 (MQ=255)
                                                                                              cagcagcTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGcc                                                                                  <  1:380974/81‑1 (MQ=255)
                                                                                              cagcagcTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGcc                                                                                  >  2:380974/1‑81 (MQ=255)
                                                                                                      tcaaaTCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGc                >  3:425177/3‑139 (MQ=255)
                                                                                                              cTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTccga                                >  4:266331/1‑113 (MQ=255)
                                                                                                              cTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTccga                                <  3:266331/115‑3 (MQ=255)
                                                                                                              cTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCTTTCGttt        >  5:25402/1‑139 (MQ=255)
                                                                                                              cTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCTTTCGttt        >  3:438252/1‑139 (MQ=255)
                                                                                                              cTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCTTTCGttt        >  8:436693/1‑139 (MQ=255)
                                                                                                              cTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCTTTCGttt        >  1:11540/1‑139 (MQ=255)
                                                                                                              cTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCGTTCGttt        >  3:33970/1‑139 (MQ=255)
                                                                                                                 atagAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCTTTCGTTTTGc     <  7:410469/136‑1 (MQ=255)
                                                                                                                    gAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTCAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTATAAGCAGCTTCGCTTTCGTTTTGccag  <  1:30528/139‑3 (MQ=255)
                                                                                                                    gAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCTTTCGTTTTGccag  <  7:158078/139‑3 (MQ=255)
                                                                                                                    gAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCTTTCGTTTTGccag  <  7:151358/139‑3 (MQ=255)
                                                                                                                    gAATAAAGATCTGCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCTTTCGTTTTGccag  <  3:211595/139‑3 (MQ=255)
                                                                                                                               |                                                                                                                               
CCTCCGGAATACCGCCAATGTTTGTTCCGATACAAGGCACTCCACAAGCCATCGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGG  >  CP006881/2175136‑2175390

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: