Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,728,699 C→T 11.5% intergenic (+51/+72) cspA‑I → / ← PP_1523 major cold shock protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,728,6990CT11.5% 130.1 / 8.6 52intergenic (+51/+72)cspA‑I/PP_1523major cold shock protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base T (3/3);  total (26/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 = 2.35e-01

GGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCC  >  NC_002947/1728556‑1728845
                                                                                                                                               |                                                                                                                                                    
gggCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGAAGCCACCACGGTggg                                                                                                                                                >  2:334950/1‑149 (MQ=255)
      cGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGGGCAGTGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTCGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGcttc                                                                                                                                          <  1:33773/150‑4 (MQ=255)
           ccTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCa                                                                                                                                   >  1:176405/1‑150 (MQ=255)
                 cGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGTCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGccc                                                                                                                             <  2:270225/150‑1 (MQ=255)
                    aGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTa                                                                                                                          <  2:546291/150‑1 (MQ=255)
                           cGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCgg                                                                                                                   <  2:392386/150‑1 (MQ=255)
                                  aGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAAcc                                                                                                            >  2:1137955/1‑150 (MQ=255)
                                        ccTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTcc                                                                                                      <  1:1137955/150‑1 (MQ=255)
                                         cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTcgcg                                                                                                                    <  2:129981/135‑1 (MQ=255)
                                         cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTcgcg                                                                                                                    >  1:129981/1‑135 (MQ=255)
                                           gCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTAc                                                                                                   <  2:661604/150‑1 (MQ=255)
                                           gCGGGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTAc                                                                                                   <  1:482021/150‑1 (MQ=255)
                                                                 gCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCCAGCCTGCCGtt                                                                             >  1:572332/1‑150 (MQ=255)
                                                                       cGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATGCGCGGGTAAACCCATTCCTACAGGAACGAGCAAGCCTGCCGTTTCatga                                                                       <  2:176405/150‑1 (MQ=255)
                                                                        gAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGa                                                                                                <  2:1151425/124‑1 (MQ=255)
                                                                        gAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGa                                                                                                >  1:1151425/1‑124 (MQ=255)
                                                                          aGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGc                                                                                                                                             <  1:772321/79‑1 (MQ=38)
                                                                          aGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGc                                                                                                                                             >  2:772321/1‑79 (MQ=38)
                                                                          aGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc                                                                                     >  1:133954/1‑133 (MQ=255)
                                                                          aGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc                                                                                     <  2:133954/133‑1 (MQ=255)
                                                                                      ccTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             <  2:517269/97‑1 (MQ=255)
                                                                                      ccTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                             >  1:517269/1‑97 (MQ=255)
                                                                                      ccTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAt                                                        >  1:488334/1‑150 (MQ=255)
                                                                                       ctctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCCGCCGTTTCATGATGTACGCCAGGTGAtt                                                       >  2:573856/1‑150 (MQ=255)
                                                                                        tctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAttt                                                      >  2:890726/1‑150 (MQ=255)
                                                                                              ccTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCtt                                                <  2:684381/150‑1 (MQ=255)
                                                                                              ccTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGGATTCGCGGGTAAACCCATTCCTACAGAAACGGGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCtt                                                >  1:684381/1‑150 (MQ=255)
                                                                                                  tGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGc                                                                                 >  1:1041489/1‑113 (MQ=255)
                                                                                                  tGGTTTCAAGACGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGc                                                                                 <  2:1041489/113‑1 (MQ=255)
                                                                                                     tttCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTc                                         >  2:108489/1‑150 (MQ=255)
                                                                                                       tCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATATCCCGTATTCGCGGGTAAACCCATTCCTACATAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCgg                                       <  1:193962/150‑1 (MQ=255)
                                                                                                         aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTc                                                     <  1:704067/134‑1 (MQ=255)
                                                                                                         aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTc                                                     >  2:704067/1‑134 (MQ=255)
                                                                                                         aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGgcg                                     >  2:178494/1‑150 (MQ=255)
                                                                                                         aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGgcg                                     <  1:890726/150‑1 (MQ=255)
                                                                                                         aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTAGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTTATGACGTACGCCCGATGATTTTCTCTTAACCGTgggac                                     >  1:249945/1‑148 (MQ=255)
                                                                                                           gCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt                                                           >  2:381962/1‑126 (MQ=255)
                                                                                                           gCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt                                                           <  1:381962/126‑1 (MQ=255)
                                                                                                                        gCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAg                      >  2:852612/1‑150 (MQ=255)
                                                                                                                              ccGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCa                >  2:456262/1‑150 (MQ=255)
                                                                                                                              ccGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCa                <  1:456262/150‑1 (MQ=255)
                                                                                                                               cgcgAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCNGGAAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCGAGGTGATTTACTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAg               >  1:831008/1‑150 (MQ=255)
                                                                                                                                gcgAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCg                             <  2:25556/135‑1 (MQ=255)
                                                                                                                                gcgAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCg                             >  1:25556/1‑135 (MQ=255)
                                                                                                                                ccgtggcagcCACCGTGGTGGCTGCTTCGCGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCa                >  1:409443/6‑150 (MQ=38)
                                                                                                                                ccgtggcagcCACCGTGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGc                     >  1:156590/6‑143 (MQ=255)
                                                                                                                                ccgtggcagcCACCGTGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGc                     <  2:156590/138‑1 (MQ=255)
                                                                                                                                  gAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCg                                        <  1:430004/122‑1 (MQ=255)
                                                                                                                                  gAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCg                                        >  2:430004/1‑122 (MQ=255)
                                                                                                                                   aaGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGAACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAtc           <  1:647202/150‑1 (MQ=255)
                                                                                                                                         ccaccaCGGTGGCTGT‑‑CACGGGCATATCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCttcttc     <  2:249945/150‑1 (MQ=255)
                                                                                                                                         ccaccaCGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCttcttc     <  1:108489/150‑1 (MQ=255)
                                                                                                                                           accacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTc                 >  2:422138/1‑136 (MQ=255)
                                                                                                                                           accacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTc                 <  1:422138/136‑1 (MQ=255)
                                                                                                                                           accacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTccc   <  2:488334/150‑1 (MQ=255)
                                                                                                                                            ccacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTcccc  <  1:191263/150‑1 (MQ=255)
                                                                                                                                            ccacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTcccc  >  2:5380/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCC  >  NC_002947/1728556‑1728845

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