Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,151,492 (G)6→7 intergenic (+88/‑123) fabG → / → acpP 3‑oxoacyl‑[acyl‑carrier‑protein] reductase/acyl carrier protein (ACP)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,151,4861.G98.1% 199.9 / ‑2.6 52intergenic (+82/‑129)fabG/acpP3‑oxoacyl‑[acyl‑carrier‑protein] reductase/acyl carrier protein (ACP)
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (34/17);  total (34/18)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

ATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGAAATTTAAGAGTATGAGCACTATCG  >  NC_000913/1151338‑1151627
                                                                                                                                                     |                                                                                                                                             
aTCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTgg                                                                                                                                             >  1:538803/1‑151 (MQ=255)
  ccGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTgggg                                                                                                                                           >  2:705779/1‑151 (MQ=255)
  ccGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTggg                                                                                                                                            <  2:912836/150‑1 (MQ=255)
    gacgaAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTgggggg                                                                                                                                         <  1:884308/151‑1 (MQ=255)
        aaGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTgggggggta                                                                                                                                      >  2:17674/1‑147 (MQ=255)
           cagcTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGc                                                                                                                                   <  2:573115/151‑1 (MQ=255)
              cTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTc                                                                                                                                >  1:524788/1‑151 (MQ=255)
              cTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTc                                                                                                                                >  2:939382/1‑151 (MQ=255)
                 aCATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCaa                                                                                                                              >  1:452069/1‑150 (MQ=255)
                   aTCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAAt                                                                                                                           >  1:26569/1‑151 (MQ=255)
                     cACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATa                                                                                                                          <  2:598237/150‑1 (MQ=255)
                       cGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACg                                                                                                                       >  1:409508/1‑151 (MQ=255)
                          gTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTaa                                                                                                                    <  2:840108/151‑1 (MQ=255)
                                  ttGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTgg                                                                                                            >  1:503224/1‑151 (MQ=255)
                                       tgtgAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTa                                                                                                          >  2:576380/1‑148 (MQ=255)
                                             cggcggGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGcc                                                                                                 >  1:613071/1‑151 (MQ=255)
                                                   gATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGAtt                                                                                           >  2:107821/1‑151 (MQ=255)
                                                   gATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGAtt                                                                                           <  1:810528/151‑1 (MQ=255)
                                                      gTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGAttt                                                                                          >  2:933747/1‑149 (MQ=255)
                                                       tACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGt                                                                                       >  1:526572/1‑151 (MQ=255)
                                                           tGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTg                                                                                     >  2:145506/1‑149 (MQ=255)
                                                              tCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAAt                                                                                >  2:507712/1‑151 (MQ=255)
                                                                 gACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAAtttt                                                                             <  2:25359/151‑1 (MQ=255)
                                                                       gATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAAc                                                                        >  1:685986/1‑150 (MQ=255)
                                                                        aTTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACAt                                                                      >  2:842154/1‑151 (MQ=255)
                                                                         tttGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACAtt                                                                     >  1:427951/1‑151 (MQ=255)
                                                                              acaAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTtat                                                                 >  1:714164/1‑150 (MQ=255)
                                                                                 aaaTGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATacac                                                             >  1:482832/1‑151 (MQ=255)
                                                                                   aTGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACt                                                            >  1:821796/1‑150 (MQ=255)
                                                                                    tGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTAc                                                          >  2:522934/1‑151 (MQ=255)
                                                                                    tGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTAc                                                          >  1:191299/1‑151 (MQ=255)
                                                                                     gCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACg                                                         >  1:593074/1‑151 (MQ=255)
                                                                                        cATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGaaa                                                      <  2:705547/151‑1 (MQ=255)
                                                                                              gcgcgcAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCAt                                                 <  1:526667/150‑1 (MQ=255)
                                                                                              gcgcgcAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATc                                                >  2:800796/1‑151 (MQ=255)
                                                                                                 cgcAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATcgcg                                             >  2:924558/1‑151 (MQ=255)
                                                                                                  gcAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATCTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGa                                            <  1:538555/151‑1 (MQ=255)
                                                                                                    aGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATcgc                                              >  2:32837/1‑147 (MQ=255)
                                                                                                    aGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGaa                                           >  2:365069/1‑150 (MQ=255)
                                                                                                     gTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAg                                         <  1:805355/151‑1 (MQ=255)
                                                                                                          cGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGt                                    >  2:562490/1‑151 (MQ=255)
                                                                                                           gTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGtt                                   >  2:275003/1‑151 (MQ=255)
                                                                                                           gTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGtt                                   <  2:538748/151‑1 (MQ=255)
                                                                                                             ggTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGttt                                  <  2:503171/150‑1 (MQ=255)
                                                                                                              gTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTg                                >  1:486689/1‑151 (MQ=255)
                                                                                                              gTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTg                                <  1:19018/151‑1 (MQ=255)
                                                                                                               ttATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTg                                >  2:970003/1‑150 (MQ=255)
                                                                                                                    aGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGaa                         <  1:931032/151‑1 (MQ=255)
                                                                                                                      taataaTTAGTGCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGaa                         >  2:923440/1‑150 (MQ=255)
                                                                                                                                 gCAAAATGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGAAATTTAAGAGTAt             >  1:344399/1‑151 (MQ=255)
                                                                                                                                   aaaaTGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGAAATTTAAGAGTa              >  1:367188/1‑148 (MQ=255)
                                                                                                                                   aaaaTGATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGAAATTTAAGAGTATGa           <  2:613009/151‑1 (MQ=255)
                                                                                                                                        gATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGAAATTTAAGAGTATGAGCACt      >  2:470304/1‑151 (MQ=255)
                                                                                                                                        gATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGAAATTTAAGAGTATGAGCACt      <  1:924641/151‑1 (MQ=255)
                                                                                                                                        gATTTGCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGAAATTTAAGAGTATGAGCACt      <  1:43069/151‑1 (MQ=255)
                                                                                                                                             gCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGAAATTTAAGAGTATGAGCACTATct  <  1:923523/150‑2 (MQ=255)
                                                                                                                                             gCGTTATTGGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGAAATTTAAGAGTATGAGCACTATCg  <  1:564365/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
ATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTTAGTTGCAAATTTTTCAACATTTTATACACTACGAAAACCATCGCGAAAGCGAGTTTTGATAGGAAATTTAAGAGTATGAGCACTATCG  >  NC_000913/1151338‑1151627

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